【正文】 。 pha=angle(H)。 mag=abs(H)。 H=(H(1:1:501))39。whole39。 N=1000。 hd=sin(wc*m)./(pi*m)。 n=[0:(M1)]。) 程序 5 function hd=ideal(wc,M)。 26 ylabel(39。頻率（單位 :pi)39。 %繪制相頻特性 axis([0 2 4 4])。相位響應 39。 plot(w,pha)。YTick39。manual39。YTickMode39。XTick39。manual39。XTickMode39。分貝數(shù) 39。)。 xlabel(39。 grid。幅度響應（單位： dB)39。 %繪制幅度響應 。) subplot(2,2,3)。 ylabel(39。n39。) %繪制實際脈沖響應 axis([0 N1 ])。 title(39。) subplot(2,2,2)。 text(N+1,0,39。w(n)39。) %繪制布萊克曼窗 axis([0 N1 0 ])。 title(39。 %頻率分辨率 Rp=(min(db(1:wp/dw+1))) %檢驗通帶波動 As=round(max(db(wst/dw+1:501))) %檢驗最小阻帶衰減 figure(1) subplot(2,2,1)。 %設計的脈沖響應應為理想脈沖響 應與窗函數(shù)乘積 [db,mag,pha,grd,w]=myfreqz(h,[1])。 %求理想脈沖響應 wdbla=(blackman(N))39。 wc=(wst+wp)/2。 %按布萊克曼窗計算所需的濾波器長度 N=N0+mod(N0+1,2)。 deltaw=wstwp。 } 程序 4 wst=*pi。 *TMCR_MGS3 =0x510。 } void TMCR_reset( void ) { ioport unsigned int *TMCR_MGS3=(unsigned int *)0x07FE。 clkmd=(unsigned int *)0x1c00。//0x2413。 *clkmd =0x2033。 } 程序 3 include void CLK_init() { ioport unsigned int *clkmd。 *sdc2 = 0X38F。 *emirst = 0。 *ce21 = 0x1fff。 *ce01 = 0X3000。 *egcr = 0x220。 //*ebsr = 0x221。 ioport unsigned int *init =(unsigned int *)0x812。 //ioport unsigned int *sdper =(unsigned int *)0x810。 //ioport unsigned int *ce33 =(unsigned int *)0x80E。 ioport unsigned int *ce31 =(unsigned int *)0x80C。 //ioport unsigned int *ce22 =(unsigned int *)0x80A。 //ioport unsigned int *ce13 =(unsigned int *)0x808。 ioport unsigned int *ce11 =(unsigned int *)0x806。 //ioport unsigned int *ce02 =(unsigned int *)0x804。 //ioport unsigned int *emibe =(unsigned int *)0x802。 ioport unsigned int *egcr =(unsigned int *)0x800。 } return(fSum)。iFIRNUMBER。 23 fSum=0。 return(fXn[0])。 fSignal2+=fStepSignal2。 fSignal1+=fStepSignal1。i ) fXn[i]=fXn[i1]。 } } } float InputWave() { for ( i=FIRNUMBER1。 nOut++。 fOutput=FIR()。 nIn++。 while ( 1 ) { fInput=InputWave()。 fStepSignal1=2*PI/30。 22 fSignal1=。 nOut=0。 int nIn,nOut。 int i。 float fStepSignal1,fStepSignal2。 float fInput,fOutput。 float fHn[FIRNUMBER]={ , , , , }。 參考文獻 [1] 程 佩青著 .數(shù)字信號處理教程 .[M]清華大學出版社 .20xx. [2] 趙順珍，馬英 .基于 DSP的 FIR數(shù)字濾波器設計與實現(xiàn) [J]. [3] 王獻峰，石東 .基于 DSP 的 FIR 數(shù)字濾波器的設計與實現(xiàn) . 天津科技大學電子信息與自動化學院，天津 300222. [4] 冉茂華 .基于 DSP的 FIR數(shù)字濾波器設計 ,武漢理工大學 .20xx. [5] 劉萬松 .基于 DSP的 FIR數(shù)字濾波器設計，貴州大學 .20xx. [6] 孫宗瀛 . DSP 原理設計與應用 [M]. 北京：清華大學出版社， 20xx. [7] 王獻峰，石東 . 基于 DSP 的 FIR 數(shù)字濾波器的設計與實現(xiàn) . 天津科技大學電子信息與自動化學院，天津 300222. [8]張雄偉 ,曹鐵勇 . DSP 芯片的原理與開發(fā)應用 . 電子工業(yè)出版社 . 20xx. 21 附錄 程序 1 include include include include define FIRNUMBER 25 define SIGNAL1F 1000 define SIGNAL2F 4500 define SAMPLEF 10000 define PI float InputWave()。對于濾波，只在軟件上進行。希望以后能對數(shù)字信號處理應用的更廣泛。而 MATLAB在編程上更靈活，功能更強大。如何通過軟件來實現(xiàn)低通濾波器設計，而之前經(jīng)常用硬件來實現(xiàn)濾波器的設計，硬件實現(xiàn)的低通濾波器濾波單一，無法改變，硬件設計一旦錯了就會造成無法挽回和經(jīng)濟的浪費，結(jié)果又得重新 開始浪費時間。 MATLAB的結(jié)果分析，在幅度響 應和頻率的圖形中，在 3dB時，對應的頻率是 ，根據(jù)計算器截止頻率為 147500HZ，相對于截止頻率高了 ，其結(jié)果不是很理想。 CCS的結(jié)果分析如下：當調(diào)用一個低頻的正弦信號和一個高頻的余弦信號，低頻信號與高頻信號的疊加就使得低頻信號的正弦上產(chǎn)生毛刺，通過 FIR濾波器的設計要求編寫的 C語言程序完成對高頻信號的過濾，最終得到光滑的低頻的正弦信號。 圖 15 INPUT 圖 16 OUTPUT 下面是 FFT的濾波結(jié)果如圖 17和輸入波形和圖 18 FFT輸出波形。 圖 14 觀察設置窗口 點擊 OK完成設置。 （ 4）下載程序：執(zhí)行 FileLoad Program，然后再打開的提示框中選擇剛剛建立 C:\ICETEKVC5509EDULab\Lab0501UseCCS\Debug\ 文件。如圖 12 圖 12 CCS編 譯圖 （ 2）現(xiàn)在就編譯：選擇菜單中的“ Project”的“ Rebuild All”進行編譯，編譯完 18 成后會在 CCS下方的提示框中顯示編譯信息，根據(jù)編譯的提示，修改程序。 （ 1） 首先新建一個工程并且命名為 FIR，此時該工程下的文件為空。 /* break point */ if ( nOut=256 ) { nOut=0。 //濾波器函數(shù)調(diào)用 fOut[nOut]=fOutput。 nIn%=256。 //輸入函數(shù)的調(diào)用 fIn[nIn]=fInput。 return(fXn[0])。 fSignal2+=fStepSignal2。 fSignal1+=fStepSignal1。i ) fXn[i]=fXn[i1]。 } 輸入子函數(shù) 17 float InputWave() { for ( i=FIRNUMBER1。i++ ) { fSum+=(fXn[i]*fHn[i])。 for ( i=0。 子函數(shù)如 float FIR() //濾波器函數(shù) { float fSum。 由于窗函數(shù)的設計是由理想的 )(ndh 和布拉克曼窗函數(shù) )(nw 相乘得到的低通濾波器函數(shù) 11,0),()()( ??? Nnnwndhnh ?，這個函數(shù)包括主函數(shù)和兩個子函數(shù)，子函數(shù)一個是由 )(ndh 和 )(nw 相乘的函數(shù)。子函數(shù)主要是根據(jù) FIR低通濾波器的算法進行編程的，子函數(shù)是由 )(ndh 和 )(n? 進行相乘得 到的，函數(shù)為fSum+=(fXn[i]*fHn[i]其函數(shù)流程圖為圖 12。 16 CCS實現(xiàn) FIR低通濾波器設計 用 C語言編程實現(xiàn) FIR低通濾波器主要是根據(jù)第二章的窗函數(shù)設計法的思想進行編程。它提供一整套的程序編制、維護、編譯、調(diào)試環(huán)境，能將匯編語言和 C 語言程序編譯連接生成 COFF (公共目標文件 )格式的可執(zhí)行文件，并能將程序下載到目標 DSP 上運行調(diào)試。 圖 11 初始化 按“空格”或任意鍵即完成初始化。 進入設置選擇 C5XX functional simulator 后選擇 Save and quit 后點擊“ NO”完成軟件的設置，設置過程如圖 10。 在使用軟件仿真方式工作時，無需連接板卡和仿真器等硬件。 CCS 提供了配置、建立、調(diào)試、跟蹤和分析程序的工具，它便 于實時、嵌入式信號處理程序的編制和測試，能夠加速開發(fā)進程，提高工作效率。安裝完畢后，在安裝“ CCS 開發(fā)軟件”目錄中的壓縮文件，“ ”中的升級程序，將 升級到 。 DSP 軟件開發(fā)環(huán)境 本次試驗的研究內(nèi)容主要就是在 DSP 芯片的集成軟件開發(fā)環(huán)境 CCS(Code Composer Studio)中調(diào)試從而實現(xiàn)數(shù)字濾波功能。與此同時， DMA 控制器可以執(zhí)行最多兩個數(shù)據(jù)周期獨立于 CPU 的活動。 CPU 支持內(nèi)部 structure 是由一個程序總線， 3 的 數(shù)據(jù)讀總線，兩條數(shù)據(jù)寫總線，并 additional buses 專用于外周和 DMA 活動。 DSP 的硬件結(jié)構(gòu)示意圖 ICETEKVC5509A 評估板接口說明實物圖如圖 8 和 ICETEKVC5509A 評估板原理框圖如圖 9。 14 此框圖是 DSP 實現(xiàn)的主要過程，因為 DSP 處理的是數(shù)字信號，所以首先應該將輸入的前置的模擬信號通過 A/D 變換器轉(zhuǎn)化為數(shù)字信號，并把它交給 DSP 處理， DSP 根據(jù)設計要求經(jīng)過編程處理（如 F