【正文】 [k(1,:),k(2,:),k(3,:),k(4,:)] h1=buffer(s3,8)39。 } IIR 低通濾波器實(shí)現(xiàn) (1)、 MATLAB 程序代碼 clear clc wp=。 } iircas51(x, h, r, dp, NBIQ, NX)。iNX。i++)x[i]=0。 // clear delay buffer (a must) for(i=0。 i5*NBIQ。i++) r[i] =0。 // error flag or index into r vector where error void main() { for (i=0。 運(yùn)行結(jié)果如圖 所示： 圖 濾波前后的波形 (二 )、 IIR 濾波器程序代碼 CCS 軟件環(huán)境下主程序代碼 define PI include include include include include short i。 DATA rtest[NX]。 圖 幅頻響應(yīng) (2)、 CCS 程序代碼 程序代碼 define NX 320 define NH 209 //define FNAME t6 define MAXERROR 10 define PI pragma DATA_SECTION(h,.coeffs) DATA h[NH]={ 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 }。39。 %80 point hamming window FIR filter, Wn= fvtool(b,1) h=buffer(round(b*32768),8)39。stop39。 wc=[wc1,wc2]。 N=N0+mod(N0+1,2) wc1=(ws1+wp1)/2。 B=ws1wp1。 wp=[wp1,wp2]。 wp1=。 運(yùn)行結(jié)果如圖 所示 圖 濾波前后的的波形 FIR 帶阻濾波器實(shí)現(xiàn) (1)、 MATLAB 程序代碼 ws1=。 DATA rtest[NX]。,h) 運(yùn)行結(jié)果： N =209 生成 文檔 3,15,9,5,0,5,10,18 4,0,12,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,302,188,94 0,102,223,392,79,0,267,74 402,286,0,314,484,98,388,0 140,761,478,241,0,271,604,1083 224,0,799,230,1297,967,0,1193 1978,439,1928,0,943,6515,5822,5561 13098,5561,5822,6515,943,0,1928,439 1978,1193,0,967,1297,230,799,0 224,1083,604,271,0,241,478,761 140,0,388,98,484,314,0,286 402,74,267,0,79,392,223,102 0,94,188,302,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,12,0 4,18,10,5,0,5,9,15 3,0,0,0,0,0,0,0 程序運(yùn)行結(jié)果如圖 所示 . 圖 幅頻特性 (2)、 CCS 程序代碼 代碼 define NX 320 define NH 209 //define FNAME t6 define MAXERROR 10 define PI pragma DATA_SECTION(h,.coeffs) DATA h[NH]={ 3,15,9,5,0,5,10,18, 4,0,12,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,302,188,94, 0,102,223,392,79,0,267,74, 402,286,0,314,484,98,388,0, 140,761,478,241,0,271,604,1083, 224,0,799,230,1297,967,0,1193, 1978,439,1928,0,943,6515,5822,5561, 13098,5561,5822,6515,943,0,1928,439, 1978,1193,0,967,1297,230,799,0, 224,1083,604,271,0,241,478,761, 140,0,388,98,484,314,0,286, 402,74,267,0,79,392,223,102, 0,94,188,302,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,12,0, 4,18,10,5,0,5,9,15, 3 }。 csvwrite(39。 b=fir1(N1,wc,hamming(N))。wc2=(ws2+wp2)/2。%計(jì)算過(guò)渡帶寬 N0=ceil(*pi/B)。ws=[ws1,ws2]。ws2=。wp2=。 運(yùn)行結(jié)果如圖 所示。 DATA rtest[NX]。 圖 幅頻響應(yīng) (2)、 CCS 程序代碼 代碼 define NX 320 define NH 209 define MAXERROR 10 define PI pragma DATA_SECTION(h,.coeffs) DATA h[NH]={ 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 }。39。 。)。 b=fir1(N1,wc,39。%計(jì)算過(guò)渡帶寬 N0=ceil(*pi/B)。 ws=。 DATA r[NX]。 pragma DATA_SECTION(db,.dbuffer) DATA db[NH+2] 。 運(yùn)行 結(jié)果如下： N=104 生成 文檔 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 圖形如圖 所示。,h) %參數(shù)表輸出到文件中 %然后只要把參數(shù)表從 中拷到代碼中，每行前面加上 39。 csvwrite(39。 %80 point hamming window FIR filter, Wn= fvtool(b,1) %Flter Visualization Tool 可心以查看幅度響應(yīng)、沖激響應(yīng)等 h=buffer(round(b*32768),8)39。%計(jì)算過(guò)渡帶寬 N=ceil(*pi/B) %按哈明窗計(jì)算濾波器長(zhǎng)度 wc=(ws+wp)/2。 ws=。//雙 MAC 運(yùn)算，速度快 while(1)。 i++) db[i] = 0。 // clear output buffer (optional) for (i=0。 iNX。i++) { x[i]=(short)(((1+sin(PI*3000*i/10000))*si