1、長沙理工大學數字信號處理課程設計報告學院計算機與通信工程專業通信工程班級_學號學生姓名 指導教師 課程成績 完成日期 課程設計成績評定學 院計算機與通信工程專 業通信工程班 級學 號學生姓名指導教師完成日期 2015年3月13日指導教師對學生在課程設計中的評價評分項目優良中及格不及格課程設計中的創造性成果學生掌握課程內容的程度課程設計完成情況課程設計動手能力文字表達學習態度規范要求課程設計論文的質量指導教師對課程設計的評定意見綜合成績指導教師簽字2015年3月15 日課程設計任務書計算機與通信工程 系通信工程 專業課程名稱數字信號處理課程設計時間 20142015學年第二學期12周學生姓名指導

2、老師題目語音信號濾波去噪一使用頻率米樣法設計的頻率米樣型FIR濾波器主要內容：用麥克風米集一段語音信號，繪制波形并觀察其頻譜特點，加入一個帶 外單頻噪聲，使用頻率米樣法設計 個滿足指標的 FIR濾波器，對該含噪語音 信號進行濾波去噪處理，畫出頻率采樣型濾波器結構圖。比較濾波前后的波形 和頻譜并進行分析，根據結果和學過的理論得出合理的結論。 與不冋濾波方法 的同學比較各種濾波方法性能的優劣。要求：（1）濾波器指標必須符合工程實際。（2）設計完后應檢查其頻率響應曲線是否滿足指標。（3 ）處理結果和分析結論應該一致，而且應符合理論（4）獨立完成課程設計并按要求編寫課程設計報告書 應當提交的文件：（1

3、 ）課程設計學年論文。（2）課程設計附件（主要是源程序）。語音信號濾波去噪一一使用頻率采樣法設計的頻率采樣型FIR濾波器學生姓名：指導老師：摘要此次課程設計是使用頻率采樣法設計一個頻率采樣型濾波器，對一段加噪的語音信號進行濾波去噪處理，并且對濾波前后的時域波形和頻譜分析濾波性能。此次課程設計仿真平臺為MATLAB7.0，開發工具是M語言編程。首先利用Windows錄音機錄 制一段語音信號，并人為加入單頻噪聲，然后對信號進行頻譜分析以確定所加噪聲頻率， 再設計濾波器進行濾波去噪處理，最后比較濾波前后的波形和頻譜并進行分析，并用 VISIO繪制濾波器的結構圖。由分析結果可知，濾波后的語音信號與原始

4、語音信號基本 一致，即設計的FIR濾波器能夠去除信號中所加單頻噪聲，即達到了設計目的。關鍵詞 MATLAB ;濾波去噪；頻率采樣法；頻率采樣型濾波器；VISIO1引言濾波去噪是信號處理中一種最基本但十分重要的技術。利用濾波可以從復雜的信 號中提取所需的信號，抑制不需要的信號。濾波器就是這樣一種可以在時域和頻域對信 號進行濾波處理的系統。通常情況下，有用信號和干擾信號是在不同頻段上的，于是通 過對濾波器的頻率特性精心設計就能達到濾波的目的。本課程設計是采用頻率采樣法設計頻率抽樣型濾波器，從而對語音信號濾波去噪。通過對比濾波前后的波形圖及回放濾波前后的語音信號，來判斷濾波器對噪聲信號確實有濾除作用

5、。1.1課程設計目的(1) 熟悉使用MATLAB ；(2) 了解FIR濾波器原理及結構；(3) 利用所學數字信號處理相關知識用 MATLAB設計一個FIR濾波器；提高自己動手的能力。(5)對加噪的語音信號進行濾波去噪處理，比較濾波前后的時域波形和頻譜并進行分 析；1.2課程設計的要求(1) 濾波器指標必須符合工程設計。(2) 設計完后應檢查其頻率響應曲線是否滿足指標。(3) 處理結果和分析結論應該一致，而且應符合理論。(4) 獨立完成課程設計并按要求編寫課程設計報告。1.3設計平臺本課程設計仿真平臺為 MATLAB7.0。MATLAB的名稱源自Matrix Laboratory ，1984年由

