語音信號的FIR濾波器處理(華南農業(yè)大學)

1、DSP課程設計報告語音信號的FIR濾波器處理姓 名：XXX XXX XXX班 級：10電信1 學 號： XXXXXXXX指導老師：徐梅宣日期： 2021.06.102021.06.21華南農業(yè)大學工程學院摘要DSP( Digital Signal Processing) 也就是我們常說的數字信號處理, 它是利用計算機或專用處理設備, 以數字形式對信號進行采集、變換、濾波、估值、增強、壓縮、識別等處理, 以得到符合人們需要的信號形式。DSP 芯片在信號處理、通信、雷達等許多領域得到廣泛的應用。Matlab是一款強大的軟件，它將數值分析、矩陣計算、科學數據可視化等諸多強大功能集成在一個易于使用的視

2、窗環(huán)境中，為科學研究、工程設計以及必須進行有效數值計算的眾多科學領域提供了一種全面的解決方案。CSS集成開發(fā)環(huán)境使用CCS內置的軟件仿真simulator對程序進行編譯，調試和運行，主要用于檢測目標程序運行的正確性和連貫性,并能通過仿真器與目標板連接，在目標板上實時觀察效果。在本次設計中，我們選擇的課題是基于DSP的語言信號的FIR濾波處理。首先利用MATLAB進行了仿真，得到濾波前后的時域波形和頻譜。然后通過調用MATLAB的分析工具FDATOOL，根據仿真結果導出了濾波器的相關參數，將原始信號數據和濾波器參數輸入CCS進行DSP編程。最后在DSP中實現了FIR帶通濾波，并通過CCS的頻譜分

3、析功能查看了最終DSP的濾波效果。關鍵詞:語音信號 DSP FIR濾波 MATLAB CCS目 錄1.前言42.設計原理42.1數字信號處理器42.2濾波器簡介52.3 fir濾波器原理62.4窗函數簡介73.FIR濾波器的MATLAB設計83.1總體方案的設計83.2語音信號的采集93.3 MATLAP的具體處理103.5 FIR濾波器的系統(tǒng)參數144.FIR濾波器的CCS設計與仿真164.1 CCS程序流程圖174.2 CCS仿真結果205.問題及解決辦法226.設計感想23參考文獻237.附錄241 前言隨著信息與數字技術的發(fā)展，數字信號處理已經成為當今極其重要而學科與技術領域之一。它在

4、通信、語音、圖像、自動控制、雷達、軍事、航空航天、醫(yī)療和家用電器等眾多領域得到了廣泛的應用。在數字信號處理的基本方法中，通常會涉及到變換、濾波、頻譜分析、調制解調和編碼解碼等處理。其中濾波是應用非常廣泛的一個環(huán)節(jié)，數字濾波器的理論和相關設計也一直都是人們研究的重點之一。FIR濾波器的是非遞歸的，穩(wěn)定性好，精度高；更重要的是，FIR濾波器在滿足幅頻響應要求的同時，可以獲得嚴格的線性相位特征。因此，它在高保真的信號處理，如數字音頻、圖像處理、數據傳輸和生物醫(yī)學等領域得到廣泛應用。在數字信號處理中，濾波占有極其重要的地位。數字濾波是語音信號處理、圖像處理、模式識別、頻譜分析等應用的基本處理算法。用D

5、SP芯片實現數字濾波除了具有穩(wěn)定性好、精確度高、不受環(huán)境影響等優(yōu)點外，還具有靈活性好等特點。本文介紹了數字濾波器的設計基礎及用窗函數法設計FIR濾波器的方法，運用MATHLAB語言實現了低通濾波器的設計并用CCS 3.3進行觀察效果。2 設計原理2.1 數字信號處理器DSP（Digital Signal Processor）是一種獨特的微處理器，是以數字信號來處理大量信息的器件。其工作原理是接收模擬信號，轉換為0或1的數字信號。再對數字信號進行修改、刪除、強化，并在其他系統(tǒng)芯片中把數字數據解譯回模擬數據或實際環(huán)境格式。它不僅具有可編程性，而且其實時運行速度可達每秒數以千萬條復雜指令程序，遠遠超

