dsp課程設(shè)計(jì)報(bào)告--fir濾波器設(shè)計(jì)

1、<p><b>　　DSP</b></p><p><b>　　課程設(shè)計(jì)報(bào)告</b></p><p>　　課程名稱： DSP原理及應(yīng)用</p><p>　　設(shè)計(jì)題目： FIR濾波器設(shè)計(jì)</p><p>　　學(xué) 院： 電氣信息學(xué)院</p><

2、p>　　專業(yè)年級(jí)： 通信工程2011級(jí)</p><p>　　指導(dǎo)教師： </p><p>　　姓 名： </p><p>　　學(xué) 號(hào)： </p><p>　　時(shí) 間： 2014.6.30-2014.7.4</p><p><b>　　目

3、 錄</b></p><p><b>　　引 言 3</b></p><p>　　第一章 課程設(shè)計(jì)要求及流程4</p><p>　　1.1FIR濾波器設(shè)計(jì)要求4</p><p>　　1.2設(shè)計(jì)流程4</p><p>　　第二章課程設(shè)計(jì)實(shí)現(xiàn)4</p>&

4、lt;p>　　2.1功能描述4</p><p>　　2.2算法研究4</p><p>　　2.3參數(shù)計(jì)算5</p><p>　　2.3.1利用Matlab計(jì)算濾波系數(shù)5</p><p>　　2.3.2編寫(xiě)濾波器輸入信號(hào)程序6</p><p>　　2.4編寫(xiě)源程序6</p>&

5、lt;p>　　2.4.1匯編源程序“.asm”文件的編寫(xiě)6</p><p>　　2.4.2鏈接命令文件“.cmd”的編寫(xiě)8</p><p><b>　　2.5調(diào)試過(guò)程8</b></p><p>　　2.5.1調(diào)試前的準(zhǔn)備8</p><p>　　2.5.2.inc和.dat的使用8</p>

6、<p>　　2.5.3編寫(xiě)及編譯程序8</p><p>　　2.5.4觀察點(diǎn)設(shè)置9</p><p>　　2.6運(yùn)行并觀察結(jié)果10</p><p>　　第三章課程設(shè)計(jì)總結(jié)12</p><p><b>　　參考文獻(xiàn)12</b></p><p><b>　　引 言

7、</b></p><p>　　隨著信息與數(shù)字技術(shù)的發(fā)展，數(shù)字信號(hào)處理已經(jīng)成為當(dāng)今極其重要而學(xué)科與技術(shù)領(lǐng)域之一。它在通信、語(yǔ)音、圖像、自動(dòng)控制、雷達(dá)、軍事、航空航天、醫(yī)療和家用電器等眾多領(lǐng)域得到了廣泛的應(yīng)用。在數(shù)字信號(hào)處理的基本方法中，通常會(huì)涉及到變換、濾波、頻譜分析、調(diào)制解調(diào)和編碼解碼等處理。其中濾波是應(yīng)用非常廣泛的一個(gè)環(huán)節(jié)，數(shù)字濾波器的理論和相關(guān)設(shè)計(jì)也一直都是人們研究的重點(diǎn)之一。FIR濾波器的是非遞

8、歸的，穩(wěn)定性好，精度高；更重要的是，F(xiàn)IR濾波器在滿足幅頻響應(yīng)要求的同時(shí)，可以獲得嚴(yán)格的線性相位特征。因此，它在高保真的信號(hào)處理，如數(shù)字音頻、圖像處理、數(shù)據(jù)傳輸和生物醫(yī)學(xué)等領(lǐng)域得到廣泛應(yīng)用。</p><p>　　在數(shù)字信號(hào)處理中，濾波占有極其重要的地位。數(shù)字濾波是語(yǔ)音信號(hào)處理、圖像處理、模式識(shí)別、頻譜分析等應(yīng)用的基本處理算法。用DSP芯片實(shí)現(xiàn)數(shù)字濾波除了具有穩(wěn)定性好、精確度高、不受環(huán)境影響等優(yōu)點(diǎn)外，還具有靈活性好

