【正文】 繪制截取后的語音信號title(39。原始語音信號39。)。xlabel(39。時間(單位：s)39。)。ylabel(39。幅度39。)。axis tightsubplot(2,2,3)。plot(t,y) % 繪制截取后的語音信號title(39。加入單頻干擾后的語音信號39。)。xlabel(39。時間(單位：s)39。)。ylabel(39。幅度39。)。axis tightX=abs(fft(x))。 Y=abs(fft(y))。 % 對原始信號和加噪信號進行fft變換,取幅度譜X=X(1:N/2)。 Y=Y(1:N/2)。 % 截取前半部分deltaf=fs/N。 % 計算頻譜的譜線間隔f=0:deltaf:fs/2deltaf。 % 計算頻譜頻率范圍subplot(2,2,2)。plot(f,abs(X)) % 繪制頻譜圖title(39。語音信號幅度譜圖39。)。xlabel(39。頻率（單位：Hz）39。)。ylabel(39。幅度譜39。) 。axis tightsubplot(2,2,4)。plot(f,abs(Y)) % 繪制頻譜圖title(39。加入干擾后的語音信號幅度譜圖39。)。xlabel(39。頻率（單位：Hz）39。)。ylabel(39。幅度譜39。) 。axis tightfp=2000。fc=2150。 fs=8000 %定義通帶和阻帶截止頻率Rp=1。As=50。 % 定義通帶波紋和阻帶衰減wp=fp/fs*2*pi。ws=fc/fs*2*pi。 %計算對應的數字頻率T=2。 %定義采樣間隔Omegap=2/T*tan(wp/2)。Omegas=2/T*tan(ws/2)。 %截止頻率預畸變[cs,ds]=afd_chb1(Omegap,Omegas,Rp,As) % 計算濾波器系統(tǒng)函數分子分母系數[b,a]=bilinear(cs,ds,1/T) % 雙線性變換得到數字濾波器系統(tǒng)函數分子分母系數[db,mag,pha,grd,w]=freqz_m(b,a)。 %驗證濾波器是否達到指定性能delta=[1,zeros(1,99)]。ha=filter(b,a,delta)。Subplot(2,2,1)。plot(w/pi,db,39。LineWidth39。,2)。axis([0,1,80,10])。grid on。% 繪制相對幅度響應曲線圖，限制橫縱坐標范圍，添加網格X_l=[0,0,wp/pi,ws/pi。1,1,wp/pi,ws/pi]。Y_l=[As,Rp,80,80。As,Rp,10,10]。 % 在wp，ws，Rp，As處畫線以更直觀判斷設計是否達標，每列參數是每個線條的端點坐標line(X_l,Y_l,39。Color39。,39。r39。,39。LineWidth39。,2,39。LineStyle39。,39。39。) % 添加線條，紅色，線寬為2title(39。濾波器幅度響應圖39。)。xlabel(39。w/pi39。)。ylabel(39。dB39。)。subplot(2,2,2)。plot(w/pi,mag,39。LineWidth39。,2)。axis([0,1,0,])。grid on。title(39。濾波器幅度響應圖39。)。xlabel(39。w/pi39。)。ylabel(39。幅度mag39。)。% 繪制絕對幅度響應曲線圖subplot(2,2,3)。plot(w/pi,pha)。axis tight。grid on。title(39。濾波器相位響應圖39。)。xlabel(39。w/pi39。)。ylabel(39。相位pha39。)。% 繪制相位特性曲線圖subplot(2,2,4)。plot(ha,39。LineWidth39。,2)。 axis([0,100,])。grid on。title(39。濾波器脈沖響應圖39。)。xlabel(39。n39。)。ylabel(39。ha39。)。% 繪制脈沖響應圖y_fil=filter(b,a,y)。 % 對信號y進行濾波處理，其中b,a為上面設計好的濾波器參數Y_fil=abs(fft(y_fil))。Y_fil=Y_fil(1:N/2)。 % 計算頻譜取前一半sound (y_fil,fs,bits)figure。subplot(3,2,1)。plot(t,x)。axis tight。grid on。title(39。原始語音信號時間x39。)。xlabel(39。時間t39。)。ylabel(39。幅度39。)。% 繪制原始錄制的語音信號時間波形圖subplot(3,2,3)。plot(t,y)。axis tight。grid on。title(39。加干擾語音信號時間y39。)。xlabel(39。時間t39。)。ylabel(39。幅度39。)。% 繪制加噪語音信號時間波形圖subplot(3,2,5)。plot(t,y_fil)。axis tight。grid on。title(39。濾波后語音信號時間y_fil39。)。xlabel(39。時間t39。)。ylabel(39。幅度39。)。% 繪制去噪語音信號時間波形圖subplot(3,2,2)。plot(f,abs(X))。axis([0 3000 0 60])。grid on。title(39。原始語音信號幅度譜X39。)。xlabel(39。頻率f39。)。ylabel(39。幅度39。)。% 繪制原始錄制的語音信號頻譜圖subplot(3,2,4)。plot(f,abs(Y))。axis([0 3000 0 60])。grid on。title(39。加干擾語音信號幅度譜Y39。)。xlabel(39。頻率f39。)。ylabel(39。幅度39。)。% 繪制加噪語音信號頻譜圖subplot(3,2,6)。plot(f,abs(Y_fil))。axis([0 3000 0 60])。grid on。title(39。濾波后語音信號幅度譜Y_fil39。)。xlabel(39。頻率f39。)。ylabel(39。幅度39。)。% 繪制去噪語音信號頻譜圖[C,B,A]=dir2par(b,a)