6、過通用微處理器，是數字化電子世界中日益重要的電腦芯片。它的強大數據處理能力和高運行速度，是最值得稱道的兩大特色。由Ti公司提供專業(yè)的開發(fā)工具CCS，自帶DSP/BIOS操作系統(tǒng)，能夠直接編寫適合DSP開發(fā)工程及文件，滿足DSP程序設計要求。由MathWorks 公司和TI 公司聯合開發(fā)的DSPMATLAB Link for CCS Development Tools（簡稱CCSLink）是MATLAB6.5 版本（Release13）中增加的一個全新的工具箱，它提供了MATLAB、CCS 和DSP 目標板的接口，利用此工具可以像操作MATLAB變量一樣來操作DSP 器件的存儲器和寄存器，使開發(fā)

7、人員在MATLAB環(huán)境下完成對DSP的操作，從而極大地提高DSP應用系統(tǒng)的開發(fā)進程。本設計主要是使用DSP對語音噪聲信號進行處理，使用CCS開發(fā)工具編寫程序，同時利用MATLAB操作DSP 器件的存儲器和寄存器。2.2 濾波器簡介濾波器，顧名思義，是對波進行過濾的器件?！安ā笔且粋€非常廣泛的物理概念，在電子技術領域，“波”被狹義地局限于特指描述各種物理量的取值隨時間起伏變化的過程。該過程通過各類傳感器的作用，被轉換為電壓或電流的時間函數，稱之為各種物理量的時間波形，或者稱之為信號。因為自變量時間是連續(xù)取值的，所以稱之為連續(xù)時間信號，又習慣地稱之為模擬信號(Analog Signal)。隨著數字

8、式電子計算機(一般簡稱計算機)技術的產生和飛速發(fā)展，為了便于計算機對信號進行處理，產生了在抽樣定理指導下將連續(xù)時間信號變換成離散時間信號的完整的理論和方法。也就是說，可以只用原模擬信號在一系列離散時間坐標點上的樣本值表達原始信號而不丟失任何信息，波、波形、信號這些概念既然表達的是客觀世界中各種物理量的變化，自然就是現代社會賴以生存的各種信息的載體。信息需要傳播，靠的就是波形信號的傳遞。信號在它的產生、轉換、傳輸的每一個環(huán)節(jié)都可能由于環(huán)境和干擾的存在而畸變，有時，甚至是在相當多的情況下，這種畸變還很嚴重，以致信號及其所攜帶的信息被深深地埋在噪聲當中了。濾波，本質上是從被噪聲畸變和污染了的信號中提

9、取原始信號所攜帶的信息的過程。濾波器特性可以用其頻率響應來描述，按其特性的不同，可以分為低通濾波器，高通濾波器，帶通濾波器和帶阻濾波器等。帶通濾波器（英語：Band-pass filter）是指能通過某一頻率范圍內的頻率分量、但將其他范圍的頻率分量衰減到極低水平的濾波器，與帶阻濾波器的概念相對。一個模擬帶通濾波器的例子是電阻-電感-電容電路(RLC circuit)。這些濾波器也可以用低通濾波器同高通濾波器組合來產生。2.3 FIR濾波器原理FIR濾波器(有限長單位沖激響應濾波器)是在數字信號處理（DSP）中經常使用的兩種基本的濾波器之一，另一個為IIR濾波器。IIR濾波器是無限沖激響應濾波器

10、。不論哪一種濾波器設計方法，都要求出濾波器的單位沖激響應h(n)，然后才能在時域中實現頻域中的濾波。在頻域，當其輸入信號為X(ej)時，如濾波器的頻率響應為H(ej)，則其輸出信號為Y(ej)=X(ej)H(ej)。在時域，設濾波器的單位沖激響應h(n)為一N點序列，即0nN-1時h(n)的值不為零，根據離散傅氏變換的性質，則可以將濾波器的輸入序列x（n）的響應y（n）表示為x（n）與h(n)的卷積和，即：這就是濾波系統(tǒng)的差分方程，它給濾波器的實現奠定了理論基礎。即求出時域的h(n)后，便可通過卷積來實現頻域的濾波。卷積和運算主要有以下幾個步驟：(1) h(n)序列N個點數值的存儲由于h(n)

