【正文】 德豐，數(shù)字信號處理與應用，北京；清華大學出版社，2022 年（3）李莉，王沛，俞玉蓮，常寧，數(shù)字信號處理原理和算法實現(xiàn)，北京；清華大學出版社，2022 年（4）萬永革，數(shù)字信號處理的 MATLAB 實現(xiàn)，北京；科學出版社，2022 年魯東大學本科畢業(yè)設(shè)計18（5）劉衛(wèi)國，MATLAB 程序設(shè)計與應用，北京；高等教育出版社，2022 年（6）楊順遼，李永全，數(shù)字信號處理實現(xiàn)與實踐，武漢；華中科技大學出版社（7）王嘉梅，基于 MATLAB 的數(shù)字信號處理與時間開發(fā)，西安；西安電子科技大學出版社，2022 年附錄 A 巴特沃斯濾波器繪制巴特沃斯模擬低通濾波器的幅頻平方響應曲線，階數(shù)分別為1,6,10,15,20 n=0::2。 %頻率點 for ii=1:5 %取 4 種濾波器 switch ii case 1,N=1。 case 2,N=6。 case 3,N=10。 case 4,N=15。 case 5,N=20。 end [z,p,k]=buttap(N)。 %設(shè)計巴特沃斯濾波器 [b,a]=zp2tf(z,p,k)。 %將零點增益形式轉(zhuǎn)換為傳遞函數(shù)形式 [H,w]=freqs(b,a,n)。 %按 n 指定的頻率點給出頻率響應 magH2=(abs(H)).^2。 %給出傳遞函數(shù)幅度平方 hold on。 plot(w,magH2)。 %繪制傳遞函數(shù)平方 end xlabel(39。w/wc39。)。 ylabel(39。︱H(jw)︱^239。)。 title(39。Butterworth 模擬低通濾波器39。)。 text(,39。n=139。) %作必要標記 text(,39。n=639。) text(,39。n=1039。) text(,39。n=1539。)魯東大學本科畢業(yè)設(shè)計19 text(,39。n=2039。) grid on。 附錄 B 切比雪夫 I 型模擬低通濾波器繪制切比雪夫 I 型模擬低通濾波器的平方幅頻響應曲線，階數(shù)為 1,4,6,10 程序如下： n=0::2。 %設(shè)置頻率點for ii=1:4 switch ii case 1,N=1。 case 2,N=4。 case 3,N=6。 case 4,N=10。 end Rp=1。 %為通帶波紋為 1dB [z,p,k]=cheb1ap(N,Rp)。 %設(shè)計切比雪夫 I 型濾波器 [b,a]=zp2tf(z,p,k)。 %轉(zhuǎn)換為傳遞函數(shù)形式 [H,w]=freqs(b,a,n)。 %求的傳遞函數(shù)的頻率特性 magH2=(abs(H)).^2。 %求的傳遞函數(shù)的幅頻響應 posplot=[39。2,2,39。 num2str(ii)]。 %將數(shù)字 ii 換為字符串，與‘2，2’合并并賦給 posplot subplot(posplot)。 plot(w,magH2)。 title([39。N=39。 num2str(N)])。%將數(shù)字 N 轉(zhuǎn)換為字符串與‘N=’合并為標題 xlabel(39。w/wc39。)。 ylabel(39。切比雪夫 I︱H(jw) ︱^239。)。 grid onend附錄 C 切比雪夫 II 型模擬低通濾波器魯東大學本科畢業(yè)設(shè)計20 繪制切比雪夫 II 型模擬低通濾波器的平方幅頻響應曲線，階數(shù)分別為 1，4,6,10程序如下： n=0::2。 %設(shè)置頻率點 for ii=1:4 switch ii case 1,N=1。 case 2,N=4。 case 3,N=6。 case 4,N=10。 end Rs=16。 [z,p,k]=cheb2ap(N,Rs)。 %設(shè)計 ChebyshevII 型模擬低通濾波器 [b,a]=zp2tf(z,p,k)。 %轉(zhuǎn)化為傳遞函數(shù) [H,w]=freqs(b,a,n)。 %求出濾波器的頻率響應 magH2=(abs(H)).^2。 %頻率響應的幅度平方 posplot=[39。2,2,39。 num2str(ii)]。 subplot(posplot)。 plot(w,magH2)。 %繪出幅頻平方響應函數(shù) title([39。N=39。 num2str(N)])。 xlabel(39。w/wc39。)。 ylabel(39。切比雪夫 II 型 ︱H(jw)︱^239。)。 grid on。 end附錄 D 例：繪制橢圓濾波器的幅頻平方響應曲線，階數(shù)分別為 1,4,6,10程序如下所示： n=0::2。 %設(shè)置頻率點for ii=1:4 switch ii case 1,N=1。 case 2,N=4。 case 3,N=6。 case 4,N=10。 end魯東大學本科畢業(yè)設(shè)計21 Rp=1。Rs=15。 %設(shè)置通帶波紋為 1dB,阻帶衰減為 15dB [z,p,k]=ellipap(N,Rp,Rs)。 %設(shè)計橢圓濾波器 [b,a]=zp2tf(z,p,k)。 %將零極點增益形式轉(zhuǎn)化為傳遞函數(shù)形式 [H,w]=freqs(b,a,n)。 %求得傳遞函數(shù)的復數(shù)頻率響應 magH2=(abs(H)).^2。 %給出幅度函數(shù)平方函數(shù) posplot=[39。2,2,39。 num2str(ii)]。 %給出繪圖區(qū)域 subplot(posplot)。 plot(w,magH2)。 ylim ([0,1]) %繪出幅度平方函數(shù) title([39。N=39。 num2str(N)])。 xlabel(39。w/wc39。)。 ylabel(39。橢圓 ︱H(jw) ︱^239。)。 grid onend 程序運行結(jié)果： 附錄 E 脈沖響應不變法設(shè)計 IIR 濾波器設(shè)計通帶截止頻率 = ，阻帶下線頻率 = ，通帶最大衰減 為 3dB，阻帶p??s??p?最小衰減 為 20dB，給定 =?STTs=Ap=3。As=20。OmegaP=*pi/Ts。OmegaS=*pi/Ts。 %模擬通帶、阻帶截止頻率[n,Wn]=buttord(OmegaP,OmegaS,Ap,As,39。s39。)。 %確定最小階數(shù) n 和反歸一化截止頻率Wnfprintf(39。濾波器的階數(shù) N=%.0f\n39。,n)。[b,a]=butter(n,Wn,39。s39。)。 %b、a 分別為模擬濾波器的分子分母按降冪排列的多項式系數(shù)[bz,az]=impinvar(b,a,1/Ts)。 %脈沖響應不變法得到數(shù)字濾波器的分子、分母系數(shù)魯東大學本科畢業(yè)設(shè)計22disp(39。分子系數(shù) b39。)。 %下面是顯示分子、分母多項式系數(shù)fprintf(39。%.4e39。,bz)。fprintf(39。\n39。)。disp(39。分母系數(shù) a39。)。 fprintf(39。%.4e39。,az)。fprintf(39。\n39。)。omega=[0::pi]。 %確定坐標軸范圍h=freqz(bz,az,omega)。 %得到模擬濾波器的單位沖擊響應系數(shù)Ampli=20*log10(abs(h)/abs(h(1)))。 %求衰減的分貝subplot(2,1,1)。plot(omega/pi,Ampli,39。k39。)。 %顯示濾波器的幅度響應xlabel(39。數(shù)字頻率 /\pi39。)。ylabel(39。幅度/dB39。)。grid。subplot(2,1,2)。theta=phasez(bz,az,omega)。 %濾波器的相位響應及坐標值plot(omega/pi,theta*360/(2*pi),39。k39。)。 %顯示濾波器的相位響應xlabel(39。數(shù)字頻率 /\pi39。)。ylabel(39。相位/度39。)。grid。附錄 F 雙線性變換法設(shè)計 II