TWI344753B - Decimation filter - Google Patents
Decimation filter Download PDFInfo
- Publication number
- TWI344753B TWI344753B TW096140522A TW96140522A TWI344753B TW I344753 B TWI344753 B TW I344753B TW 096140522 A TW096140522 A TW 096140522A TW 96140522 A TW96140522 A TW 96140522A TW I344753 B TWI344753 B TW I344753B
- Authority
- TW
- Taiwan
- Prior art keywords
- filter
- coefficients
- coefficient
- computing devices
- memory
- Prior art date
Links
- 230000015654 memory Effects 0.000 claims description 44
- 238000006243 chemical reaction Methods 0.000 claims description 27
- 238000005070 sampling Methods 0.000 description 19
- 238000010586 diagram Methods 0.000 description 7
- 239000004744 fabric Substances 0.000 description 7
- 238000001914 filtration Methods 0.000 description 5
- 238000000034 method Methods 0.000 description 5
- 230000006870 function Effects 0.000 description 4
- 239000003638 chemical reducing agent Substances 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 230000009467 reduction Effects 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 235000011389 fruit/vegetable juice Nutrition 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 241000282320 Panthera leo Species 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000001934 delay Effects 0.000 description 1
- 230000035784 germination Effects 0.000 description 1
- 238000000227 grinding Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 238000010079 rubber tapping Methods 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
- H03H17/00—Networks using digital techniques
- H03H17/02—Frequency selective networks
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
- H03H17/00—Networks using digital techniques
- H03H17/02—Frequency selective networks
- H03H17/06—Non-recursive filters
- H03H17/0621—Non-recursive filters with input-sampling frequency and output-delivery frequency which differ, e.g. extrapolation; Anti-aliasing
- H03H17/0635—Non-recursive filters with input-sampling frequency and output-delivery frequency which differ, e.g. extrapolation; Anti-aliasing characterized by the ratio between the input-sampling and output-delivery frequencies
- H03H17/065—Non-recursive filters with input-sampling frequency and output-delivery frequency which differ, e.g. extrapolation; Anti-aliasing characterized by the ratio between the input-sampling and output-delivery frequencies the ratio being integer
- H03H17/0664—Non-recursive filters with input-sampling frequency and output-delivery frequency which differ, e.g. extrapolation; Anti-aliasing characterized by the ratio between the input-sampling and output-delivery frequencies the ratio being integer where the output-delivery frequency is lower than the input sampling frequency, i.e. decimation
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Mathematical Physics (AREA)
- Complex Calculations (AREA)
- Analogue/Digital Conversion (AREA)
Description
1344753 九、發明說明: (優先權之主張) 本申請案主張於2006年η月9日向曰本專利局所提 出之日本專利申請案第2006-303461號之優先權。該優先 權申請案之整個内容併入本案作為參考。 【發明所屬之技術領域】 本發明係關於降頻濾波器,其將訊號頻率以預定之比 率轉換至較低頻率。 •【先前技術】 於訊號處理中,有時需要改變以預定的取樣頻率所取 樣之訊號的取樣率。-升高頻率之系統稱之為升頻器 (interpolator ) ’而降低頻率之系統稱之為降頻器 (decimator)。本發明係相關於降頻器。 於降頻器中,對資料貫施濾、波處理,然後資料被降頻 乂實現所希望之處理。於減少取樣(d〇wn_sampi丨叫)至較低 籲頻率的情形中,當包含高於較低頻率之1/2 (尼奎斯特頻 率(Nyquist freqUenCy))的成分時,會產生重疊失真 (foldover distortion)(頻疊成分(aliasing component))。因此,係執行濾波處理以便將頻寬限制至 尼奎斯特頻率或更低。用於濾波處理之濾波器稱之為降頻 遽波器(decimation filter)。 當降頻器依照該原理而實現降頻時,能轉換取樣率而 不需使用類比訊號。然而,有譬如處理負擔增加之問題, 這是因為於減少取樣之前必須以高取樣率執行缝處理之 319700 5 1344753 •故,而不能使用理想的濾波器。因此,係施行計算處理, 同時將轉移函數(transfer functi〇n) H(z)分成對應於 ,頻比率Μ的複數個濾波器係數群。此濾、波器組構稱之為 多相組構(p〇lyphase c〇nfigurati〇n),而具有該組構之濾 波器稱之為多相濾波器(polyphase fi Iter)。 藉由從原來的濾波器係數中於每一個M數抽取係數而 • 2得於多相組構中之濾波器係數。當輸入資料至濾波器之 刖施行減少取樣時,濾波器能以由該減少取樣所引起之低 籲操作速度進行操作。因此,能組構有效的濾波器。 第4圖為顯示降頻器例子之圖示,該降頻器由多相組 構所實現。參照圖式,由輸入侧以較高取樣頻率Fs輸 入之訊號係從輸出側51以較低取樣頻率Fd輸出。延遲元 件(delay element)52係為將訊號延遲l/Fs週期(對應於 一個頻率)的移位暫存器(shift register)。減少取樣器 (down-sampler)53 (於圖式中由「n>」表示)係減少取樣 籲輸入訊號至降頻比率Μ。藉由令轉移函數Η (ζ)經過多相 分解(polyphase decomposition),而獲得多相濾波器54 之Do ( z )…Dd ( z),而該多相濾波器54之Do ( z )…Dh (z)係構成一群濾波器,其中濾波器係數係於每一 μ數被 選擇。加法器55係將該等濾波之訊號相加於彼此。 舉例而言,假定其分接數(tap number)為η、且以 1 : 3執行降頻的濾波器之轉移函數η ( ζ )係由下列公式所 表示(“*表示乘法符號,,): 6 319700 1344753 H(z) = h0.+ h,z-1 + h2*z,2 + h3*z-3 + h4*z-< + …+ hn-2*zn-2 +</ RTI> </ RTI> <RTIgt; </ RTI> <RTIgt; </ RTI> <RTIgt; </ RTI> <RTIgt; </ RTI> <RTIgt; </ RTI> </ RTI> </ RTI> </ RTI> </ RTI> </ RTI> <RTIgt; The entire content of this priority application is incorporated into this case for reference. BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a down-conversion filter that converts a signal frequency to a lower frequency at a predetermined ratio. • [Prior Art] In signal processing, it is sometimes necessary to change the sampling rate of the signal sampled at a predetermined sampling frequency. - A system that raises the frequency is called an interpolator. The system that reduces the frequency is called a decimator. The invention is related to a frequency reducer. In the down-converter, filtering and wave processing are performed on the data, and then the data is down-converted to achieve the desired processing. In the case of reducing the sampling (d〇wn_sampi squeaking) to the lower frequency, when the component is higher than 1/2 (Nyquist freqUenCy) of the lower frequency, overlapping distortion occurs ( Foldover distortion) (aliasing component). Therefore, filtering processing is performed to limit the bandwidth to the Nyquist frequency or lower. The filter used for the filtering process is called a decimation filter. When the downconverter is down-converted according to this principle, the sampling rate can be converted without using an analog signal. However, there is a problem such as an increase in the processing load, because the slit processing must be performed at a high sampling rate before the sampling is reduced 319700 5 1344753, and the ideal filter cannot be used. Therefore, the calculation process is performed, and the transfer function (transfer functi〇n) H(z) is divided into a plurality of filter coefficient groups corresponding to the frequency ratio Μ. This filter and waver structure is called a polyphase fabric (p〇lyphase c〇nfigurati〇n), and the filter having this structure is called a polyphase fiiter. The filter coefficients in the polyphase fabric are obtained by decimation of the coefficients from each of the M coefficients of the original filter coefficients. When the input data is applied to the filter to reduce the sampling, the filter can operate at the low operating speed caused by the reduced sampling. Therefore, an effective filter can be constructed. Figure 4 is a diagram showing an example of a downconverter implemented by a polyphase fabric. Referring to the drawing, the signal input from the input side at the higher sampling frequency Fs is output from the output side 51 at the lower sampling frequency Fd. The delay element 52 is a shift register that delays the signal by a period of 1/Fs (corresponding to a frequency). The down-sampler 53 (indicated by "n>" in the figure) is to reduce the sampling and call the input signal to the down-conversion ratio Μ. Do (z)...Dd(z) of the polyphase filter 54 is obtained by causing the transfer function Η(ζ) to undergo polyphase decomposition, and the Do(z)...Dh of the polyphase filter 54 is obtained. (z) constitutes a group of filters in which the filter coefficients are selected for each μ number. The adder 55 adds the filtered signals to each other. For example, a transfer function η ( ζ ) of a filter whose filter is assumed to have a tap number of η and is down-converted by 1:3 is represented by the following formula ("* indicates a multiplication symbol,"): 6 319700 1344753 H(z) = h0.+ h,z-1 + h2*z,2 + h3*z-3 + h4*z-< + ...+ hn-2*zn-2 +
Wz11-1 (a) 則’形成三個多相濾波器54或形成為j)D ( Z )至〇2 ( Z ), 並由下列公式所表示(於此情況分接數是3之倍數): D〇(z) = h〇 + h3*z_l + h6*z'2 +...+ h„-3*z'n/3+1 (b)Wz11-1 (a) then 'forms three polyphase filters 54 or is formed as j) D ( Z ) to 〇 2 ( Z ) and is represented by the following formula (in this case, the number of taps is a multiple of 3) : D〇(z) = h〇+ h3*z_l + h6*z'2 +...+ h„-3*z'n/3+1 (b)
Di(z) = hi + h4*z_1 + h7*z'2 +...+ hn-2*z'n/3+1 (c) d2(z) = h2 + hs^z-*1 + hB*z'2 +...+ hft_i*z"n/3+1 (d) 以此方式,多相濾波器54於減少取樣後能夠以所獲 知的處理率執行濾波計算。因此,能夠減少處理負擔。諸 輸出由加法器55加在一起,且最後的減少取樣訊號係從輸 出侧51輸出。於多相濾波器54之計算中,可使用ρ IR濾 皮器(有限脈衝反應濾'波器(Finite Impulse Response er) )、11R ;慮波器(無限脈衝反應濾波器(I n ^ i n i te impulseResponseFilter))、或 FFT (快速傅利葉轉換) 計算。 • [非專利參考多比率訊號處理(Multi rate signalDi(z) = hi + h4*z_1 + h7*z'2 +...+ hn-2*z'n/3+1 (c) d2(z) = h2 + hs^z-*1 + hB *z'2 +...+ hft_i*z"n/3+1 (d) In this manner, the polyphase filter 54 can perform the filtering calculation at the learned processing rate after the sampling is reduced. Therefore, the processing load can be reduced. The outputs are added together by an adder 55, and the last downsampled signal is output from the output side 51. In the calculation of the polyphase filter 54, a ρ IR filter (Finite Impulse Response er), 11R, and a filter (infinite impulse response filter (I n ^ ini te) can be used. impulseResponseFilter)), or FFT (Fast Fourier Transform) calculation. • [Non-patent reference multi-rate signal processing (Multi rate signal
Processing),,(SH0K0D0 有限公司,ΚΙγΑ,, 尤其第4章。 然而,於具有上述相關技術組構之降頻器中,當將 執行:頻比率大於100之大降頻時,需要正比於降;比 之大置之延遲元件52 ’而增加電路尺寸(當如上述例子 頻比率為3時’所需要的延遲元件數為3)。再者 =相遽波器54之計算尺寸。於如此情況,係使用—種 將降頻器分成數級之同時執行降頻的方法。 319700 7 相較之下,於譬如儀器 沒右宣止4 ~ 儎态的應用中,有降頻比率很大而 及有事先決定的情況。於此 之組構並非總是適當的。因&,:難成多個級 現所希望料。 料供-種料㈣實 【發明内容】Processing),,(SH0K0D0 Co., Ltd., ΚΙγΑ,, especially Chapter 4. However, in the downconverter having the above-mentioned related art fabric, when performing a frequency reduction ratio greater than 100, it is necessary to be proportional to the drop. Increasing the circuit size (when the frequency ratio is 3 as in the above example), the number of delay elements required is 3). Again = the calculated size of the phase chopper 54. It is a method of performing frequency reduction while dividing the frequency reducer into several stages. 319700 7 In contrast, in applications where the instrument does not have a 4~ state, there is a large frequency reduction ratio. There is a situation in advance. The organization of this is not always appropriate. Because &,: It is difficult to achieve multiple levels of what is expected. Material supply - seed material (four) real [invention content]
替之例示實施例提供—種降頻濾波器,其將輸入 降:轉換至較低頻率,並具有能夠容易處理任意的 降頻比率而無需增加硬體的濾波器組構。 依照本發明之—個或多個實施例,降頻濾波器包括: 艘個計算裝置,各具有乘法器(multiplier)和累加器 (aCCUD1Ulat〇r);複數個係數記憶體(coefficient memory) ’其係儲存濾波器係數並分別對應於計算裝置;以 及選擇器’其與時脈訊號同步地依序選擇性地輸出該複數 個計算裝置之輸出;其中,當降頻比率為η時,被依序移 位達η個濾、波器係數之錢器係數係從該複數個係數記憶Instead of the exemplary embodiment, a down-converting filter is provided that converts the input down: to a lower frequency and has a filter fabric that can easily handle any down-conversion ratio without adding hardware. In accordance with one or more embodiments of the present invention, the down-converting filter includes: a plurality of computing devices each having a multiplier and an accumulator (aCCUD1Ulat〇r); a plurality of coefficient memories 'its And storing the filter coefficients respectively corresponding to the computing device; and the selector 'selectively outputs the output of the plurality of computing devices in synchronization with the clock signal; wherein, when the down-conversion ratio is η, the sequence is sequentially The money coefficient of shifting up to n filters and wave coefficients is recovered from the plurality of coefficients
體被讀出,並乘上該計算裝置之乘法器中之訊號,而相乘 ,之結果係累加於該累加器中以待輸出。依照此組構,能 提供能夠容易處理任意的降頻比率之攄波器組構而無需增 加硬體。 複數個係數記憶體可具有:環式記憶體 memory),其儲存所有先前已計算過的濾波器係數,並且從 該環式記憶體以與該時脈訊號同步之方式依序循環地讀出 該濾波器係數;以及複數個移位暫存器,其梯級連接 reascade-connected)至環式記憶體;並且,該等移位暫存 319700 8 丄 • Γ、有‘降頻比率為11時能夠儲存η個濾波器係數的容 環式記憶體讀出之遽波器係數係儲存入該複數個 ^立暫存器以及從該複數個移位暫存器讀取出,同時依序 ^位該濾、波器係數。依照該組構,能夠大量地減少於該 4係數記憶體中所需之儲存容量。 較^的情況是’計算裝置之數目等於或大於藉由將遽 ^係數之數目除以該降頻比率所獲得之值。或者較佳的 情兄是儲存於係數記憶體中之遽波器係數之數目等於或小 於藉由將計算裝置之數目乘上該降頻比率所獲得之值。依 照該組構,能約有效地使用所提供之計算裝置之數目而執 行計算。 本發明能夠提供-種降_波器,其將輸人訊號之頻 率轉換至較低頻率,並具有能夠容易處理任意的降頻比率 而無需增加硬體之濾波器組構。 由下列之詳細說明、所附的圖式和申請專利範圍,其 籲他的特徵和優點可以很清楚。 【實施方式】 _現將說明本發明之降㈣波器之實施例。帛!圖為顯 不降頻濾波器之實施例之組構之圖示,而第2圖為顯示取 樣訊號和濾波器係數之間的關係之圖示。本發明之降頻濾 波器非為多相濾波器。於下列實施例中所表示之係數、除 數(division number)、其他特定的數目等僅係例示用來輔 助了解本發明者,除非有特別說明否則該等數值並不限制 本發明。 319700 9 丄 JH·件/:):) 於第1圖中戶斤千夕ι _ _ _ β不之降頻濾波器中,從輸入側10以較 南取樣頻率Fs郎·&λ a φ , ± 厅輸入之訊號係從輸出侧11以較低取樣頻 率F d輪出。降萌、、磨、由 味领慮波裔包括:複數個計算裝置20,各由 乘法器21和累加51 99。α 、斋22 ν且構而成;環式記憶體3〇,其為係 =體之例子;複數個移位暫存器3卜其為係數記憶體 •^歹,以及選擇器13 ’其依序地選擇計算裝置2〇之輸 出以將其輸出。 %式記憶體30能夠儲存先前已計算之所有的遽波器 糸、> 而該等係數能夠以與時脈訊號同步之方式依序循環 地被讀出°移位暫存器31係以梯級(caseade)方式(階層 、)連接至環式記憶體30,俾將從環式記憶體30讀出 、’《皮器係數儲存入該複數個移位暫存器,以及從該複數 移暫存益予以讀取出,同時依序地移位該濾波器係數。 以下將詳細說明環式記憶體30和移位暫存器31之操 作。虽從%式記憶體30讀出濾波器係數而輸入對應之計算 。 時濾波益係數之資料係同時被送至第一移位暫存 器3:1。、於環式記憶體30巾,次一個待讀取之遽波器係數 ^位址被移位。被送出之濾波器係數係儲存於各移位暫存 斋31中並同時依序被移位,而在當移位暫存器31之容量 被填滿時的時點,則執行從移位暫存器31讀取至計算裝置 20的處理。與此讀取同時地,濾波器係數之資料係被送至 二個移位暫存器3卜重複此程序,以使得濾波器係數之 資料在移位暫#|| 31巾以有如濾波器健被往前推動之 方式移動。 319700 10 於降頻比率為n之情況中,環式記憶體30係儲存n * k in個濾波器係數。此移位暫存器μ具有能 個遽波器係數的容量(深度)。係數k表示計算裝置2〇 、、、’且疋決疋濾波器分接長度(tap length)之係 數並且為例如大約24之值。係數出為整數,其不少於◦',、 以及為了防止產生超額之乘法器21,較佳為不大於〇二卜 免、J佳匱况疋’計算裝置20之數目等於或大於藉由將The body is read and multiplied by the signal in the multiplier of the computing device, and multiplied, and the result is accumulated in the accumulator for output. According to this configuration, it is possible to provide a chopper configuration that can easily handle an arbitrary down-conversion ratio without adding hardware. The plurality of coefficient memories may have: a memory (memory) storing all previously calculated filter coefficients, and sequentially reading the loop memory from the loop memory in synchronization with the clock signal Filter coefficients; and a plurality of shift registers whose steps are connected to reascade-connected) to the ring memory; and, the shift registers 319700 8 丄• Γ, can be stored when the down-conversion ratio is 11. The chopper coefficients of the n-th filter coefficients are stored in the plurality of registers and read from the plurality of shift registers, and the filter is sequentially , wave coefficient. According to this configuration, the required storage capacity in the 4-factor memory can be largely reduced. The case of ^ is that the number of computing devices is equal to or greater than the value obtained by dividing the number of 遽^ coefficients by the down-conversion ratio. Or preferably, the number of chopper coefficients stored in the coefficient memory is equal to or less than the value obtained by multiplying the number of computing devices by the down-conversion ratio. According to this configuration, the calculation can be performed approximately efficiently using the number of computing devices provided. The present invention is capable of providing a down-waveform that converts the frequency of the input signal to a lower frequency and has a filter fabric that can easily handle an arbitrary down-conversion ratio without adding hardware. The features and advantages of the invention are apparent from the following detailed description, appended claims and claims. [Embodiment] An embodiment of the falling (four) wave device of the present invention will now be described. silk! The figure shows an illustration of the configuration of an embodiment of the reduced down filter, and Fig. 2 is a diagram showing the relationship between the sampled signal and the filter coefficients. The down-converting filter of the present invention is not a polyphase filter. The figures, the division numbers, the other specific numbers, and the like, which are represented in the following examples, are merely intended to be illustrative of the present invention and are not intended to limit the invention unless otherwise specified. 319700 9 丄JH·pieces/:):) In the down-converter filter of Figure 1 in Figure 1, from the input side 10 to the south sampling frequency Fs Lang·&λ a φ The signal input to the hall is rotated from the output side 11 at a lower sampling frequency Fd. The germination, grinding, and taste of the wave include: a plurality of computing devices 20, each of which is multiplier 21 and accumulates 51 99. α, 斋 22 ν and constructed; ring memory 3〇, which is an example of a system = body; a plurality of shift registers 3 are coefficient memory • ^ 歹, and selector 13 ' The output of the computing device 2 is sequentially selected to output it. The %-type memory 30 is capable of storing all of the previously calculated choppers, > and the coefficients can be sequentially cyclically read out in synchronization with the clock signal. The shift register 31 is stepped. The (caseade) mode (hierarchy) is connected to the ring memory 30, and the 俾 will be read from the ring memory 30, 'the skin coefficient is stored in the plurality of shift registers, and the temporary number is temporarily stored. It is read out and the filter coefficients are shifted sequentially. The operation of the ring memory 30 and the shift register 31 will be described in detail below. The filter coefficients are read from the % memory 30 and the corresponding calculation is input. The data of the time filter benefit coefficient is simultaneously sent to the first shift register 3:1. In the ring memory 30 towel, the next chopper coefficient to be read ^ address is shifted. The filter coefficients sent out are stored in each shift temporary storage 31 and are sequentially shifted at the same time, and when the capacity of the shift register 31 is filled, the shift from the shift is performed. The device 31 reads the processing to the computing device 20. Simultaneously with this reading, the data of the filter coefficients is sent to the two shift registers 3 to repeat the procedure, so that the data of the filter coefficients are shifted in the temporary #|| Moved in a way that is pushed forward. 319700 10 In the case where the down-conversion ratio is n, the ring memory 30 stores n * k in filter coefficients. This shift register μ has a capacity (depth) capable of chopper coefficients. The coefficient k represents the calculation means 2 〇 , , , ' and the coefficient of the filter tap length and is, for example, a value of about 24. The coefficient is an integer, which is not less than ◦', and in order to prevent the generation of the excess multiplier 21, preferably not more than 〇2, the number of computing devices 20 is equal to or greater than
濾波器係數之數目除以該降頻比率所獲得之值。換言之, 較佳情況是’待儲存^環式記憶體3G中之濾波器係數之數 目專^或小於藉由將計算裝X 20之數目乘上該降頻比率 所獲%•之值。依照此組構,能有效地使用計算裝置所提供 之數目而執行計算。 y、 輸入側1 0輸入之訊號係被供應至該複數個計算裝 置2〇»各者之乘法11 21之一個輸入。環式記憶冑30和移位 暫存盗31之輸出係分別連接至乘法器21之另—輸入,俾 執行:入訊號與各個濾波器係數的相乘。相乘結果由各個 、W 2累加直到元成濾波器分接長度之一系列計算為 止,並且當連接選擇器13時,輸出至輸出侧u。 一現將參照帛2圖說明降頻濾、波器之操作之例子。舉例 :言’茲假設降頻比率是3,係數k為24,係數!„為〇。 ”Ί /慮波器分接長度為72。因此,環式記憶體3〇具有 用於儲存72個分接之遽波器係數之容量,而移位暫存器 31具有,於儲存3個分接之濾波器係數之容量。於第^圖 十开裝置20係由從左邊開始依序由2〇_ι、2〇-2、 319700 11 1344753 20-3、…、20-24所表示。於第2圖中,滤波器係數之72 個分接係由h0、hi、h2、…、h71所表示,而於取樣點之 輸入訊號係由d0、dl、h2、…所表示。 於上述狀況下,計算裝置20-1、20-2、20-3、20-24 分別執行由下列公式(e) 、(f) 、(g) 、(h)所表示 •之計算: (20-1) d〇*h〇 + dx*hi + d2*h2 + d3*h3 +...+ d69*h6? + d7〇*h7fr + dn*h71 (e) φ (20-2) d3*h〇 + d4*hi + d5*h2 + d«*h3 +...+ d72*h69 + d73*h7〇 + d74*h7i (f) (20-3) d6*h〇 + d7*hi + d8*h2 + d9*h3 +...+ d75*h69 + d76*h7〇 + dr?*h7i (g) (20-4) d69*h〇 + d7〇*hi + d7i*h2 + d72*h3 +...+ d138*h69 + di39*h7〇 + di4〇*h?i (h) 在由計算裝置20-1從開始計算經過n個時脈(於此實施例 中為3個時脈)後,計算裝置20-2開始計算。同樣地,計 算裝置20於每η個時脈開始操作。 * 於計算裝置20-1中’當由公式(e)所表示之計算結 束時,重設累加器22,然後於時序d72接續開始由公式(i) .所表示之濾波器計算。於計算裝置20-2至20-24中’同樣 地,處理接續的資料,藉此能夠進行連續的頻率降頻。 ^2〇-1) d72*li〇 + ^73*hi + d74*h2 + d75*h3 +...+ d14a*h69 + di42*h7〇 + di43*hn ⑴ 在當計算裝置20_1結束72個分接之計算時之時間 热,選擇器!3連接至計算裝置2G-1。然後,以1/3取樣 12 319700 ^^4/53 之比率’計算&置2G被依序地切換 管輸出累加之資料。也就是說,在當公式⑷之: ::束時之時間點’選擇計算裝置2(M之輸出。在從該時 .·溫過3個時脈後,遠擇計算裝置2〇 —2之輸出,並且之 ^依序地選擇剩餘之計算裝置,直到計算裝i 2㈣為The number of filter coefficients divided by the value obtained by the down-conversion ratio. In other words, it is preferable that the number of filter coefficients in the to-be-stored memory 3G is smaller or smaller than the value obtained by multiplying the number of the calculations X 20 by the down-conversion ratio. According to this configuration, the calculation can be performed efficiently using the number provided by the computing device. y. The signal input to the input side 10 is supplied to one of the plurality of multiplications 11 21 of the plurality of computing devices 2». The output of the ring memory cassette 30 and the shift register 31 are respectively connected to the other input of the multiplier 21, and the execution is performed by multiplying the input signal by each filter coefficient. The multiplication result is calculated by the sum of each of W 2 until the tapping length of the elementary filter is calculated, and when the selector 13 is connected, it is output to the output side u. An example of the operation of the down-converting filter and the waver will be described with reference to FIG. Example: Words assume that the down-conversion ratio is 3, the coefficient k is 24, and the coefficient! „为〇.”Ί/wave filter tap length is 72. Therefore, the ring memory 3A has a capacity for storing 72 tapped chopper coefficients, and the shift register 31 has a capacity for storing 3 tapped filter coefficients. In the figure, the ten-opening device 20 is represented by 2〇_ι, 2〇-2, 319700 11 1344753 20-3, ..., 20-24 in order from the left. In Fig. 2, the 72 tap coefficients of the filter coefficients are represented by h0, hi, h2, ..., h71, and the input signals at the sampling points are represented by d0, dl, h2, .... Under the above conditions, the computing devices 20-1, 20-2, 20-3, and 20-24 respectively perform calculations represented by the following formulas (e), (f), (g), and (h): (20) -1) d〇*h〇+ dx*hi + d2*h2 + d3*h3 +...+ d69*h6? + d7〇*h7fr + dn*h71 (e) φ (20-2) d3*h 〇+ d4*hi + d5*h2 + d«*h3 +...+ d72*h69 + d73*h7〇+ d74*h7i (f) (20-3) d6*h〇+ d7*hi + d8* H2 + d9*h3 +...+ d75*h69 + d76*h7〇+ dr?*h7i (g) (20-4) d69*h〇+ d7〇*hi + d7i*h2 + d72*h3 +. ..+d138*h69 + di39*h7〇+ di4〇*h?i (h) After n calculations (three clocks in this embodiment) from the start of calculation by the computing device 20-1, Computing device 20-2 begins the calculation. Similarly, computing device 20 begins operation every n clock cycles. * In the computing device 20-1, when the calculation represented by the formula (e) is completed, the accumulator 22 is reset, and then the filter calculation represented by the formula (i) is successively started at the timing d72. Similarly, in the computing devices 20-2 to 20-24, the connected data is processed, whereby continuous frequency down-conversion can be performed. ^2〇-1) d72*li〇+ ^73*hi + d74*h2 + d75*h3 +...+ d14a*h69 + di42*h7〇+ di43*hn (1) At the end of the computing device 20_1, 72 points The time of the calculation is hot, the selector! 3 is connected to the computing device 2G-1. Then, the ratio of 1/3 sampling 12 319700 ^^4/53 'calculated & 2G is sequentially switched to the output of the tube output. That is to say, at the time point of the formula (4): :: beam time 'select the output of the computing device 2 (M.) After the temperature has passed 3 hours, the remote computing device 2〇2 Output, and select the remaining computing devices sequentially until the calculation of i 2 (four) is
之^後,再選擇計算裝置20-1之輸出。在計算裝置20-24 之汁算(h)結束後經過3個時脈後之時間點結 置20-1之次一個計算。 四於如此組構之降頻濾波器中,藉由僅調整係數記憶體 式記憶體30和移位暫存器31)之深度(容量)而能 夠奋·易組構任意降頻比率之濾波器。於相關技術例子中, 必須改變延遲元件52和多相濾波器54之數目,而此改變 係=及電路之大規模修改。因此,很難實現任意的降頻滤 f态。依照本發明,相較之下’能提供其能夠容易處理任 忍的降頻比率而無需增加硬體的濾波器組構。 尤其疋,該係數記憶體係由一個環式記憶體3〇和複 數個移位暫存器31所組構,而且該濾波器係數資料在被依 序移位的同時被儲存入該係數記憶體而輸入至後續的級 (=tage)。因此,各係數記憶體不需要令所有的濾波器係數 負料對應於汁算裝置,而能夠大量地減少於該等係數記憶 體中所需之儲存容量。 舉例而言,上述降頻濾波器能夠由Fpga組構。近年 來的FPGA包含1〇〇個或更多組之此種乘法器21和記憶體 (是式§己憶體30和移位暫存器31),並因此適宜用來實 13 319700 1344753 現此組構。當濾波器係藉由使用FPGA而組構時,計算裝置 2〇之數目和記憶體之容量(深度)能以動態方式適當地設 定。 計算裝置20之數目係對應於濾波器分接長度。能藉 由增加或減少計算裝置20之數目而調整所需之濾波器特 性。在截止頻率區(cut-off frequency region)需要陡山肖 之暫態特性的情況中,計算裝置20之數目會增加。相比之 下’在濾波器狀態能放鬆至允許程度之情況中,能夠減少 。十鼻裝置2 0之數目,俾能夠組構能經濟地實現之降頻濾波 器。 /’、’ 計算裝置20總是以不產生浪費的時間週期的方式來 進行操作,並能夠以較具有一般多核心組構之電腦執行軟 ^處理情況之計算處理速度為高之計算處理速度來實施計 算處理。於譬如1C測試器或記憶體測試器之儀器中,在輸 入訊號必須依照將被測量之目標以不同的比率降頻之情 況,尤其疋,濾波器係依上述方式而組構,藉此能夠實 以即時操作之高速和高度準確之濾波器。 [其他實施例] 、於上述實施例中’係數記憶體係由—個環式記憶體 和複數個移位暫存11 31所組構。可取而代之者為,所有 糸數-己隐體可由%式記憶體所組構,由此實現本發明。 第3圖為顯示降頻濾波器之另—個組構例子 =圖中’未設置移位暫存器31,而所有的係數記^ 係由壤式記3〇所組構。各環式記·%係儲存^ 14 獅。; 1344753 所計算的所有濾波器係數。在降頻比率為η之情況中,彼 此移位η個濾波器係數之濾波器係數係依序地從各環式記 * 憶體30讀出。於各計算裝置2〇中,山 上丄 1开衣罝π甲所讀出之濾波器係數 和訊號係由乘法器21彼此相乘’並累加於累加器22中, 然後輸出,由此能執行如上述實施例中的相同計算處理。 當降頻比率為1時,本發明之降頻遽波器作用為-般 •的FIR濾波器,而因此亦能用作為即時操作之高速和高度 準確的FIR濾波器。 籲耗本發明已參考所附圖式參照較佳實施例而加以 說明,但是本發明當然並不受該等實施例之限制。很顯然 的熟悉此項技術者在所附申請專利範圍内可作各種的修飾 和變化,而這些修飾和變化將落在本發明之技術範圍内。 本發明能用作為降頻滤波器,該降頻據波器將訊號之 頻率以一定的比率轉換至較低頻率。 【圖式簡單說明】 # 第1圖為顯示降頻濾波器之實施例之組構之圖示。 一第2圖為顯示取樣訊號和濾波器係數之間關係之圖 〇 第3圖為顯示降頻濾波器之另一個組構例子之圖示。 -第4 ®為顯示由多相組構所實現之降頻器例子之圖 不 〇 【主要元件符號說明】 輸入侧 輸出侧 - 319700 15 11 1344753 13 選擇器 20 計算裝置 計算裝置 20-1 、 20-2 、 20-3 、…、20-24 21 乘法器 22 累加器 30 環式記憶體 31 移位暫存器 50 輸入側 51 輸出侧 52 延遲元件 53 減少取樣器 54 多相濾波器 55 加法器 h0、hi、h2.....h71 分接 dO、dl、h2、…輸入訊號 Fs 取樣頻率After that, the output of the computing device 20-1 is selected. The next calculation of the time point 20-1 after the passage of 3 clocks after the end of the juice calculation (h) of the computing device 20-24. In the down-converting filter thus constructed, it is possible to construct a filter having an arbitrary down-conversion ratio by adjusting only the depth (capacity) of the coefficient memory type memory 30 and the shift register 31). In the related art example, the number of delay elements 52 and polyphase filters 54 must be changed, and this change is a large-scale modification of the circuit. Therefore, it is difficult to achieve arbitrary down-conversion filtering. According to the present invention, it is possible to provide a filter structure which can easily handle the tolerable down-conversion ratio without adding hardware. In particular, the coefficient memory system is constructed by a ring memory 3〇 and a plurality of shift registers 31, and the filter coefficient data is stored into the coefficient memory while being sequentially shifted. Enter to the next level (=tage). Therefore, each coefficient memory does not need to have all of the filter coefficient negatives correspond to the juice calculation device, and can greatly reduce the required storage capacity in the coefficient memories. For example, the down-converting filter described above can be constructed by Fpga. In recent years, FPGAs have included one or more such multipliers 21 and memories (the § 己 体 30 and the shift register 31), and thus are suitable for real 13 319700 1344753. Fabrication. When the filter is constructed by using an FPGA, the number of computing devices 2 and the capacity (depth) of the memory can be appropriately set in a dynamic manner. The number of computing devices 20 corresponds to the filter tap length. The required filter characteristics can be adjusted by increasing or decreasing the number of computing devices 20. In the case where the cut-off frequency region requires a transient characteristic of the steep mountain, the number of computing devices 20 increases. In contrast, in the case where the filter state can be relaxed to the allowable level, it can be reduced. With the number of 10 nose devices, 俾 can construct a down-converting filter that can be economically realized. /', 'The computing device 20 always operates in a manner that does not generate a wasted time period, and can calculate the processing speed at a higher processing speed than a computer having a general multi-core configuration. Implement calculation processing. In the instrument such as the 1C tester or the memory tester, the input signal must be down-converted according to the target to be measured at different ratios. In particular, the filter is constructed in the above manner. High speed and highly accurate filter for instant operation. [Other Embodiments] In the above embodiment, the coefficient memory system is constructed by a ring memory and a plurality of shift registers 11 31 . Alternatively, all of the parameters - the hidden bodies can be constructed from % memory, thereby implementing the present invention. Fig. 3 is a diagram showing another configuration example of the down-converting filter. In the figure, the shift register 31 is not provided, and all the coefficient records are constructed by the magnetic type. Each ring type·% is stored ^ 14 lion. ; 1344753 All filter coefficients calculated. In the case where the down-conversion ratio is η, the filter coefficients which are shifted by n filter coefficients are sequentially read from the respective ring-shaped memories 30. In each of the computing devices 2, the filter coefficients and signals read by the 丄1 开 甲 甲 are multiplied by the multiplier 21 and accumulate in the accumulator 22, and then output, thereby being able to perform The same calculation process in the above embodiment. When the down-conversion ratio is 1, the down-converter chopper of the present invention functions as a general-purpose FIR filter, and thus can be used as a high-speed and highly accurate FIR filter for immediate operation. The invention has been described with reference to the preferred embodiments thereof, but the invention is of course not limited by the embodiments. It is obvious that those skilled in the art can make various modifications and changes within the scope of the appended claims, and such modifications and variations will fall within the technical scope of the present invention. The present invention can be used as a down-converting filter that converts the frequency of the signal to a lower frequency at a certain ratio. BRIEF DESCRIPTION OF THE DRAWINGS # Fig. 1 is a diagram showing the configuration of an embodiment of a down-converting filter. Figure 2 is a diagram showing the relationship between the sampling signal and the filter coefficients. Figure 3 is a diagram showing another example of the configuration of the down-converting filter. - 4th is a diagram showing an example of a downconverter implemented by a polyphase fabric. [Main component symbol description] Input side output side - 319700 15 11 1344753 13 Selector 20 Calculation device calculation device 20-1, 20 -2, 20-3, ..., 20-24 21 Multiplier 22 Accumulator 30 Ring Memory 31 Shift Register 50 Input Side 51 Output Side 52 Delay Element 53 Sample Reducer 54 Polyphase Filter 55 Adder H0, hi, h2.....h71 tap dO, dl, h2, ... input signal Fs sampling frequency
Fd 取樣頻率 16 319700Fd sampling frequency 16 319700
Claims (1)
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2006303461A JP4957194B2 (en) | 2006-11-09 | 2006-11-09 | Decimation filter |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| TW200828794A TW200828794A (en) | 2008-07-01 |
| TWI344753B true TWI344753B (en) | 2011-07-01 |
Family
ID=39370460
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| TW096140522A TWI344753B (en) | 2006-11-09 | 2007-10-29 | Decimation filter |
Country Status (4)
| Country | Link |
|---|---|
| US (1) | US8131790B2 (en) |
| JP (1) | JP4957194B2 (en) |
| KR (1) | KR100893740B1 (en) |
| TW (1) | TWI344753B (en) |
Families Citing this family (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP5522893B2 (en) * | 2007-10-02 | 2014-06-18 | キヤノン株式会社 | Image processing apparatus, image processing method, and program |
| TWI390841B (en) | 2009-10-14 | 2013-03-21 | Novatek Microelectronics Corp | Multi-rate filter bank |
| KR101317180B1 (en) * | 2009-12-16 | 2013-10-15 | 한국전자통신연구원 | Second Order Sinc Decimation filter |
| GB2487361A (en) * | 2011-01-17 | 2012-07-25 | Sony Corp | A digital interpolating FIR filter using fewer multipliers |
| JP5812774B2 (en) * | 2011-09-08 | 2015-11-17 | ルネサスエレクトロニクス株式会社 | Semiconductor device |
| KR101613521B1 (en) * | 2015-02-17 | 2016-05-02 | 서강대학교산학협력단 | Filter Assembly for medical image signal and dynamic decimation method using thereof |
| US10879877B1 (en) | 2018-09-28 | 2020-12-29 | The Mitre Corporation | Systems and method for a low power correlator architecture using distributed arithmetic |
| US10410700B1 (en) | 2018-09-28 | 2019-09-10 | The Mitre Corporation | Systems and method for a low-power correlator architecture using shifting coefficients |
| US11979131B2 (en) | 2021-05-19 | 2024-05-07 | Hughes Network Systems, LLC. | High-rate decimation filter with low hardware complexity |
| WO2022245531A1 (en) * | 2021-05-19 | 2022-11-24 | Hughes Network Systems, Llc | High-rate decimation filter with low hardware complexity |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE69024213T2 (en) * | 1990-09-18 | 1996-07-04 | Alcatel Nv | Multi-channel clock frequency reducer |
| JP3243831B2 (en) * | 1992-03-31 | 2002-01-07 | ソニー株式会社 | FIR type filter |
| JPH0677769A (en) | 1992-08-27 | 1994-03-18 | Fujitsu Ltd | Fir digital filter |
| US6427158B1 (en) * | 2000-12-14 | 2002-07-30 | Texas Instruments Incorporated | FIR decimation filter and method |
-
2006
- 2006-11-09 JP JP2006303461A patent/JP4957194B2/en not_active Expired - Fee Related
-
2007
- 2007-10-29 TW TW096140522A patent/TWI344753B/en not_active IP Right Cessation
- 2007-10-30 US US11/927,927 patent/US8131790B2/en not_active Expired - Fee Related
- 2007-11-08 KR KR1020070113824A patent/KR100893740B1/en not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JP2008124593A (en) | 2008-05-29 |
| US20080114821A1 (en) | 2008-05-15 |
| US8131790B2 (en) | 2012-03-06 |
| JP4957194B2 (en) | 2012-06-20 |
| KR100893740B1 (en) | 2009-04-17 |
| TW200828794A (en) | 2008-07-01 |
| KR20080042729A (en) | 2008-05-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| TWI344753B (en) | Decimation filter | |
| US7409417B2 (en) | Polyphase filter with optimized silicon area | |
| US7262716B2 (en) | Asynchronous sample rate converter and method | |
| US8631059B2 (en) | Methods or structures for reconstruction of substantially uniform samples from substantially nonuniform samples | |
| US5515402A (en) | Quadrature filter with real conversion | |
| US20130110898A1 (en) | Apparatus for signal processing | |
| JPH07500461A (en) | Folding and decimation filter architecture | |
| US6532273B1 (en) | Efficient polyphase decimation filter | |
| EP1372264A2 (en) | Low power decimation system and method of deriving same | |
| CN107979354B (en) | Multiphase device for sample rate conversion and method for manufacturing and testing the same | |
| TWI539742B (en) | Digital filter-decimator-tuner | |
| TW200414675A (en) | Digital filter design method and device, digital filter design program, and digital filter | |
| US7908306B1 (en) | SRC with multiple sets of filter coefficients in memory and a high order coefficient interpolator | |
| JP2008021119A (en) | Digital filter and image processor using the same | |
| KR100905153B1 (en) | Filters, Interpolation Filters, and Decimation Filters for Digital Up-Down Converters | |
| JP6311601B2 (en) | Multistage filter processing apparatus and method | |
| JP2008219560A (en) | Decimation filter | |
| JP4242362B2 (en) | Signal processing apparatus having post-processing | |
| CN120074513B (en) | Anti-aliasing downsampling circuit, circuit module, chip and signal processing device | |
| TWI482427B (en) | Time division decimation filter bank and time division decimation filtering method | |
| JP4413413B2 (en) | Semiconductor memory device and digital film | |
| JP2004128858A (en) | Fir digital filter | |
| EP1609241B1 (en) | Systems and methods for implementing a sample rate converter using hardware and software to maximize speed and flexibility | |
| US20100091829A1 (en) | Optimised architecture for a downsampling fir filter | |
| CN102201606B (en) | Time-sharing frequency-reducing filter and time-sharing frequency-reducing filtering method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| MM4A | Annulment or lapse of patent due to non-payment of fees |