【正文】 B, A 變?yōu)榻频葍r的數(shù)字濾波器的系數(shù)Bd, Ad.（5）把( 2) , ( 3) , ( 4)合為一步的數(shù)字濾波器設計函數(shù)butter(N, w c, ? ftype! )。 [ B, A] = butter( N, w c)設計低通或帶通數(shù)字濾波器系數(shù)B, A (當為帶通濾波器時, 第( 1)類函數(shù)由wp = [ w p1, w p2 ] 會自動生成w c = [ w 1,w2] )。[ B, A] = butter( N, w c, ? h igh! )設計高通數(shù)字濾波器系數(shù)B, A。[ B, A] = butter( N, w c, ? stop! )設計帶阻數(shù)字濾波器系數(shù)B, A。由此可見, 通常情況下, 當采用的是雙線性變換法, 可用( 1)和( 5)兩類函數(shù)設計數(shù)字濾波器。 使用MATLAB設計程序用雙線性變換法設計數(shù)字低通濾波器t= 0: 0. 001: 1。X= sin( 2* p *i 100* t) + sin( 2* p*i 200*t) 。 輸入信號subplot( 2, 3, 1)plot( X ( 1: 50) )。 g rid。 t itle ( ‘輸入信號!’)。x labe l( t! )。 y labe l( f( t) ! )Xa= fft( X, 512) 。 f= 1000* ( 0: 256) /512。 對序列進行快速傅立葉變換subplot( 2, 3, 2)plot( ,f Xa( 1: 257) ) 。 g rid。 tit le( |Ha( jf) |! )x labe l( f( Hz) ! )。 y labe l( 幅度! ) Wp= 0. 2。 W s= 0. 4。 Rp= 3。 Rs= 20輸入數(shù)字低通濾波器的指標[ N, w c] = buttord(W p, W s, Rp, R s)計算數(shù)字濾波器的階數(shù)N 和頻率參數(shù)w c[ B, A] = butter( N, w c)求低通數(shù)字濾波器系數(shù)B, Aw d= [ 0: 118]* p i/512設定數(shù)字歸一化頻率[ h, w ] = freqz( B, A, wd)計算數(shù)字濾波器頻率響應subplot( 2, 3, 3)plot( w d /p,i 20* log10( abs( h) ) ) 。 grid。 xlabe l( ‘ f(H z)’)。 y labe l( ‘幅度dB’ )fv too l( B, A)[H, T] = impz( B, A )s= conv(H, X )用時域卷積求濾波后的信號subplot( 2, 3, 4)plot( s( 1: 50) )。 grid。 x labe l( ‘t’ )。 y labe l( ‘ s( t)‘ )對程序的說明: 程序有四個部分,第一部分是產(chǎn)生需要濾波的信號, 畫出其時域圖。第二部分是對信號作頻譜分析, 并作出頻譜圖。第三部分是根據(jù)濾波器指標產(chǎn)生數(shù)字低通濾波器, 并作出其幅頻特性。第四部分是進行濾波, 作出濾波后信號的時域圖及幅頻特性圖這里使用的是巴特沃思濾波器,濾波范圍為低通,選取的抽樣頻率20000hz,通帶拐角頻率為2000hz,幅頻特性,想頻特性和歸一化處理后的幅頻特性及想頻特性: define IIR16_COEFF { 0,3794,0,307,307, 2159,7894,668,1335,668, 3483,8904,541,1082,541, 6131,10923,13364,26729,13364} define IIR16_ISF 4398define IIR16_NBIQ 4define IIR16_QFMAT 13 程序運行過程及結果1. 運行過程啟動CCS, 新建工程文件, 文件名為iir. asm。選擇Pro ject菜單中的Add F ille to Pro ject選項,將編好的匯編源程序添加到工程文件中。新建配置文件, 配置文件編輯工具同時會生成iir. cmd文件, 然后將它們添加到工程文件中。選擇project bu ild all進行編譯、連接。編譯無錯后, 生成iir. out文件, 選擇F ile / load prog ram 將其裝入DSP的片內(nèi)或外部擴展寄存器(確定裝入地址有CMD 文件) , 然后設置指針及斷點。選擇V iew G raph T ime /Frequency, 將S tartAddress 改為0x3100, 將 D isplay Data Size 設置為256, 將 DSP Data Type 改為 16 bit signedinteger 。選擇Debug Run運行程序。這樣,將在圖形窗口中顯示從0x3220 (信號輸入緩沖器)開始的128個點的16位有符號整數(shù)。將鼠標移至圖上, 點擊右鍵選擇propert ies, 重新打開一個窗口, 將其中的?? Signal T ime 改為 FFTMagnitude , 顯示頻域波形。再打開一個窗口。顯示從地址0x3100 (濾波信號輸出緩沖)開始的128點的16位有符號整數(shù)的濾波輸出。程序運行結果如下：濾波前CCS中的數(shù)據(jù)時域波濾波前ccs的數(shù)據(jù)頻域波形如上濾波后ccs中的時域波形濾波后ccs中的頻域波形參考文獻 （漢科技大學電子技術學院）[J]2. 趙曦，劉睿芝，（陜西理工學院物理系）[J]3. 艷,（宿州學院）[J]4．姜乃卓，都思丹IIR 數(shù)字濾波器的系數(shù)量化效應分析P811 ( 南京大學電子科學與工程系) [J] 5．周耀輝，王蕓波，朱維新，張玉仲，武煥舟. IIR 數(shù)字濾波器設計P129131（1. 中國石化管道儲運公司儀長輸油處，湖北武漢430077；2. 深圳山東核電工程有限責任公司，山東濟南250100）[J]6. 戴朝華 , 陳維榮 ,　朱云芳 , 鄭永康 ,　李　 ( 西南交通大學電氣工程學院)[J] IIR數(shù)字濾波器優(yōu)化及FPGA仿真驗證（西南交通大學）[R]，（科學出版社）[M]（清華大學出版社）[M]，黃繼業(yè)，（先電子科技大學出版社）[M]（山東大學）[R]. 基于DSP 的數(shù)字濾波器(IIR)的實現(xiàn)P26 ( 吉林省公主嶺市廣播電視局)[J]　星,唐　寧,嚴　舒,楊　白. 基于FPGA的I IR數(shù)字濾波器的設計與仿真P304307 (桂林電子科技大學信息與通信學院)[J]（西南交通大學）[R],彭杰. 基于FPGA 的IIR 數(shù)字濾波器的設計與實現(xiàn)P105106[J] 華 ,孫運強 ,姚愛琴. 基于FPGA的IIR數(shù)字濾波器的實現(xiàn)（1 中北大學 信息與通信工程學院 山西 太原 030051；2 中北大學 儀器科學與動態(tài)測試教育部重點實驗室 山西 太原 030051）[J]17. 梁二虎，劉文怡，張文棟基于FPGA 的IIR 數(shù)字濾波器硬件模塊的設計P218219 (電子測試技術國家重點實驗室)[J]18. 閻威,劉鳳新基于FPGA 的高精度浮點IIR 濾波器設計P203205（北京化工大學）[J]19. 田永毅，倪宏斌. 基于M atlab的帶通IIR數(shù)字濾波器設計與仿真( 1. 陜西工業(yè)職業(yè)技術學院電氣工程學院, 陜西咸陽712000。2. 中國通信建設第二工程局有限公司, 陜西西安710119)[J]20. 鄭 群, 徐 南,王麗華. 基于多項式優(yōu)化的IIR濾波器設計P45(秦皇島職業(yè)技術學院)[J]