【正文】 第 43 頁。濾波后的信號頻譜39。plot(f(1:600),abs(X1(1:600)))。)。title(39。subplot(2,1,1)。X1=fft(X,1200)。 X=conv(hn,x)。wc=[(wls+wlp)/2/pi,(wus+wup)/2/pi]。 N0=ceil(11*pi/Bt)。wus=7000*2*pi/FS。 wlp=600*2*pi/FS,wup=6000*2*pi/FS。加入噪聲后的信號頻譜39。plot(f(1:600),abs(y1(1:600)))。)。title(39。subplot(2,1,1)。sound(x1,FS,bits)y1=fft(x1,1200)。)。title(39。zs1=fft(zs0,1200)。噪聲信號波形39。subplot(2,1,1)plot(zs0)。zs0=*cos(2*pi*100*t/22050)。)。title(39。 subplot(2,1,2)。)y=fft(x,3260)。sound(x,FS,bits)title(39。subplot(2,1,1)。x=x(:,1)。C:\WINDOWS\Media\39。附錄F：blackman帶通濾波器程序Fs=22050。濾波后的信號頻譜39。plot(f(1:600),abs(X1(1:600)))。)。title(39。subplot(2,1,1)。X1=fft(X,1200)。 X=conv(hn,x)。wc=(wp+ws)/2/pi。 N0=ceil(11*pi/Bt)。ws=2*pi*fc/FS。fp=3000,fc=3500。加入噪聲后的信號頻譜39。plot(f(1:600),abs(y1(1:600)))。)。title(39。subplot(2,1,1)。sound(x1,FS,bits)y1=fft(x1,1200)。)。title(39。zs1=fft(zs,1200)。噪聲信號波形39。subplot(2,1,1)plot(zs0) 。zs0=*cos(2*pi*10000*t/22050000)。t=0:length(x)1。語音信號頻譜圖39。plot(f(1:1630),abs(y(1:1630)))。f=(FS/1630)*[1:1630]。語音信號時域波形圖39。sound(x,FS,bits)。subplot(2,1,1)。x=x(:,1)。C:\WINDOWS\Media\39。附錄E：blackman低通濾波器程序Fs=22050。濾波后的信號頻譜39。plot(f(1:600),abs(X1(1:600)))。)。title(39。subplot(211)。X1=fft(X,1200)。 X=conv(hn,x)。high39。wc=(wp+ws)/2/pi。 N0=ceil(11*pi/Bt)。ws=2*pi*fc/FS。fp=600,fc=400。加入噪聲后的信號頻譜39。plot(f(1:600),abs(y1(1:600)))。)。title(39。subplot(2,1,1)。sound(x1,FS,bits)y1=fft(x1,1200)。)。 title(39。zs1=fft(zs0,1200)。噪聲信號波形39。subplot(2,1,1)plot(zs0)。zs0=*cos(2*pi*100*t/22050)。)。title(39。 subplot(2,1,2)。)y=fft(x,3260)。 title(39。plot(x)。figure(1)。)。[x,FS,bits]=wavread(39。sound(yd,FS,bits)。濾波后的信號頻譜39。plot(f(1:600),abs(ydd(1:600)))。ydd=fft(yd,1200)。濾波后的信號波形39。plot(yd)。 figure(5)。 grid on。耗損（dB）39。)。 xlabel(39。巴特沃斯帶通濾波器39。 plot(w,abs(h))。 figure(4)。)。 [B,A]=butter(N,wc,39。s39。 was=2*tan(ws/2)。Fs1=1。 ws=2*pi*fs/Fs。rs=10。Fs=22050。 fp=[600,6000]。加入噪聲后的信號頻譜39。 plot(f(1:600),abs(y1(1:600)))。)。 title(39。 subplot(2,1,1)。 y1=fft(x1,1200)。 x1=x+zs039。噪聲信號頻譜39。plot(f(1:600),abs(zs1(1:600)))。 sound(zs,FS,bits)。)。 title(39。 subplot(2,1,1)。 zs0=*cos(2*pi*100*t/22050)。)。 title(39。 subplot(2,1,2)。 y=fft(x,3260)。語音信號時域波形圖39。 sound(x,FS,bits)。 subplot(2,1,1)。 x=x(:,1)。C:\WINDOWS\Media\39。附錄C：巴特沃斯帶通濾波器程序Fs=22050。)。 title(39。 subplot(2,1,2)。)。 title(39。 subplot(2,1,1)。 yd=filter(Bz,Az,x1)。)。ylabel(39。頻率HZ39。)。 title(39。 [h,w]=freqz(Bz,Az,512,Fs1*22050)。 [Bz,Az]=bilinear(B,A,Fs1)。s39。high39。)。 [N,wc]=buttord(wap,was,rp,rs,39。 wap=2*tan(wp/2)。 T=1。 wp=2*pi*fp/Fs。 rp=1。fs=400。)。 title(39。 subplot(2,1,2)。加入噪聲后的信號波形39。plot(x1)。 figure(3)。 sound(x1,FS,bits)。)。 title(39。 subplot(2,1,2)。 zs1=fft(zs0,1200)。噪聲信號波形39。 plot(zs0)。 figure(2)。 t=0:length(x)1。語音信號頻譜圖39。 plot(f(1:1630),abs(y(1:1630)))。 f=(FS/1630)*[1:1630]。語音信號時域波形圖39。sound(x,FS,bits)。 subplot(2,1,1)。 x=x(:,1)。C:\WINDOWS\Media\39。附錄B：巴特沃斯高通濾波器程序Fs=22050。)。 title(39。subplot(2,1,2)。)。title(39。subplot(2,1,1)。yd=filter(Bz,Az,x1)。)。ylabel(39。頻率(HZ)39。)。title(39。[h,w]=freqz(Bz,Az,512,Fs1*22050)。[Bz,Az]=bilinear(B,A,Fs1)。s39。)。[N,wc]=buttord(wap,was,rp,rs,39。wap=2*tan(wp/2)。ws=2*pi*fs/FS。rs=10。 Fs=22050。)。title(39。subplot(2,1,2)。加入噪聲后的信號波形39。plot(x1)。figure(3)。sound(x1,FS,bits)。)。 title(39。 sound(x,FS,bits)。)。title(39。figure(2)。zs=*cos(2*pi*10000*t/22050)。)。title(39。 subplot(2,1,2)。)y=fft(x,3260)。title(39。subplot(2,1,1)plot(x)。x=x(:,1)。C:\WINDOWS\Media\39。（5）《數(shù)字信號處理》，北京：科學(xué)出版社（6）《數(shù)字信號處理——理論、算法與實現(xiàn)（第二版）》胡廣書編著，北京：電子工業(yè)出版社附錄A：巴特沃斯低通濾波器程序Fs=22050。 這次課程設(shè)計的每個實驗細節(jié)