【正文】 x=wavread(39。D:\Backup\我的文檔 \My Files\39。)。 sound(x)。 fs=100。N=128。 y=fft(x,N)。 magy=abs(y)。 f=(0:length(y)1)39。*fs/length(y)。 subplot(221)。plot(f,magy)。 xlabel(39。頻率（ Hz） 39。)。ylabel(39。幅值 39。)。 title(39。N=128(a)39。)。grid subplot(222)。plot(f(1:N/2),magy(1:N/2))。 xlabel(39。頻率（ Hz） 39。)。ylabel(39。幅值 39。)。 title(39。N=128(b)39。)。grid fs=100。N=1024。 y=fft(x,N)。 magy=abs(y)。 f=(0:length(y)1)39。*fs/length(y)。 subplot(223)。plot(f,magy)。 xlabel(39。頻率（ Hz） 39。)。ylabel(39。幅值 39。)。 title(39。N=1024(c)39。)。grid subplot(224)。plot(f(1:N/2),magy(1:N/2))。 xlabel(39。頻率（ Hz） ?)。ylabel(39。幅值 39。)。 title(39。N=1024(d)39。)。grid 運行結(jié)果如圖 42 所示： 湖北師范學院 2020屆教育信息與技術學院學士學位論文（設計） 11 圖 42 語音信號定點分析 (3)N 階高通濾波器 在這里，以 5 階為例，其中 wc 為其 3dB 邊緣頻率 ,程序設計如下： x=wavread(?D:\Backup\我的文檔 \My Files\?)。 sound(x)。 N=5。wc=。 [b,a]=butter(N,wc,?high?)。 X=fft(x)。 subplot(321)。plot(x)。title(?濾波前信號的波形 ?)。 ylabel(?amplitude?)。xlabel(?time(ms)?)。 subplot(322)。 plot(abs(x))。axis([0,6144,0,1])。 title(?濾波前信號的頻譜 ?)。 ylabel(?amplitude?)。xlabel(?frequency(hz)?)。 y=filter(b,a,x)。 Y=fft(y)。 subplot(323)。plot(y)。title(?IIR 濾波后信號的波形 ?)。 ylabel(?amplitude?)。xlabel(?time(ms)?)。 subplot(324)。 plot(abs(y))。axis([0,6144,0,1])。 湖北師范學院 2020屆教育信息與技術學院學士學位論文（設計） 12 title(?IIR 濾波后信號的頻譜 ?)。 ylabel(?amplitude?)。xlabel(?frequency(hz)?)。 z=fftfilt(b,x)。 Z=fft(z)。 subplot(325)。 plot(z)。 title(?FIR 濾波后信號的波形 ?)。 ylabel(?amplitude?)。xlabel(?time(ms)?)。 subplot(326)。 plot(abs(z))。axis([0,6144,0,1])。 title(?FIR 濾波后信號的頻譜 ?)。 ylabel(?amplitude?)。xlabel(?frequency(hz)?)。 得到結(jié)果如圖 43 所示： 圖 43 N階高通濾波器 湖北師范學院 2020屆教育信息與技術學院學士學位論文（設計） 13 (4)N 階低通濾波器 在這里，同樣以 5 階為例，其中 wc 為其 3dB 邊緣頻率，程序設計如下： x=wavread(39。D:\Backup\我的文檔 \My Files\39。)。 sound(x)。N=5。wc=。 [b,a]=butter(N,wc)。 X=fft(x)。 subplot(321)。plot(x)。 title(39。濾波前信號的波形 39。)。 ylabel(39。amplitude39。)。xlabel(39。time(ms)39。)。 subplot(322)。plot(abs(x))。axis([0,6144,0,1])。 title(39。濾波前信號的頻譜 39。)。 ylabel(39。amplitude39。)。xlabel(39。frequency(hz)39。)。 y=filter(b,a,x)。 Y=fft(y)。 subplot(323)。plot(y)。 title(39。IIR 濾波后信號的波形 39。)。 ylabel(39。amplitude39。)。xlabel(39。time(ms)39。)。 subplot(324)。plot(abs(y))。axis([0,6144,0,1])。 title(39。IIR 濾波后信號的頻譜 39。)。 ylabel(39。amplitude39。)。xlabel(39。frequency(hz)39。)。 z=fftfilt(b,x)。 Z=fft(z)。subplot(325)。plot(z)。 title(39。FIR 濾波后信號的波形 39。)。 ylabel(39。amplitude39。)。xlabel(39。time(ms)39。)。 subplot(326)。plot(abs(z))。axis([0,6144,0,1])。 title(39。FIR 濾波后信號的頻譜 39。)。 ylabel(39。amplitude39。)。xlabel(39。frequency(hz)39。)。 得到結(jié)果如圖 44 所示： 湖北師范學院 2020屆教育信息與技術學院學士學位論文（設計） 14 圖 44 N階低通濾波器 (5)2N 階帶通濾波器 2N 階帶通濾波器的設計（在這里，以 10 階為例，其中 wc 為其 3dB 邊緣頻率， wc=[w1,w2],w1 wc w2），程序設計如下： x=wavread(?D:\Backup\我的文檔 \My Files\?)。 sound(x)。 N=5。wc=[,]。 [b,a]=butter(N,wc)。 X=fft(x)。 subplot(321)。plot(x)。title(?濾波前 信號的波形 ?)。 ylabel(?amplitude?)。xlabel(?time(ms)?)。 subplot(322)。plot(abs(x))。axis([0,6144,0,1])。title(?濾波前信號的頻譜 ?)。 湖北師范學院 2020屆教育信息與技術學院學士學位論文（設計） 15 ylabel(?amplitude?)。xlabel(?frequency(hz)?)。 y=filter(b,a,x)。Y=fft(y)。 subplot(323)。plot(y)。title(?IIR 濾波后信號的波形 ?)。 ylabel(?amplitude?)。xlabel(?time(ms)?)。 subplot(324)。plot(abs(y))。axis([0,6144,0,1])。title(?IIR 濾波后信號的頻譜 ?)。 ylabel(?amplitude?)。xlabel(?frequency(hz)?)。 z=fftfilt(b,x)。 Z=fft(z)。 subplot(325)。plot(z)。title(?FIR 濾波后信號的波形 ?)。 ylabel(?amplitude?)。xlabel(?time(ms)?)。 subplot(326)。plot(abs(z))。axis([0,6144,0,1])。title(?FIR 濾波后信號的頻譜 ?)。 ylabel(?amplitude?)。xlabel(?frequency(hz)?)。 得到結(jié)果如圖 45 所示： 圖 45 帶通濾波器 (6) 2N 階帶阻濾波器 2N 階帶阻濾波器的設計（在這里，以 10 階為例，其中 wc 為其 3dB 邊緣頻率， wc=[w1,w2],w1 wc w2），程序設計如下： 湖北師范學院 2020屆教育信息與技術學院學士學位論文（設計） 16 x=wavread(39。D:\Backup\我的文檔 \My Files\39。)。 sound(x)。 N=5。wc=[,]。 [b,a]=butter(N,wc,39。stop39。)。 X=fft(x)。 subplot(321)。plot(x)。 title(39。濾波前信號的波形 39。)。 ylabel(39。amplitude39。)。xlabel(39。time(ms)39。)。 subplot(322)。plot(abs(x))。axis([0,6144,0,1])。 tit