9、等特點(diǎn)。</p><p>　　本次課程設(shè)計(jì)的是FIR低通濾波器，首先運(yùn)用MATLAB產(chǎn)生濾波系數(shù)，再用DSP仿真軟件CCS 進(jìn)行仿真，觀察輸入輸出，看是否符合低通濾波器的功能。</p><p>　　第一章 課程設(shè)計(jì)要求及流程</p><p>　　1.1 FIR濾波器設(shè)計(jì)要求</p><p>　　設(shè)計(jì)一個(gè)FIR低通濾波器,通帶邊界頻率為1500

10、Hz，通帶波紋小于1dB；阻帶邊界頻率為2000Hz，阻帶衰減大于40dB；采樣頻率為10000Hz，濾波器階數(shù)為16。FIR濾波器的設(shè)計(jì)也可以用MATLAB窗函數(shù)法進(jìn)行。</p><p><b>　　1.2設(shè)計(jì)流程</b></p><p><b>　　設(shè)計(jì)流程如圖1。</b></p><p><b>　　圖1

11、設(shè)計(jì)流程圖</b></p><p>　　第二章 課程設(shè)計(jì)實(shí)現(xiàn)</p><p><b>　　2.1功能描述</b></p><p>　　低通濾波器是允許低于截至頻率的信號(hào)經(jīng)過(guò)， 但高于截止頻率的信號(hào)不能經(jīng)過(guò)的電子濾波安裝。</p><p><b>　　2.2算法研究</b></p&g

12、t;<p><b>　　FIR的原理及公式</b></p><p><b>　　FIR設(shè)計(jì)原理：</b></p><p>　　如果一個(gè)FIR濾波有一個(gè)沖激響應(yīng)，h(0),h(1),…,h(N-1),和x(n)描繪輸入的時(shí)常濾波n，輸出濾波y(n)的n給出以下方程式：</p><p>　　Y(n)=h(0)x(

13、n)+h(1)x(n-1)+h(2)x(n-2)+…+h(N-1)x[n-(n-1)]</p><p>　　N階有限沖激響應(yīng)濾波器（FIR）公式：</p><p>　　Y(n)= N=0,1,2,3…</p><p><b>　　2.3參數(shù)計(jì)算</b></p><p>　　2.3.1利用Matlab計(jì)算濾波系數(shù)<

14、;/p><p>　　濾波器階數(shù)為16階，截止頻率1500HZ，采樣頻率10000HZ。</p><p>　　2.3.2編寫(xiě)濾波器輸入信號(hào)程序</p><p>　　因?yàn)榻刂诡l率為1500Hz，所以輸入信號(hào)應(yīng)該至少包含兩種頻率成分的正弦信號(hào)，一種信號(hào)頻率小于1500Hz，一種信號(hào)頻率大于1500Hz。利用Matlab生成濾波器輸入信號(hào)程序“.inc”文件和“.dat”文件。

15、 </p><p><b>　　2.4編寫(xiě)源程序</b></p><p>　　2.4.1匯編源程序“.asm”文件的編寫(xiě)</p><p>　　.globalstart,fir</p><p><b>　　.mmregs</b></p><p>　　COFF_FIR_STAR

16、T: .sect "coff_fir"</p><p>　　.include"0127.inc"</p><p>　　K_FIR_BFFR .set20</p><p>　　d_data_buffer .usect"fir_bfr",64</p><p>　　FIR_DP

17、.usect"fir_vars",0</p><p>　　d_filin.usect"fir_vars",1</p><p>　　output.usect"fir_vars",1</p><p>　　input.usect "fir_vars",1</p&g

18、t;<p>　　d_filout.usect"fir_vars",100h</p><p>　　stacksize .set 256</p><p>　　stack .usect "fir_vars",stacksize</p><p>　　.asg AR4,FIR_DATA_P</p>&l

