【正文】 濾波后被污染的信號頻譜39。 subplot(2,1,2)。title(39。sound(M)。title(39。)。 plot([bz,az])。,blackman(N))。 N=N0+mod(N0+1,2)。 wp=2*pi*fp/FS。被噪聲污染的信號頻譜39。 subplot(2,1,2)。 plot(y) title(39。y=x+zs39。axis([0,100,0,60])。)。 figure(2)。)。plot(f(1:1630),abs(y(1:1630)))。語音信號時域波形圖39。 subplot(2,1,1)。d:\media\39。濾波后被污染的信號頻譜39。 subplot(2,1,2)。title(39。sound(M)。帶通損耗函數(shù)曲線39。)。 plot([bz,az])。 [bz,az]=fir1(N1,wc,blackman(N))。 Bt=wlpwls。wlp=100*2*pi/FS。plot(f(1:1000),abs(y1(1:1000)))。被噪聲污染的信號波形39。%sound(x,FS,bits)。title(39。)。 figure(2)。)。plot(f(1:1630),abs(y(1:1630)))。語音信號時域波形圖39。 subplot(2,1,1)。d:\media\39。濾波后被污染的信號頻譜39。M1=fft(M,2000)。plot(M)。M=conv([bz,az],y)。axis([0,1,400,2])。h(n)波形39。subplot(2,1,1)。 wc=(wp+ws)/2/pi。 ws=2*pi*fs/FS。)。 subplot(2,1,2)。 plot(y) title(39。y=x+zs39。plot(f(1:600),abs(zs1(1:600)))。噪聲信號波形39。 zs0=*cos(2*pi*10000*t/22050000)。語音信號頻譜圖39。 f=(FS/1630)*[1:1630]。 %sound(x,FS,bits)。x=x(:,1)。附錄D Blackman窗低通濾波器仿真程序：clf。plot(f(1:1000),abs(M1(1:1000)))。濾波后被污染的信號波形39。 %回放污染語音 figure(5)。freqz(num,den,128,fs)。[at1,bt1,ct1,dt1]=lp2hp(a,b,c,d,wc)。s39。rp=3。被噪聲污染的信號頻譜39。 subplot(2,1,2)。 plot(y) title(39。y=x+zs39。axis([0,50,0,60])。)。 figure(2)。)。plot(f(1:1630),abs(y(1:1630)))。語音信號時域波形圖39。 subplot(2,1,1)。d:\media\39。濾波后被污染的信號頻譜39。 subplot(2,1,2)。title(39。%sound(M)。)。ylabel(39。plot(fk/1000,20*log(abs(Hk)))。fk=0:2000/512:2000。)。Rp=1。被噪聲污染的信號頻譜39。y1=fft(y,2000)。 subplot(2,1,1)。)。 subplot(2,1,2)。 plot(zs0) title(39。 zs=*cos(2*pi*10000*t/22050)。title(39。 f=(FS/1630)*[1:1630]。 %sound(x,FS,bits)。x=x(:,1)。附錄B 切比雪夫帶通濾波器程序：clf。%axis([0,1000,0,150])。)。subplot(2,1,1)。[Bz,Az]=impinvar(B,A)。)title(39。頻率（Hz）39。Hk=freqs(B,A,wk)。s39。[N1,wp1]=cheb1ord(wp,ws,Rp,As,39。wp=2*pi*500。plot(f(1:1000),abs(y1(1:1000)))。被噪聲污染的信號波形39。%sound(y,FS,bits)。title(39。)。 %回放噪音 figure(2)。 t=0:length(x)1。axis([0,500,0,300])。)y=fft(x,3260)。 plot(x)。)。懂得了所有的事都得自己親手去操作才會在腦海中留下深刻的印象，這次課程設(shè)計讓我可以更加熟練地操作MATLAB軟件，也了解了各種數(shù)字濾波器的功能的應(yīng)用，也加深了對MATLAB認識和理解。5心得體會在這次為期一周的語音信號的數(shù)字濾波處理課程設(shè)計中，剛開始覺得時間有點緊，對書本知識有點生疏，但真正做起來時，發(fā)現(xiàn)前期做的比較慢，后期越來越熟練，做的也越來越快。 Blackman窗高通濾波器分析圖Blackman窗帶通濾波器可以實現(xiàn)對低頻或高頻噪聲進行濾除，將混合信號恢復(fù)至原始信號，為驗證此結(jié)論的正確性，本次課程設(shè)計中我們添加噪聲信號，高頻噪聲zs0=*cos(2*pi*100*t/22050)于原始信號上，并使用自己設(shè)計的Blackman窗濾波器對其進行處理，仿真結(jié)果證明所設(shè)計濾波器可實現(xiàn)要求。 切比雪夫低通濾波器分析圖 切比雪夫高通濾波器可以實現(xiàn)對低頻噪聲進行濾除，將混合信號恢復(fù)至原始信號，為驗證此結(jié)論的正確性，本次課程設(shè)計中我們添加噪聲信號，低頻噪聲zs0=*cos(2*pi*100*t/22050)于原始信號上，并使用自己設(shè)計的切比雪夫濾波器對其進行處理，仿真結(jié)果證明所設(shè)計濾波器可實現(xiàn)要求。title(39。M1=fft(M,2000)。plot(M)。M=filter(bz,az,y)。axis([300,1400,300,100])。h(n)波形39。subplot(2,1,1)。high39。 N0=ceil(11*pi/Bt)。fs=400。,blackman(N))。 N=N0+mod(N0+1,2)。 wp=2*pi*fp/FS。axis([300,1400,300,100])。h(n)波形39。subplot(2,1,1)。濾波后被污染的信號頻譜39。 subplot(2,1,2)。title(39。sound(M,FS,bits)。帶通損耗函數(shù)曲線39。)。 plot([bz,az])。 [bz,az]=fir1(N1,wc,blackman(N))。 Bt=wlpwls。 給被污染的信號加窗：wlp=100*2*pi/FS。 plot(w/pi*FS,20*log(abs(h)))。title(39。figure(4)。 N=N0+mod(N0+1,2)。 wls=30*2*pi/FS。濾波后被污染的信號頻譜39。M1=fft(M,2000)。plot(M)。M=filter(bz,az,y)。axis([0,3500,400,2])。h(n)波形39。subplot(2,1,1)。 wc=(wp+ws)/2/pi。 ws=2*pi*fs/FS。)。 plot(w/pi*FS,20*log(abs(h)))。title(39。figure(4)。 N=N0+mod(N0+1,2)。 wp=2*pi*fp/FS。濾波后被污染的信號頻譜39。M1=fft(M,2000)。plot(M)。M=filter(num,den,y)。%雙線性變換法[num,den]=ss2tf(at2,bt2,ct2,dt2)。[z,p,k]=cheb2ap(n,rs)。fs=1000。 給被污染的信號加高通數(shù)字濾波器wp=2*pi*350。%雙線性變換法[num,den]=ss2tf(at2,bt2,ct2,dt2)。[z,p,k]=cheb2ap(n,rs)。fs=1000。： 設(shè)計數(shù)字高通濾波器 設(shè)計指標(biāo)：先將高通濾波器的指標(biāo)裝換為對應(yīng)的低通數(shù)字濾波氣的指標(biāo)，設(shè)計相應(yīng)的低通系統(tǒng)函數(shù)，對低通系統(tǒng)函數(shù)進行頻率變換，得到高通數(shù)字濾波器，本次取高通濾波器通帶邊界頻率350Hz,阻帶邊界截止頻率150Hz MATLAB程序：wp=2*pi*350。plot(f(1:1000),abs(M1(1:1000)))。濾波后被污染的信號波形39。 %回放污染語音 figure(5)。axis([0,160,5])。幅度（dB）39。xlabel(39。wk=2*pi*fk。[B,A]=cheby1(N1,Rp,wp1,39。As=20。幅度（dB）39。xlabel(39。fk=0:1500/512:1500。)。Rp=1。濾波后被污染的信號頻譜39。M1=fft(M,2000)。plot(M)。 %脈沖響應(yīng)不變法M=filter(Bz,Az,y)。濾波器的頻譜分析39。)。figure(4)。)。s39。ws=2*pi*2000。)。Hk=freqs(B1,A1,wk)。s39。[N1,wp1]=cheb1ord(wp,ws,Rp,As,39。 MATLAB程序如下：wp=2*pi*500。被噪聲污染的信號頻譜39。y1=fft(y,2000)。 subplot(2,1,1)。噪聲信號頻譜39。zs1=fft(zs,1200)。 subplot(2,1,1)。 運行程序得到語音信號時域和頻譜圖形如圖所示:，頻譜分析：MATLAB程序如下：t=0:length(x)1。plot(f(1:1630),abs(y(1:1630)))。語音信號時域波形圖39。 subplot(2,1,1)。d:\media\39。3設(shè)計題第二部分（數(shù)字濾波器設(shè)計）：① 利用Windows下的錄音機或其他軟件，進行語音信號的采集（*.wav）；② 語音信號的頻譜分析，畫出采樣后語音信號的時域波形和頻譜圖；③ 產(chǎn)生噪聲信號并加到語音信號中，得到被污染的語音信號，并回放語音信號；④ 污染信號的頻譜分析，畫出被污染的語音信號時域波形和頻譜； ⑤ 根據(jù)有關(guān)的頻譜特性，采用間接法設(shè)計IIR數(shù)字濾波器，并畫出相應(yīng)濾波器的幅頻、相頻圖（設(shè)計3個IIR濾波器）a. 模擬濾波器類型：切比雪夫濾波器（低通、帶通、高通）b．總體要求：Matlab原程序+仿真波形+技術(shù)指標(biāo) ⑥ 根據(jù)有關(guān)的頻譜特性，采用直接法設(shè)計FIR數(shù)字濾波器，并畫出相應(yīng)濾波器的幅頻、相頻圖（設(shè)計3個FIR濾波器）a. 濾波器類型：Blackman窗（低通、帶通、高通）b．總體要求：Matlab原程序+仿真波形+技術(shù)指標(biāo)+窗函數(shù) ⑦ 用自己設(shè)計的這些濾波器分別對被不同噪聲污染的信號進行濾波；⑧ 分析得到信號的頻譜，畫出濾波后信號的時域波形和頻譜，并對濾波前后的信號進行對比，分析信號的變化；⑨ 回放語音信號。plot(abs(y2))。 n=0: