【文章內容簡介】 log(abs(Hk)))。 xlabel(39。頻率（ Hz） 39。)。ylabel(39。幅度（ dB）39。) axis([0,160,5]) 給被污染的信號加帶通數(shù)字濾波器 MATLAB 程序： wp=2*pi*[100,500]。ws=2*pi*[10,800]。Rp=1。As=20。 [N1,wp1]=cheb1ord(wp,ws,Rp,As,39。s39。)。 [B,A]=cheby1(N1,Rp,wp1,39。s39。)。 fk=0:2020/512:2020。wk=2*pi*fk。 Hk=freqs(B,A,wk)。 figure(4)。 plot(fk/1000,20*log(abs(Hk)))。 xlabel(39。頻率（ Hz） 39。)。ylabel(39。幅度（ dB）39。) title(39。帶通濾波器的頻譜分析 39。)。 axis([0,160,5])。 [Bz,Az]=impinvar(B,A,Fs)。 M=filter(Bz,Az,y)。 sound(M)。 %回放污染語音 figure(5)。 subplot(2,1,1)。plot(M)。 title(39。濾波后被污染的信號波形 39。)。 M1=fft(M,2020)。 subplot(2,1,2)。 plot(f(1:1000),abs(M1(1:1000)))。 title(39。濾波后被污染的信號頻譜 39。)。 一語音信號的數(shù)字濾波處理 第 23 頁 運行結果： 一語音信號的數(shù)字濾波處理 第 24 頁 設計數(shù)字高通濾波器 設計指標： 先將高通濾波器的指標裝換為對應的低通數(shù)字濾波氣的指標，設計相應的低通系統(tǒng)函數(shù)，對低通系統(tǒng)函數(shù)進行頻率變換，得到高通數(shù)字濾波器，本次取高通濾波器通帶邊界頻率 350Hz,阻帶邊界截止頻率150Hz MATLAB 程序： wp=2*pi*350。 ws=2*pi*150。 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)。 一語音信號的數(shù)字濾波處理 第 25 頁 figure(4)。 freqz(num,den,128,fs)。 axis([150,400,60,0]) grid on。 給被污染的信號加高通數(shù)字濾波器 wp=2*pi*350。 ws=2*pi*150。 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([150,400,60,0]) grid on。 M=filter(num,den,y)。 sound(M,FS,bits)。 %回放污染語音 figure(5)。 subplot(2,1,1)。plot(M)。 title(39。濾波后被污染的信號波形 39。)。 M1=fft(M,2020)。 subplot(2,1,2)。 plot(f(1:1000),abs(M1(1:1000)))。 axis([0,1400,0,150]) title(39。濾波后被污染的信號頻譜 39。)。 程序運行結果： 一語音信號的數(shù)字濾波處理 第 26 頁 一語音信號的數(shù)字濾波處理 第 27 頁 采用直接法設計 FIR 數(shù)字濾波器 設計 Blackman 窗 低通： 設計指標： 選取低通濾波器指標，通帶邊界截止頻率 500Hz，阻帶邊界截止頻率 1000Hz Matlab 程序： 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(4)。 subplot(2,1,1)。 plot([bz,az])。 axis([0,40,])。 title(39。h(n)波形 39。)。 subplot(2,1,2)。 plot(w/pi*FS,20*log(abs(h)))。 axis([0,3500,400,2])。 title(39。低通 損耗函數(shù)曲線 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)。 一語音信號的數(shù)字濾波處理 第 28 頁 figure(4)。 subplot(2,1,1)。 plot([bz,az])。 axis([0,40,])。 title(39。h(n)波形 39。)。 subplot(2,1,2)。 plot(w/pi*FS,20*log(abs(h)))。 axis([0,3500,400,2])。 title(39。低通 損耗函數(shù)曲線 39。)。 M=filter(bz,az,y)。 sound(M,FS,bits)。 %播放加窗后語音 figure(5)。 subplot(2,1,1)。plot(M)。 title(39。濾波后被污染的信號波形 39。)。 M1=fft(M,2020)。 subplot(2,1,2)。 plot(f(1:500),abs(M1(1:500)))。 title(39。濾波后被污染的信號頻譜 39。)。 程序運行結果： 一語音信號的數(shù)字濾波處理 第 29 頁 一語音信號的數(shù)字濾波處理 第 30 頁 設計 Blackman 窗 帶通 設計指標： 選取帶 通指標參數(shù)：通帶截止頻率 100Hz， 500Hz，阻帶截止頻率 30Hz,800Hz Matlab 程序： wlp=100*2*pi/FS。wup=500*2*pi/FS。 wls=30*2*pi/FS。wus=1500*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*FS,20*log(abs(h)))。 title(39。帶通 損耗函數(shù)曲線 39。)。 給被污染的信號加窗： wlp=100*2*pi/FS。wup=500*2*pi/FS。 wls=30*2*pi/FS。wus=1500*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*FS,20*log(abs(h)))。 title(39。帶通 損耗函數(shù)曲線 39。)。 %axis([0,1,600,50])。 一語音信號的數(shù)字濾波處理 第 31 頁 M=filter(bz,az,y)。 sound(M,FS,bits)。 %播放加窗后語音 figure(5)。 subplot(2,1,1)。plot(M)。 title(39。濾波后被污染的信號波形 39。)。 M1=fft(M,2020)。 subplot(2,1,2)。 plot(f(1:1000),abs(M1(1:1000)))。 axis([0,2020,0,150])。 title(39。濾波后被污染的信號頻譜 39。)。 程序運行結果： 一語音信號的數(shù)字濾波處理 第 32 頁 一語音信號的數(shù)字濾波處理 第 33 頁 設計 Blackman 窗 高通： 設計指標： 選取帶通指標參數(shù)：通帶截止頻率 600Hz 阻帶截止頻率 400Hz Matlab 程序： [h,w]=freqz(bz,az)。 figure(4)。 subplot(2,1,1)。 plot([bz,az])。 axis([0,70,1])。 title(39。h(n)波形 39。)。 subplot(2,1,2)。 plot(w/pi*FS,20*log(abs(h)))。 axis([300,1400,300,100])。 title(39。高通 損耗函 fp=600。fs=400。 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))。 數(shù)曲線 39。)。 給被污染的信號加窗： fp=600。fs=400。 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)。 一語