6、美國Mathworks 公司推向市場。它是一種科學計算軟件，專門以矩陣的形 式處理數據。MATLAB將高性能的數值計算和可視化集成在一起，并提供了大量的內 置函數，從而被廣泛地應用于科學計算、控制系統、信息處理等領域的分析、仿真和設 計工作。1993年MathWorks 公司從加拿大滑鐵盧大學購得 MAPLE軟件的使用權，從而以MAPLE為“引擎”開發了符號數學工具箱(Symbolic Math Toolbox ) 22設計原理用windows自帶的錄音機錄制一段語音信號，繪制波形并觀察其頻譜，給定相應 技術指標，用頻率采樣法設計的一個滿足指標的頻率采樣型FIR濾波器，對該語音信號進行濾波去噪

7、處理，比較濾波前后的波形和頻譜并進行分析。2.1 FIR濾波器的設計FIR(Finite Impulse Response)濾波器3:有限長單位沖激響應濾波器，又稱為非遞歸型濾波器，是數字信號處理系統中最基本的元件，它可以在保證任意幅頻特性的同 時具有嚴格的線性相頻特性，同時其單位抽樣響應是有限長的，因而濾波器是穩定的系 統。因此，FIR濾波器在通信、圖像處理、模式識別等領域都有著廣泛的應用。有限長單位沖激響應(FIR)濾波器有以下特點：(1) 系統的單位沖激響應h (n)在有限個n值處不為零(2) 系統函數H(z)在|z|>0處收斂，極點全部在z = 0處(因果系統)(3) 結構上主要

8、是非遞歸結構，沒有輸出到輸入的反饋，但有些結構中(例如頻率 抽樣結構)也包含有反饋的遞歸部分。2.2頻率采樣型結構把一個有限長序列(長度為N點)的z變換H (z)在單位圓上作N等分抽樣，就得 到H (k)，其主值序列就等于h (n)的離散傅里葉變換H (k)。那里也說到用H (k)表示的 H (z)的內插公式為H(z)(1N、1 N 1 H (k)N k o1 Wn z 1(2.1)這個公式就為FIR濾波器提供了另外一種結構，這種結構由兩部分級聯組成。1 N 1H(z) Hc(z)-Hk(z)(2.2)N k o其中級聯的第一部分為梳狀濾波器，其結構如下圖所示：Hc(z)(1 z N)(2.3

9、)X(n)仁Y(n)m-Z -N圖2-1梳狀濾波器結構圖第二部分由N個諧振器組成的諧振柜。它是由N個一階網絡并聯組成，而這每一個一階網絡都是一個諧振器Hk(z)1 WNkz1(2.4)其結構如下圖所示：H(k)Hk(z>圖2-2 一階諧振器頻率抽樣型結構特點：2(1) 它的系數H(k)直接就是濾波器在Wk k處的頻率響應。因此，控制濾波器的N頻率響應是很直接的。(2) 結構有兩個主要缺點：(a) 所有的相乘系數及 H(k)都是復數，應將它們先化成二階的實數，這樣乘起來較復雜，增加乘法次數，存儲量。(b) 所有諧振器的極點都是在單位圓上，由WNk決定考慮到系數量化的影響，當系數量化時，極點

10、會移動，有些極點就不能被梳狀濾波器的零點所抵消。(零點由延時單元決 定，不受量化的影響)系統就不穩定了。(3) 將一階網絡合并為二階網絡(a)第k和第N-k個諧振器合并為一個實系數的二階網絡，因為h(n)是實數，它的DFT也是圓周共軛對稱的(2.5)H(k) H*(N k) k 1,2,3丄 N 1因此，可以將第k和第N-k個諧振器合并為一個二階網絡Hk(z)H(k)1 rWNkz 1H (N k)H (k)(NT)1k11 Wnz 1 Wn z*H (k)1 rWNk*z 1H(k)1 z 1Wn k* WNkr r2WNkWNk*z210k 1kz1 z 12r cos(2 k)N其中：0

11、k 2ReH( k), 1k2rReH( k)W,(2.6)(b)第k和第N-k個諧振器合并為一個二階網絡的極點在單位圓內，而不是在單位圓上，因而從頻率響應的幾何解釋可知，它相當于一個有限Q的諧振器。其諧振頻率Wk 為：丨 .f1 z1丁 圖2-3二階網絡結構圖除了共軛復根外，還有實根。當N=偶數時，有一對實根,H°(z)0它們分別為嚴和H1 rz 1(2.7)當N=奇數時，只有一個實根z=r（k=O ）,即只有HO（z）.（c）修正頻率抽樣結構流圖（N=偶數）x(n)圖2-4修正頻率抽樣結構流圖（N=偶數）H(z)N N、(1 r z )1 H(0)N1 rz 1H(j)N/211

12、 rz 1Ok11kzz12rco畤k)(2.8)修正頻率抽樣結構流圖（N=奇數）x(n)7y(n),v -y-rz2rcos()"nJ1-1-rz .Ar2rcos().V/(0J圖2-5修正頻率抽樣結構流圖（N=奇數）N N、1H(z) (1 r z)1H(0)1 rz 1N/2 11Ok 1kz2k 1 1 z 12r cos(k)N(2.9)2.3頻率采樣法設計思路:這種設計方法是從頻域進行設計的一種方法，首先給定一個希望逼近的頻率響應。H(k)Hd(ej )| 2k 0,1,2,., N 1(2.10)kN知道H(k)后,由IDFT定義，可以用這N個采樣值H(k)來惟一確定

13、有限長序列h(n)， 即N 11nkh(n)H(k)WNn 0,1,2,., N 1(2.11)N k oN 1H (z)h(n)z n(2.12)n 0N 1H (ejw ) h(n )e jwn(2.13)n 0內插公式H(ej )N 12H (k)kk 0N(2.14)式中,0(3 )是內插函數()sin( N / 2) j (n 1)/2eN sin(/2)(2.15)四種線性相位的FIR濾波器如下表2-1所示表2-1四種線性相位的FIR濾波器類型甲®)1型及關于疋 亦偶對稱®厘型刃山慶于=Ck Qjt偶對茹*關于tu=ic奇對稱2型仙戶也Mid). M旗迫做血關于

