dsp課程設(shè)計--有限沖擊響應(yīng)fir濾波器的設(shè)計

1、<p>　　課程名稱： DSP原理及應(yīng)用課程設(shè)計 </p><p>　　實驗項目： 有限沖擊響應(yīng)FIR濾波器的設(shè)計 </p><p>　　實驗地點： </p><p>　　2012年 12月17日</p><p><b>

2、　　一、設(shè)計目的</b></p><p>　　1使學(xué)生加深對DSP芯片TMS320C55x的結(jié)構(gòu)、工作原理的理解；2了解FIR 的原理和特性；3熟悉設(shè)計FIR 數(shù)字濾波器的原理和方法；4 學(xué)習(xí)FIR 濾波器的DSP 實現(xiàn)原理；5 學(xué)習(xí)使用ccs 的波形觀察窗口觀察輸入、輸出信號波形和頻譜變化情況。</p><p><b>　　二、設(shè)計原理</b>&

3、lt;/p><p>　　濾波器就是在時間域或頻域內(nèi)，對已知激勵產(chǎn)生規(guī)定響應(yīng)的網(wǎng)絡(luò)，使其能夠從信號中提取有用的信號，抑制并衰減不需要的信號，濾波器的設(shè)計實質(zhì)上就是對提出的要求給出相應(yīng)的性能指標(biāo)，再通過計算，使物理可實現(xiàn)的實際濾波器響應(yīng)特性逼近給出的頻率響應(yīng)特性。</p><p>　　FIR數(shù)字濾波器是一種非遞歸系統(tǒng)，其傳遞函數(shù)為：</p><p><b>　　=

4、</b></p><p>　　由此可得到系統(tǒng)的差分方程為：</p><p>　　則，F(xiàn)IR濾波器的結(jié)構(gòu)如下圖：</p><p>　　其沖激響應(yīng)是有限長序列，它 濾波器系數(shù)向量,N 為FIR濾波器的階數(shù)。</p><p>　　在數(shù)字信號處理應(yīng)用中往往需要設(shè)計線性相位的濾波器，F(xiàn)IR濾波器在保證幅度特性滿足技術(shù)要求的同時，很容易做到嚴(yán)

5、格的線性相位特性為了使濾波器滿足線性相位條件，要求其單位脈沖響應(yīng)為實序列，且滿足偶對稱或奇對稱條件，即或.這樣，當(dāng)N為偶數(shù)時，偶對稱線性相位FIR濾波器的差分方程表達(dá)為</p><p>　　由上可見，F(xiàn)IR濾波器不斷地對輸入樣本延時后，再做乘法累加運算，將濾波器結(jié)果輸出。因此，F(xiàn)IR實際上是一種乘法累加運算。而對于線性相位FIR濾波器系數(shù)的對稱特性，可以采用結(jié)構(gòu)精簡的FIR結(jié)構(gòu)將乘法器數(shù)目減少一半。</p&

6、gt;<p>　　應(yīng)用MATLAB設(shè)計FIR濾波器的主要任務(wù)就是根據(jù)給定的性能指標(biāo)設(shè)計一個，使其逼近這一指標(biāo)，進(jìn)而計算并確定濾波器的系數(shù)，再將所設(shè)計濾波器的幅頻響應(yīng)、相頻響應(yīng)曲線作為輸出，與設(shè)計要求進(jìn)行比較，對設(shè)計的濾波器進(jìn)行優(yōu)化，設(shè)計完成之后將得到的FIR濾波器的單位沖激響應(yīng)序列的各個參數(shù)值。</p><p>　　三、總體方案設(shè)計1 啟動CCS ，在CCS 中建立一個匯編源文件、建立一個C 源文

7、件和一個命令文件，并將這三個文件添加到工程，再編譯并裝載程序2 設(shè)置波形時域觀察窗口，得到其濾波前后波形變化圖；3 設(shè)置頻域觀察窗口，得到其濾波前后頻譜變化圖。</p><p><b>　　四、源程序</b></p><p><b>　　1、新建工程文件</b></p><p>　　打開Setup CCS選擇c55芯片

