【正文】 次實習中選擇方法一。 方法二：獨立編寫 C 語言。其調(diào)用形式有： ushort oflag = iircas5 (DATA *x, DATA *h, DATA *r, DATA *dbuffer, short nbiq,ushort nx),其中 h 的存放系數(shù)順序為 h = a11 a21 b21 b01 b11 ... a1i a2i b2i b0i b1i。 2）、高階 IIR 濾波器的實現(xiàn) 一個高階 IIR 濾波器可以由若干個二階基本節(jié)相級聯(lián)或并聯(lián)構成。其方法是將高階 IIR 濾波器分解成若干個二階基本節(jié)，然后進行級聯(lián)來實現(xiàn)。 圖 4 i 個二階基本節(jié) IIR 濾波器程序流程圖 開始初始化DSP產(chǎn)生輸入數(shù)據(jù)讀入新數(shù)據(jù)并存放到棧頂濾波處理保存濾結果到輸出單元待濾波的樣本點數(shù)處理完成結束否是 圖 5 程序流程圖 五、 程序代碼及運行結果 (一 )、 FIR 濾波器程序代碼 CCS 軟件環(huán)境下主程序代碼如下 include include include include short i。 short eflag2= PASS。 DATA *dbptr = amp。 void main() { int i。 for(i=0。i++)x[i]=0。iNX。 } for (i=0。 i++) r[i] = 0。 iNH+2。 // clear delay buffer (a must) fir2(x, h, r, dbptr, NX, NH)。 } FIR 低通濾波器實現(xiàn) （ 1） MATLAB 程序代碼 clear all clc wp=。 B=wswp。 b=fir1(N1,wc)。 %參數(shù) *2^15,把小數(shù)點移最高位后面后取整，然后排成 8個 %一行便于輸出。39。 .word 39。 圖 幅頻響應 （ 2）、 CCS 軟件程序代碼 代碼 define NX 320 define NH 104 define MAXERROR 10 define PI pragma DATA_SECTION(h,.coeffs) DATA h[NH]={ 11,12,12,13,15,16,18,20, 23,26,29,33,37,42,47,52, 59,65,73,81,89,99,109,120, 132,145,158,173,189,207,226,247, 269,294,321,352,385,423,466,515, 572,639,719,816,938,1095,1307,1612, 2087,2936,4911,14758,14758,4911,2936,2087, 1612,1307,1095,938,816,719,639,572, 515,466,423,385,352,321,294,269, 247,226,207,189,173,158,145,132, 120,109,99,89,81,73,65,59, 52,47,42,37,33,29,26,23, 20,18,16,15,13,12,12,11 }。 DATA rtest[NX]。 運行結果如圖 所示： 圖 a 濾波前 x 的波形 圖 b 濾波后 r 的波形 圖 濾波前后的波形 FIR 高通濾波器實現(xiàn) (1)、 MATLAB 程序代碼 clear all clc wp=。 B=wpws。 N=N0+mod(N0+1,2) wc=(ws+wp)/2。high39。 fvtool(b,1) %Flter Visualization Tool 可心以查看幅 度響應、沖激響應等 h=buffer(round(b*32768),8)39。 csvwrite(39。,h) 運行結果： N=209 生成 文檔 0,6,8,6,0,6,9,7 0,8,11,8,0,9,14,10 0,12,18,13,0,15,23,17 0,19,29,22,0,24,36,27 0,30,45,34,0,37,56,41 0,46,68,50,0,55,82,61 0,67,99,73,0,80,118,88 0,96,141,104,0,114,168,124 0,135,200,148,0,161,239,177 0,193,286,212,0,234,347,258 0,287,428,320,0,360,541,408 0,469,715,548,0,654,1021,806 0,1043,1725,1468,0,2454,5211,7375 24578,7375,5211,2454,0,1468,1725,1043 0,806,1021,654,0,548,715,469 0,408,541,360,0,320,428,287 0,258,347,234,0,212,286,193 0,177,239,161,0,148,200,135 0,124,168,114,0,104,141,96 0,88,118,80,0,73,99,67 0,61,82,55,0,50,68,46 0,41,56,37,0,34,45,30 0,27,36,24,0,22,29,19 0,17,23,15,0,13,18,12 0,10,14,9,0,8,11,8 0,7,9,6,0,6,8,6 0,0,0,0,0,0,0,0 圖形如圖 所示。 pragma DATA_SECTION(db,.dbuffer) DATA db[NH+2] 。 DATA r[NX]。 圖 a 濾波前 x 的波形 圖 b 濾波后 r 的波形 圖 濾波前后的波形 FIR 帶通濾波器實現(xiàn) (1)、 MATLAB 程序代碼 wp1=。 ws1=。 wp=[wp1,wp2]。 B=wp1ws1。 N=N0+mod(N0+1,2) wc1=(ws1+wp1)/2。 wc=[wc1,wc2]。 fvtool(b,1) h=buffer(round(b*32768),8)39。39。 pragma DATA_SECTION(db,.dbuffer) DATA db[NH+2] 。 DATA r[NX]。ws2=。wp2=。ws=[ws1,ws2]。%計算過渡帶寬 N0=ceil(*pi/B)。wc2=(ws2+wp2)/2。 b=fir1(N1,wc,39。)。 csvwrite(39。,h) 運行結果如下： N=209 生成的 文檔 3,15,9,5,0,5,10,18 4,0,13,4,19,14,0,16 24,5,20,0,7,40,25,13 0,14,32,57,12,0,40,11 62,45,0,50,77,16,62,0 22,121,76,38,0,42,92,162 33,0,111,31,167,119,0,130 199,40,158,0,56,303,188,94 0,102,224,393,79,0,267,74 402,287,0,315,485,98,388,0 140,762,478,241,0,272,605,1085 224,0,800,231,1298,969,0,1194 1981,439,1930,0,945,6523,5829,5569 19672,5569,5829,6523,945,0,1930,439 1981,1194,0,969,1298,231,800,0 224,1085,605,272,0,241,478,762 140,0,388,98,485,315,0,287 402,74,267,0,79,393,224,102 0,94,188,303,56,0,158,40 199,130,0,119,167,31,111,0 33,162,92,42,0,38,76,121 22,0,62,16,77,50,0,45 62,11,40,0,12,57,32,14 0,13,25,40,7,0,20,5 24,16,0,14,19,4,13,0 4,18,10,5,0,5,9,15 3,0,0,0,0,0,0,0 程序運行結果如圖 1 所示。 pragma DATA_SECTION(db,.dbuffer) DATA db[NH+2] 。 DATA r[NX]。 short eflag = PASS。iNX。 for (i=0。 i++) dbuffer[i] = 0。iNX。 for (i=0。i++) { x[i]=(short)((1+sin(PI*3000*i/12020))*sin(PI*4000*i/12020)*32768/3) 。 while(1)。 ws=。 csvwrite(39。,h1) 程序運行結果如下： 生成 文檔 8192,8192,0,6022,0,8192,12450,8192 7816,5556,8192,7398,8192,2514,7031,8192 5145,8192,180,7881,0,0,0,0 程序運行結果如圖 所示 . 圖 幅頻和相頻特性 (2)、將生成的 文檔粘貼到頭文件中的 DATA h[5*NBIQ]{ 8192,8192,0,6022,0,8192,12450,8192， 7816,5556,8192,7398,8192,2514,7031,8192， 5145,8192,180,7881 }中 ，然后將 NBIQ 的值改為 4。 圖 濾波前后的波形 IIR 高通濾波器實現(xiàn) (1)、 MATLAB 程序代碼 clear clc ws= wp= %[n,wn]=cheb1ord(wp,ws,3,80) %[b,a]=cheby1(n,3,wn,39。) [n,wc]=ellipord(wp,ws,3,80) [b,a]=ellip(n,3,80,wc,39。) freqz(b,a,512,10000) [s, g]=tf2sos(b,a) %s2=s/20*32768 s1=round(s/4*32768) k=[s1(:,1:3),s1(:,5:6)] s3=[k(1,:),k(2,:),k(3,:)] h1=buffer(s3,8)39。39。 程序運行結果如圖 所示。 ws=[,]。bandpass39。 csvwrite(39。,h1) 程序運行結果如下 生成 文檔 8192,14384,8192,5594,6913,8192,16161,8192 9581,7105,8192,7753,8192,1969,7378,8192 15192,8192,12046,7650,8192,5014,8192,150 7940,8192,14644,8192,13081,8042,0,0 程序運行結果如圖 所示 . 圖 幅頻和相頻特性 (2)、將生成的 文檔粘貼到頭文件中的 DATA h[5*NBIQ]{ 8192,14384,8192,5594,6913,8192,16161,8192 9581,7105,8192,7753,8192,1969,7378,8192 15192,8192,12046,7650,8192,5014,8192,150 7940,8192,14644,8192,13081,8042 }中 ，然后將 NBIQ 的值改為 5。 圖 濾波前后的波形 IIR 帶阻濾波器實現(xiàn) (1)、 MATLAB 程序代碼 clear clc ws=[,]。 %[n,wn]=cheb1ord(wp,ws,3,80) %[b,a]=cheby1(n,3,wn,39。) [n,wc]=ellipord(wp,ws,3,80) [b,a]=ellip(n,