11、是根據濾波性能要求已經設計好的有限長單位沖激響應，故其N個點的數值是已知的，因此可以存放在ROM或RAM當中，且對應著N個不同的地址，便于尋址。(2) 輸入序列x（n）的移位寄存輸入序列x（n）是不斷變化的，因此只能對其進行移位寄存，寄存器的個數為N，即N個寄存器中分別存放著x（n）、x（n-1）x（n-N+1），它們都隨著n的變化而變化。(3) 乘法器用以完成兩個數值的乘法，即h(m)x（n-m），也就是將存儲器中N地址所對應的N個固定數值h(m)分別與N個移位寄存器中的不斷變化的N個變化數值x（n-m）相乘。(4) 累加器用以實現N個乘積的累加，即將當前x（n）所對應的N個乘積進行累加，所

12、得到的和就是y（n）。當濾波器的下一個輸入值即x（n +1）到來時，累加器清零，并重新將下一組x（n +1）所對應的N個乘積進行累加，所得到的和就是y（n +1）。2.4 窗函數簡介數字信號處理的主要數學工具是傅里葉變換。而傅里葉變換是研究整個時間域和頻率域的關系。不過，當運用計算機實現工程測試信號處理時，不可能對無限長的信號進行測量和運算，而是取其有限的時間片段進行分析。做法是從信號中截取一個時間片段，然后用觀察的信號時間片段進行周期延拓處理，得到虛擬的無限長的信號，然后就可以對信號進行傅里葉變換、相關分析等數學處理。無線長的信號被截斷以后，其頻譜發(fā)生了畸變，原來集中在f(0)處的能量被分散

13、到兩個較寬的頻帶中去了（這種現象稱之為頻譜能量泄漏）。為了減少頻譜能量泄漏，可采用不同的截取函數對信號進行截短，截斷函數稱為窗函數，簡稱為窗。信號截短以后產生的能量泄漏現象是必然的，因為窗函數w(t)是一個頻帶無限的函數，所以即使原信號x(t)是有限帶寬信號，而在截短以后也必然成為無限帶寬的函數，即信號在頻域的能量與分布被擴展了。又從采樣定理可知，無論采樣頻率多高，只要信號一經截短，就不可避免地引起混疊，因此信號截短必然導致一些誤差。泄漏與窗函數頻譜的兩側旁瓣有關，如果兩側瓣的高度趨于零，而使能量相對集中在主瓣，就可以較為接近于真實的頻譜，為此，在時間域中可采用不同的窗函數來截短信號。3 FI

14、R濾波器的Matlab設計3.1. 總體方案設計本次課程設計所采用的主要軟件為Matlab。設計中用到的函數主要有wavread,sound,fft,subplot,plot,fir1,freqz,filter。課程設計利用窗函數法設計FIR數字濾波器，并利用MATLAB作為輔助工具仿真濾波前后的時域波形和頻譜。對于DSP的處理，首先，通過Matlab給語音信號添加高頻噪音，即將原始信號調制到某一高頻上（本設計用10000hz），再將其與原始語音信號相疊加，最后將這個疊加后的信號轉化為數據（matlab編程），并存入后綴為dat的文件中。之后根據仿真結果使用MATLAB的設計分析工具FDATO

