【正文】 +zs39。sound(x,FS,bits)。 %回放污染語音 figure(3)。 subplot(2,1,1)。 plot(y) title(39。被噪聲污染的信號波形39。)。y1=fft(y,2000)。 subplot(2,1,2)。plot(f(1:1000),abs(y1(1:1000)))。axis([0,200,0,200])。title(39。被噪聲污染的信號頻譜39。)。wp=2*pi*50。ws=2*pi*25。rp=3。rs=10。fs=1000。[n,wc]=cheb2ord(wp,ws,rp,rs,39。s39。)。[z,p,k]=cheb2ap(n,rs)。[a,b,c,d]=zp2ss(z,p,k)。[at1,bt1,ct1,dt1]=lp2hp(a,b,c,d,wc)。[at2,bt2,ct2,dt2]=bilinear(at1,bt1,ct1,dt1,fs)。%雙線性變換法[num,den]=ss2tf(at2,bt2,ct2,dt2)。figure(4)。freqz(num,den,128,fs)。axis([25,200,30,0])grid on。M=filter(num,den,y)。sound(M)。 %回放污染語音 figure(5)。subplot(2,1,1)。plot(M)。title(39。濾波后被污染的信號波形39。)。M1=fft(M,2000)。 subplot(2,1,2)。plot(f(1:1000),abs(M1(1:1000)))。axis([0,500,0,40])title(39。濾波后被污染的信號頻譜39。)。附錄D Blackman窗低通濾波器仿真程序：clf。[x,FS,bits]=wavread(39。d:\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。)。y=x+zs39。%sound(x,FS,bits)。 %回放污染語音 figure(3)。 subplot(2,1,1)。 plot(y) title(39。被噪聲污染的信號波形39。)。y1=fft(y,2000)。 subplot(2,1,2)。plot(f(1:1000),abs(y1(1:1000)))。title(39。被噪聲污染的信號頻譜39。)。fp=500。fs=1000。 wp=2*pi*fp/FS。 ws=2*pi*fs/FS。 Bt=wswp。 N0=ceil(11*pi/Bt)。 N=N0+mod(N0+1,2)。 wc=(wp+ws)/2/pi。 [bz,az]=fir1(N1,wc,blackman(N))。 [h,w]=freqz(bz,az)。figure()。subplot(2,1,1)。 plot([bz,az])。axis([0,40,])。title(39。h(n)波形39。)。subplot(2,1,2)。 plot(w/pi,20*log(abs(h)))。axis([0,1,400,2])。title(39。低通損耗函數(shù)曲線39。)。M=conv([bz,az],y)。%sound(M)。 %播放加窗后語音 figure(5)。subplot(2,1,1)。plot(M)。title(39。濾波后被污染的信號波形39。)。M1=fft(M,2000)。 subplot(2,1,2)。plot(f(1:500),abs(M1(1:500)))。title(39。濾波后被污染的信號頻譜39。)。附錄E Blackman窗帶通濾波器仿真程序：clf。[x,FS,bits]=wavread(39。d:\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)))。axis([0,500,0,300])。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。)。y=x+zs39。%sound(x,FS,bits)。 %回放污染語音 figure(3)。 subplot(2,1,1)。 plot(y) title(39。被噪聲污染的信號波形39。)。y1=fft(y,2000)。 subplot(2,1,2)。plot(f(1:1000),abs(y1(1:1000)))。title(39。被噪聲污染的信號頻譜39。)。wlp=100*2*pi/FS。wup=500*2*pi/FS。 wls=30*2*pi/FS。wus=800*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]。 [bz,az]=fir1(N1,wc,blackman(N))。 [h,w]=freqz(bz,az)。figure(4)。subplot(2,1,1)。 plot([bz,az])。axis([100,200,])。title(39。h(n)波形39。)。subplot(2,1,2)。 plot(w/pi,20*log(abs(h)))。title(39。帶通損耗函數(shù)曲線39。)。axis([0,1,600,50])。M=conv([bz,az],y)。sound(M)。 %播放加窗后語音 figure(5)。subplot(2,1,1)。plot(M)。title(39。濾波后被污染的信號波形39。)。M1=fft(M,2000)。 subplot(2,1,2)。plot(f(1:1000),abs(M1(1:1000)))。axis([0,500,0,150])。title(39。濾波后被污染的信號頻譜39。)。附錄F Blackman窗高通濾波器仿真程序：clf。[x,FS,bits]=wavread(39。d:\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)))。axis([0,500,0,300])。title(39。語音信號頻譜圖39。)。 t=0:length(x)1。 zs=*cos(2*pi*100*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)))。axis([0,100,0,60])。title(39。噪聲信號頻譜39。)。y=x+zs39。%sound(x,FS,bits)。 %回放污染語音 figure(3)。 subplot(2,1,1)。 plot(y) title(39。被噪聲污染的信號波形39。)。y1=fft(y,2000)。 subplot(2,1,2)。plot(f(1:1000),abs(y1(1:1000)))。axis([0,500,0,150])。title(39。被噪聲污染的信號頻譜39。)。fp=50。fs=25。 wp=2*pi*fp/FS。 ws=2*pi*fs/FS。 Bt=wpws。 N0=ceil(11*pi/Bt)。 N=N0+mod(N0+1,2)。 wc=(wp+ws)/2/pi。 [bz,az]=fir1(N1,wc,39。high39。,blackman(N))。 [h,w]=freqz(bz,az)。figure(4)。subplot(2,1,1)。 plot([bz,az])。axis([400,500,1])。title(39。h(n)波形39。)。subplot(2,1,2)。 plot(w/pi,20*log(abs(h)))。axis([0,200,2])。title(39。高通損耗函數(shù)曲線39。)。M=conv([bz,az],y)。sound(M)。 %播放加窗后語音 figure(5)。subplot(2,1,1)。plot(M)。title(39。濾波后被污染的信號波形39。)。M1=fft(M,2000)。 subplot(2,1,2)。plot(f(1:1000),abs(M1(1:1000)))。axis([0,500,0,150])。title(39。濾波后被污染的信號頻譜39。)。