8、，打開CCS軟件，選擇project/new，命名為myproject。</p><p><b>　　2、添加文件到工程</b></p><p>　?。?）建立文件 選擇file/new，匯編源程序，如下：</p><p>　　.global _fir,_init,_B,_outdata</p><p><b&g

9、t;　　_fir</b></p><p><b>　　bset frct</b></p><p>　　amov #_B,xdp</p><p>　　mov #_B,cdp</p><p>　　mov t0,ac0</p><p>　　sub #1,ac0</p><

10、p>　　mov ac0,mmap(csr)</p><p>　　add ac0,ar0</p><p>　　mov #0,ac0</p><p><b>　　rpt csr</b></p><p>　　macmz *ar0-,*cdp+,ac0</p><p>　　mov hi(ac0),t

11、0</p><p><b>　　ret</b></p><p><b>　　_init</b></p><p>　　mov mmap(t0),ac0</p><p>　　sub #1,ac0</p><p>　　mov ac0,ar7</p><p>

12、　　rptz ac0,ar7</p><p>　　mov ac0,*ar0+</p><p><b>　　ret</b></p><p><b>　　_outdata</b></p><p>　　mov t1,ac0</p><p>　　sub #2,ac0</p>

13、;<p>　　mov ac0,mmap(csr)</p><p>　　add ac0,ar0</p><p><b>　　rpt csr</b></p><p>　　delay *ar0-</p><p><b>　　mar *ar0+</b></p><p>

14、　　mov t0,*ar0</p><p><b>　　ret</b></p><p><b>　　c文件</b></p><p>　　#include "s.h"</p><p>　　#include "math.h"</p><p>

15、;　　#define signal_1_f 200</p><p>　　#define signal_2_f 620</p><p>　　#define signal_sample_f 2000</p><p>　　#define pi 3.1415926</p><p>　　#define coff_L 23</p><

16、p>　　#define bufer_L 256</p><p>　　int data_in[bufer_L];</p><p>　　int out[bufer_L] ;</p><p>　　int firout;</p><p>　　int x[coff_L+1];</p><p><b>　　int

17、k=0;</b></p><p>　　int bufer=bufer_L;</p><p>　　extern int fir(int *,int);</p><p>　　extern int init(int *,int);</p><p>　　extern int outdata(int *,int,int);</p>

18、;<p>　　void inputwave();</p><p>　　void main()</p><p><b>　　{</b></p><p>　　inputwave();</p><p>　　init(x,BL);</p><p><b>　　while(1)<

19、;/b></p><p><b>　　{</b></p><p>　　x[0]=data_in[k];</p><p>　　firout=fir(x,BL);</p><p>　　outdata(out,firout,bufer);</p><p><b>　　k++;</b

20、></p><p>　　if(k>=bufer_L)</p><p><b>　　{</b></p><p><b>　　k=0;</b></p><p><b>　　}</b></p><p><b>　　}</b>&

21、lt;/p><p><b>　　}</b></p><p>　　void inputwave()</p><p><b>　　{</b></p><p>　　float wt1;</p><p>　　float wt2;</p><p><b>

22、　　int i;</b></p><p>　　for(i=0;i<=bufer_L;i++)</p><p><b>　　{</b></p><p>　　wt1=2*pi*i*signal_1_f;</p><p>　　wt1=wt1/signal_sample_f;</p><p&

23、gt;　　wt2=2*pi*i*signal_2_f;</p><p>　　wt2=wt2/signal_sample_f;</p><p>　　data_in[i]=(cos(wt1)+cos(wt2))/2*32768;</p><p><b>　　}</b></p><p><b>　　}</b&g

24、t;</p><p><b>　　S.h文件</b></p><p>　　int BL=23;</p><p>　　int B[23]={-417,342,1961,1692,-762,-1003,1676,1115,-3200,-1159,10335,17562,10335,-1159,-3200,1115,1676,-1003,-762,1

25、692,1961,342,-417};</p><p><b>　?。?）保存文件</b></p><p>　　（3）添加文件 右鍵點擊建立的工程，選擇Add file to project，在選中建立的文件。</p><p>　　3、匯編源程序的鏈接命令文件：</p><p>　　-stack 0x500</p

