【正文】 [alpha*(2*pi)/M*k1,alpha*(2*pi)/M*(Mk2)]。H=Hrs.*exp(j*angH)。h=real(ifft(H,M))。%求FIR單位沖擊響應[db,mag,pha,grd,w]=freqz_m(h,1)。%幅頻特性曲線[Hr,ww,a,L]=hr_type2(h)。%畫圖―――――――――――――――――――――――――――――subplot(2,2,1)。stem(l,h)。title(39。濾波器單位沖擊響應39。)。xlabel(39。n39。)。ylabel(39。y(n)39。)。axis([0 40 1 1])。grid on。subplot(2,2,2)。plot(w/pi,pha)。title(39。濾波器相位響應39。)。xlabel(39。w/pi39。)。ylabel(39。相位pha39。)。axis([0 1 4 4])。grid on。subplot(2,2,3)。plot(w/pi,db)。title(39。濾波器幅度響應39。)。xlabel(39。w/pi39。)。ylabel(39。dB39。)。axis([0 1 50 20])。grid on。%畫線X_l=[0,0,wp/pi,ws/pi。1,1,wp/pi,ws/pi]。Y_l=[As,Rp,50,50。As,Rp,20,20]。 % 在wp，ws，Rp，As處畫線以更直觀判斷設(shè)計是否達標，每列參數(shù)是每個線條的端點坐標line(X_l,Y_l,39。Color39。,39。r39。,39。LineWidth39。,2,39。LineStyle39。,39。39。) % 添加線寬為2的紅色線條subplot(2,2,4)。plot(w/pi,mag)。title(39。濾波器幅度響應39。)。xlabel(39。w/pi39。)。ylabel(39。幅度mag39。)。axis([0 1 0 ])。grid on。%―――――――――――――――――――――――――――――――%進行濾波處理z=filter(h,1,y)。 % 用設(shè)計好的濾波器對x1進行濾波Z=fft(z)。Z=Z(1:length(Z)/2)。 % 計算頻譜取前一半%―――――――――――――――――――――――――――――――%比較濾波前后的性能%在同一張大圖里分別繪制原始信號x，加噪信號y，濾波去噪信號z的時域波形和頻譜，以便比較和分析。figure(2)。subplot(321)。plot(t,x)。axis tight。grid on。title(39。原始語音信號時間x39。)。xlabel(39。時間t39。)。ylabel(39。幅度39。)。% 繪制原始錄制的語音信號時間波形圖subplot(323)。plot(t,y)。axis tight。grid on。title(39。加干擾語音信號時間y39。)。xlabel(39。時間t39。)。ylabel(39。幅度39。)。% 繪制加噪語音信號時間波形圖subplot(325)。plot(t,z)。axis tight。grid on。title(39。濾波后語音信號時間z39。)。xlabel(39。時間t39。)。ylabel(39。幅度39。)。% 繪制去噪語音信號時間波形圖subplot(322)。plot(f,abs(X))。axis tight。grid on。title(39。原始語音信號幅度譜X39。)。xlabel(39。頻率f39。)。ylabel(39。幅度39。)。% 繪制原始錄制的語音信號頻譜圖subplot(324)。plot(f,abs(Y))。axis tight。grid on。title(39。加干擾語音信號幅度譜Y39。)。xlabel(39。頻率f39。)。ylabel(39。幅度39。)。axis([0 4000 0 200])。% 繪制加噪語音信號頻譜圖subplot(326)。plot(f,abs(Z))。axis tight。grid on。title(39。濾波后語音信號幅度譜Z39。)。xlabel(39。頻率f39。)。ylabel(39。幅度39。)。% 繪制去噪語音信號頻譜圖sound(z,fs,bits)。%聽取濾波后的聲音信號figure(3)。subplot(2,2,1)。plot(t,x)。xlabel(39。時間（單位：S）39。)。ylabel(39。幅度39。)。title(39。原始語音信號39。)。axis([0 1 1])。grid on。subplot(2,2,2)。plot(f,X)。xlabel(39。頻率（單位：Hz）39。)。ylabel(39。幅度譜39。)。title(39。語音信號幅度譜圖39。)。axis([0 4000 0 200])。grid on。subplot(2,2,3)。plot(t,y)。xlabel(39。時間（單位：S）39。)。ylabel(39。幅度39。)。title(39。加入單頻干擾后的語音信號39。)。axis([0 1 1])。grid on。subplot(2,2,4)。plot(f,Y)。xlabel(39。頻率（單位：Hz）39。)。ylabel(39。幅度譜39。)。title(39。加入干擾后的語音信號幅度譜圖39。)。axis([0 4000 0 200])。grid on。%―――――――――――――――――――――――――――――――[C,B,A]=dir2fs(h)%調(diào)用函數(shù)，得出濾波器參數(shù)wavwrite(z,fs,bits,39。濾波后的語音信號z39。)%生成濾波后的wav聲音文件21