14、=0、JL亦奇時稱V- 1J74型辺)關于-0、加苛對稱關于厲=討禺對稱2 2逼近誤差及其改進措施這樣設計所得到的頻響H(ej 3)與要求的理想頻響Hd(ej 3)會有一定的差別(1)采樣點上濾波器的實際頻率響應是嚴格地和理想頻率響應數值相等的。但是在采樣點之間的頻響則是由各采樣點的加權內插函數的延伸疊加而成的,因而有一定的逼近誤差，誤差大小取決于理想頻率響應曲線形狀。(3) 理想頻率響應特性變化越平緩，則內插值越接近理想值，逼近誤差越小。(4) 如果采樣點之間的理想頻率特性變化越陡，則內插值與理想值的誤差就越大，因而在理想頻率特性的不連續點附近，就會產生肩峰和起伏。濾波器性能的改善：(1)

15、增加過渡帶采樣點，它可以大大減少振蕩，阻帶衰減也可以得到進一步改善。一般一點到二點的過渡帶采樣即可得到滿意的結果。(2) 增加采樣點密度，過渡帶的寬度與采樣點數N成反比。但N值意味著或長度的增加，濾波器運算量必然增大4。3設計步驟3.1設計流程圖錄制語音信號ht.wav。并用wavread函數采集加入2500Hz的噪聲。進行快速傅里葉變換，繪制加噪前后的頻譜對比圖，和時域對比圖設定濾波器的性能指標回放語音信號圖3-1流程圖3.2錄制語音信號使用電腦錄音機錄制一段語音信號屬性設置為8000Hz ,8位單聲道的音頻格式。如圖3.1-2所示:用頻率采樣法設計頻率采樣型的FIR濾波器丿是否、符合標準/

16、VN用設計的濾波器進行濾波處理， 回放濾波后的信號，并保存1比較濾波前后語音信言號的波形和頻譜圖3-2語音信號設置然后在MATLAB軟件平臺下，利用函數 wavread對語音信號進行采樣，源程序為：x,fs,bits=wavread('ht.wav')，記住采樣頻率和采樣點數，matlab 實現得：fs=8000 ；bits=8。3.3語音信號的頻譜分析在MATLAB中編寫加噪之后時域，頻域對比圖的源程序如下：x,fs,bits=wavread('D:ht.wav');%輸入參數為文件的全路徑和文件名，輸出的第一個參數是每個樣本的值，fs是生成該波形文件時的采樣

