【正文】 參考文獻 [1]、高西全 ,丁玉美 .數(shù)字信號處理 (第三版 ).西安電子科技大學出版社 .2021 年 8 月 [2]、劉樹棠 [譯 ].數(shù)字信號處理 (MATLAB 版 )(第二版 ).西安交通大學出版社 .2021 年 1 月 。從圖 2 26看來， IIR 濾波器的的相位是非 線性的，其處理后的信號存在嚴重的相位失真。 (3)兩種濾波器的損耗函數(shù)、相頻特性的比較圖見圖 2 26： 圖 25 損耗函數(shù)比較圖 圖 26 相頻特性比較圖 (4)IIR 濾波器的階數(shù)： N1=3 FIR 濾波器的階數(shù)： N=24 三、 比較結(jié)果分析： 在相同的濾波器幅頻特性指標下， FIR 濾波器的階數(shù)都要0 0 . 1 0 . 2 0 . 3 0 . 4 0 . 5 0 . 6 0 . 7 0 . 8 0 . 9 1 1 2 0 1 0 0 8 0 6 0 4 0 2 0020w/ ?幅度/dB損耗函數(shù)曲線 IIRF IR0 0 . 1 0 . 2 0 . 3 0 . 4 0 . 5 0 . 6 0 . 7 0 . 8 0 . 9 11 0 . 8 0 . 6 0 . 4 0 . 200 . 20 . 40 . 60 . 81w/ ?相位/?相頻特性曲線 IIRF I R比 IIR 高 5~10 倍，設計成本較高，信號延時比較大。FIR39。IIR39。)。 title(39。相位 /\pi39。)。grid on hold off xlabel(39。linewidth39。:39。 figure(2) plot(wk1/pi,pah/pi)。FIR39。IIR39。)。) title(39。ylabel(39。w/\pi39。,3)。,39。grid on hold on plot(wk/pi,20*log10(abs(Hk(k+1))),39。pah=angle(Hk1)。) [Hk1,wk1]=freqz(B,A)。 0 0 . 1 0 . 2 0 . 3 0 . 4 0 . 5 0 . 6 0 . 7 0 . 8 0 . 9 1 1 2 0 1 0 0 8 0 6 0 4 0 2 0020w/ ?幅度/dB損耗函數(shù)曲線 IIRF IR0 1 2 3 . 4 0 . 5 0 . 6 0 . 7 0 . 8 0 . 9 11 0 . 8 0 . 6 0 . 4 0 . 200 . 20 . 40 . 60 . 81w/ ?相位/?相頻特性曲線 IIRF IR[B,A]=ellip(N1,Rp,As,wp1,39。As=50。ws2=2*fs/Fs。 wk=(2*pi/M)*k。 Hk=fft(hn,M)。 hn=real(ifft(Hk))。 thetak=pi*(N1)*(0:N1)/N。 Ak(Np+2)=T。Ns=N2*Np1。 N=N+mod(N+1,2)。 m=1。 datB=wpws。 wp=2*pi*fp/Fs。 fp=600。 (2)兩種濾波器的損耗函數(shù)、相頻特性的比較圖見圖 2 24： 圖 23 損耗函數(shù)比較圖 圖 24 相頻特性比較圖 (3)IIR 濾波器的階數(shù)： N1=3 FIR 濾波器的階數(shù)： N=17 二、高通濾波器部分： (1)matlab 代碼： T=。FIR39。IIR39。)。 title(39。相位 /\pi39。)。grid on hold off xlabel(39。linewidth39。:39。 figure(2) plot(wk1/pi,pah/pi)。FIR39。IIR39。)。) title(39。ylabel(39。w/\pi39。,3)。,39。grid on hold on plot(wk/pi,20*log10(abs(Hk(k+1))),39。pah=angle(Hk1)。 [B,A]=ellip(N1,Rp,As,wp1) [Hk1,wk1]=freqz(B,A)。As=50。ws2=2*fs/Fs。 wk=(2*pi/M)*k。 Hk=fft(hn,M)。 hn=fir1(N,wc,kaiser(N+1,alph))。 N=ceil((rs8)/)。 Bt=wswp。 wp=2*pi*fp/Fs。fs=600。 一、低通濾波器部分： (1)matlab 代碼： Fs=3800。畫出所設計的濾波器頻率特性圖、信號時域圖；給出濾波器設計的 MATLAB 代碼與濾波器實現(xiàn)的代碼；選擇合適的信號采樣周期 T。)。ylabel(39。頻率 /Hz39。)。 title(39。 figure(4) plot(fs,abs(S))。 S=fft(y,N)。n39。)。 ylabel(39。信號 y(n)39。 figure(3) plot(n,y)。幅值 39。)。 xlabel(39。原信號幅度頻譜 (采樣點數(shù)為 500)39。 axis([0,1500,0,180])。fs=n/(N*T)。)。xlabel(39。幅值 39。)。 title(39。m1=x(m)。m3=m2。m5=m4。 m7=m6。m9=m8。m11=m10。m13=m12。m15=m14。 for m=1:length(x) y(m)=*x(m)m1**+m3*+m4*+m5*98... m6**+m8***+m11*+m12*+... m13****。 m15=0。m13=0。m11=0。m9=0。m7=0。m5=0。m3=0。 m1=0。t=n*T。 f=2800。 (2)數(shù)字濾波器的損耗函數(shù)和相頻特性分別如圖 1 18所示： 0 0 . 1 0 . 2 0 . 3 0 . 4 0 . 5 0 . 6 0 . 7 0 . 8 0 . 9 1 8 0 7 0 6 0 5 0 4 0 3 0 2 0 1 00w/ ?幅度/dB損耗函數(shù)0 0 . 1 0 . 2 0 . 3 0 . 4 0 . 5 0 . 6 0 . 7 0 . 8 0 . 9 11 0 . 8 0 . 6 0 . 4 0 . 200 . 20 . 40 . 60 . 81w/ ?相位/?相頻特性曲線 圖 17 損 耗函數(shù)曲線 圖 18 相頻特性曲線 二、 按直接型網(wǎng)絡結(jié)構(gòu)編程 編寫濾波程序 : (1)matlab 代碼： N=500。相頻特性曲線 39。)。ylabel(39。w/\pi39。grid on figure(3) plot(wk/pi,angle(Hk(k+1))/pi)。損耗函數(shù) 39。)。ylabe