15、OL產生濾波器參數。將混疊后的信號數據和濾波器參數輸入CCS進行DSP編程，在DSP中實現FIR帶通濾波，并使用CCS的頻譜分析功能，查看最終DSP的濾波效果。3.2 語音信號的采集使用PC 機上的聲卡和WINDOWS 操作系統(tǒng)可以進行數字信號的采集。將話筒插入計算機的語音輸入插口上,啟動MOOo錄音專家。按下錄音按鈕,接著對話筒錄音,錄音結束后屏幕左側將顯示所錄聲音的長度。點擊放音按鈕,可以實現所錄音的重現。如圖3.1所示，以文件名“voice”保存入matlab當前工作路徑中?？梢钥吹?文件存儲器的后綴默認為.wav ,這是WINDOWS 操作系統(tǒng)規(guī)定的聲音文件存的標準。3.3 matla

16、b具體處理初始化讀入語音信號原始音頻抽樣加噪聲并帶通濾波導出混合信號供DSP處理3.3.1 原始音頻抽樣Fs=10000;bits=16;T=1;ts=1/Fs;N=38;x,Fs,bits=wavread('voice.wav'); %x=x(:,1); m=length(x);%sound(x,Fs,bits);y=fft(x,m); %對錄制的信號進行傅立葉變換f=(Fs/m)*1:m;t=1:m/Fs;figure(1);subplot(211);plot(t,x);title('原始信號波形');xlabel('time(s)');su

