【正文】 )% sound(y,fs,bits)。 title(39。)。 title(39。 %對輸入信號進行濾波 y4=fft(y,32768)。grid on。相位 /^o39。)。 subplot(2,1,2),plot(f,180/pi*unwrap(angle(H))) xlabel(39。)。ylabel(39。頻率 /Hz39。 %用 Nn 點繪出頻率特性 subplot(2,1,1),plot(f,20*log10(abs(H)))。)。 %求得數(shù)字濾波器的最小階數(shù)和截止頻率（歸一化頻率） 17 [b,a]=butter(N,Wn,39。 %通帶波紋和阻帶衰減 Nn=128。 %根據(jù)采樣頻率將濾波器邊界頻率進行轉(zhuǎn)換 Rp=1。 %采樣頻率 wp=50*2/Fs。加單頻余弦信號的語音信號頻譜 39。 subplot(2,1,2) plot(abs(y3))。加單頻余弦信號的語音信號時域波形 39。 %sound(x3,fs,bits)。 d=[*cos(2*pi*40*t)]39。) N=length(x1)1。 title(39。)。 title(39。 y1=fft(x1,32768)。)。 高通濾波 [x1,fs,bits]=wavread(39。加單頻余弦信號濾波后的語音信號頻譜 39。 subplot(2,1,2) plot(abs(y4))。加單頻余弦信號濾波后的語音信號時域波形 39。 subplot(2,1,1) plot(y)。 figure(4) y=filter(b,a,x3)。)。ylabel(39。頻率 /Hz39。grid on。振幅 /dB39。)。 xlabel(39。 %設計 Butterworth 低通數(shù)字濾波器 figure(3) [H,f]=freqz(b,a,Nn,Fs)。 %顯示濾波器頻率特性的數(shù)據(jù)長度 [N,Wn]=buttord(wp,ws,Rp,Rs)。Rs=20。ws=4300*2/Fs。) Fs=fs。 title(39。) y3=fft(x3,32768)。 figure(2) subplot(2,1,1) plot(x3) title(39。 x3=x1+d。 t=0:1/fs:N/fs。原始語音信號頻譜 39。 subplot(2,1,2) plot(abs(y1))。原始語音信號時域波形 39。 figure(1) subplot(2,1,1) plot(x1)。 %sound(x1,fs,bits)。i:/好男人 _new39。)% %sound(y,fs,bits)。 title(39。)。 title(39。 %對輸入信號進行濾波 y4=fft(y,32768)。grid on。相位 /^o39。)。 subplot(2,1,2),plot(f,180/pi*unwrap(angle(H))) xlabel(39。)。ylabel(39。頻率 /Hz39。 %求得濾波器的頻率特性 subplot(2,1,1),plot(f,20*log10(abs(H)))。)。%求得數(shù)字濾波器的最小階數(shù)和歸一化截止頻率 [b,a]=cheby1(N,Rp,Wn,39。Nn=128。 %阻帶邊界頻率（歸一化頻率） Rp=1。 %采樣頻率 wp=[4500 5500]*2/Fs。加單頻余弦信號的語音信號頻譜 39。 subplot(2,1,2) plot(abs(y3))。加單頻余弦信號的語音信號時域波形 39。 %sound(x3,fs,bits)。 d=[*cos(2*pi*5000*t)]39。) N=length(x1)1。 title(39。)。 title(39。 y1=fft(x1,32768)。)。 14 6. 附件 程序源代碼： 帶阻濾波 [x1,fs,bits]=wavread(39。雖然中間遇到很多問題，但通過自己的仔細排查，老師的幫助，同學之間的互相討論基本上解決了。有互聯(lián)網(wǎng) 這 個幫手，對于不懂的可以從網(wǎng)上搜索，也可以查找很多有用的函數(shù)、源代碼及其它有用的資料。我這次做的題目是 采用 IIR 濾波器的完全設計函數(shù)來設計各型 IIR 濾波器（ Butterworth 型與切比雪夫 I 型、切比雪夫Ⅱ型和橢圓型的低通、高通、帶通、帶阻中的至少 3 種類型） ，繪 出 濾波器的頻域響應； 并用這些數(shù)字濾波器對含噪語音信號分別進行濾波處理，比較不同方法下設計出來的數(shù)字濾波器的濾波效果，并從理論上進行分析。如下圖所示： 13 圖 41 程序調(diào)試錯誤圖示 后來經(jīng)過調(diào)試才知道程序中有多處應用了 sound 函數(shù)，才導致這樣的錯誤，在注釋掉多余的 sound 函數(shù)，只剩下最后一個后，程序正常工作，而且聽到了濾波后干凈的歌聲。)中有個 stop 表示帶阻，但首先一直沒加，后來看以前做過的上機實驗才知道寫錯了。 在設計帶阻濾波器時，語句 [b,a]=cheby1(N,Rp,Wn,39。 0 50 100 150 200 250 300 350 400 450 500 8 0 6 0 4 0 2 00頻率 / H z振幅/dB0 50 100 150 200 250 300 350 400 450 500 4 0 0 3 0 0 2 0 0 1 0 00頻率 / H z相位/o 圖 35 高通濾波器性能指標 12 0 0 . 5 1 1 . 5 2 2 . 5 3 3 . 5x 1 051 0 . 500 . 51加單頻余弦信號濾波后的語音信號時域波形0 0 . 5 1 1 . 5 2 2 . 5 3 3 . 5x 1 04050100150200加單頻余弦信號濾波后的語音信號頻譜 圖 36 加單頻余弦信號濾波后 波形 4. 程序調(diào)試 在設計低通濾波器時，首先通帶邊界頻率 wp 和阻帶邊界頻率 ws 取值不當，分別取 3500*2/Fs 和 4000*2/Fs，使濾波效果太差，后來經(jīng)過不斷調(diào)試才取值為4000*2/Fs 和 4300*2/Fs，使濾波效果很理想。加單頻余弦信號濾波后的語音信號頻譜 39。 subplot(2,1,2) plot(abs(y4))。加單頻余弦信號濾波后的語音信號時域波形 39。 subplot(2,1,1) plot(y)。 figure(4) 11 y=filter(b,a,x3)。)。ylabel(39。頻率 /Hz39。grid on。振幅 /dB39。)。 xlabel(39。 %設計 Butterworth 高通數(shù)字濾波器 figure(3) [H,f]=freqz(b,a,Nn,Fs)。high39。 %顯示濾波器頻率特性的數(shù)據(jù)長度 [N,Wn]=buttord(wp,ws,Rp,Rs)。Rs=20。ws=20*2/Fs。 0 2021 4000 6000