26、><p>　　-sysstack 0x500</p><p>　　-heap 0x1000</p><p><b>　　-c</b></p><p><b>　　-u _Reset</b></p><p>　　-l rts55.lib</p><p>&l

27、t;b>　　MEMORY</b></p><p><b>　　{</b></p><p><b>　　PAGE 0:</b></p><p>　　RAM(RWIX): origin=0x000100, length=0x01ff00</p><p>　　ROM(RIX): or

28、igin=0x020100, length=0x01ff00</p><p>　　VECS(RIX): origin=0xffff00, length=0x000200</p><p><b>　　PAGE 2:</b></p><p>　　IOPORT(RWI):origin=0x000000, length=0x020000<

29、/p><p><b>　　}</b></p><p><b>　　SECTIONS</b></p><p><b>　　{</b></p><p>　　.text >ROM PAGE 0</p><p>　　.data >ROM PAGE 0&l

30、t;/p><p>　　.bss >RAM PAGE 0</p><p>　　.const >RAM PAGE 0</p><p>　　.sysmem >RAM PAGE 0</p><p>　　.stack >RAM PAGE 0</p><p>　　.sysstack >RAM PAGE

31、0</p><p>　　.switch >RAM PAGE 0</p><p>　　.cinit >RAM PAGE 0</p><p>　　.pinit >RAM PAGE 0</p><p>　　.vectors >VECS PAGE 0</p><p>　　.ioport >IOP

32、ORT PAGE 2</p><p><b>　　}</b></p><p><b>　　4、加載可執(zhí)行文件</b></p><p>　　選擇工程，點擊file/load program，選擇.out文件</p><p>　　5、運行 執(zhí)行.out文件</p><p>　　6

33、、設(shè)置波形時域/頻域觀察窗口，得到其濾波前后波形變化圖；</p><p><b>　　五、實驗結(jié)果及分析</b></p><p>　　Acquisition buffer size：256</p><p>　　Display data size：250</p><p>　　DSP data type：16-bits si

34、gned integer</p><p><b>　　時域輸入輸出波形：</b></p><p>　　由時域圖可得：濾波前的時域圖是帶很多毛刺不規(guī)則的正弦波，而濾波后的波形是平滑的等幅正弦波且波形振幅相對濾波前的最大振幅有所減小。</p><p><b>　　頻域輸入輸出波形：</b></p><p&g

35、t;　　由頻域圖可得：濾波前的頻域圖會出現(xiàn)兩個帶尖峰的圖像,而濾波后的波形只會出現(xiàn)一個帶尖峰圖像，而且是出現(xiàn)在低頻部分，說明信號經(jīng)過濾波器后，將信號中的高頻成分濾除掉了。</p><p><b>　　六、設(shè)計總結(jié) </b></p><p>　　通過這次基于DSP的數(shù)字濾波器設(shè)計，我初步學(xué)會了應(yīng)用CCS進(jìn)行DSP設(shè)計的實踐操作。熟悉了DSP集成開發(fā)環(huán)境（CCS），掌握了

36、FIR數(shù)字濾波器原理和方法。通過對系統(tǒng)的調(diào)試，對線性相位FIR數(shù)字濾波器特性及不同階數(shù)對濾波器特性的影響進(jìn)行了研究。通過觀察時域和頻域圖，對FIR數(shù)字濾波器特性進(jìn)行了驗證,在應(yīng)用DSP對數(shù)字濾波器實現(xiàn)的過程中，體會到了DSP在數(shù)字濾波實現(xiàn)中的優(yōu)缺點。用DSP 實現(xiàn)FIR數(shù)字濾波,具有穩(wěn)定性好、準(zhǔn)確度高、靈活性好、不受環(huán)境影響等優(yōu)點。在這次的課程設(shè)計中學(xué)會了如何和老師、同學(xué)一起交流、探討問題，在學(xué)習(xí)的過程中同時培養(yǎng)了和同學(xué)之間的合作能力，