DSP高通濾波器課程設(shè)計報告

1、而船石油火毒DSP課程設(shè)計報告題目：FIR高通濾波器設(shè)計姓名學(xué)號教學(xué)院系專業(yè)年級指導(dǎo)教師目錄一、設(shè)計題目1二、設(shè)計目標(biāo)1三、算法研究與參數(shù)計算11、FIR的原理和參數(shù)生成公式12、利用MATLA時算濾波系數(shù)13、輸入信號參數(shù)計算2四、編寫源程序3五、調(diào)試過程41、調(diào)試前準(zhǔn)備42、MATLAB勺使用53、編寫及編譯程序54、設(shè)置斷點和探針55、打開觀察窗口6六、實驗結(jié)果及分析61、輸入信號的時域波形和頻域波形62、輸出信號的時域波形和頻域波形6七、設(shè)計心得7一、設(shè)計題目FIR高通濾波器設(shè)計二、設(shè)計目標(biāo)設(shè)計一個FIR高通濾波器，通帶邊界頻率為6000Hz,采樣頻率為20000HzFIR濾波器的設(shè)計

2、用MATLA函函數(shù)法進行。三、算法研究與參數(shù)計算1、FIR的原理和參數(shù)生成公式ojn階有限沖激啊而通波器cm)公式；/(，)=(無*)+-(H-1+燈)公式】TNR12.FIR儀計原理；根據(jù)索嬖h是偶對稱為T褐化運算產(chǎn)生如F計算方法如果一個網(wǎng)代濾波有一個沖激唁皮.h,h(1)和苒描說地人的時需濾波m輜出漁波式n)的n給出以卜.方程式；y(n)-h(O)x(n)+h(l)x(n-lhh(2)x(a-2>.+h(N1聞Xm切公式17圖3-12、利用MATLA升算濾波系數(shù)在MATLABF面輸入圖3-2所示程序，可得到濾波系數(shù)并生成INC文件。y=roond(firl(64j6000/20000

3、*2,!high1八3#焚)；,Eid=fopen(T0126.inc1jTwr，);fprintf(fid,',-wordayJ;fclose(fid)圖3-2輸入freqz(y,1,512),MATLA中顯示局通濾波器的濾波特性曲線。如圖3-3所示圖3-33、輸入信號參數(shù)計算MATLA計輸入圖3-4中所示程序，包含兩種頻率成分的正弦信號，一種信號頻率1000Hz,一種信號6000H4CoiriVTi&rLdlWizidow»i=D：l:255;y=romd(sin(2pi*i*1000/20000+sin(Z*pi*i*6000/20000)*32768/2);f

4、idfopenC1207030126.daf/v)；fprintf(fid/1651工DI0n);Eptirrtf(£l4jy);fclose(fid)圖3-4其頻譜特性曲線如圖3-5。五三名mu>3一落*ufc附而丫H2印系> >lr«杳：口I圖3-5四、編寫源程序參考資料，編寫匯編語言源程序：HIGHPASS.set1;ifyouwanttouse,pleasesetthevalueto1.globalstart,fir.mmregsCOFF_FIR_START:.sect"coff_fir".copy"0126.inc&

5、quot;K_FIR_BFFR.set64d_data_buffer.usect"fir_bfr",64FIR_DP.usect"fir_vars",0d_filin.usect"fir_vars",1output.usect"fir_vars",1input.usect"fir_vars",1d_filout.usect"fir_vars",100hstacksize.set256stack.usect"fir_vars",stacksize.asgAR

6、4,FIR_DATA_P.asgAR6,INBUF_P.asgAR7,OUTBUF_P.asgAR3,OUTBUF.asgAR2,INBUF.sect"fir_prog"nopstart:stm#stack+stacksize,SPLD#FIR_DP,DPSTM#d_data_buffer,FIR_DATA_PRPTZA,#K_FIR_BFFR-1STLA,*FIR_DATA_P+STM#d_filin,INBUF_PSTM#d_filout,OUTBUF_PSTM#output,OUTBUFSTM#input,INBUFSTM#100h,BKfir_loop:NOP;Ad

7、dBreakpoint&porbepointLD*INBUF_P,ASTLA,*INBUFCALLfirSTHA,*OUTBUF_P+%STHA,*OUTBUFmain_end:bfir_loopfir:;SSBXSXM;SSBXFRCTSTM#d_data_buffer,FIR_DATA_PSTLA,*FIR_DATA_PSTM#(d_data_buffer+K_FIR_BFFR-1),FIR_DATA_Pfir_task:RPTZA,#K_FIR_BFFR-1MACD*FIR_DATA_P-,COFF_FIR_START,ARET.end五、調(diào)試過程1、調(diào)試前準(zhǔn)備1）啟動SETU所

8、選擇芯片。2）在CCS勺安裝目錄下的myproject的子目錄下創(chuàng)建新文件夾，命名120703012&2、MATLAB勺使用使用MATLABk成0126.Inc文件和1207030126.dat文件，保存在1207030126文件夾中。如圖5-1:irectorvLV-Pron-amFilesKMATLrBTlVworl:囪于地I國0,AlLFil電工L|Fi"Typ.SimL”ModiMalSl207030126.datDATFile2KB201S7-91107Ol26jncjINCFile1KB20157-911:072J2J_/_圖5-13、編寫及編譯程序1）新建pro

9、ject并命名120703012&堂 Pj1«CEL £門”-Projects日圖 1207Q012C,vjt國 oia »d：I nsMioe j皿 Lj Gf ma叵 |Ifl 0L26, inrI Libruri«sE L_ Source2）在project中添加0126.asm文件和0126.cmd文件。如圖5-2。3）編譯并鏈接程序。4）設(shè)置工程文件，選擇菜單“project”的“buildoptions”項。圖5-25）加載程序。File-loadprogram4、設(shè)置斷點和探針在標(biāo)號"fir_loop”下的“NOP語句處設(shè)

10、置軟件斷點和探針選擇“filefile I/O;按照實驗程序操作5、打開觀察窗口選擇菜單“view”、“Graph"、“time/Freauency”。按要求設(shè)置。六、實驗結(jié)果及分析1、輸入信號的時域波形和頻域波形輸入信號的時域波形如圖6-1所示，頻域波形如圖6-2所示圖6-1圖6-2這個信號由1000Hz和6000Hz的信號疊加而得2、輸出信號的時域波形和頻域波形輸出信號的時域波形如圖6-3所示，頻域波形如圖6-4所示圖6-3圖6-4從輸入信號波形和輸出信號波形對比看，低頻信號（1000Hz）不能順利通過濾波器，被削減到幾乎為0,而高頻信號（6000HZ）可以順利通過。因此這是一個性能良好的高通濾波器。七、設(shè)計心得本次我的課程設(shè)計是基于窗函數(shù)法的FIR數(shù)字高通濾波器，由于數(shù)字信號處理課程是上學(xué)期學(xué)的，再加上自己對MATLA歆件操作不熟練。導(dǎo)致前期用了很長一段時間查閱相關(guān)書籍，重新學(xué)習(xí)了數(shù)字濾波器的相關(guān)知識和MATLA歆件