|
累積梳狀（CIC）濾波器分析與設(shè)計 1、累積梳狀（CIC）濾波器的分析 所謂累積梳狀濾波器，是指該濾波器的沖激響應(yīng)具有如下形式：  （1） 式中 為梳狀濾波器的系數(shù)長度（后面將會看到這里的 也就是抽取因子）。根據(jù)Z變換的定義，濾波器的Z變換為：  （2） 式中，  （3）  （4） 其實(shí)現(xiàn)框圖如圖1所示： 可見，CIC濾波器是由兩部分組成：累積器 和梳狀濾波器 的級聯(lián)，這就是為什么稱之為累積梳狀濾波器的原因。下面分析一下梳狀濾波器的幅頻特性。 把 代入可得 的頻率響應(yīng)為：  （5） 其幅頻特性為：  （6） 若設(shè)N＝7，就可以得到如圖2所示的相應(yīng)的頻譜特性曲線： 由圖2可以清楚地看到： 的形狀猶如一把梳子，故把其形象地稱之為梳狀濾波器。同樣可以求得累積器 的頻率響應(yīng)為：  （7） 故CIC濾波器的總頻率響應(yīng)為：  （8） 式中， 為抽樣函數(shù)，且 ，所以CIC濾波器在 處的幅度值為 ，即：  （9） CIC濾波器的幅頻特性如圖3所示： 在 區(qū)間上稱 的區(qū)間為CIC濾波器的主瓣，而其它區(qū)間稱為旁瓣。由圖3可知，在 區(qū)間上隨著頻率的增大，旁瓣電平不斷減小，其中第一旁瓣電平為：  （10） 比如N＝7，則第一旁瓣與主瓣的電平差值為：16.9dB。不過，當(dāng) 的時候，有 ，所以第一旁瓣電平 為：  （11） 因此，旁瓣與主瓣的差值 （用dB數(shù)表示）為：  （12） 可見，單級CIC濾波器的旁瓣電平是比較大的，只比主瓣低13.46dB，這也就意味著阻帶衰減差，一般很難直接滿足實(shí)用要求。為了減低旁瓣電平，自然會想到的方法是采用多級CIC濾波器級聯(lián)的辦法。設(shè)用Q級CIC級聯(lián)，那么總的頻率響應(yīng)為：  （13） 同理可求得Q級CIC濾波器的旁瓣抑制為：  （14） 比如當(dāng) 時，主瓣與旁瓣的差值為：67.3dB，這樣的阻帶衰減基本上能滿足實(shí)際要求。但在實(shí)際的抽樣率變換系統(tǒng)中，CIC濾波器旁瓣區(qū)域往往作為不確定 帶來處理，也就是說在這些旁瓣區(qū)域不會有信號頻譜（鏡像或混疊頻譜），因此在CIC濾波器設(shè)計中所要考慮的重要指標(biāo)是抗混疊問題。這個問題放在具體的抽取率變換濾波器的分析中加以討論。 2、累積梳狀（CIC）抽樣率變換濾波器的分析 抽樣率變換濾波器包括抽取濾波器和內(nèi)插濾波器。根據(jù)多采樣率數(shù)字信號處理理論，內(nèi)插濾波器和抽取濾波器是對偶關(guān)系，也就是說抽取濾波器通過轉(zhuǎn)置可以得到內(nèi)插濾波器，反之亦然。因此，這里主要針對抽取濾波器展開分析，其結(jié)果也同樣適用于內(nèi)插濾波器。   根據(jù)多采樣率數(shù)字信號處理理論，抽取濾波器的結(jié)構(gòu)框圖如下：  現(xiàn)若假設(shè)用N級CIC濾波器來代替 ，每一級的濾波器系數(shù)長度為R，每一級的差分延遲為M，抽取數(shù)為R，那么可以得到如圖5所示的CIC抽取濾波器結(jié)構(gòu)圖：
據(jù)CIC濾波器的傳遞函數(shù)，圖5所示的抽取濾波器結(jié)構(gòu)圖和圖6所示的抽取濾波器結(jié)構(gòu)圖是等價的。  根據(jù)分支運(yùn)算的換位和電路的恒等關(guān)系（具體見附錄1），可以得到如圖7所示的CIC抽取濾波器結(jié)構(gòu)圖： 據(jù)圖6、圖7可以得到相對于高采樣率 的系統(tǒng)函數(shù)為：  （15） 從式（15）可以看到，N級CIC濾波器在功能上相當(dāng)于N級完全相同的FIR濾波器的級聯(lián)。如果按傳統(tǒng)的FIR濾波器方式實(shí)現(xiàn)，那么N級FIR濾波器的每一級都需要RM個存儲單元和一個累加器，但如果用CIC方式實(shí)現(xiàn)，那么N級CIC濾波器的每一級只需要M個存儲單元。 CIC濾波器幅頻特性具有帶帶的低通特性，這可以從圖3中看出來，下面具體來分析CIC抽取濾波器的幅頻特性。設(shè)：  （16） 其中 是相對于低采樣率 的歸一化頻率，把（16）式代入（15）式得到CIC濾波器的幅頻特性為：  （17） 如果抽取率R足夠大，那么有 ，從而有：   （18）  從上式可以看到，CIC的增益可達(dá)（RM）N, 差分延遲M影響著零點(diǎn)分布，也就是說差分延遲M可以影響幅頻特性，一般M＝1或2。若設(shè)N＝5、M＝1及R＝4，那么可以得到如圖8所示的幅頻特性曲線：
對于CIC抽取濾波器來說，那些零點(diǎn)附近的區(qū)域?qū)徽郫B到通帶而引起混疊誤差，具體而言這些混疊帶為：  （18） 其中 ，而 1，2，┅， ，這里 是指不大于x的最大整數(shù)。在實(shí)際設(shè)計時，混疊誤差是以所有混疊帶中的最大混疊誤差來衡量的。如設(shè)在圖8中的 為通帶截止頻率，那么最大的混疊誤差的位置為：  （19） 而在任意頻率點(diǎn) （ ），相對于最大值的衰減（以dB表示）為：  （20） 對于通帶截止頻率點(diǎn)的衰減，不一定需要知道 和 ，其實(shí)只要知道 即可，但對于阻帶截止頻率點(diǎn)的衰減，如果只知道 的條件是不夠的，應(yīng)該同時需要知道和的值。據(jù)上述分析，就可以得到如表1所示的不同通帶截止頻率點(diǎn)的衰減值和如表2所示的不同通帶截止頻率下的阻帶截止頻率點(diǎn)衰減值。 相對帶寬和差分延遲乘積（  ） | 在  處的通帶衰減（dB） 級數(shù)（N） | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | |
差分延遲 | 相對帶寬 | 在  處的阻帶衰減（dB） 級數(shù)（N） | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | |
3、CIC抽樣率變換濾波器的設(shè)計 上面已經(jīng)提到過，根據(jù)多采樣率數(shù)字信號處理理論，內(nèi)插濾波器和抽取濾波器是對偶關(guān)系。所以不管是內(nèi)插濾波器還是抽取濾波器，只要設(shè)計了一種濾波器，那么另一種濾波器也相當(dāng)于得到了設(shè)計。因此，這里只研究抽取濾波器的設(shè)計問題。 根據(jù)第2部分中對抽取濾波器的分析及自己的具體設(shè)計要求，我們可以確定所要選用的CIC濾波器類型。比如：我們要求把帶寬為30kHz、采樣率為6MHz的信號降低為采樣率為240kHz的信號，而要求通帶衰減最大不得超過3dB，阻帶衰減不得低于60dB。據(jù)此條件參考表1和2可知，用級數(shù)N＝4、差分延遲M＝1的CIC濾波器可以滿足要求。但對于CIC濾波器的設(shè)計來說，那是遠(yuǎn)遠(yuǎn)不夠的，這是因?yàn)镃IC是數(shù)字濾波器，所以必須考慮在濾波過程中為了不產(chǎn)生溢出而需要的最大位數(shù)和每級可以舍去的位數(shù)。下面就來分析這兩個問題。 在圖7中從第j級開始到最后一級的總的系統(tǒng)函數(shù)為：   （21）其中，   （22）這組公式的獲得參見附錄2。 現(xiàn)把濾波器濾波過程中可能增加的最大值定義為最大幅度輸入信號情況下最大的輸出幅度。若假設(shè)信號輸入位數(shù)為Bin(滿刻度輸入)，即輸入信號為2Bin，那么根據(jù)：  （23） 故增加的值為：  （24） 由（15）式和（21）式得到：  （25） 因正系數(shù)的多項(xiàng)式乘積的展開式依然是正系數(shù)的多項(xiàng)式，因此有 ，這樣（24）式變?yōu)椋?/font>  （26） 另外，由（25）式，令 得到：  （27） 故有：  （28） 從而在輸入信號位數(shù)為 的情況下，濾波輸出的最大可能位數(shù)為：  （29） 這里之所以需要減1，那是因?yàn)槟莻�符號數(shù)相乘結(jié)果中會有兩位相同的符號位，所以需要去掉一位。另外， 是指不小于x的最小整數(shù)。  不僅是濾波器輸出的最大可能的位數(shù)，也是每一級濾波器的最大可能位數(shù)。但在許多實(shí)際情況中，太大會使寄存器位數(shù)過寬，這是很不經(jīng)濟(jì)的。為了能解決這個問題，在每一級濾波器中引入四舍五入或截取的方法來減少位數(shù)。這自然需要研究一下因四舍五入或截取所引起的濾波器輸出端誤差，具體的方法就是通過單獨(dú)研究每一級誤差源（指每一級濾波器截取或四舍五入引起的誤差）的均值和方差，利用統(tǒng)計特性研究濾波器輸出端總的均值和方差。根據(jù)實(shí)際實(shí)現(xiàn)時的濾波器結(jié)構(gòu)（下面介紹），可以認(rèn)為誤差源有 個： 個濾波器級因四舍五入或截取所產(chǎn)生的誤差源和一個輸出寄存器因四舍五入或截取所產(chǎn)生的誤差源。 一般我們認(rèn)為截取比四舍五入的效果要好，但在CIC濾波器中，除了第一級和第2N＋1級誤差源以外，其它級誤差源的輸出誤差統(tǒng)計結(jié)果不管是四舍五入還是截取都是相同的（見下面的分析）。如果在實(shí)際設(shè)計中，第一級采用全精度運(yùn)算（也就是以 位進(jìn)行運(yùn)算，不截取也不舍入），那么對于設(shè)計者來講只需要在最后一級（輸出寄存器誤差源）考慮采用截取還是四舍五入。 一般認(rèn)為每一級誤差源都是與輸入信號和其它誤差源不相關(guān)的白噪聲，并且服從均勻概率分布。對于第j級誤差源，我們可以得到誤差函數(shù)為：  （30） 這里 是第j級被拋掉的位數(shù)，因此可以得到第j級誤差源的均值 和方差 表達(dá)式：  （31）  （32）  據(jù)式（21），對于第k個系數(shù)，對應(yīng)的誤差均值和方差為： 和 ，其隨機(jī)過程模型如下圖所示： 因此，第j級誤差源在輸出端所產(chǎn)生的總的均值 為：  （33） 其中，  （34） 同樣，第j級誤差源在輸出端所產(chǎn)生的總的方差 為：  （35） 其中，  （36） 另外，可以獲證（我還未證明之）：  （37） 從這里可以看出，對于選用截取還是舍入除了對第一級和最后一級有不同的影響外，對其它級都不會有影響。這樣可以得到因截取或舍入所引起的在輸出端所產(chǎn)生的總均值 和方差 為：  （38）  （39） 現(xiàn)在，我們可以根據(jù)上述分析結(jié)果，在給定合適的誤差限制條件下反過來確定每一級濾波器需要拋棄的位數(shù)。由于方差在所有的誤差源中都受截取或舍入的影響，而均值只有對第一和最后一級產(chǎn)生影響，所以這里選用方差作為唯一的設(shè)計參數(shù)。 如果現(xiàn)在假定在輸出寄存器中所保留的位數(shù)為 ，那么最后需要被拋棄的最低位數(shù)為：  （40） 這里只所以需要加1，那是因?yàn)楸粧仐壍臄?shù)是符號數(shù)，所以除了拋棄的數(shù)據(jù)位，還需要指定拋棄數(shù)據(jù)的符號。 為了得到前2N級濾波器中的每一級所要拋棄的位數(shù) ，現(xiàn)提出以下的設(shè)計準(zhǔn)則：前2N個誤差源在輸出端所產(chǎn)生的方差小于或等于最后一級誤差源的方差，同時假定誤差等同地分布在這些源中。這樣就可以得到如下的不等式：  （41） 而由式（30）、（32）和（35），可以得到：  （42） 根據(jù)式（41）和（42）得到：  （43） 以2為底求對數(shù)得到：  （44） 也就是：   （45） 具體的計算步驟如下： - 根據(jù)已知的R、M和N的值及式（22）求
 ； - 據(jù)式（36）求
 ； - 據(jù)式（40）、（30）、（32）和（35）式求
 ； - 據(jù)（45）式求
 。
