【文章內(nèi)容簡介】 11,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 }。 pragma DATA_SECTION(db,.dbuffer) DATA db[NH+2] 。 DATA rtest[NX]。 DATA r[NX]。 運(yùn)行結(jié)果如圖 所示 圖 濾波前后的的波形 FIR 帶阻濾波器實(shí)現(xiàn) (1)、 MATLAB 程序代碼 ws1=。ws2=。 wp1=。wp2=。 wp=[wp1,wp2]。ws=[ws1,ws2]。 B=ws1wp1。%計算過渡帶寬 N0=ceil(*pi/B)。 N=N0+mod(N0+1,2) wc1=(ws1+wp1)/2。wc2=(ws2+wp2)/2。 wc=[wc1,wc2]。 b=fir1(N1,wc,39。stop39。)。 %80 point hamming window FIR filter, Wn= fvtool(b,1) h=buffer(round(b*32768),8)39。 csvwrite(39。39。,h) 運(yùn)行結(jié)果如下： 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 程序運(yùn)行結(jié)果如圖 1 所示。 圖 幅頻響應(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 }。 pragma DATA_SECTION(db,.dbuffer) DATA db[NH+2] 。 DATA rtest[NX]。 DATA r[NX]。 運(yùn)行結(jié)果如圖 所示： 圖 濾波前后的波形 (二 )、 IIR 濾波器程序代碼 CCS 軟件環(huán)境下主程序代碼 define PI include include include include include short i。 short eflag = PASS。 // error flag or index into r vector where error void main() { for (i=0。iNX。i++) r[i] =0。 for (i=0。 i5*NBIQ。 i++) dbuffer[i] = 0。 // clear delay buffer (a must) for(i=0。iNX。i++)x[i]=0。 for (i=0。iNX。i++) { x[i]=(short)((1+sin(PI*3000*i/12020))*sin(PI*4000*i/12020)*32768/3) 。 } iircas51(x, h, r, dp, NBIQ, NX)。 while(1)。 } IIR 低通濾波器實(shí)現(xiàn) (1)、 MATLAB 程序代碼 clear clc wp=。 ws=。 [n,wc]=ellipord(wp,ws,3,80) [b,a]=ellip(n,3,80,wc) freqz(b,a,512,10000) [s, g]=tf2sos(b,a) s1=round(s/4*32768) k=[s1(:,1:3),s1(:,5:6)] s3=[k(1,:),k(2,:),k(3,:),k(4,:)] h1=buffer(s3,8)39。 csvwrite(39。39。,h1) 程序運(yùn)行結(jié)果如下： 生成 文檔 8192,8192,0,6022,0,8192,12450,8192 7816,5556,8192,7398,8192,2514,7031,8192 5145,8192,180,7881,0,0,0,0 程序運(yùn)行結(jié)果如圖 所示 . 圖 幅頻和相頻特性 (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。 . 運(yùn)行結(jié)果如圖 所示。 圖 濾波前后的波形 IIR 高通濾波器實(shí)現(xiàn) (1)、 MATLAB 程序代碼 clear clc ws= wp= %[n,wn]=cheb1ord(wp,ws,3,80) %[b,a]=cheby1(n,3,wn,39。high39。) [n,wc]=ellipord(wp,ws,3,80) [b,a]=ellip(n,3,80,wc,39。high39。) 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。 csvwrite(39。39。,h1) 運(yùn)行結(jié)果如下： 生成 文檔 8192,16277,8192,237,2948,8192,15678,8192 10398,6810,8192,15204,8192,12912,7894,0 運(yùn)行結(jié)果如圖 1 所示 圖 幅頻和相頻特性 (2)、將生成的 文檔粘貼到頭文件中的 DATA h[5*NBIQ]{ 8192,16277,8192,237,2948,8192,15678,8192, 10398,6810,8192,15204,8192,12912,7894 }中 ，然后將 NBIQ 的值改為 3。 程序運(yùn)行結(jié)果如圖 所示。 圖 濾波前后的波形圖 IIR 帶通濾波器實(shí)現(xiàn) (1)、 MATLAB 程序代碼 clear clc wp=[,]。 ws=[,]。 [n,wc]=ellipord(wp,ws,3,80) [b,a]=ellip(n,3,80,wc,39。bandpass39。) 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=[