19、t;p>　　.asgAR6,INBUF_P</p><p>　　.asgAR7,OUTBUF_P</p><p>　　.asgAR3,OUTBUF</p><p>　　.asgAR2,INBUF</p><p>　　.sect"fir_prog" </p><p><b&

20、gt;　　nop</b></p><p><b>　　start:</b></p><p>　　stm #stack+stacksize,SP</p><p>　　LD#FIR_DP,DP</p><p>　　STM#d_data_buffer,FIR_DATA_P;</p><

21、;p>　　RPTZ A,#K_FIR_BFFR-1;RC=31</p><p>　　STLA,*FIR_DATA_P+;</p><p>　　STM#d_filin,INBUF_P</p><p>　　STM#d_filout,OUTBUF_P</p><p>　　STM#output,OUTBUF</p&g

22、t;<p>　　STM#input,INBUF</p><p>　　STM#100h,BK</p><p>　　fir_loop: </p><p><b>　　NOP </b></p><p>　　LD*INBUF_P,A</p><p>　　STL A

23、,*INBUF</p><p><b>　　CALLfir</b></p><p>　　STHA,*OUTBUF_P+%</p><p>　　STHA,*OUTBUF</p><p><b>　　main_end:</b></p><p>　　bfir_loop

24、</p><p><b>　　fir:</b></p><p>　　STM#d_data_buffer,FIR_DATA_P</p><p>　　STLA,*FIR_DATA_P</p><p>　　STM#(d_data_buffer+K_FIR_BFFR-1),FIR_DATA_P </p&g

25、t;<p>　　fir_task: </p><p>　　RPTZA,#K_FIR_BFFR-1</p><p>　　MACD*FIR_DATA_P-,COFF_FIR_START,A</p><p><b>　　RET</b></p><p><b>　　.end</b>

26、</p><p>　　2.4.2鏈接命令文件“.cmd”的編寫(xiě)</p><p><b>　　MEMORY</b></p><p><b>　　{</b></p><p><b>　　PAGE 0:</b></p><p>　　PROG:o= 100h,

27、l= 2000h </p><p><b>　　PAGE 1:</b></p><p>　　DATA1:o= 2600h, l= 1000h</p><p>　　DATA2:o= 2100h, l= 100h</p><p>　　DATA3:o= 2200h,l= 100h</p><p&

28、gt;　　DATA4:o= 2300h,l= 100h</p><p>　　DATA5:o= 2400h,l= 100h</p><p>　　DATA6:o= 2500h,l= 100h</p><p><b>　　}</b></p><p><b>　　SECTIONS</b>&l

29、t;/p><p><b>　　{</b></p><p>　　coff_fir : {}> PROGPAGE 0</p><p>　　fir_prog : {}> PROGPAGE 0</p><p>　　fir_vars : {}> DATA1PAGE 1</p><

30、p>　　fir_coff : {}> DATA2PAGE 1</p><p>　　fir_bfr : {}> DATA3PAGE 1</p><p><b>　　}</b></p><p><b>　　2.5 調(diào)試過(guò)程</b></p><p>　　2.5.1調(diào)試前的

31、準(zhǔn)備</p><p>　　1）啟動(dòng)SETUP并選擇芯片C5402 simulator。</p><p>　　2）在CCS的安裝目錄myproject子目錄下創(chuàng)建一個(gè)文件夾0127。</p><p>　　2.5.2 .inc和.dat的使用</p><p>　　將MATLAB生成的0127.inc文件和0127.dat文件復(fù)制到0127文件夾中

32、。</p><p>　　2.5.3編寫(xiě)及編譯程序</p><p>　　1）創(chuàng)建工程文件：在Project菜單中選擇 New項(xiàng)，在Project中輸入0127并保存。</p><p>　　2）編寫(xiě)匯編源程序和鏈接程序：在File菜單中選中New-Source File，編寫(xiě)源程序，然后以.asm格式保存在0127文件夾里；再在File菜單中選中New-Source F

