TW201616810A - 有限脈衝響應濾波器與濾波方法 - Google Patents

有限脈衝響應濾波器與濾波方法 Download PDF

Info

Publication number
TW201616810A
TW201616810A TW103137590A TW103137590A TW201616810A TW 201616810 A TW201616810 A TW 201616810A TW 103137590 A TW103137590 A TW 103137590A TW 103137590 A TW103137590 A TW 103137590A TW 201616810 A TW201616810 A TW 201616810A
Authority
TW
Taiwan
Prior art keywords
coefficients
input values
impulse response
finite impulse
input
Prior art date
Application number
TW103137590A
Other languages
English (en)
Other versions
TWI566523B (zh
Inventor
呂明和
孫際恬
Original Assignee
財團法人工業技術研究院
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 財團法人工業技術研究院 filed Critical 財團法人工業技術研究院
Priority to TW103137590A priority Critical patent/TWI566523B/zh
Priority to US14/598,242 priority patent/US20160126933A1/en
Publication of TW201616810A publication Critical patent/TW201616810A/zh
Application granted granted Critical
Publication of TWI566523B publication Critical patent/TWI566523B/zh

Links

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
    • H03H17/06—Non-recursive filters
    • 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/0223—Computation saving measures; Accelerating measures
    • H03H2017/0245—Measures to reduce power consumption
    • 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/0223—Computation saving measures; Accelerating measures
    • H03H2017/0247—Parallel structures using a slower clock

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Mathematical Physics (AREA)
  • Complex Calculations (AREA)

Abstract

一種有限脈衝響應濾波器及其對應的濾波方法。上述有限脈衝響應濾波器接收輸入序列。輸入序列包括多個輸入值。上述有限脈衝響應濾波器包括至少一第一加法器、至少一乘法器、以及第二加法器。每一上述第一加法器以平行方式同時進行多個加法運算。每一上述加法運算輸出上述多個輸入值其中兩個之和數。每一上述乘法器以平行方式同時進行多個乘法運算。每一上述乘法運算輸出上述多個和數其中之一與有限脈衝響應濾波器的多個係數其中之一的乘積。第二加法器輸出上述多個乘積之總和。

Description

有限脈衝響應濾波器與濾波方法
本揭露是有關於一種有限脈衝響應(FIR：finite impulse response)濾波器，且特別是有關於一種平行有限脈衝響應濾波器與其相對應的濾波方法。
有限脈衝響應濾波器通常用在無線通訊系統的傳送端，用來調整待傳送的訊號的頻譜(spectrum)，使訊號符合規格所需要的頻譜遮罩(spectrum mask)。
近年來由於通訊技術的發展，從無線區域網路(WLAN：wireless local area network)到第四代(4G)技術，以至於目前的第五代(5G)技術，通訊技術越來越複雜多元。通訊系統的功耗、傳輸速度、以及硬體面積等議題也更受重視。
本揭露提供一種有限脈衝響應濾波器與其相對應的濾波方法，可減少通訊系統的功耗與硬體面積，並提高通訊系統的吞 吐率(throughput)。
本揭露的有限脈衝響應濾波器接收輸入序列。輸入序列 包括多個輸入值。此有限脈衝響應濾波器包括至少一第一加法器、至少一乘法器、以及第二加法器。每一上述第一加法器以平行方式同時進行多個加法運算。每一上述加法運算輸出上述多個輸入值其中兩個之和數。乘法器耦接第一加法器。每一上述乘法器以平行方式同時進行多個乘法運算。每一上述乘法運算輸出上述多個和數其中之一與有限脈衝響應濾波器的多個係數其中之一的乘積。第二加法器耦接乘法器，輸出上述多個乘積之總和。
本揭露的濾波方法包括以下步驟：接收輸入序列，其中輸入序列包括多個輸入值；在多個時脈週期的每一時脈週期，以平行方式同時進行多個加法運算，其中每一上述加法運算輸出上述多個輸入值其中兩個之和數；在每一上述時脈週期，以平行方式同時進行多個乘法運算，其中每一上述乘法運算輸出上述多個和數其中之一與多個係數其中之一的乘積；以及輸出上述多個乘積之總和。
上述的有限脈衝響應濾波器與濾波方法能藉由平行架構減少功耗並提高吞吐率，能藉由關閉一部分乘法運算以進一步降低功耗，也能藉由簡化乘法運算來減少硬體面積。
為讓本揭露的上述特徵能更明顯易懂，下文特舉實施例，並配合所附圖式作詳細說明如下。
100、500、600、700‧‧‧有限脈衝響應濾波器
110、510、610、710‧‧‧延遲鏈
111~113、131~137‧‧‧延遲器
121~127、150~160、204、302、521~527、621~627、721~727‧‧‧加法器
141~147‧‧‧乘法器
201~203、301‧‧‧移位器
170‧‧‧輸入端
s n ‧‧‧輸入值的和數
X n ~X n+13‧‧‧輸入值的批次
圖1是依照本揭露的一實施例的一種有限脈衝響應濾波器的示意圖。
圖2是依照本揭露的一實施例的乘法器的部分示意圖。
圖3是依照本揭露的另一實施例的乘法器的部分示意圖。
圖4是802.11p通訊標準的傳輸頻譜遮罩的示意圖。
圖5至圖7是依照本揭露的一實施例的多個有限脈衝響應濾波器的示意圖。
有限脈衝響應濾波器(以下簡稱FIR濾波器)可用下面的公式(1)表示。
公式(1)其中的x(n)是FIR濾波器的輸入值，y(n)是FIR濾波器的輸出值，n的數值範圍是0至無限大，h(i)是FIR濾波器的係數，N是h(i)的數量。FIR濾波器的輸出值y(n)就是輸入值x(n)至x(n-(N-1))和係數h(0)至h(N-1)的卷積(convolution)。
FIR濾波器的係數h(i)是首尾對稱的，也就是說，係數h(i)符合下面的公式(2)。根據公式(2)可將公式(1)簡化而得到下面的公式(3)。
h(i)=h(N-1-i)……………………………………………………(2)
上面的公式(3)假設N是奇數。如果N是偶數，則公式(3)應改為下面的公式(4)。
圖1是依照本揭露的一實施例的一種FIR濾波器100的示意圖。FIR濾波器100是根據公式(3)設計的實體數位電路。FIR濾波器100採用四路平行架構，有51個係數，也就是說N等於51。FIR濾波器100自輸入端170接收一個輸入序列，此輸入序列包括多個輸入值。FIR濾波器100包括延遲鏈110、加法器121~127、延遲器131~137、乘法器141~147、以及加法器150。加法器121~127耦接延遲鏈110。延遲器131~137分別耦接加法器121~127。乘法器141~147分別耦接延遲器131~137。加法器150耦接乘法器141~147。
延遲鏈110接收輸入序列，並依照輸入值在輸入序列的順序，將輸入值組成多個批次(batch)X n 至X n +13，其中每一個批次包括4個輸入值。舉例而言，批次X n+1的四個輸入值分別表示為X n+1,1、X n+1,2、X n+1,3和X n+1,4，依此類推。延遲鏈110可包括串列耦接的至少一個延遲器，例如延遲器111~113。其中的第一個延遲器直接自輸入序列逐一接收批次X n 至X n +13。其餘的每一個延遲器自前一個延遲器逐一接收批次X n 至X n +13。每一個延遲器將所接收的批次延遲一段預設時間，然後輸出延遲後的批次。上述 的預設時間可以是一個時脈訊號的一個週期。FIR濾波器100的每一個延遲器都可以接收此時脈訊號，做為延遲的依據。
每一個加法器121~127以平行方式同時進行多個加法運算，其中每一個加法運算輸出輸入序列其中兩個輸入值之和數。每一個加法器121~127都自延遲鏈110的延遲器所輸出的批次直接取得輸入值。每一個乘法器141~147以平行方式同時進行多個乘法運算，其中每一個乘法運算輸出上述多個和數其中之一與FIR濾波器100的多個係數其中之一的乘積。下面的表1列出加法器121~127進行的加法運算以及乘法器141~147進行的乘法運算。
表1的加法運算其中的運算元(operand)X n+1,3至X n+13,1相當於公式(3)的輸入值x(n)至x(n-(N-1))。s n,1至s n,26分別是同一列的加法運算所產生的和數。例如s n,24=X n+8,4+X n+7,4，依此類推。h 1至h 26相當於公式(3)的係數h(i)。
由表1可以看出，每一個加法器121~127最多可同時進行四個加法運算。每一個乘法器141~147最多可同時進行四個乘法運算。輸入值X n+7,2是整個輸入序列的中點。對於每一個加法運算，產生其和數的兩個輸入值的輸入順序分別位在輸入序列的中 點之前和中點之後，而且上述兩個輸入值在輸入序列的位置對於中點X n+7,2是對稱的。如果FIR濾波器100的係數數量為偶數，則輸入序列的中點是在最中間的兩個輸入值之間。從公式(3)、(4)就可以看出上述的輸入值的位置對稱關係。另外從表1也可以看出，每一個加法器121~127分別使用輸入值其中的兩組連續輸入值以同時進行多個加法運算。例如加法器125使用的第一組連續輸入值是X n+9,1至X n+9,4，第二組連續輸入值是X n+5,3至X n+6,4。
延遲器131~137的作用是讓每一個加法器121~127和對應的乘法器141~147在不同的時脈週期運作。例如加法器123在某一個時脈週期以平行方式同時計算出四個和數s n,9至s n,12，然後經過延遲器133的延遲，乘法器143在下一個時脈週期取得這四個和數s n,9至s n,12，以同時進行四個乘法運算。
加法器150包括多個加法器151~160以及多個延遲器。加法器151將乘法器141以平行方式計算產生的四個乘積加總之後輸出。加法器152將乘法器142以平行方式計算產生的四個乘積加總之後輸出，依此類推。加法器158將加法器151~154的輸出數值加總之後輸出。加法器159將加法器155~157的輸出數值加總之後輸出。加法器160將加法器158和159的輸出數值相加之後輸出。因此，加法器150的最終輸出就是乘法器141~147計算產生的所有乘積的總和，相當於公式(3)的y(n)。加法器150其中的多個延遲器是為了在連續的兩級加法器之間加入一個時脈週期的緩衝。
圖1的加法器150只是一個範例。在另一實施例中，可以改變加法器150的架構，只要能將乘法器141~147計算產生的所有乘積加總輸出即可。
由表1可知，加法器121~126各自進行四個加法運算，加法器127進行兩個加法運算。乘法器141~146各自進行四個乘法運算，乘法器147進行兩個乘法運算。因此，和一般的非平行FIR濾波器相比，FIR濾波器100可達到近乎四倍的吞吐率。如果是在同樣吞吐率的需求下，可以降低操作頻率，以降低功率需求。例如FIR濾波器100採用四路平行架構，只需要四分之一的操作頻率，可因此大幅降低功耗。
圖1的FIR濾波器100的係數數量為奇數。在此技術領域具有通常知識者應當知曉，只要對FIR濾波器100做小幅修改，就能將係數數量改為偶數。
圖1的FIR濾波器100採用四路平行架構。在另一實施例中，FIR濾波器100可以是任意的L路平行架構，其中L為大於一的預設整數。對於採用L路平行架構的FIR濾波器100而言，延遲鏈110提供的每一個輸入批次包括L個輸入值，每一個加法器121~127最多可同時進行L個加法運算，而且每一個乘法器141~147最多可同時進行L個乘法運算。
圖1的FIR濾波器100具有51個係數，相當於公式(3)其中的N等於51的情況。在另一實施例中，FIR濾波器100的係數數量N可以改變。在這樣的實施例中，延遲鏈110其中的延遲 器的數量、加法器121~127的數量、延遲器131~137的數量、乘法器141~147的數量、以及加法器150的架構都可以調整，以配合不同的N值。一般的規則是，延遲鏈110的延遲器的數量、加法器121~127的數量、延遲器131~137的數量、以及乘法器141~147的數量都和N成正比。因此FIR濾波器100的架構可適應任意的N值。
FIR濾波器100的平行化架構會增加硬體的數量與面積。為了降低硬體面積，可以簡化係數h 1~h 26。舉例而言，假設每一個係數h 1~h 26都是十位元的預設常數，則乘法器141~147進行的每一個乘法運算需要λ次的移位(shift)與加法運算，λ為對應的係數的非零位元的數量，λ的上限為10。如果能將每一個係數h 1~h 26都簡化為只有兩個或三個非零位元，就能大幅簡化乘法運算與其對應的硬體面積。
如上所述，表1的係數h 1~h 26相當於公式(3)、(4)的係數h(i)。以下就用h(i)表示係數h 1~h 26。在一個實施例中，可用公式(5)為每一個原始係數h(i)計算一個相對應的簡化係數。
公式(5)之中，λ i 等於2或3，c k,i 等於-1、0或1，g k,i 為大於或等於0而且小於原始係數h(i)的位元數的整數。上述的c k,i 和g k,i 為透過塔布搜尋法(tap search)在時域(time domain)和頻域(frequency domain)搜尋所得的最佳參數。本實施例的塔布搜尋法 在後面有進一步說明。在表1的每一個乘法運算中，用簡化係數取代原始係數h(i)就能簡化乘法運算與其相對應的硬體電路。
圖2是依照本揭露的一實施例的乘法器141~147的部分示意圖。以乘法器141為例。對於乘法器141執行的每一個乘法運算，乘法器141可包括一個如圖2所示的乘法電路以執行此乘法運算。其餘的每一個乘法器142~147也是如此。此乘法電路對應的λ i 等於3。因此圖2的乘法電路包括三個移位器201~203以及一個加法器204。加法器204耦接移位器201~203。
移位器201~203接收對應的乘法運算的和數s n 。移位器201~203分別對應於對應此乘法電路的係數h(i)的參數g 1,i 、g 2,i 和g 3,i 。移位器201將和數s n 移位g 1,i 次之後輸出，相當於將和數s n 乘上。移位器202將和數s n 移位g 2,i 次之後輸出，相當於將和數s n 乘上。移位器203將和數s n 移位g 3,i 次之後輸出，相當於將和數s n 乘上。加法器204根據對應的係數h(i)的參數c k,i 將移位器201~203的輸出相加或相減，以產生公式(5)的簡化係數。
如果，在另一實施例中，圖2的乘法電路對應的λ i 等於2，則可以省略移位器203。
圖3是依照本揭露的另一實施例的乘法器141~147的部分示意圖。以乘法器141為例。對於乘法器141執行的每一個乘法運算，乘法器141可包括一個如圖3所示的乘法電路。其餘的每一個乘法器142~147也是如此。圖3的乘法電路包括移位器301 以及加法器302。加法器302耦接移位器301。
移位器301接收此乘法電路對應的乘法運算的和數s n 。 移位器301在一個時脈訊號的第k個週期將和數s n 移位g k,i 次之後輸出，相當於將和數s n 乘上。加法器302可累加移位器301的k次輸出，以產生公式(5)的簡化係數。
以下說明在本揭露的一個實施例中如何透過塔布搜尋法來搜尋最佳的參數c k,i 和g k,i 。首先，依照公式(6)在時域為每一個原始係數h(i)搜尋相對應的參數c k,i 和g k,i 。
公式(6)之中，N為本實施例的FIR濾波器100的係數的數量，本實施例的N為奇數。G是一個有多個可能數值的參數。例如可以定義一個等差數列，而G可以是該等差數列中的任何一個數值。例如可將0.5到1的範圍分為500等分，每個等分的長度為(1-0.5)/500=0.001。上述的等差數列可以是這500個等分的501個端點，0.5和1是這501個端點其中兩個，而G可以是這501個端點其中任何一個。Q( )是一個量化函數(quantization function)，可將一個實數映射至一個定義域(domain)D之中的最接近該實數的一個元素(element)。下面的公式(7)是定義域D的定義。
公式(7)之中的β表示定義域D的元素，R是所有實數的集合。公式(7)的β和公式(5)的簡化係數有相似的定義，因此公 式(7)的λ也就是公式(5)的λ i ，λ等於2或3，c k 等於-1、0或1，g k 為大於或等於0而且小於原始係數h(i)的位元數的整數。定義域D就是可用形式表示的所有實數組成的集合。
公式(6)相當於在計算所有原始係數h(i)和所有簡化係數之間的誤差值e q (G)。公式(6)其中的簡化係數的定義如同公式(5)。對於每一個簡化係數，其對應的參數c k,i 和g k,i 各有多個可能的數值。參數G也有多個可能的數值。公式(6)其中的(N-1)/2+1個參數c k,i 、(N-1)/2+1個參數g k,i 、以及一個參數G的多個可能數值的每一種組合可計算出一個對應的誤差值e q (G)。將全部的組合所得的誤差值e q (G)排序，可取得這些誤差值其中的最小誤差值e q(min)，然後可選出這些誤差值其中小於M*e q(min)的多個誤差值e q 。這些被選出的誤差值e q 也包括最小誤差值e q(min)。M是一個預設倍數，本實施例的M等於5。在另一個實施例中，M可以是其他的大於一的整數。
由公式(6)可看出每一個被選出的誤差值e q 對應多個參數c k,i 和多個參數g k,i 。用這些參數c k,i 和g k,i 計算產生的簡化係數取代原始係數h(i)，可計算出FIR濾波器100的一個頻率響應(frequency response)。因此每一個被選出的誤差值e q 對應一個頻率響應。接下來可進行頻域的參數搜尋，也就是將每一個被選出的誤差值e q 對應的頻率響應和FIR濾波器100的原始頻率響應相比較，如此可找出最近似原始頻率響應的頻率響應，也可以找出這 個最近似的頻率響應所對應的誤差值e q 。這個最近似的誤差值e q 所對應的多個參數c k,i 和多個參數g k,i 就是公式(5)所採用的最佳參數。
現在已經有多種方法可以比較兩個頻率響應是否相似，上述的頻域搜尋可以採用其中任何一種方法。例如可以在頻域計算每一個被選出的誤差值e q 對應的頻率響應在FIR濾波器100的通過頻帶(pass band)的平均值(mean)，並在頻域計算FIR濾波器100的原始頻率響應在同一頻帶的平均值。比較以上平均值就能決定哪一個被選出的誤差值e q 對應的頻率響應最近似原始頻率響應。
前述的公式(6)適用於係數數量N為奇數的情況。如果FIR濾波器100的係數數量N為偶數，則可用下面的公式(8)取代公式(6)。
FIR濾波器100可以視應用的需要而關閉一部分乘法運算，使被關閉的乘法運算的輸出為零。如此用同一個FIR濾波器就能滿足多種頻譜遮罩，也能減少不必要的功耗。
更詳細的說，可將FIR濾波器100的係數h(i)編號為0至N-1，也就是h(0)至h(N-1)。可將係數h(i)分為兩個集合S1和S2。集合S1包括係數h(i)其中的第j個係數至第N-1-j個係數，也就是h(j)至h(N-1-j)，j為正整數而且小於N/2。另一個集合S2則 包括其餘的係數h(i)。FIR濾波器100可以關閉集合S2之中的係數所對應的乘法運算，使被關閉的乘法運算的輸出為零。如圖2與圖3的實施例所述，每一個乘法運算各有一個對應的乘法電路。關閉乘法運算，也就是關閉對應的乘法電路。
舉例來說，電機電子工程師學會(IEEE：Institute of Electrical and Electronics Engineers)的802.11p通訊標準的專用短距離通訊(DSRC：Dedicated Short Range Communications)系統的每一個裝置等級(device class)都有一個對應的傳輸頻譜遮罩。圖4繪示操作在5.9 DSRC頻譜的裝置等級A~D的傳輸頻譜遮罩，其中每一張遮罩圖的縱軸為功率衰減(power attenuation)，橫軸為偏移頻率(offset frequency)。
以本揭露的某一個實施例的FIR濾波器100為例，假設係數h(i)的數量N為71。如圖4所示，裝置等級A和裝置等級B需要將操作頻帶以外的功率抑制到大約-20dBr。此時集合S1只需要包括h(i)中間的23個係數，也就是h(24)至h(46)。考慮到FIR濾波器100的係數h(i)的首尾對稱性，只需要用到對應h(24)至h(46)的12個乘法運算。其餘的乘法運算所對應的乘法電路可以關閉。
同理，裝置等級C需要將操作頻帶以外的功率抑制到大約-30dBr。此時集合S1需要包括h(i)中間的39個係數，也就是h(16)至h(54)，只需要用到對應h(16)至h(54)的20個乘法運算。其餘的乘法運算所對應的乘法電路可以關閉。
裝置等級D需要將操作頻帶以外的功率抑制到大約-45 dBr。此時需要用到所有係數，需要全部的36個乘法運算。每一個乘法電路都必須啟動。
裝置等級A和B使用的乘法電路數量只有裝置等級D的三分之一。也就是說在裝置等級A及B時可以關掉FIR濾波器100的三分之二的乘法電路，避免不必要的功耗。為了適用於上述的關閉一部分乘法電路的操作，FIR濾波器100可以是使用窗口演算法(windowing algorithm)設計所得的濾波器。
在本揭露的另一個實施例中，可以用多個類似FIR濾波器100的平行FIR濾波器的組合達到更高的平行度。例如可將圖1的FIR濾波器100、圖5的FIR濾波器500、圖6的FIR濾波器600、以及圖7的FIR濾波器700這四個FIR濾波器組成一個具有更高的平行度的FIR濾波器。
圖5是本實施例的FIR濾波器500的示意圖。圖5僅繪示FIR濾波器500的延遲鏈510和加法器521~527。FIR濾波器500的其餘部分和圖1的FIR濾波器100相同。圖6是本實施例的FIR濾波器600的示意圖。圖6僅繪示FIR濾波器600的延遲鏈610和加法器621~627。FIR濾波器600的其餘部分和圖1的FIR濾波器100相同。圖7是本實施例的FIR濾波器700的示意圖。圖7僅繪示FIR濾波器700的延遲鏈710和加法器721~727。FIR濾波器700的其餘部分和圖1的FIR濾波器100相同。加法器521~527、621~627、以及721~727皆可同時進行多個加法運算。下面的表2列示加法器521~527、621~627、以及721~727所進行 的加法運算。
從表1可以看出FIR濾波器100計算的是輸入值X n+1,3至X n+13,1和係數h 1至h 51的卷積。因為h 1至h 51的首尾對稱，實際上FIR濾波器100只用到係數h 1至h 26。從表2可以看出，FIR濾波器500計算輸入值X n+1,4至X n+13,2和係數h 1至h 51的卷積，FIR濾波器600計算輸入值X n,1至X n+13,3和係數h 1至h 51的卷積，FIR濾波器700則計算輸入值X n,2至X n+13,4和係數h 1至h 51的卷積。如此就是四個FIR濾波器同時計算四個不同的卷積。FIR濾波器100、500、600和700的組合可以在每個時脈週期以平行方式同時進行十六個加法運算與十六個乘法運算，可將吞吐率提高為一般非平行FIR濾波器的十六倍。
為了更清楚地說明上述的平行FIR濾波器，可將輸入序列中的輸入值連續編號。例如批次X 1包括輸入值x(1)~x(4)，批次X 2包括輸入值x(5)~x(8)，依此類推。下面的表3列示僅使用FIR濾波器100的情況下，在四個時脈週期中的每一時脈週期所計算 的卷積，以及這些卷積和公式(3)的輸出值y(n)的關係。後面的時脈週期可依此類推。下面的表4則列示使用FIR濾波器100、500、600和700組合成平行FIR濾波器的情況下，在四個時脈週期中的每一時脈週期所計算的卷積，以及這些卷積和公式(3)的輸出值y(n)的關係。後面的時脈週期可依此類推。
從表3可以看出，若僅使用FIR濾波器100，則可在每個時脈週期接收一個輸入值x(n)並計算一個輸出值y(n)。從表4可以看出，若使用包括FIR濾波器100、500、600和700的組合式平行FIR濾波器，則FIR濾波器100、500、600或700皆可在每個時脈週期各自接收一個輸入值x(n)並計算一個輸出值y(n)。於是整個組合式平行FIR濾波器可在每個時脈週期接收四個輸入值x(n)並計算四個輸出值y(n)。在另一實施例中，可沿用上述規則組合任意數量的FIR濾波器以達成更低或更高的平行度。
本揭露的一實施例提供一種濾波方法。圖1的FIR濾波器100也可視為此濾波方法的流程示意圖。首先，延遲鏈110自輸入端170接收包括多個輸入值的輸入序列。然後在多個時脈週期的每一個時脈週期以平行方式同時進行多個加法運算。例如加法器127在某一個時脈週期以平行方式同時進行兩個加法運算， 加法器126在下一個時脈週期以平行方式同時進行另外四個加法運算，加法器125在更下一個時脈週期以平行方式同時進行另外四個加法運算，依此類推。在每一個時脈週期，只有加法器121~127其中之一在進行加法運算。然後在多個時脈週期的每一個時脈週期以平行方式同時進行多個乘法運算。例如乘法器147在某一個時脈週期以平行方式同時進行兩個乘法運算，乘法器146在下一個時脈週期以平行方式同時進行另外四個乘法運算，乘法器145在更下一個時脈週期以平行方式同時進行另外四個乘法運算，依此類推。在每一個時脈週期，只有乘法器141~147其中之一在進行乘法運算。最後加法器150輸出乘法器141~147輸出的所有乘積之總和。至於此濾波方法的技術細節，在前面的實施例中已經說明過，就不再贅述。在另一個實施例中，本揭露的濾波方法可以如同圖5至圖7的實施例，以同時計算多個輸出值y(n)的方式提高吞吐率。
總之，上述的FIR濾波器可降低通訊系統的傳送端的操作頻率，進而降低功率消耗。上述的FIR濾波器可用加法器跟移位器取代乘法器，可大幅節省乘法電路的硬體面積。上述的FIR濾波器可動態關閉一部分乘法電路以降低功耗，而且只用一個FIR濾波器即可滿足多重頻譜遮罩的需求。
雖然本揭露已以實施例揭露如上，然其並非用以限定本揭露，任何所屬技術領域中具有通常知識者，在不脫離本揭露的精神和範圍內，當可作些許的更動與潤飾，故本揭露的保護範圍 當視後附的申請專利範圍所界定者為準。
100‧‧‧有限脈衝響應濾波器
110‧‧‧延遲鏈
111~113、131~137‧‧‧延遲器
121~127、150~160‧‧‧加法器
141~147‧‧‧乘法器
170‧‧‧輸入端
X n ~X n+13‧‧‧輸入值的批次

Claims (19)

  1. 一種有限脈衝響應濾波器，接收一輸入序列，該輸入序列包括多個輸入值，該有限脈衝響應濾波器包括：至少一第一加法器，每一上述第一加法器以平行方式同時進行多個加法運算，每一上述加法運算輸出上述多個輸入值其中兩個之和數；至少一乘法器，耦接上述至少一第一加法器，每一上述乘法器以平行方式同時進行多個乘法運算，每一上述乘法運算輸出上述多個和數其中之一與該有限脈衝響應濾波器的多個係數其中之一的乘積；以及一第二加法器，耦接上述至少一乘法器，輸出上述多個乘積之總和。
  2. 如申請專利範圍第1項所述的有限脈衝響應濾波器，更包括：一延遲鏈，耦接上述至少一第一加法器，接收該輸入序列，並依照上述多個輸入值在該輸入序列的順序將上述多個輸入值組成多個批次，其中每一上述批次包括L個輸入值，L為大於一的整數，上述至少一第一加法器自上述多個批次取得上述多個輸入值。
  3. 如申請專利範圍第2項所述的有限脈衝響應濾波器，其中L為每一上述第一加法器同時進行的上述加法運算的數量上限，也是每一上述乘法器同時進行的上述乘法運算的數量上限。
  4. 如申請專利範圍第2項所述的有限脈衝響應濾波器，其中該延遲鏈包括：串列耦接的至少一延遲器，其中的第一個延遲器直接自該輸入序列逐一接收上述多個批次，其餘的每一上述延遲器自前一個延遲器逐一接收上述多個批次，每一上述延遲器將所接收的該批次延遲一段預設時間然後輸出該批次。
  5. 如申請專利範圍第1項所述的有限脈衝響應濾波器，其中對於每一上述加法運算，產生該和數的上述兩個輸入值的輸入順序分別在該輸入序列的中點之前和該中點之後，而且上述兩個輸入值在該輸入序列的位置對於該中點是對稱的。
  6. 如申請專利範圍第1項所述的有限脈衝響應濾波器，其中每一上述第一加法器使用上述多個輸入值其中的兩組連續輸入值進行該第一加法器的上述多個加法運算。
  7. 如申請專利範圍第1項所述的有限脈衝響應濾波器，其中上述多個係數的數量為N而且被編號為0至N-1，N為正整數，上述多個係數被分為一第一集合和一第二集合，該第一集合包括上述多個係數中的第j個係數至第N-1-j個係數，該第二集合包括其餘的係數，j為正整數而且小於N/2，該有限脈衝響應濾波器關閉該第二集合中的係數所對應的上述乘法運算，使被關閉的上述乘法運算的輸出為零。
  8. 如申請專利範圍第1項所述的有限脈衝響應濾波器，其中 每一上述乘法運算中的該係數被簡化為，λ等於2或3，c k 等於-1、0或1，g k 為大於或等於0而且小於該係數的位元數的整數。
  9. 如申請專利範圍第8項所述的有限脈衝響應濾波器，其中上述的c k 和g k 為透過塔布搜尋法在時域和頻域搜尋所得。
  10. 如申請專利範圍第8項所述的有限脈衝響應濾波器，其中對於每一上述乘法運算，每一上述乘法器包括：多個移位器，每一上述移位器對應上述g k 其中之一，將該乘法運算對應的該和數移位g k 次之後輸出，相當於將該和數乘上；以及一第三加法器，耦接上述多個移位器，根據上述c k 將上述多個移位器的輸出相加或相減，以產生簡化後的該係數。
  11. 如申請專利範圍第8項所述的有限脈衝響應濾波器，其中對於每一上述乘法運算，每一上述乘法器包括：一移位器，在一時脈訊號的第k個週期將該和數移位g k 次之後輸出，相當於將該和數乘上；以及一第三加法器，耦接該移位器，累加該移位器的k次輸出，以產生簡化後的該係數。
  12. 一種濾波方法，包括：接收一輸入序列，其中該輸入序列包括多個輸入值；在多個時脈週期的每一時脈週期，以平行方式同時進行多個 加法運算，其中每一上述加法運算輸出上述多個輸入值其中兩個之和數；在每一上述時脈週期，以平行方式同時進行多個乘法運算，其中每一上述乘法運算輸出上述多個和數其中之一與多個係數其中之一的乘積；以及輸出上述多個乘積之總和。
  13. 如申請專利範圍第12項所述的濾波方法，更包括：依照上述多個輸入值在該輸入序列的順序，將上述多個輸入值組成多個批次，其中每一上述批次包括L個輸入值，L為大於一的整數，上述多個加法運算自上述多個批次取得上述多個輸入值。
  14. 如申請專利範圍第13項所述的濾波方法，其中L為每一上述時脈週期之中同時進行的上述加法運算的數量上限，也是每一上述時脈週期之中同時進行的上述乘法運算的數量上限。
  15. 如申請專利範圍第12項所述的濾波方法，其中對於每一上述加法運算，產生該和數的上述兩個輸入值的輸入順序分別在該輸入序列的中點之前和該中點之後，而且上述兩個輸入值在該輸入序列的位置對於該中點是對稱的。
  16. 如申請專利範圍第12項所述的濾波方法，更包括：在每一上述時脈週期之中，使用上述多個輸入值其中的兩組連續輸入值進行上述多個加法運算。
  17. 如申請專利範圍第12項所述的濾波方法，其中上述多個 係數的數量為N而且被編號為0至N-1，N為正整數，上述多個係數被分為一第一集合和一第二集合，該第一集合包括上述多個係數中的第j個係數至第N-1-j個係數，該第二集合包括其餘的係數，j為正整數而且小於N/2，而且該濾波方法更包括：關閉該第二集合中的係數所對應的上述乘法運算，使被關閉的上述乘法運算的輸出為零。
  18. 如申請專利範圍第12項所述的濾波方法，其中每一上述乘法運算中的該係數被簡化為，λ等於2或3，c k 等於-1、0或1，g k 為大於或等於0而且小於該係數的位元數的整數。
  19. 如申請專利範圍第18項所述的濾波方法，其中上述的c k 和g k 為透過塔布搜尋法在時域和頻域搜尋所得。
TW103137590A 2014-10-30 2014-10-30 有限脈衝響應濾波器與濾波方法 TWI566523B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
TW103137590A TWI566523B (zh) 2014-10-30 2014-10-30 有限脈衝響應濾波器與濾波方法
US14/598,242 US20160126933A1 (en) 2014-10-30 2015-01-16 Finite impulse response filter and filtering method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
TW103137590A TWI566523B (zh) 2014-10-30 2014-10-30 有限脈衝響應濾波器與濾波方法

Publications (2)

Publication Number Publication Date
TW201616810A true TW201616810A (zh) 2016-05-01
TWI566523B TWI566523B (zh) 2017-01-11

Family

ID=55853799

Family Applications (1)

Application Number Title Priority Date Filing Date
TW103137590A TWI566523B (zh) 2014-10-30 2014-10-30 有限脈衝響應濾波器與濾波方法

Country Status (2)

Country Link
US (1) US20160126933A1 (zh)
TW (1) TWI566523B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109525214A (zh) * 2018-10-11 2019-03-26 上海辛格林纳新时达电机有限公司 Fir滤波器的滤波方法、fir滤波器及伺服驱动器

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5258939A (en) * 1991-10-10 1993-11-02 Harris Corporation Fold and decimate filter architecture
US5479363A (en) * 1993-04-30 1995-12-26 The Regents Of The University Of California Programmable digital signal processor using switchable unit-delays for optimal hardware allocation
US6889238B2 (en) * 2001-06-25 2005-05-03 Lockheed Martin Corporation Parallel decimator adaptive filter and method for all-rate gigabit-per-second modems
US7353243B2 (en) * 2002-10-22 2008-04-01 Nvidia Corporation Reconfigurable filter node for an adaptive computing machine
US8725785B1 (en) * 2004-08-09 2014-05-13 L-3 Communications Corp. Parallel infinite impulse response filter
US20060218213A1 (en) * 2005-03-24 2006-09-28 Shereef Shehata Optimization technique for FIR and IIR filter design

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109525214A (zh) * 2018-10-11 2019-03-26 上海辛格林纳新时达电机有限公司 Fir滤波器的滤波方法、fir滤波器及伺服驱动器
CN109525214B (zh) * 2018-10-11 2022-10-14 上海辛格林纳新时达电机有限公司 Fir滤波器的滤波方法、fir滤波器及伺服驱动器

Also Published As

Publication number Publication date
TWI566523B (zh) 2017-01-11
US20160126933A1 (en) 2016-05-05

Similar Documents

Publication Publication Date Title
Fernandez-Vazquez et al. Maximally flat CIC compensation filter: Design and multiplierless implementation
Kumar et al. Generalized rational sampling rate conversion polyphase FIR filter
Goel et al. Design of FIR filter using FCSD representation
TWI566523B (zh) 有限脈衝響應濾波器與濾波方法
Vinod et al. An efficient coefficient-partitioning algorithm for realizing low-complexity digital filters
Eshtawie et al. An algorithm proposed for FIR Filter coefficients representation
Kumar et al. FPGA Implementation of Systolic FIR Filter Using Single-Channel Method
JPH05259813A (ja) ディジタルフィルタ
KR100924934B1 (ko) 폴라 모듈레이터를 위한 멀티플라이어가 없는 인터폴레이션기술
Reddy et al. Shift add approach based implementation of RNS-FIR filter using modified product encoder
Dai et al. FPGA Realization of Hardware-Flexible Parallel Structure FIR Filters Using Combined Systolic Arrays
Macpherson et al. Low FPGA area multiplier blocks for full parallel FIR filters
Meyer-Baese et al. Infinite impulse response (IIR) digital filters
Ramstad Digital two-rate IIR and hybrid IIR/FIR filters for sampling rate conversion
Chodoker et al. Multiple Constant Multiplication Technique for Configurable Finite Impulse Response Filter Design
Skulina et al. Computational Efficiency of Interpolated Band-Stop Filters for Even Spectral Bands
Umasankar et al. Design and analysis of various slice reduction algorithm for low power and area efficient FIR filter
Liu et al. A pipelined sign-error LMS adaptive filter architecture with low computational complexity
Chulet et al. FIR filter designing using wallace multiplier
Mathur et al. High throughput arbitrary sample rate converter for software radios
Meyer-Base et al. Optimal algebraic integer implementation with application to complex frequency sampling filters
Sarband et al. Obtaining minimum depth sum of products from multiple constant multiplication
Mishra et al. Design of sixth order CIC compensation filter
Kotha et al. An RNS based reconfigurable FIR filter design using shift and add approach
Priyadharshini et al. Retracted: Implementation of Distributed Arithmetic based Sum-of-Products