17、率，bits是波形文件每樣本的編碼位數。sound(x,fs,bits); %按指定的采樣率和每樣本編碼位數回放N=length(x); %計算信號x的長度fn=2500;%單頻噪聲頻率，此參數可改t=0:1/fs:(N-1)/fs; %計算時間范圍，樣本數除以采樣頻率x=x' y=x+0.1*si n(fn*2*pi*t); % 加噪sound(y,fs,bits); %應該可以明顯聽出有尖銳的單頻嘯叫聲X=abs(fft(x); Y=abs(fft(y); %對原始信號和加噪信號進行fft變換，取幅度譜X=X(1:N/2); Y=Y(1:N/2);%截取前半部分deltaf=fs/

18、N;%計算頻譜的譜線間隔f=0:deltaf:fs/2-deltaf;%計算頻譜頻率范圍運行結果如下圖所示:圖3-3加噪前后時域頻域對比圖由圖3.3-1可以看出，在頻域為2500Hz處加入了一個單頻噪聲，而加入噪聲之后, 時域的波形出現了明顯失真，通過聽取原聲音信號 x和加噪聲音信號y，可以明顯聽到 y語音信號中有一明顯尖銳噪聲。3.4濾波器設計設計一個帶阻濾波器，將單頻信號濾出去，源程序如下所示：fpd=2400;fsd=2450;fsu=2550;fpu=2600;Rp=0.5;As=45;% 帶阻濾波器設計指標fcd=(fpd+fsd)/2;fcu=(fpu+fsu)/2;df=mi n

19、( (fsd-fpd),(fpu-fsu);%計算上下邊帶中心頻率fcd fcu，和頻率間隔dfwp=fcd/fs*2*pi;ws=fcu/fs*2*pi;dw=df/fs*2*pi;% 將 Hz 為單位的模擬頻率換算為rad為單位的數字頻率M=41;alpha=(M-1)/2;%設置抽樣點數l=0:M-1;wl=2*pi*l/M;T1=0.6025,T2=0.127% 過渡帶的抽樣值NP仁fix(M*fpd/8000);NS仁fix(M*fsd/8000);NS2=ceil(M*fsu/8000);NP2=ceil(M*fpu/8000);%各邊界頻率的抽樣值 1Hrs=o nes(1,NP

20、1-2),0.6025,0.127,zeros(1,NP2-NP1+1),0.127,0.6025, on es(1,M-2 *NP2-5),0.6025,0.127,zeros(1,NP2-NP1+1),0.127,0.6025,o nes(1,NP1-3);Hdr=1,1,0,0,1,1;% 帶阻濾波器的參數設置wdl 仁 fpd*2/fs;wdl2=fsd*2/fs;wdl3=fsu*2/fs;wdl4=fpu*2/fs;%rad wdl=0,wdl1,wdl2,wdl3,wdl4,1;k1=0:floor(M-1)/2);k2=floor(M-1)/2)+1:M-1;an gH=-al

