【文章內(nèi)容簡(jiǎn)介】 fft(M,2000)。 subplot(2,1,2)。plot(f(1:1000),abs(M1(1:1000)))。title(39。濾波后被污染的信號(hào)頻譜39。)。： 設(shè)計(jì)數(shù)字高通濾波器 設(shè)計(jì)指標(biāo)：先將高通濾波器的指標(biāo)裝換為對(duì)應(yīng)的低通數(shù)字濾波氣的指標(biāo)，設(shè)計(jì)相應(yīng)的低通系統(tǒng)函數(shù)，對(duì)低通系統(tǒng)函數(shù)進(jìn)行頻率變換，得到高通數(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。 給被污染的信號(hào)加高通數(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)。 %回放污染語(yǔ)音 figure(5)。subplot(2,1,1)。plot(M)。title(39。濾波后被污染的信號(hào)波形39。)。M1=fft(M,2000)。 subplot(2,1,2)。plot(f(1:1000),abs(M1(1:1000)))。axis([0,1400,0,150])title(39。濾波后被污染的信號(hào)頻譜39。)。 程序運(yùn)行結(jié)果：： 設(shè)計(jì)指標(biāo)：選取低通濾波器指標(biāo)，通帶邊界截止頻率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。)。 給被污染的信號(hào)加窗：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)。 %播放加窗后語(yǔ)音 figure(5)。subplot(2,1,1)。plot(M)。title(39。濾波后被污染的信號(hào)波形39。)。M1=fft(M,2000)。 subplot(2,1,2)。plot(f(1:500),abs(M1(1:500)))。title(39。濾波后被污染的信號(hào)頻譜39。)。 程序運(yùn)行結(jié)果： 設(shè)計(jì)指標(biāo)：選取帶通指標(biāo)參數(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。)。 給被污染的信號(hào)加窗：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)。 %播放加窗后語(yǔ)音 figure(5)。subplot(2,1,1)。plot(M)。title(39。濾波后被污染的信號(hào)波形39。)。M1=fft(M,2000)。 subplot(2,1,2)。plot(f(1:1000),abs(M1(1:1000)))。axis([0,2000,0,150])。title(39。濾波后被污染的信號(hào)頻譜39。)。 程序運(yùn)行結(jié)果：： 設(shè)計(jì)指標(biāo)：選取帶通指標(biāo)參數(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。)。 給被污染的信號(hào)加窗：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)。 %播放加窗后語(yǔ)音 figure(5)。subplot(2,1,1)。plot(M)。title(39。濾波后被污染的信號(hào)波形39。)。M1=fft(M,2000)。 subplot(2,1,2)。plot(f(1:1000),abs(M1(1:1000)))。axis([0,1400,0,400])。title(39。濾波后被污染的信號(hào)頻譜39。)。 程序運(yùn)行結(jié)果：5 針對(duì)不同噪聲信號(hào)濾波及分析切比雪夫低通濾波器可以實(shí)現(xiàn)對(duì)高頻噪聲進(jìn)行濾除，將混合信號(hào)恢復(fù)至原始信號(hào)，為驗(yàn)證此結(jié)論的正確性，本次課程設(shè)計(jì)中我們添加噪聲信號(hào)，高頻噪聲zs=*cos(2*pi*10000*t/22050)于原始信號(hào)上，并使用自己設(shè)計(jì)的切比雪夫?yàn)V波器對(duì)其進(jìn)行處理，仿真結(jié)果證明所設(shè)計(jì)濾波器可實(shí)現(xiàn)要求。 切比雪夫低通濾波器分析圖 切比雪夫高通濾波器可以實(shí)現(xiàn)對(duì)低頻噪聲進(jìn)行濾除，將混合信號(hào)恢復(fù)至原始信號(hào)，為驗(yàn)證此結(jié)論的正確性，本次課程設(shè)計(jì)中我們添加噪聲信號(hào)，低頻噪聲zs0=*cos(2*pi*100*t/22050)于原始信號(hào)上，并使用自己設(shè)計(jì)的切比雪夫?yàn)V波器對(duì)