【正文】 DSP課程設(shè)計實驗報告 自動化0804 姚笑菲 08212116DSP課程設(shè)計實驗報告語音噪聲濾波院（系）：電子信息工程學(xué)院 班級： 自動化0804 姓名： 姚笑菲 學(xué)號： 08212116 指導(dǎo)教師： 楊恒 目錄一、實驗背景 3二、設(shè)計要求 3三、設(shè)計思路 3四、算法原理 4直接存儲器訪問DMA 4A/D和D/A轉(zhuǎn)換器 6FIR濾波器算法 6LMS算法 7利用DSP產(chǎn)生回波 8利用自適應(yīng)濾波實現(xiàn)回波對消 9五、用C語言實現(xiàn)LMS算法 11六、用MATLAB設(shè)計FIR濾波器 15七、程序設(shè)計 17主程序（） 17中斷服務(wù)程序（） 22鏈接命令文件（） 24中斷向量表（） 26八、程序調(diào)試 30九、結(jié)果分析 31十、感想 40一、實驗背景我們的日常生活中離不開語音通信，如接打電話、收聽音樂等。語音通信的過程主要分為三個部分，即發(fā)送部分——將音頻信號轉(zhuǎn)化為電信號，經(jīng)發(fā)送端設(shè)備變換為適合傳輸?shù)男问桨l(fā)送出去；傳輸部分——通過傳輸信道將電信號進行傳輸；接收部分——經(jīng)接收端設(shè)備恢復(fù)出原來的語音信號，經(jīng)耳機或者喇叭轉(zhuǎn)換為接收者可以聽到的聲音信號。但是由于環(huán)境的原因，我們采集到語音信號經(jīng)常含有不同程度的噪聲。與語音通信系統(tǒng)的三個部分相對應(yīng)，語音通信系統(tǒng)中的噪聲一般也來自三個方面：一是信號發(fā)送端空間環(huán)境中的音頻噪聲信號經(jīng)麥克風(fēng)變換為電信號之后，與有用信號其同傳遞到接收端；二是信號接收端空間環(huán)境中的音頻噪聲對信號接收者的影響；三是信號處理設(shè)備產(chǎn)生的電噪聲及傳輸信道中的電噪聲。當(dāng)語音信號受到背景噪聲干擾而使語音通信質(zhì)量變得不可接受時，要對語音信號中的噪聲進行濾除，也就是本次實驗要完成的語音噪聲濾波。二、設(shè)計要求基本部分：（1） 對DMA進行初始化；（2） 對A/D、D/A進行初始化；（3） 編寫DMA中斷服務(wù)程序和濾波算法程序，實現(xiàn)語音信號的實時濾波；發(fā)揮部分：（1） 使用DSP產(chǎn)生帶回波的語音信號；（2） 利用自適應(yīng)濾波實現(xiàn)語音信號的回波對消。（3） 比較采用不同窗和階數(shù)時濾波器的濾波效果；（4） 測試所設(shè)計濾波器的幅頻特性和相頻特性，并與MATLAB下的設(shè)計結(jié)果進行比較。三、設(shè)計思路首先利用DSP的DMA方式對外部含噪聲的語音信號進行實時采集，語音信號先經(jīng)過A/D轉(zhuǎn)換為數(shù)字信號，利用MCBSP的接收寄存器接收數(shù)據(jù)。編寫濾波算法程序，或調(diào)用DSPLIB中的濾波函數(shù)，對信號進行濾波。濾波后的數(shù)據(jù)利用DMA方式送到D/A轉(zhuǎn)換器轉(zhuǎn)換為模擬信號。設(shè)計流程框圖如下：上電復(fù)位A/D和D/A初始化DMA通道初始化DMA方式接收A/D數(shù)據(jù)等待DMA接收中斷信號濾波信號經(jīng)D/A輸出主程序中斷服務(wù)程序模式選擇模式0普通濾波模式1自適應(yīng)濾波加有噪聲色語音信號由DSK板的MIC端口輸入，經(jīng)A/D轉(zhuǎn)換器轉(zhuǎn)換為數(shù)字信號后到達DSP的DMA通道2，DSP接收到信號后調(diào)用中斷服務(wù)程序進行濾波，可選擇濾波模式進行不同種類的濾波，濾波后的信號在CCS上顯示，或通過DMA通道3將信號輸出到D/A轉(zhuǎn)換器，轉(zhuǎn)換為模擬信號后從SPEAKER端口輸出。四、算法原理直接存儲器訪問DMA直接存儲器訪問（Direct Memory Access）是C54x DSP非常重要的片上外設(shè)，DMA控制器可以在不影響CPU的情況下完成數(shù)據(jù)的傳輸，因此數(shù)據(jù)傳輸速度快，在要求信號實時采集和處理的系統(tǒng)中常采用DMA方式進行信號采集和傳輸。C5402有6個可獨立編程的DMA通道，每個DMA通道受各自的5個16位寄存器控制：源地址寄存器DMSRC、目的地址寄存器DMDST、單元計數(shù)寄存器DMCTR、同步事件和計數(shù)寄存器DMSFC、發(fā)送模式控制寄存器DMMCR。目的地址寄存器DMDST規(guī)定DMA要傳送數(shù)據(jù)目的地址和首地址。單元計數(shù)寄存器DMCTR規(guī)定DMA傳送數(shù)據(jù)的個數(shù)為DMCTR寄存器的值加1。同步事件和幀計數(shù)寄存器DMSFC規(guī)定DMA傳送數(shù)據(jù)的同步事件類型和傳送一塊數(shù)據(jù)所含幀信號的個數(shù)。發(fā)送模式控制寄存器DMMCR規(guī)定DMA通道的傳輸模式，當(dāng)DMA工作在自動初始化模式時，CPU在一個DMA事件完成后自動裝載下一個DMA初始化設(shè)置并繼續(xù)進行數(shù)據(jù)傳送。DINM位和IMOD位設(shè)置DMA中斷產(chǎn)生方式。CTMOD位設(shè)置發(fā)送計數(shù)模式控制，CTMOD位為0時，DMA工作在多幀模式，CTMOD位為1時，DMA工作在自動緩沖ABU模式。SIND位和DIND位用來設(shè)置源地址和目的地址的自動調(diào)整方式。此外，DMA的6個通道還受通道優(yōu)先級和使能控制寄存器DMPREC控制。本實驗利用DMA通道2與MCBSP1通道結(jié)合來讀取A/D轉(zhuǎn)換器的數(shù)據(jù)，利用DMA通道3與MCBSP1通道結(jié)合來將處理后的數(shù)據(jù)發(fā)送至D/A轉(zhuǎn)換器。選擇MCBSP1通道的接收寄存器DRR11（41H）為DMA傳送數(shù)據(jù)的首地址，并選擇DMA源地址工作在不調(diào)整方式，選擇MCBSP1接收事件為DMA同步事件，以實現(xiàn)DMA和MCBSP的結(jié)合。A/D轉(zhuǎn)換器的數(shù)據(jù)按MCBSP1的設(shè)置被傳送到DSP內(nèi)部接收寄存器DRR11中，再由DMA通道2將DRR11中的數(shù)據(jù)讀到指定數(shù)據(jù)存儲區(qū)INPBUFFER中完成數(shù)據(jù)采集。DMA在傳送外部來的數(shù)據(jù)時不會影響CPU的正常運行，當(dāng)DMA通道2采集完一組數(shù)據(jù)后產(chǎn)生一個DMA中斷事件中斷CPU，來通知CPU對其進行相應(yīng)的處理，此時DMA通道2可以按照設(shè)定繼續(xù)采集下一組數(shù)據(jù)，實現(xiàn)了數(shù)據(jù)采集與CPU處理的并行操作。當(dāng)一組數(shù)據(jù)處理完成后，將數(shù)據(jù)存放在數(shù)據(jù)存儲區(qū)OUTBUFFER中，選擇OUTBUFFER為DMA通道3的首地址，并選擇源工作地址工作在不調(diào)整方式。選擇MCBSP1通道的發(fā)送寄存器DXR11（43H）為DMA通道3傳送數(shù)據(jù)的目的地址，并選擇DMA源地址工作在不調(diào)整方式，選擇MCBSP1接收事件為DMA同步事件，以實現(xiàn)DMA和MCBSP的結(jié)合。在CCS集成開發(fā)環(huán)境中。在這兩個頭文件中，定義了DMA的寄存器資源、使用方法和DMA寄存器的地址和基本訪問方式，以及寄存器的各個比特域和訪問方法，因此需要進行C源程序的開始處包含這兩個頭文件。輸入信號經(jīng)A/D轉(zhuǎn)換后，利用DMA通道2進行傳輸進入DSP處理，而濾波后的信號輸出利用了DMA通道3進行傳輸進入D/A轉(zhuǎn)換器等輸出，因此需要對DMA通道2和通道3進行初始化配置。詳細的DMA初始化程序請見后面的程序清單。A/D和D/A轉(zhuǎn)換器A/D模數(shù)轉(zhuǎn)換器是將輸入電壓信號轉(zhuǎn)換為輸出的數(shù)字信號。由于數(shù)字信號本身不具有實際意義，僅僅表示一個相對大小，所以任何一個模數(shù)轉(zhuǎn)換器都需要一個參考模擬量作為轉(zhuǎn)換的標準，而輸出的數(shù)字量則表示輸入信號相對于參考信號的大小。A/D轉(zhuǎn)換器能夠?qū)⒔邮盏恼Z音信號（模擬信號）轉(zhuǎn)化為數(shù)字信號，供CPU進行處理。D/A轉(zhuǎn)換器是將數(shù)字信號轉(zhuǎn)換為模擬信號的電路。實現(xiàn)原理就是將輸入的二進制代碼按其權(quán)值的大小轉(zhuǎn)換成相應(yīng)的模擬量，然后將所有的模擬量相加，使其與數(shù)字量成正比。D/A轉(zhuǎn)換器能夠?qū)?shù)字信號轉(zhuǎn)換為語音信號送SPEAKER端口輸出。要使用A/D和D/A轉(zhuǎn)換器，必須首先對A/D和D/A轉(zhuǎn)換器進行初始化設(shè)置，即設(shè)置A/D轉(zhuǎn)換器的工作模式（15+1bit模式，16bit模式）、輸入增益（0dB,6dB,12dB）以及抽樣頻率（8000Hz,16000Hz）等。FIR濾波器算法數(shù)字濾波器原理一般具有如下差分方程:y(n)=∑akx(nk)+ ∑bky(nk) （1）式中x(n)為輸人序列，y(n)為輸出序列，ak和bk為濾波器系數(shù)，,N是濾波器階數(shù)。當(dāng)所有的bk均為零時，則有y(n)= ∑akx(nk) （2）(2)式是FIR濾波器的差分方程，其一般形式為y(n)= ∑h(k)x(nk) （3）對(3)式進行z變換，整理后可得FIR濾波器的傳遞函數(shù):H(z)= ∑h(k)zk （4）FIR(Finite Impulse Response)濾波器，即有限長單位沖激響應(yīng)濾波器，是指系統(tǒng)的單位脈沖響應(yīng)h[k]僅在有限范圍內(nèi)有非零值的濾波器。FIR數(shù)字濾波器的設(shè)計方法主要有窗函數(shù)法和頻率抽樣設(shè)計法，其中窗函數(shù)法是基本而有效的設(shè)計方法。下面為FIR濾波器的算法實現(xiàn)方法。FIR濾波為有限沖擊響應(yīng)濾波，其濾波結(jié)構(gòu)是一個分節(jié)的延時線，每節(jié)的輸出加權(quán)累加,得到濾波器的輸出。由前面的分析知，F(xiàn)IR濾波器數(shù)學(xué)上可表示為:y(n)= ∑h(k)x(nk) （5）式（5）中 x[n]為最近(t=nT)的輸人信號，x[nk]是延時了k個取樣周期的輸人信號，h[k]是第k個延時的加權(quán)值(即濾波器系數(shù))，可由MATLAB設(shè)計出來，y[n]是時刻t=nT時濾波器的輸出信號，N是濾波器的階數(shù)也稱濾波器的抽頭數(shù)，為實整數(shù)，且須滿足N*t1/fs，其中fs為采樣頻率。因為該濾波器的沖激響應(yīng)在N個周期后變?yōu)?也就是每次乘加運算都要用到之前N個數(shù)。濾波器系數(shù)可通過MATLAB設(shè)計得到。 LMS算法LMS(leastmeansquare)算法，即最小均方算法，是由Widrow和Hof在1960年創(chuàng)建的，直到現(xiàn)在它還是應(yīng)用最廣泛的自適應(yīng)濾波算法。LMS算法是隨機梯度算法族中的一員。該算法在隨機輸入維納濾波器遞推計算中使用確定性梯度。LMS算法的一個顯著特點就是它的簡單性。此外，它不需要計算有關(guān)的相關(guān)函數(shù)，也不需要矩陣求逆運算。事實上，正是因為LMS算法的簡單性，使得它成為其他自適應(yīng)濾波算法的參照標準。這兩個過程一起工作組成一個反饋環(huán)，如下圖所示。首先有一個橫向濾波器(圍繞它構(gòu)造LMS算法)，該部件的作用在于完成濾波過程。其次，有一個對橫向濾波器抽頭權(quán)值進行自適應(yīng)控制過程的算法。橫向濾波器各部分的細節(jié)如下圖所示。抽頭輸入u(k),u(k1),……,u(kM+1)，其中M+1是延遲單元的個數(shù)，這些輸入張成一個多維空間。相應(yīng)的h0(k), h1 (k), ...， hm1（k)為抽頭權(quán)向量h(k)的元素。通過LMS算法計算這個向量所得的值表示一個估計，當(dāng)?shù)螖?shù)趨于無窮時，該估計的期望值可能接近維納最優(yōu)解,在濾波過程中，期望響應(yīng)d(k)與抽頭輸入向量z(k)一道參與處理。在這種情況下，給定一個輸入，橫向濾波器產(chǎn)生一個輸入，橫向濾波器產(chǎn)生一個輸出y(k)作為期望響應(yīng)d(k)的估計。因此，我們可把估計誤差e(k)定義為期望響應(yīng)與實際濾波器輸出之差，估計誤差e(k)與抽頭輸入向量z(k)都被加到自適應(yīng)部分，因此圍繞抽頭權(quán)值的反饋環(huán)是閉環(huán)的。利用DSP產(chǎn)生回波在打電話的過程中，常常能聽見自己說話的聲音在話筒里重復(fù)，實際上就是自己的聲音泄露到了接收系統(tǒng)中。在一些特定的系統(tǒng)中出現(xiàn)這種現(xiàn)象將會對信號的正常收發(fā)產(chǎn)生不利影響。因此，回聲在上述系統(tǒng)中會嚴重影響了通話的清晰度，多點回聲甚至?xí)鹜ㄓ嵕W(wǎng)之間聲音的振蕩。