21、pha*(2*pi)/M*k1,alpha*(2*pi)/M*(M-k2);H=Hrs.*exp(j*a ngH);h=real(ifft(H,M);db,mag,pha,grd,w=freqz_m(h,1);Hr,ww,a,L=hr_type2(h);程序運行結果如下圖3.3-2所示:M ig j re 1Fil Edit Vi曹訓n榊rt 口Qk 總衣hfiop' WwHfow 屮*P."：' X固 i 旦遞液器覃怔沖擊響N圖3-4濾波器參數圖由圖3-4可以看出，濾波器的衰減大于設定值As=45,滿足性能指標，濾波器的衰減可由增加過度帶寬來得到。上圖中的抽樣值M二

22、41，因此阻帶的帶寬比較寬。可以通過調大M值來使得阻帶的帶寬變窄，但這樣會使得濾波器的階數相應的增加，使得 其結構圖標的復雜，為了是濾波器能達到相應的性能要求，又使得其結構變得簡單，最 終取M值為41 03.5信號濾波處理原程序如下所示：%進行濾波處理z=filter(h,1,y); %用設計好的濾波器對y進行濾波Z=fft(z);Z=Z(1:length(Z)/2);% 計算頻譜取前一半程序運行結果如下圖3.3-3所示:圖3-5濾波前后時域，頻域對比圖頻譜取前一半由上圖3-5可以看出，在濾波之后時域圖能得到恢復，頻域圖中的單頻噪聲信號也的到濾除，說明了設計的濾波器能濾除加入的噪聲信號，因此說

23、達到了 設計的要求。3.6結果分析開始通過分析決定設計一個帶阻的濾波器來濾除加入的單頻噪聲，根據噪聲的頻率來設計阻帶的范圍。在采集語音信號后，按照步驟用頻率采樣法設計頻率采樣型濾波器。由圖3-4可知,設計的濾波器達到要求。我們觀察到圖3-5濾波前后語音信號的波形對比圖，發現時域波形中加干擾噪聲后有明顯的變化，不過經過濾波后幾乎沒有變化，說明設計的濾波 器達到要求。再通過聽取原始語音信號，加噪信號 y,和濾波之后的信號z。對比之后, 發現濾波器確實濾除了噪聲。從理想的角度考慮，該帶阻濾波器的阻帶帶寬應該可以變得更窄， 讓濾波效果更好, 但是這樣的采樣值會變得非常大，考慮到實際的情況，通過權衡之后

24、，決定犧牲帶寬來 使得濾波器的階數降低，因此在上圖3-5中我們可以看出，在噪聲頻譜左右兩邊的信號 也被濾除了。3.6濾波器結構設計通過計算得出濾波器的系數如下圖所示:B -2. QQdQ1, 0000-E97661. OOOQ亂 gotof10.W11. 000 Q-1.S0681.00002.嘰 9083-0,98831. 0000-1.79231. 00002.0000-0-97370. 9737昇 0000-L«36»1.00002. ODOQ0. 95347. 9534】* 0000-U44101,00002. 0000-0.»2750.A275L 000

25、0-1. 212&1,00002, 00000, 8»«27” B962.0000-0.95541. 00002. ODCO-0, 8&960,95961. 0000-0.4760l.ooao2- 00000. 8179-0.町"h 0D00-CL 38081. 0000-0*77151,77161” 00000,07661. 0000隊 2540叭 72057+ 7205L 00000.22941. 00000 . 0000663Q. 6653.00000.53001.0000a. oata0. 3E62» »«h

26、0D000.81B11, 0000此 oodo(L 00000.»1351,0000Loan1.0000仇 26400. 117+ "酣1.00001.3307E_ 0000U2050-0.40110” 40911. 00001.54301.00002. 00000.-0.33M1.00<0L719I1.00002. 0000一0. 2«500. 2臨益0I.OUOQL>5501.00002.00(10-0*10041. 00001.94741.00002-0000-0. 11470, IU71. 00001.99H1.00D01. 0000S 283

27、-0,0383L 0000-LOQOO0圖3-6濾波器系數根據以上求出的系數，利用 VISIO軟件來繪制濾波器的結構圖如下圖所示:Z1Z1-11-1Z1-11.6359Z1-10.9275 Az1z1才-0.89621.2125-1.1*z1-10.3808 z1.1-1z1z1-0.6653 0.25400.665307-1z1L-1.7191z1-1z111-11.0000z1-0.7205z1 -0.9534 I-11.9068Z1 -0.98830.98832.0000-0.7715 2.00000.95342.00001.44104-1-0.26502.0000-0.4091 0.2

28、540 b 0.9971 i0.89622.00000.2650z1>1.97661.7923Z1 0.97370.33801.2050-0.8596 2.0000Z10.8596-0.9737 2.0000-0.8179J-0.99712.0000>0.81792.0000Z1 0.4091 1z1-0.92752.00000.72051.2050-0.11472.0000/ 0.1147：0.9554-12.0000 0.0383z1Z-0.0383I -1.9474z10.19042.0000-0.1904 Ay（n）-0.2294-1Z1-1.5430 z1 -033801

29、-1-1.8550-1Z1F1圖3-7濾波器的結構圖4出現的問題及解決辦法在這次課程設計當中，由于基礎不扎實，出現了很多問題，即有MATLAB軟件操作不當，也有知識掌握程度不夠出現的各種問題。1，在一開始錄制語音信號時，沒有修改參數值，導致頻率抽樣過高，后在老師的 提醒之下將頻率修改為8000hz。2，在調用音樂文件時沒有將文件放在 MATLAB的工作文件夾下面，到時文件找不 到，還有程序中用到的各種函數都沒有放入正確的位置。3，在繪制加噪前后頻率對比圖時，留白過多，對比不明顯，通過axis函數對橫縱坐標進行限定4，在使用visio軟件繪制濾波器的結構圖時，沒有合理規劃整個頁面，導致重新繪 制

30、。5，在使用visio繪制結構圖時，前面濾波器參數設置出現錯誤，所以濾波器的系數 發生了改變，只能重新修改參數。6，在設計濾波器的過程中，沒有準確理解每一條指令代表的含義，導致程序前后 不對應，出現很多錯誤，通過 matlab中的錯誤提示，準確找到錯誤的那一行代碼，進 行修改。7,在設計濾波器的過程當中，濾波器的衰減小于開始所設置的值，通過犧牲過渡 帶和調節過渡帶的采樣值，即T1和T2來使得衰減大于所設定值As。8，頻譜圖沒有進行去前半段的處理。9，最后聽取濾波后聲音，沒有將其保存，在請教了同學之后，學會了如何保存濾 波后的聲音文件。10,在設計濾波器的過程當中，抽樣點數過高，是的濾波器的階數

31、很高，導致繪圖 困難。通過降低抽樣的點數，來平衡濾波器階數和阻帶帶寬之間的矛盾。5結束語這是第三次進行課程設計了，在前面的課程設計當中我們學習到了如何使用 matlab，所以對于matlab軟件的使用并沒有那么陌生了，盡管如此，在使用 matlab 的過程當中還是出現了很多錯誤，比如說忘記添加函數文件，參數前后不對應之類的低 級錯誤。在這次課程設計當中，還學習到了另外一個繪圖軟件，VISIO，相對而言，是比較容易操作的，在這次課程設計中主要用來繪制各類結構圖，以及流程圖。此次課程設計，讓我更深入的了解到頻率抽樣法以及頻率抽樣型的濾波器，由開始 的無從下手，再翻閱書上的例子，首先決定使用低通的濾

32、波器將噪聲濾除，但由于是加 入的是一個單頻的信號噪聲，就決定使用帶阻濾波器，在網上查找相關資料時，大多是 低通濾波器，沒有帶阻濾波器的例子，所以就通過書籍信號與系統上高通的例題來 進行修改，當然這過程不是一蹴而就的。經過了一次又一次的錯誤，才得出來最后的模 型。還有這也是我和其他同學一起討論出來的。這次課程設計讓我對濾波器的類型有了一個更加完整了解，在設計中也使我對一些概念有了更深刻的認識。例如：在濾波器分類方面，我深刻的了解了低通，高通濾波器 與帶通，帶阻濾波器的特性區別。還有在課程設計中每一次的數據輸入都有其重要意義, 用MATLAB編譯程序時，可以根據濾波器指標的要求實時知道對濾波器的影

33、響。通過 一次次的調試和權衡使濾波器的性能達到最佳。課程設計不僅要求對濾波器理論的研究，更重要的是培養一種遇到問題解決問題 的思維。因為有了這次課程設計，我懂得了書本知識只是實際應用的理論指導。如果僅 僅只學習書本知識，不去在實踐中運用，那只是停留在只知其表，不知其因的層面。比 如在數學計算上，可以將噪聲完全濾除，而在這次設計中，若要完全濾除噪音，濾波器 的階數就會增高，在現實生活中是很難實現的，所以噪聲是不能完全濾除的。課程設計結束了，我相信這次課程設計對今后的學習是很有幫助的，它讓我將理論 更好地和實踐相結合，提高了動手的能力，也填補了自己學習上的一些不足。這次課程設計能取得成功，不僅僅是

34、我一個人的努力的結果，更離不開指導老師與 同學的幫助，在此向老師和同學們表示衷心的感謝。參考文獻1 吳鎮揚，數字信號處理M.高等教育出版社，20042 張圣勤.MATLAB7.0 實用教程M.北京：機械工程出版社，20063 程佩青數字信號處理教程M 北京：清華大學出版社，20024 高西全，丁玉美.數字信號處理M.第三版.西安：西安科大出版社，1994附錄：源程序%程序名稱：htl.m%程序功能：設計一頻率采樣法設計的頻率采樣型的濾波器，然后對一加單憑頻噪聲信 號去噪%程序作者：賀婷%最后修改日期：2015-3-10x,fs,bits二wavread('D:ht.wav');

35、%輸入參數為文件的全路徑和文件名，輸出的第一個參數是每個樣本的值，fs是生成該波形文件時的采樣率，bits是波形文件每樣本的編碼位數。 sound(x,fs,bits);%按指定的采樣率和每樣本編碼位數回放N=length(x); %計算信號x的長度fn=2500;%單頻噪聲頻率，此參數可改t=0:1/fs:(N-1)/fs; % 計算時間范圍，樣本數除以采樣頻率x=x' y=x+0.1*si n(fn*2*pi*t); % 添加噪聲sound(y,fs,bits); %可以明顯聽出有尖銳的單頻嘯叫聲X=abs(fft(x);Y二abs(fft(y);%對原始信號和加噪信號進行fft變

36、換,取幅度譜X=X(1:N/2); Y=Y (1:N/2);% 截取前半部分deltaf=fs/N;%計算頻譜的譜線間隔f=0:deltaf:fs/2-deltaf;%計算頻譜頻率范圍%設計濾波器fpd=2400;fsd=2450;fsu=2550;fpu=2600;Rp=0.5;As=45;% 帶阻濾波器設計指標fcd=(fpd+fsd)/2;fcu=(fpu+fsu)/2;df二min(fsd-fpd),(fpu-fsu);% 計算上下邊帶中心頻率fcd fcu，和頻率間隔dfwp=fcd/fs*2*pi;ws=fcu/fs*2*pi;dw=df/fs*2*pi;% 將 Hz 為單位的模擬

37、頻率換算為rad為單位的數字頻率M=41;alpha=(M-1)/2;%設置抽樣點數l=0:M-1;T1=0.6025,T2=0.127% 過渡帶的抽樣值NP1= fix(M*fpd/8000);NS 1= fix(M*fsd/8000);NS2二ceil(M*fsu/8000);NP2二ceil(M*fpu/8000);%各邊界頻率的抽樣值1Hrs=o nes(1,NP1-2),0.6025,0.127,zeros(1,NP2-NP1+1),0.127,0.6025,on es(1,M-2*NP2-5),0.6025,0.127,zeros(1,NP2-NP1+1),0.127,0.6025

38、,on es(1,NP1-3);k1=0:floor(M-1)/2);k2=floor(M-1)/2)+1:M-1;an gH=-alpha*(2*pi)/M*k1,alpha*(2*pi)/M*(M-k2);H=Hrs.*exp(j*a ngH);h=real(ifft(H,M);% 求 FIR 單位沖擊響應db,mag,pha,grd,w=freqz_m(h,1);% 幅頻特性曲線Hr,ww,a,L=hr_type2(h);%畫圖subplot(2,2,1);stem(l,h);title(' 濾波器 單位沖 擊響應');xlabel(' n');ylabe

39、l('y( n)');axis(O 40 -1 1);grid on; subplot(2,2,2);plot(w/pi,pha);title('濾 波器相 位響 應');xlabel('w/pi');ylabel('相位 pha');axis(0 1-4 4);grid on;subplot(2,2,3);plot(w/pi,db);title('濾 波 器 幅 度 響 應');xlabel('w/pi');ylabel('dB');axis(O 1 -50 20);grid on

40、;%畫線X_l二0,0,wp/pi,ws/pi;1,1,wp/pi,ws/pi;Y_l二-As,-Rp,-50,-50;-As,-Rp,20,20;%在wp , ws , Rp , As處畫線以更直觀判斷設計是否達標，每列參數是 每個線條的端點坐標line(X_l,Y_l,'Color','r','LineWidth',2,'LineStyle','-')% 添加線寬為 2的紅色線條subplot(2,2,4);plot(w/pi,mag);title('濾 波器幅 度響應');xlabel('w/pi');ylabel(' 幅度 mag');axis(0 1 0 1.2);grid on;%進行濾波處理z=filter(h,1,y); % 用設計好的濾波器對x1進行濾波Z=fft(z);Z=Z(1:length(Z)/2);% 計算頻譜取前一半%比較濾波前后的性能%在同一張大圖里分別繪制原始信號x，加噪信號y，濾波去噪信號z的時域波形和頻譜，以便比較和分析。figure(2);subplot(321);plot(t,x);axis tight;grid on;title('原始語音信號時間x