【正文】 =freqz(hn,1,15000)。 figure(3)freqz(hn,1,15000)。xlabel(39。頻率(HZ)39。)。ylabel(39。幅度(db)39。)。title(39。高通漢寧窗數(shù)字濾波器幅頻特性39。)。 x=fftfilt(hn,y2)。X=fft(x,15000)。figure(4)。subplot(2,1,1)。plot(abs(Y2(1:7500)))。 title(39。濾波前信號頻譜39。)。subplot(2,1,2)。plot(abs(X(1:7500)))。title(39。濾波后信號頻譜39。)。 figure(5)。subplot(2,1,1)。plot(y2)。 語音信號的數(shù)字濾波處理第 30 頁 共 36 頁title(39。濾波前信號波形39。)。subplot(2,1,2)。plot(x)。title(39。濾波后信號波形39。)。axis([0,18000,])。sound(x,fs,bits)。5. 切比雪夫Ⅱ型帶通濾波器[y,fs,bits]=wavread(39。e:\39。,[1 16384])。f=fs*(0:7500)/15000。t=0:1/16384:(size(y)1)/16384。 %將所加噪聲信號的點數(shù)調(diào)整到與原始信號相同Au=。d=[Au*cos(2*pi*5000*t)]39。y2=y+d。sound(y2,16384)。 %播放加噪聲后的語音信號Y=fft(y,15000)。f=fs*(0:7500)/15000。Y2=fft(y2,15000)。 figure(1)subplot(2,1,1)。plot(y)。title(39。原始語音信號39。)。subplot(2,1,2)。plot(y2)。title(39。加噪后的信號39。)。xlabel(39。時間 n39。)。ylabel(39。幅值39。)。 figure(2)subplot(2,1,1)。plot(abs(Y(1:7500)))。title(39。原始語音信號頻譜39。)。xlabel(39。Hz39。)。ylabel(39。幅值39。)。subplot(2,1,2)。plot(abs(Y2(1:7500)))。title(39。加噪后的信號頻譜39。)。axis([0 8000 0 100])。xlabel(39。Hz39。)。ylabel(39。幅值39。)。 %濾波器設計fss=15000wp=[600*2*pi/fss,4000*2*pi/fss]。Rp=1。Rs=100。ws=[300*2*pi/fss,7000*2*pi/fss]。wo=sqrt((600*2*pi/fss)*(4000*2*pi/fss))。Bw=4000*2*pi/fss600*2*pi/fss。 語音信號的數(shù)字濾波處理第 31 頁 共 36 頁[n,wn]=cheb2ord(wp,ws,1,100,39。s39。)。[z,p,k]=cheb2ap(n,99)。[b0,a0]=zp2tf(z,p,k)。[bt,at]=lp2bp(b0,a0,wo,Bw)。[bz,az]=bilinear(bt,at,1)。[h,w]=freqz(bz,az)。f=w*fss/(2*pi)。 figure(3)plot(f,20*log10(abs(h)))。xlabel(39。頻率(HZ)39。)。ylabel(39。幅度(db)39。)。title(39。帶通切比雪夫數(shù)字濾波器幅頻特性39。)。axis([0,12022,150,10])。 x=filter(bz,az,y2)。X=fft(x,15000)。figure(4)。subplot(2,1,1)。plot(abs(Y2(1:7500)))。 title(39。濾波前信號頻譜39。)。subplot(2,1,2)。plot(abs(X(1:7500)))。 title(39。濾波后信號頻譜39。)。 figure(5)。subplot(2,1,1)。plot(y2)。title(39。濾波前信號波形39。)。subplot(2,1,2)。plot(x)。 title(39。濾波后信號波形39。)。axis([0,18000,])。sound(x,fs,bits)。 6. 切比雪夫Ⅱ型低通濾波器[y,fs,bits]=wavread(39。e:\39。,[1 16384])。f=fs*(0:7500)/15000。t=0:1/16384:(size(y)1)/16384。 %將所加噪聲信號的點數(shù)調(diào)整到與原始信號相同Au=。d=[Au*cos(2*pi*5000*t)]39。y2=y+d。sound(y2,16384)。 %播放加噪聲后的語音信號 Y=fft(y,15000)。f=fs*(0:7500)/15000。Y2=fft(y2,15000)。 語音信號的數(shù)字濾波處理第 32 頁 共 36 頁figure(1)subplot(2,1,1)。plot(y)。title(39。原始語音信號39。)。subplot(2,1,2)。plot(y2)。title(39。加噪后的信號39。)。xlabel(39。時間 n39。)。ylabel(39。幅值39。)。 figure(2)subplot(2,1,1)。plot(abs(Y(1:7500)))。title(39。原始語音信號頻譜39。)。xlabel(39。Hz39。)。ylabel(39。幅值39。)。subplot(2,1,2)。plot(abs(Y2(1:7500)))。title(39。加噪后的信號頻譜39。)。axis([0 8000 0 100])。xlabel(39。Hz39。)。ylabel(39。幅值39。)。%濾波器設計 wp=4000*2*pi/fs。Rp=1。Rs=100。ws=6000*2*pi/fs。[n,wn]=cheb2ord(wp,ws,1,100,39。s39。)。[z,p,k]=cheb2ap(n,99)。[b0,a0]=zp2tf(z,p,k)。[bt,at]=lp2lp(b0,a0,ws)。[bz,az]=bilinear(bt,at,1)。[h,w]=freqz(bz,az)。f=w*fs/(2*pi)。 figure(3)plot(f,20*log10(abs(h)))。 title(39。低通切比雪夫數(shù)字濾波器幅頻特性39。)。xlabel(39。頻率(HZ)39。)。ylabel(39。幅度(db)39。)。axis([0 13000 150 5])。 x=filter(bz,az,y2)。X=fft(x,15000)。figure(4)。subplot(2,1,1)。plot(abs(Y2(1:7500)))。 title(39。濾波前信號頻譜39。)。 語音信號的數(shù)字濾波處理第 33 頁 共 36 頁subplot(2,1,2)。plot(abs(X(1:7500)))。title(39。濾波后信號頻譜39。)。 figure(5)。subplot(2,1,1)。plot(y2)。title(39。濾波前信號波形39。)。subplot(2,1,2)。plot(x)。 title(39。濾波后信號波形39。)。axis([0,18000,])。sound(x,fs,bits)。7. 切比雪夫Ⅱ型高通濾波器[y,fs,bits]=wavread(39。e:\39。,[1 16384])。f=fs*(0:7500)/15000。t=0:1/16384:(size(y)1)/16384。 %將所加噪聲信號的點數(shù)調(diào)整到與原始信號相同Au=。d=[Au*cos(2*pi*5000*t)]39。y2=y+d。sound(y2,16384)。 %播放加噪聲后的語音信號 Y=fft(y,15000)。f=fs*(0:7500)/15000。Y2=fft(y2,15000)。 figure(1)subplot(2,1,1)。plot(y)。title(39。原始語音信號39。)。subplot(2,1,2)。plot(y2)。title(39。加噪后的信號39。)。xlabel(39。時間 n39。)。ylabel(39。幅值39。)。 figure(2)subplot(2,1,1)。plot(abs(Y(1:7500)))。title(39。原始語音信號頻譜39。)。xlabel(39。Hz39。)。ylabel(39。幅值39。)。subplot(2,1,2)。plot(abs(Y2(1:7500)))。title(39。加噪后的信號頻譜39。)。axis([0 8000 0 100])。xlabel(39。Hz39。)。ylabel(39。幅值39。)。%濾波器設計 wp=8000*2*pi/fs。Rp=1。Rs=100。ws=6000*2*pi/fs。[n,wn]=cheb2ord(wp,ws,1,100,39。s39。)。 語音信號的數(shù)字濾波處理第 34 頁 共 36 頁[z,p,k]=cheb2ap(n,99)。[b0,a0]=zp2tf(z,p,k)。[bt,at]=lp2hp(b0,a0,ws)。[bz,az]=bilinear(bt,at,1)。[h,w]=freqz(bz,az)。f=w*fs/(2*pi)。 figure(3)plot(f,20*log10(abs(h)))。xlabel(39。頻率(HZ)39。)。ylabel(39。幅度(db)39。)。title(39。高通切比雪夫數(shù)字濾波器幅頻特性39。)。axis([0,12022,150,10])。 x=filter(bz,az,y2)。X=fft(x,15000)。figure(4)。subplot(2,1,1)。plot(abs(Y2(1:7500)))。 title(39。濾波前信號頻譜39。)。subplot(2,1,2)。plot(abs(X(1:7500)))。 title(39。濾波后信號頻譜39。)。 figure(5)。subplot(2,1,1)。plot(y2)。title(39。濾波前信號波形39。)。subplot(2,1,2)。plot(x)。 title(39。濾波后信號波形39。)。axis([0,18000,])。sound(x,fs,bits)。