到此為止，我們已經(jīng)講述了CIC濾波器的選取依據(jù)和每一級位數(shù)的設(shè)計過程，下面通過一設(shè)計實(shí)例來詳細(xì)介紹設(shè)計過程和實(shí)現(xiàn)方法。 4、設(shè)計實(shí)例及實(shí)現(xiàn)方法 這里講述的例子還是要求把帶寬為30kHz、采樣率為6MHz的信號降低為采樣率為240kHz的信號，而要求通帶衰減最大不得超過3dB，阻帶衰減不得低于60dB，同時要求輸入和輸出位數(shù) 。因此，我們可以知道抽取因子R＝25，而相對于低采樣率的通帶截止頻率 。根據(jù)表1和表2可以知道，N＝4級和差分延遲M＝1的CIC濾波器已經(jīng)能夠滿足要求。為了簡單起見，這里每一級濾波器都采用截取方式，這樣我們據(jù)（29）式得到 ，據(jù)（40）式得到 。然后根據(jù)上面的分析結(jié)果，可以從理論上推得2N濾波器級所要被拋棄的位數(shù)分別為：1，6，9，13，14，15，16和17。如果假設(shè)小數(shù)點(diǎn)在輸出寄存器最低位的右邊，這樣就可以得到均值 以及標(biāo)準(zhǔn)偏差 。 從設(shè)計的角度來看，可以說工作已經(jīng)做完了，但就具體的實(shí)現(xiàn)而言，還需要對上述結(jié)果進(jìn)行修正。比如在硬件上我們現(xiàn)在用4比特的基本單元來實(shí)現(xiàn)抽取濾波器，其具體結(jié)構(gòu)如圖10所示。這樣為了盡量減小截尾誤差但又要防止溢出，我們假定除了第一級累積級外一律使保留位數(shù)盡量靠近4的倍數(shù)，但又不達(dá)到4的倍數(shù)，具體過程如圖11所示。從而得到2N濾波器級每級被拋棄的位數(shù)分別為：0，3，7，11，11，15，15和15。此時誤差均值變?yōu)?img id="aimg_t1jY1" onclick="zoom(this, this.src, 0, 0, 0)" class="zoom" width="104" height="23" src="http://c.51hei.com/a/huq/a/a/7/36/36.143.jpg" border="0" alt="" />（因?yàn)?img id="aimg_d3ann" onclick="zoom(this, this.src, 0, 0, 0)" class="zoom" width="45" height="21" src="http://c.51hei.com/a/huq/a/a/7/36/36.144.jpg" border="0" alt="" />）以及標(biāo)準(zhǔn)偏差變?yōu)?img id="aimg_EEQ1y" onclick="zoom(this, this.src, 0, 0, 0)" class="zoom" width="103" height="23" src="http://c.51hei.com/a/huq/a/a/7/36/36.145.jpg" border="0" alt="" />。 從而可以得到如圖12所示的最終實(shí)現(xiàn)結(jié)構(gòu)： 這里累積單元中的最低級因?yàn)闆]有進(jìn)位，而且4bit累積單元內(nèi)沒有反向器，故累積單元的最后一級的進(jìn)位端接低電平（即0）；而對于梳狀單元則不同，那時因?yàn)?bit梳狀單元的內(nèi)部包含有方向器，我們知道對于補(bǔ)碼形式的數(shù)據(jù)其所對應(yīng)的反向數(shù)據(jù)為各位取反并加一。比如4位補(bǔ)碼數(shù)據(jù)0111（即7），那么其反向數(shù)據(jù)－7的補(bǔ)碼形式為1001，顯然是0111各位取反并加一。故梳狀單元的最后一級進(jìn)位端接高電平（即1）。 附 錄 - 分支運(yùn)算的換位和電路的恒等關(guān)系
 先以抽取濾波器的單級梳狀級為例： 依次類推就可以得到如圖7所示的抽取濾波器N級梳狀級結(jié)構(gòu)框圖。 2、CIC抽取濾波器的系統(tǒng)函數(shù)推導(dǎo) CIC抽取濾波器從j級開始到最后一級的系統(tǒng)函數(shù)包括兩種情況，一種是 的情況，另一種是 的情況。下面分別分析這兩種情況下的結(jié)果。 當(dāng)時，其系統(tǒng)函數(shù)為：  （A1） 也可以寫為：  （A2） （A2）式的展開多項(xiàng)式中 最高次冪的階數(shù)為：  （A3） 因此，式（A1）可變?yōu)椋?/font>  （A4） 其中為多項(xiàng)式系數(shù)。另外，我們還可以把（A1）式作二項(xiàng)式展開得到：  （A5） 如果以兩多項(xiàng)式的交叉乘積項(xiàng)來表示，那么（A5）式變?yōu)椋?/font>  （A6） 現(xiàn)令 ，則可以得到 ，從而我們可以得到：   （A7）據(jù)（A4）和（A7）式可知，當(dāng)時，有：   （A8）同理，當(dāng) 時，可以得到：  （A9）  ，  （A10） 參 考 文 獻(xiàn) 【1】E．B．Hogenauer, “An Economical Class of Digital Filters for Decimation and Interpolation,” IEEE Trans on ASSP,VOL.ASSP-29, NO.2,April 1981. 【2】楊小牛，樓才義，徐建良著. 軟件無線電原理與應(yīng)用. 電子工業(yè)出版社，2001。 【3】R. E. 克勞切 L. R. 拉賓納著，酆廣增譯，多采樣率數(shù)字信號處理， 北京：人民郵電出版社，1988。 【4】張健，向敬成編著. 軟件無線電技術(shù)導(dǎo)論. 電子科技大學(xué)出版社, 2001。 【5】丁玉美，高西全編著. 數(shù)字信號處理（第二版）. 西安電子科技大學(xué)出版社，2001。
完整的Word格式文檔51黑下載地址：
累積梳狀(CIC)濾波器分析與設(shè)計.doc
(632 KB, 下載次數(shù): 16)
2018-4-20 13:32 上傳
點(diǎn)擊文件名下載附件
下載積分: 黑幣 -5
| 
QQ好友和群
QQ空間
騰訊微博
騰訊朋友
收藏
淘帖
頂
踩
