【文章內容簡介】 fft(M,2000)。 subplot(2,1,2)。plot(f(1:1000),abs(M1(1:1000)))。title(39。濾波后被污染的信號頻譜39。)。： 設計數(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)。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,2000)。 subplot(2,1,2)。plot(f(1:1000),abs(M1(1:1000)))。axis([0,1400,0,150])title(39。濾波后被污染的信號頻譜39。)。 程序運行結果：： 設計指標：選取低通濾波器指標，通帶邊界截止頻率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)。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,2000)。 subplot(2,1,2)。plot(f(1:500),abs(M1(1:500)))。title(39。濾波后被污染的信號頻譜39。)。 程序運行結果： 設計指標：選取帶通指標參數(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])。M=filter(bz,az,y)。sound(M,FS,bits)。 %播放加窗后語音 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,2000,0,150])。title(39。濾波后被污染的信號頻譜39。)。 程序運行結果：： 設計指標：選取帶通指標參數(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)。 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。高通損耗函數(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,2000)。 subplot(2,1,2)。plot(f(1:1000),abs(M1(1:1000)))。axis([0,1400,0,400])。title(39。濾波后被污染的信號頻譜39。)。 程序運行結果：5 針對不同噪聲信號濾波及分析切比雪夫低通濾波器可以實現(xiàn)對高頻噪聲進行濾除，將混合信號恢復至原始信號，為驗證此結論的正確性，本次課程設計中我們添加噪聲信號，高頻噪聲zs=*cos(2*pi*10000*t/22050)于原始信號上，并使用自己設計的切比雪夫濾波器對其進行處理，仿真結果證明所設計濾波器可實現(xiàn)要求。 切比雪夫低通濾波器分析圖 切比雪夫高通濾波器可以實現(xiàn)對低頻噪聲進行濾除，將混合信號恢復至原始信號，為驗證此結論的正確性，本次課程設計中我們添加噪聲信號，低頻噪聲zs0=*cos(2*pi*100*t/22050)于原始信號上，并使用自己設計的切比雪夫濾波器對