17、bplot(212);plot(f,abs(y);title('原始信號頻譜');xlabel('frequency(hz)');noise4=0.01*sin(2*pi*5000*t)+0.01*sin(2*pi*200*t); 帶通噪聲%3.3.2 加噪聲并帶通濾波fpl4=500;fph4=2000;fsl4=700;fsh4=1800;wp4=2*fpl4/Fs 2*fph4/Fs;ws4=2*fsl4/Fs 2*fsh4/Fs;X4=x+noise4'%加了帶通噪聲的信號%sound(X4,Fs,bits);%回放噪聲信號Y41=fft(X4,

18、m);%對加了帶通噪聲的信號進行傅立葉變換figure(2);subplot(211);plot(t,X4);title('加了帶通噪聲的信號時域圖');xlabel('time(s)');subplot(212);plot(f,abs(Y41);title('加了帶通噪聲的信號頻譜圖');xlabel('frequency(hz)');figure(3);wc4=(wp4+ws4)/2;b4=fir1(N,wc4,kaiser(N+1);%調用kaiser計算帶通FIRDF的b4（n）hn4 wp4=freqz(b4,1,512

19、);subplot(211);plot(wp4/pi,20*log10(abs(hn4);title('帶通濾波器幅頻');xlabel('w');subplot(212);plot(wp4/pi,unwrap(angle(hn4);title('帶通濾波器相頻');xlabel('w');figure(4);x4get=filter(b4,1,X4);Y4get=fft(x4get,m);%sound(x4get,Fs,bits);subplot(211);plot(t,x4get);title('加了帶通噪聲的信號濾波

20、后時域圖');xlabel('time(s)');subplot(212);plot(f,abs(Y4get);title('加了帶通噪聲的信號濾波后頻譜圖');xlabel('frequency(hz)');3.3.3 導出混合信號供DSP處理xout=X4/max(X4);xto_ccs=round(32767*xout);%由于TMS32054X是16位定點處理器，在進行匯編程序設計時，FIR濾波器系數需要采用Q15格式。 myfile=fopen('fir1.dat','w');%新建fprintf

21、(myfile,'1651 1 e5 1 12cn');fprintf(myfile,'%dn',xto_ccs);%把調制信號+原始信號輸出fir.dat文件fclose(myfile);3.5 FIR濾波器的系統(tǒng)參數在本設計中，使用MATLAB模擬產生合成信號，然后利用CCS進行濾波。因此設計一個FIR帶通濾波器，其參數為：采樣頻率： Fs=44100Hz通帶1：Fp=700Hz阻帶1： Fc=500Hz通帶2： Fp=1800Hz阻帶2： Fp=2000Hz通帶最大衰減： 1dB阻帶最少衰減： 100dB濾波器級數： N=38（由MATLAB根據前述參數求

22、得。）其中，采樣頻率由y,fs,bits=wavread()求出。其中fs為采樣頻率。執(zhí)行程序可得fs為22050HZ。濾波器級數由%據*；c*c，（帶寬）c，（窗口長度）（），（階數）；算出階數為38。濾波器系數使用matlab濾波器的設計/分析工具FDATOOL設計FIR濾波器。步驟如下圖所示：圖 1 Matlab 聲音信號調制示意圖4 FIR濾波器的檢驗圖 2 Matlab濾波器設計在選項中選擇或輸入濾波器參數，然后點擊”Design Filter”按鈕完成濾波器設計。設計成功后的結果如圖下圖所示。在fdatool中，選擇Targets->Code Composer StudioI

23、DE。在出現的對話框中選擇輸出文件類型為C.header file,輸出系數類型為signed 16-bit integer，如圖1113。點擊Generate按鈕，選擇路徑，即可輸出前一步設計出的FIR濾波器系數表。4 FIR濾波器的CCS設計與仿真CCS是一種針對TMS320系列DSP的集成開發(fā)環(huán)境，在windows操作系統(tǒng)下，采用圖形接口界面，提供環(huán)境配置、源文件編輯、程序調試、跟蹤和分析等工具，提供了配置、建立、調試、跟蹤和分析程序的工具，從而完成編輯、編譯、鏈接、調試和數據分析等工作。在CCS中編寫匯編語言程序，進行調試，實現FIR濾波的功能。本課程設計使用CCS開發(fā)應用程序的一般步

24、驟為： （1）創(chuàng)建一個工程項目文件fir2.h，用匯編語言編寫處理主程序fir2.asm。另外根據板上的存儲器配置方式,編寫存儲器配置文件fir d文件,將MATLAB生成的fir2.h和indata.dat文件拷貝到工程文件夾下。 （2）編輯各類文件。使用CCS提供的集成編輯環(huán)境，對頭文件、鏈接命令文件和源程序進行編輯。 （3）對工程項目進行編譯，生成可執(zhí)行文件。如出現語法錯誤，將在構建窗口中顯示錯誤信息。用戶可以根據顯示信息定位錯誤信息，更改錯誤。 （4）下載程序、輸入數據，執(zhí)行程序，對結果和數據進行分析和算法評估。利用CCS提供的探測點、圖形顯示等工具，對運行結果、輸出數據進行分析，評估

25、算法性能。4.1 CCS程序流程圖圖 3 CCS匯編程序流程圖開始DSP初始化將濾波系數分別存放在存儲單元bn中從第L個數開始讀入抽樣值抽樣值存入在存儲單元xn中將bn中和xn中對應的值相乘累加bn*xn+ACC=>ACCL=L+1輸出結果4.1.1 程序分析 FIR濾波器的DSP實現主要有線性緩沖區(qū)法和循環(huán)緩沖區(qū)法兩種。 線性緩沖區(qū)法又稱延遲線法。其特點： （1）對于N級的FIR濾波器，在數據存儲器中開辟一個N單元的緩沖區(qū)（滑窗），用來存放最新的N個輸入樣本； （2）從最老樣本開始取數，每取一個樣本后，將此樣本向下移位； （3）讀完最后一個樣本后，輸入最新樣本存入緩沖區(qū)的頂部。 循環(huán)緩

26、沖區(qū)法的特點如下： （1）對于N級FIR濾波器，在數據存儲器中開辟一個N單元的緩沖區(qū)（滑窗），用來存放最新的N個輸入樣本； （2）從最新樣本開始取數； （3）讀完最后一個樣本（最老樣本）后，輸入最新樣本來代替最老樣本，而其他數據位置不變； （4）用片內BK（循環(huán)緩沖區(qū)長度）寄存器對緩沖區(qū)進行間接尋址，使循環(huán)緩沖區(qū)地址首尾相鄰。 本次設計的FIR濾波器所采用的就是循環(huán)緩沖區(qū)法。對DSP進行初始化，并設置好存儲單元。為數字濾波作準備，將預先設計好的有N個抽頭的FIR數字濾波器的沖激響應序列h(n)中的N個數值放入存儲單元bn中，從indata數據存儲段中讀入抽樣值，將抽樣值逐一放入存儲單元xn中。

27、主要循環(huán)程序如下：STM #-1,AR0RPTZ A,# ;將ACC累加器清零，循環(huán)次MAC *AR2+0%,*AR3+0%,ASTH A,*(yn) ;保存輸出 MVKD *(yn),*AR5+將累加器清零，將存儲單元bn與xn值對應相乘并累加到ACC中，循環(huán)次將ACC中的值保存輸出。再從indata數據存儲段中讀入一個新的數據，代替原來的最老樣本 MVDD *AR4+,*AR2+0% ,重復循環(huán)直到數據全部處理完。最后將處理完的數據輸入到outdata數據存儲段中。4.1.2 數據輸入MATLAB輔助DSP 實現FIR濾波器濾波 ,其總體過程是為DSP 中編寫處理程序提供濾波器系數。在MA

28、TLAB中利用濾波器設計、分析工具( FDATOOL) ,根據指定的濾波器性能快速設計一個FIR ,然后把濾波器系數以頭文件形式導入CCS 中,頭文件中包含MATLAB輔助DSP 實現FIR 數字濾波器的階數和系數數組,在MATLAB中調試、運行DSP 程序并顯示、分析處理后的數據。打開CCS IDE 中file的data加載初始化data數據，在Address中輸入數據存放的起始地址indata，并設置數據棧長度。4.1.3 運行程序仿真運行程序，查看輸入輸出波形，修改相應參數進行調試。在View的Graph中單擊Time/frequency出現graph property dialog框。

29、將顯示類型，圖形名稱，起始地址，抽樣點數，數據類型等分別進行設置，輸出各種波形。輸出濾波前的信號時域波形。首先，設置如下圖所示，輸出濾波前的信號時域波形。圖 4 濾波前語音信號的時域Graph Property Dialog設置參數圖將Display Type改為FFT Magnitude 就可以看到濾波前信號的頻譜，將Start Address改為outdata就可以看到濾波后信號的時域波形和頻譜。4.2 CCS仿真結果圖 5 濾波前語音信號的時域波形圖圖 6 濾波前信號頻譜波形圖4.2.2 FIR濾波器濾波后的信號圖 7 濾波后信號時域波形圖圖 8 濾波后信號頻譜圖圖17波前的信號，由高頻

30、信號和原信號疊加而成的波形較為雜亂，加噪后的信號幅值較大。圖19波后的信號，從圖上可以看到信號幅度較低，波形較為清晰。對比圖18圖20可以看到加躁后信號幅值較大，信號的能量集中在高頻段，濾波后，高頻成分大部分被濾去，信號的能量主要集中在低頻段。5 問題及解決方法1.為什么存入DAT文件的信號數據要乘上32768？由于TMS32054X是16位定點處理器，在進行匯編程序設計時，FIR濾波器系數需要采用Q15格式，即必須將上述系數據轉化位Q15定點式。 2.matlab編譯錯誤出現如下錯誤？The carrier frequency must be less than half the sampl

31、ing frequency.故降低載波頻率3. 在View的Graph中單擊Time/frequency出現graph property dialog框，未修改抽樣點數，顯示的圖形出現差異.4.在CCS編譯時，發(fā)現fir.h文件老是出錯，最后發(fā)現濾波系數要以整數體現，而且不能加分號。6 設計感想通過本次課程設計，我們發(fā)現在課本上學到的知識實在是不能使我們有能力駕馭這個課程設計。一開始決定做這個課題，我們就去查找了大量資料，通過查找資料我們都清楚了大概要做些什么。首先要進行MATLAB仿真，并通過MATLAB設計的FIR濾波器輸出濾波器系數。接著便是通過CCS編程，并在DSP上運行。我們覺得思路

32、是很清晰的，但是實際操作起來卻困難重重。原始信號有了，我們也得了教訓，便開始計劃開來。于是，翻書看完有關FIR濾波器的相關介紹，重點看了凱塞窗。然后查資料找出與本次設計有關的MATLAB知識。做完這些自然是開始設計，總體框架，具體步驟，都在腦袋瓜里轉了幾轉。動工，這時才深刻體會到“紙上學來終覺淺，用到實際才叫真！”開始設計從噪聲選取，到低通，高通，帶阻，帶通濾波器的設計，其中是錯漏百出。當然，作為老手自然不會再去犯什么程序設計上的低級錯誤，主要還是對濾波器的認知不夠深刻。有的只是照抄書本，最后是連自己都為自己所犯的錯誤驚愕不已。其中許多錯誤差點讓自己從頭來過，一度讓自己有放棄的打算。最后還是堅

33、持了下來，自己大嘆不容易。設計的時間過了，但它帶來的影響我想是不滅的！要有一個良好的就業(yè)前景。擁有扎實的理論知識，認真細致的處事態(tài)度，隨機應變的實際應用方法，不拘一格的思路都將使得自己能有好的未來。參 考 文 獻1 鄒彥，DSP原理及應用（修訂版），北京：電子工業(yè)出版社，20212 王宏. MATLAB6.5及其在信號處理中的應用.M. 北京:清華大學出版社,20043 高西全,丁玉美,闊永紅.數字信號處理-原理、實現及應用.M.北京:電子工業(yè)出版社,20064 胡廣書編著，數字信號處理理論、算法與實現.北京：清華大學出版社，2005；附 錄.Matlab源程序%原始音頻抽樣%Fs=10000

34、;bits=16;T=1;ts=1/Fs;N=38;x,Fs,bits=wavread('voice.wav'); %利用Windows下的錄音機，錄制一段自己的語音信號，時間為1sx=x(:,1); m=length(x);%sound(x,Fs,bits);y=fft(x,m); %對錄制的信號進行傅立葉變換f=(Fs/m)*1:m;t=1:m/Fs;figure(1);subplot(211);plot(t,x);title('原始信號波形');xlabel('time(s)');subplot(212);plot(f,abs(y);tit

35、le('原始信號頻譜');xlabel('frequency(hz)');noise4=0.01*sin(2*pi*5000*t)+0.01*sin(2*pi*200*t); %帶通噪聲%帶通濾波%fpl4=500;fph4=2000;fsl4=700;fsh4=1800;wp4=2*fpl4/Fs 2*fph4/Fs;ws4=2*fsl4/Fs 2*fsh4/Fs;X4=x+noise4'%加了帶通噪聲的信號%sound(X4,Fs,bits);%回放噪聲信號Y41=fft(X4,m);%對加了帶通噪聲的信號進行傅立葉變換figure(2);subpl

36、ot(211);plot(t,X4);title('加了帶通噪聲的信號時域圖');xlabel('time(s)');subplot(212);plot(f,abs(Y41);title('加了帶通噪聲的信號頻譜圖');xlabel('frequency(hz)');figure(3);wc4=(wp4+ws4)/2;b4=fir1(N,wc4,kaiser(N+1);%調用kaiser計算帶通FIRDF的b4（n）hn4 wp4=freqz(b4,1,512);subplot(211);plot(wp4/pi,20*log10(

37、abs(hn4);title('帶通濾波器幅頻');xlabel('w');subplot(212);plot(wp4/pi,unwrap(angle(hn4);title('帶通濾波器相頻');xlabel('w');figure(4);x4get=filter(b4,1,X4);Y4get=fft(x4get,m);%sound(x4get,Fs,bits);subplot(211);plot(t,x4get);title('加了帶通噪聲的信號濾波后時域圖');xlabel('time(s)');subplot(212);plot(f,abs(Y4get);title('加了帶通噪聲的信號濾波后頻譜圖');xlabel('frequency(hz)');%導出混合信號供DSP處理%xout=X4/max(X4);xto_ccs=round(32767*xout);%由于TMS32054X是16位定點處理器，在進行匯編程序設計時，FIR濾波器系數需要采用Q15格式。 myfile=fopen('fir1.dat','w')