【正文】 Fs1=1。 wap=2*tan(wp/2)。 was=2*tan(ws/2)。 [N,wc]=buttord(wap,was,rp,rs,39。s39。)。 [B,A]=butter(N,wc,39。s39。)。 [Bz,Az]=bilinear(B,A,Fs1)。 figure(4)。 [h,w]=freqz(Bz,Az,512,Fs1*22050)。 plot(w,abs(h))。 title(39。巴特沃斯帶通濾波器39。)。 xlabel(39。頻率HZ39。)。ylabel(39。耗損（dB）39。)。 grid on。 yd=filter(Bz,Az,x1)。 figure(5)。 subplot(2,1,1)。plot(yd)。 title(39。濾波后的信號波形39。)。ydd=fft(yd,1200)。 subplot(2,1,2)。plot(f(1:600),abs(ydd(1:600)))。title(39。濾波后的信號頻譜39。)。sound(yd,FS,bits)。附錄D：blackman高通濾波器程序Fs=22050。[x,FS,bits]=wavread(39。C:\WINDOWS\Media\39。)。x=x(:,1)。figure(1)。subplot(2,1,1)。plot(x)。sound(x,FS,bits)。 title(39。語音信號時域波形圖39。)y=fft(x,3260)。f=(FS/1630)*[1:1630]。 subplot(2,1,2)。plot(f(1:1630),abs(y(1:1630)))。title(39。語音信號頻譜圖39。)。t=0:length(x)1。zs0=*cos(2*pi*100*t/22050)。figure(2)。subplot(2,1,1)plot(zs0)。title(39。噪聲信號波形39。)。zs1=fft(zs0,1200)。 subplot(2,1,2)plot(f(1:600),abs(zs1(1:600)))。 title(39。噪聲信號頻譜39。)。x1=x+zs039。sound(x1,FS,bits)y1=fft(x1,1200)。figure(3)。subplot(2,1,1)。plot(x1)。title(39。加入噪聲后的信號波形39。)。subplot(2,1,2)。plot(f(1:600),abs(y1(1:600)))。title(39。加入噪聲后的信號頻譜39。)。fp=600,fc=400。 wp=2*pi*fp/FS。ws=2*pi*fc/FS。Bt=wpws。 N0=ceil(11*pi/Bt)。 N=N0+mod(N0+1,2)。wc=(wp+ws)/2/pi。 hn=fir1(N1,wc,39。high39。,blackman(N))。 X=conv(hn,x)。 sound(X,FS,bits)。X1=fft(X,1200)。figure(4)。subplot(211)。plot(X)。title(39。濾波后的信號波形39。)。subplot(212)。plot(f(1:600),abs(X1(1:600)))。 title(39。濾波后的信號頻譜39。)sound(yd,FS,bits)。附錄E：blackman低通濾波器程序Fs=22050。[x,FS,bits]=wavread(39。C:\WINDOWS\Media\39。)。x=x(:,1)。figure(1)。subplot(2,1,1)。plot(x)。sound(x,FS,bits)。 title(39。語音信號時域波形圖39。)y=fft(x,3260)。f=(FS/1630)*[1:1630]。 subplot(2,1,2)。plot(f(1:1630),abs(y(1:1630)))。title(39。語音信號頻譜圖39。)。t=0:length(x)1。zs=*cos(2*pi*10000*t/22050)。zs0=*cos(2*pi*10000*t/22050000)。figure(2)。subplot(2,1,1)plot(zs0) 。title(39。噪聲信號波形39。)。zs1=fft(zs,1200)。subplot(2,1,2)plot(f(1:600),abs(zs1(1:600)))。title(39。噪聲信號頻譜39。)。x1=x+zs39。sound(x1,FS,bits)y1=fft(x1,1200)。figure(3)。subplot(2,1,1)。plot(x1)。title(39。加入噪聲后的信號波形39。)。subplot(2,1,2)。plot(f(1:600),abs(y1(1:600)))。title(39。加入噪聲后的信號頻譜39。)。fp=3000,fc=3500。 wp=2*pi*fp/FS。ws=2*pi*fc/FS。Bt=wswp。 N0=ceil(11*pi/Bt)。 N=N0+mod(N0+1,2)。wc=(wp+ws)/2/pi。 hn=fir1(N1,wc,blackman(N))。 X=conv(hn,x)。 sound(X,FS,bits)。X1=fft(X,1200)。figure(4)。subplot(2,1,1)。plot(X)。title(39。濾波后的信號波形39。)。subplot(2,1,2)。plot(f(1:600),abs(X1(1:600)))。 title(39。濾波后的信號頻譜39。)sound(yd,FS,bits)。附錄F：blackman帶通濾波器程序Fs=22050。[x,FS,bits]=wavread(39。C:\WINDOWS\Media\39。)。x=x(:,1)。figure(1)。subplot(2,1,1)。plot(x)。sound(x,FS,bits)title(39。語音信號時域波形圖39。)y=fft(x,3260)。f=(FS/1630)*[1:1630]。 subplot(2,1,2)。plot(f(1:1630),abs(y(1:1630)))。title(39。語音信號頻譜圖39。)。t=0:length(x)1。zs0=*cos(2*pi*100*t/22050)。figure(2)。subplot(2,1,1)plot(zs0)。title(39。噪聲信號波形39。)。zs1=fft(zs0,1200)。subplot(2,1,2)plot(f(1:600),abs(zs1(1:600)))。title(39。噪聲信號頻譜39。)。x1=x+zs039。sound(x1,FS,bits)y1=fft(x1,1200)。figure(3)。subplot(2,1,1)。plot(x1)。title(39。加入噪聲后的信號波形39。)。subplot(2,1,2)。plot(f(1:600),abs(y1(1:600)))。title(39。加入噪聲后的信號頻譜39。)。 wlp=600*2*pi/FS,wup=6000*2*pi/FS。wls=400*2*pi/FS。wus=7000*2*pi/FS。 Bt=wlpwls。 N0=ceil(11*pi/Bt)。 N=N0+mod(N0+1,2)。wc=[(wls+wlp)/2/pi,(wus+wup)/2/pi]。 hn=fir1(N1,wc,blackman(N))。 X=conv(hn,x)。 sound(X,FS,bits)。X1=fft(X,1200)。figure(4)。subplot(2,1,1)。plot(X)。title(39。濾波后的信號波形39。)。subplot(2,1,2)。plot(f(1:600),abs(X1(1:600)))。 title(39。濾波后的信號頻譜39。)sound(yd,FS,bits)。第 43 頁