33、ile，編寫(xiě)鏈接程序，然后以.cmd格式保存在0127文件夾里。</p><p>　　3）在工程文件中添加程序文件：在Project菜單中選中Add File to Project…，將以上程序加載到工程目錄下。</p><p>　　4）生成和運(yùn)行程序：在Project菜單中選中Rebuild All對(duì)工程進(jìn)行編譯、匯編和鏈接（如果有錯(cuò)就修改，然后再編譯鏈接直到無(wú)誤為止）；在File菜單中

34、選中Load Program，在當(dāng)前目錄的Debug目錄下選擇0127.out并打開(kāi)，將Build生成的程序加載到DSP中；在Debug菜單中選中Run，運(yùn)行該程序。</p><p>　　2.5.4觀察點(diǎn)設(shè)置</p><p>　　1）設(shè)置斷點(diǎn)：在.asm文件中標(biāo)號(hào)為fir_loop下面的Nop語(yǔ)句設(shè)置軟件斷點(diǎn)（Toggle breakpoint）探針（ToggleProbe point）。

35、</p><p>　　2）在File菜單中選中File I/O，將會(huì)看到File I/O對(duì)話框，按照?qǐng)D2操作然后選中Add Probe Point…再按圖3操作。</p><p><b>　　圖2</b></p><p><b>　　圖3</b></p><p><b>　　3）打開(kāi)觀察窗

36、口 </b></p><p>　　選擇菜單“View”、“Graph”、“Time/Frequency…”按照?qǐng)D4所示依次操作。</p><p><b>　　圖 4</b></p><p>　　2.6運(yùn)行并觀察結(jié)果</p><p>　　輸入時(shí)域和頻域波形如圖5，輸出時(shí)域和頻域波形如圖6。從圖形可以看出，高于1

37、500Hz的頻率成分被濾去，低于1500Hz的頻率成分通過(guò)，符合截止頻率為1500Hz的低通濾波器的特性，達(dá)到了設(shè)計(jì)要求。</p><p><b>　　圖 5</b></p><p><b>　　圖6</b></p><p>　　第三章 課程設(shè)計(jì)總結(jié)</p><p>　　在本次課程設(shè)計(jì)中，我了解到了

38、數(shù)字濾波是信號(hào)處理技術(shù)中的重要部分，重溫了數(shù)字濾波器的基本理論知識(shí)以及實(shí)現(xiàn)它的方法。復(fù)習(xí)了MATLAB的相關(guān)知識(shí)，并且利用MATLAB與CCS設(shè)計(jì)了低通濾波器，更加熟悉CCS環(huán)境下的程序開(kāi)發(fā)方法、調(diào)試工具的使用等。使我將課堂知識(shí)與實(shí)際運(yùn)用結(jié)合起來(lái)，理論聯(lián)系實(shí)際。在這過(guò)程中，我也遇到了很多問(wèn)題，總是無(wú)法正確的顯示輸入輸出波形。在老師和同學(xué)的幫助下，最終找到了問(wèn)題的所在，原來(lái)是濾波器的階數(shù)設(shè)置的不合適。經(jīng)過(guò)多番修改，終于成功完成了本次課程設(shè)

39、計(jì)。</p><p>　　綜合運(yùn)用模擬電子、數(shù)字電子和DSP基本原理等課程中所學(xué)的理論知識(shí)，合理利用MATLAB去獨(dú)立完成一個(gè)項(xiàng)目的設(shè)計(jì)，還需要我不斷地努力、總結(jié)。本次課程設(shè)計(jì)對(duì)我還是有一定的積極意義的，我將以更加飽滿的熱情運(yùn)用所學(xué)知識(shí)投入到今后的工作和生活中。</p><p><b>　　參考文獻(xiàn)</b></p><p>　　《DSP原理及應(yīng)