DSP濾波器的設(shè)計(jì)

1、目 錄1引言12 matlab及ccs軟件介紹2 2.1 matlab的介紹2 2.2 ccs的介紹23 fir 數(shù)字低通濾波器的設(shè)計(jì)43.1 fir濾波器參數(shù)設(shè)定4 3.1.1 輸入方波matlab設(shè)計(jì)5 3.1.2 正弦疊加信號matlab設(shè)計(jì)63.2 fir濾波器的設(shè)計(jì)6 3.3 fir濾波器性能驗(yàn)證9 3.3.1 matlab的仿真驗(yàn)證9 3.3.2 ccs的仿真結(jié)果104 iir 數(shù)字低通濾波器的設(shè)計(jì)13 4.1iir濾波器設(shè)計(jì)流程圖13 4.3 iir濾波器性能驗(yàn)證14 4.3.1 matlab的仿真驗(yàn)證14 4.3.2 ccs的仿真結(jié)果155設(shè)計(jì)總結(jié)18參考文獻(xiàn)19附錄201引言

2、當(dāng)今，數(shù)字信號處理(dsp：digtal signal processing)技術(shù)正飛速發(fā)展，它不但自成一門學(xué)科，更是以不同形式影響和滲透到其他學(xué)科：它與國民經(jīng)濟(jì)息息相關(guān)，與國防建設(shè)緊密相連；它影響或改變著我們的生產(chǎn)、生活方式，因此受到人們普遍的關(guān)注。數(shù)字濾波技術(shù)是數(shù)字信號分析、處理技術(shù)的重要分支2-3。無論是信號的獲取、傳輸，還是信號的處理和交換都離不開濾波技術(shù)，它對信號安全可靠和有效靈活地傳輸是至關(guān)重要的。在所有的電子系統(tǒng)中，使用最多技術(shù)最復(fù)雜的要算數(shù)字濾波器了。數(shù)字濾波器的優(yōu)劣直接決定產(chǎn)品的優(yōu)劣。在數(shù)字信號處理中，濾波占有極其重要的地位。數(shù)字濾波是語音信號處理、圖像處理、模式識(shí)別、頻譜分

3、析等應(yīng)用的基本處理算法。用dsp芯片實(shí)現(xiàn)數(shù)字濾波除了具有穩(wěn)定性好、精確度高、不受環(huán)境影響等優(yōu)點(diǎn)外，還具有靈活性好等特點(diǎn)。本設(shè)計(jì)的題目是低通濾波器的設(shè)計(jì)，主要針對dsp芯片進(jìn)行程序設(shè)計(jì)，要求采用匯編語言編寫程序。通過此次課程設(shè)計(jì)，熟悉dsp芯片的硬件結(jié)構(gòu)，加深對dsp芯片的認(rèn)識(shí)；掌握對dsp芯片編程方法，掌握間接尋址方法的使用，重點(diǎn)掌握算法的設(shè)計(jì)。要求采用匯編語言針對dsp芯片設(shè)計(jì)低通濾波器：1.給出低通濾波器技術(shù)指標(biāo)，在matlab中首先設(shè)計(jì)符合要求的濾波器，并得到濾波器參數(shù)，并得到該濾波器的幅頻特性和相頻特性；2.fir低通濾波器設(shè)計(jì)，給出程序流程圖，在ccs軟件環(huán)境下編寫程序；3.分別用方

4、波和正弦波驗(yàn)證fir低通濾波器的設(shè)計(jì)；4.iir低通濾波器設(shè)計(jì)，給出程序流程圖，在ccs軟件環(huán)境下編寫程序；5. 分別用方波和正弦波驗(yàn)證iir低通濾波器設(shè)計(jì)；6.與matlab中設(shè)計(jì)的濾波器，進(jìn)行對比研究。2 matlab及ccs軟件介紹2.1 matlab的介紹matlab是由美國mathworks公司發(fā)布的主要面對科學(xué)計(jì)算、可視化以及交互式程序設(shè)計(jì)的高科技計(jì)算環(huán)境。它將數(shù)值分析、矩陣計(jì)算、科學(xué)數(shù)據(jù)可視化以及非線性動(dòng)態(tài)系統(tǒng)的建模和仿真等諸多強(qiáng)大功能集成在一個(gè)易于使用的視窗環(huán)境中，為科學(xué)研究、工程設(shè)計(jì)以及必須進(jìn)行有效數(shù)值計(jì)算的眾多科學(xué)領(lǐng)域提供了一種全面的解決方案，代表了當(dāng)今國際科學(xué)計(jì)算軟件的先

5、進(jìn)水平。主要應(yīng)用于工程計(jì)算、信號處理與通訊、圖像處理、信號檢測、金融建模設(shè)計(jì)與分析、控制系統(tǒng)設(shè)計(jì)以及計(jì)算生物學(xué)等眾多應(yīng)用領(lǐng)域。matlab是矩陣實(shí)驗(yàn)室（matrixlaboratory）之意。除具備卓越的數(shù)值計(jì)算能力外，它還提供了專業(yè)水平的符號計(jì)算，文字處理，可視化建模仿真和實(shí)時(shí)控制等功能。 matlab的基本數(shù)據(jù)單位是矩陣，它的指令表達(dá)式與數(shù)學(xué),工程中常用的形式十分相似,故用matlab來解算問題要比用c,fortran等語言完相同的事情簡捷得多。 當(dāng)前流行的matlab 5.3/simulink 3.0包括擁有數(shù)百個(gè)內(nèi)部函數(shù)的主包和三十幾種工具包(toolbox)。工具包又可以分為功能性工

6、具包和學(xué)科工具包。功能工具包用來擴(kuò)充matlab的符號計(jì)算,可視化建模仿真，文字處理及實(shí)時(shí)控制等功能。學(xué)科工具包是專業(yè)性比較強(qiáng)的工具包,控制工具包,信號處理工具包,通信工具包等都屬于此類。除內(nèi)部函數(shù)外,所有matlab主包文件和各種工具包都是可讀可修改的文件,用戶通過對源程序的修改或加入自己編寫程序構(gòu)造新的專用工具包。2.2 ccs的介紹ccs軟件簡介ccs（code composer studio）代碼調(diào)試器是一種合成開發(fā)環(huán)境。 它是一種針對標(biāo)準(zhǔn)tms320調(diào)試器接口的交互式方法。ccs目前有ccs1.1， ccs1.2 和ccs2.0等三個(gè)不同時(shí)期的版本，又有cc2000 ( 針對c2xx

7、 )，ccs5000 ( 針對c54xx )和ccs6000 (針對c6x )三個(gè)不同的型號。我們所使用的是ccs50002.2的版本。 ccs5000具有以下特性： ti編譯器的完全集成的環(huán)境： ccs5000目標(biāo)管理系統(tǒng)，內(nèi)建編輯器，所有的調(diào)試和分析能力集成在一個(gè)windows環(huán)境中。 對c和dsp匯編文件的目標(biāo)管理： 目標(biāo)編輯器保持對所有文件及相關(guān)內(nèi)容的跟蹤。它只對最近一次編譯中改變過的文件重新編譯，以節(jié)省編譯時(shí)間。 高集成的編輯器調(diào)整c和dsp匯編代碼： ccs5000的內(nèi)建編輯器支持c和匯編文件的動(dòng)態(tài)語法加亮顯示。使用戶能很容易地閱讀代碼和當(dāng)場發(fā)現(xiàn)語法錯(cuò)誤。 編輯和調(diào)試時(shí)的后臺(tái)編輯：

8、 用戶在使用編譯器和匯編器時(shí)沒有必要退出系統(tǒng)到dos環(huán)境中，因?yàn)閏cs5000會(huì)自動(dòng)將這些工具裝載在它的環(huán)境中。在其窗口中，錯(cuò)誤會(huì)加亮顯示 只要雙擊錯(cuò)誤就可以直接到達(dá)出錯(cuò)處。 在含有浮點(diǎn)并行調(diào)試管理器(pdm)的原有的ms窗口下支持多處理器ccs5000在 windows95和windows-me中支持多處理。pdm允許將命令傳播給所有的或所選擇的處理器。 在任何算法點(diǎn)觀察信號的圖形窗口探針：圖形顯示窗口使用戶能夠觀察時(shí)域或頻域內(nèi)的信號。對于頻域圖，fft在主機(jī)內(nèi)執(zhí)行，這樣就可以觀察所感興趣的部分而無須改變它的dsp代碼。圖顯示也可以同探針連接，當(dāng)前顯示窗口被更新時(shí)，探針被指定，這樣當(dāng)代碼執(zhí)行

9、到達(dá)該點(diǎn)時(shí)，就可以迅速地觀察到信號。 文件探針在算法處通過文件提取或加入信號或數(shù)據(jù)： ccs5000允許用戶從pc機(jī)讀或?qū)懶盘柫?。而不是?shí)時(shí)的讀信號，這就可以用已知的例子來仿真算法。 圖形分析： ccs5000的圖形分析能力在其環(huán)境中是集成的。 在后臺(tái)(系統(tǒng)命令)執(zhí)行用戶的dos程序： 用戶可以執(zhí)行ccs5000中的dos程序，并將其輸出以流水方式送到ccs5000的輸出窗口。且允許用戶將應(yīng)用程序集成到ccs5000。 技術(shù)狀態(tài)觀察窗口： ccs5000的可視窗口允許用戶進(jìn)入c表達(dá)式及相關(guān)變量。結(jié)構(gòu)、數(shù)組、指針都能很簡單地遞歸擴(kuò)展和減少，以便進(jìn)入復(fù)雜結(jié)構(gòu)。 代數(shù)分解窗口： 允許用戶選擇查看寫成

10、代數(shù)表達(dá)式的c格式，提高可讀性 。 目標(biāo)dsp上的幫助： dsp結(jié)構(gòu)和寄存器上的在線幫助可以使用戶不必查看技術(shù)手冊。 3 fir 數(shù)字低通濾波器的設(shè)計(jì)3.1 fir濾波器參數(shù)設(shè)定 fir濾波器的參數(shù)用matlab來設(shè)計(jì)，其程序?yàn)椋簄=40; %濾波器的階數(shù)wn=0.3 %濾波器的截止頻率 b=fir1(n,wn)freqz(b,512,1000) %濾波器的幅頻特性和相頻特性 其運(yùn)行結(jié)果為：wn = 0.3000 b = -0.0000 -0.0012 -0.0017 -0.0008 0.0020 0.0046 0.0036 -0.0025 -0.0101 -0.0110 0.0000 0.0

11、175 0.0259 0.0105 -0.0253 -0.0552 -0.0428 0.0312 0.1483 0.2566 0.3006 0.2566 0.1483 0.0312 -0.0428 -0.0552 -0.0253 0.0105 0.0259 0.0175 0.0000 -0.0110 -0.0101 -0.0025 0.0036 0.0046 0.0020 -0.0008 -0.0017 -0.0012 -0.0000 濾波器的幅頻特性和相頻特性如3-1所示:圖2-1 濾波器的相頻特性圖3-1 fir濾波器的幅相頻特性3.1.1 輸入方波matlab設(shè)計(jì) 方波的matlab程序

12、為：t=0:0.25:50;x=square(pi*t);plot(x);xout=x/max(x); %歸一化 xto_ccs=round(32767*xout) fid=fopen(square.dat,w); %打開文件 fprintf(fid,1651 1 229 1 201n); %輸出文件頭 fprintf(fid,%dn,xto_ccs); %輸出 fclose(fid); %關(guān)閉文件 輸出方波如圖3-2所示：圖3-2 方波信號3.1.2 正弦疊加信號matlab設(shè)計(jì) 正弦疊加信號的matlab程序?yàn)椋篺11=500; f12=9000; fs=20000; %采樣頻率 n=30

13、0 %數(shù)據(jù)個(gè)數(shù) t=1/fs; %采樣周期 n=-n/2:n/2; x11=0.5*sin(2*pi*f11*n*t); x12=0.5*sin(2*pi*f12*n*t); x_base=(x11+x12); plot(x_base); %正弦疊加信號的波形yff=abs(fft(x_base) df=n*(fs/n) xout=x_base/max(x_base); %歸一化 xto_ccs=round(32767*xout) fid=fopen(input.dat,w); %打開文件 fprintf(fid,1651 1 229 1 301n); %輸出文件頭 fprintf(fid,%

14、dn,xto_ccs); %輸出 fclose(fid); 輸出正弦疊加信號的波形如圖3-3所示：圖:3-3 正弦疊加信號3.2 fir濾波器的設(shè)計(jì) fir結(jié)構(gòu)中的延遲z-1 在程序中實(shí)現(xiàn)是數(shù)據(jù)在存儲(chǔ)器中的滑動(dòng)，這對一個(gè)l階的fir數(shù)字濾波器來說，就是要求在存儲(chǔ)器中開辟一個(gè)l個(gè)樣本并進(jìn)行乘累加運(yùn)算。這在dsp中可以用兩種方法來實(shí)現(xiàn)該緩沖區(qū)。對于線性緩沖區(qū)來說，每當(dāng)dsp讀入一個(gè)樣本后都將此樣本向后移動(dòng)，讀完最后一個(gè)樣本后，最前面的樣本（最老的）被推出緩沖區(qū)，并輸入最新的樣本至緩沖區(qū)頂部。線性緩沖區(qū)的存儲(chǔ)器延時(shí)可以由具有數(shù)據(jù)移動(dòng)的dmov指令來實(shí)現(xiàn)，它需要將緩沖區(qū)置于內(nèi)部daram塊中；對于循環(huán)

15、緩沖區(qū)，每當(dāng)輸入新樣本時(shí)，以新樣本改寫區(qū)中最老的樣本，而緩沖區(qū)中的其他樣本不需要其他樣本移動(dòng)。循環(huán)緩沖區(qū)可采用循環(huán)尋址技術(shù)來實(shí)現(xiàn)，它不需要必須定位在daram塊中，但需要適當(dāng)?shù)木彌_區(qū)對準(zhǔn)以滿足循環(huán)尋址的要求。本次課設(shè)用的是循環(huán)尋址，下面主要介紹一下循環(huán)尋址的方法。循環(huán)緩沖實(shí)現(xiàn)延時(shí)的特點(diǎn)是：1 實(shí)現(xiàn)n階fir數(shù)字濾波器時(shí)，需要在數(shù)據(jù)存儲(chǔ)器中開辟一個(gè)稱為滑窗的n個(gè)單元的緩沖區(qū)，滑窗中存放最新的n個(gè)樣本。2 當(dāng)每次輸入新樣本時(shí)以新的樣本改寫滑窗中最老的樣本，而滑窗其他的數(shù)據(jù)不需要移動(dòng)。3 利用片內(nèi)bk（循環(huán)緩沖區(qū)長度）寄存器對滑窗進(jìn)行間接尋址，循環(huán)緩沖區(qū)地址首位相鄰。 其程序流程圖3-4所示：定義4

16、1位空間li，shuai；200位空間in將li值給in空間 將li空間的首41個(gè)數(shù)據(jù)給in；并分別 將ar3，ar4指向buf和li的空間尾地址設(shè)置循環(huán)緩沖區(qū)長度bk為41位，ar0賦值為-1，并將ar5指向xn空間首地址將li和shuai中的數(shù)據(jù)進(jìn)行乘累加運(yùn)算，結(jié)果給out空間，ar5指向地址自加1將ar2指向的in空間的下一個(gè)數(shù)據(jù)給ar3指向的內(nèi)容結(jié)束開始圖3-4 程序流程圖3.3 fir濾波器性能驗(yàn)證3.3.1 matlab的仿真驗(yàn)證 matlab中輸入為方波信號信號和疊加正弦信號時(shí)的輸入輸出信號波形如圖3-5所示：圖3-5 方波及正弦波濾波后波形 由 matlabd所得仿真圖可以看出

17、，方波信號經(jīng)過此濾波器后輸出信號為單一頻率的正弦波；而疊加的正弦信號通過此低通濾波器后，將高頻的正弦波信號濾掉，只留下了頻率最低的單一正弦波信號。因此可以看出，所設(shè)計(jì)的低通濾波器在 matlab 中能夠很好地滿足設(shè)計(jì)要求，實(shí)現(xiàn)低頻通過，高頻過濾的作用。3.3.2 ccs的仿真結(jié)果ccs中方波信號和濾波后波形如圖3-6所示：圖3-6 方波波形及濾波后波形圖8 方波波形及濾波后波形ccs中方波信號的幅頻特性如圖3-7所示：圖3-7 方波信號的幅頻特性ccs中方波信號濾波后的幅頻特性如圖3-8所示圖3-8 方波信號濾波的幅頻特性ccs中正弦疊加信號與濾波后波形如圖3-9所示：圖3-9 正弦疊加信號波

18、形及濾波后波形ccs中正弦疊加信號的幅頻特性如圖3-10所示：圖3-10正弦疊加信號的幅頻、相頻特性ccs中正弦疊加信號濾波后的幅頻特性如圖3-11所示：圖3-11 正弦疊加信號濾波后的幅頻特性 通過上圖可看出當(dāng)輸入信號為方波時(shí)，輸出信號基本為正弦波；當(dāng)輸入信號為疊加正弦波時(shí)，輸出為低頻正弦波。這說明濾波器起到了較好的濾波效果，符合設(shè)計(jì)要求。通過matlab仿真和ccs仿真濾波效果比較，可以看出用兩種軟件所得的仿真圖基本一致，均證明了此濾波器達(dá)到了設(shè)計(jì)要求。同時(shí)也可以看出：dsp程序的濾波效果要好于matlab的濾波效果。 4 fir 數(shù)字低通濾波器的設(shè)計(jì)4.1iir濾波器設(shè)計(jì)流程圖由前面ii

19、r濾波器的設(shè)計(jì)原理確定運(yùn)用雙線性變換法設(shè)計(jì)iir濾波器的流程圖如圖4-1所示：數(shù)字濾波器指標(biāo)模擬濾波器指標(biāo)設(shè)計(jì)模擬濾波器對應(yīng)的數(shù)字濾波器器雙線性變換開始結(jié)束圖4-1 iir設(shè)計(jì)流程圖 iir濾波器設(shè)計(jì)方法有兩類，經(jīng)常用到的一類設(shè)計(jì)方法是借助于模擬濾波器的設(shè)計(jì)方法進(jìn)行的。其設(shè)計(jì)思路是:先設(shè)計(jì)模擬濾波器得到傳輸函數(shù)ha(s)，然后將ha(s)按某種方法轉(zhuǎn)換成數(shù)字濾波器的系統(tǒng)函數(shù)h (z)。這一類方法是基于模擬濾波器的設(shè)計(jì)方法相對比較成熟，它不僅有完整的設(shè)計(jì)公式，也有完整的圖表供查閱。更可以直接調(diào)用matlab中的對應(yīng)的函數(shù)進(jìn)行設(shè)計(jì)。另一類是直接在頻域或者時(shí)域中進(jìn)行設(shè)計(jì)的，設(shè)計(jì)時(shí)必須用計(jì)算機(jī)作輔助設(shè)

20、計(jì)，直接調(diào)用matlab中的一些程序或者函數(shù)可以很方便地設(shè)計(jì)出所需要的濾波器。我們采取直接調(diào)用的方法。4.2 iir濾波器參數(shù)設(shè)計(jì) fir濾波器的參數(shù)用matlab來設(shè)計(jì)，其程序?yàn)椋篺s=1000; fp=100; fs=300; ap=3; as=20;wp=fp/(fs/2);%計(jì)算歸一化角頻率ws=fs/(fs/2);n,wc=buttord(wp,ws,ap,as);%計(jì)算階數(shù)和截止頻率b,a=butter(n,wc);%計(jì)算h(z)分子、分母多項(xiàng)式系數(shù)h,f=freqz(b,a,500,fs);%計(jì)算h(z)的幅頻響應(yīng),freqz(b,a,計(jì)算點(diǎn)數(shù),采樣速率) 濾波器的幅頻特性和相頻

21、特性如圖4-2所示：圖4-2 iir濾波器的幅相頻特性4.3 fir濾波器性能驗(yàn)證4.3.1 matlab的仿真驗(yàn)證matlab中輸入為正弦疊加信號信號時(shí)的輸入輸出信號波形如圖4-3所示：圖4-3正弦疊加信號濾波后波形matlab中輸入信號為方波信號時(shí)的輸入輸出信號波形如圖4-4所示：4-4方波信號信號時(shí)的輸入輸出信號圖4.3.2 ccs的仿真結(jié)果ccs中方波信號和濾波后波形如圖4-5所示：圖4-5 方波波形及濾波后波形ccs中方波信號的幅頻特性如圖4-6所示:圖4-6 輸入方波信號的幅頻特性ccs中方波信號濾波后的幅頻特性如圖4-7所示：圖4-7 輸入方波信號的幅頻特性ccs中正弦疊加信號與

22、濾波后波形如圖4-8所示：圖4-8 正弦疊加信號的輸出波形ccs中正弦疊加信號的幅頻特性如圖4-9所示：圖4-9 正弦疊加信號的副頻特性ccs中正弦疊加信號濾波后的幅頻特性如圖4-10所示：圖4-10 正弦疊加信號濾波后的副頻特性 通過以上各圖可看出在matlab中輸入信號為方波時(shí)，然后經(jīng)過所設(shè)計(jì)的低通濾波器輸出信號理論上應(yīng)為正弦波，但是由于采用二階低通濾波器及其他原因，濾波效果并不是特別好，當(dāng)在matlab中輸入信號為疊加正弦波時(shí)，輸出為低頻正弦波，濾波效果還好，而在ccs中方波和正弦波都得到了較好的濾波效果，通過比較看出ccs的仿真效果要比matlab的仿真效果要好，總之所設(shè)計(jì)的iir低通

23、濾波器起到了應(yīng)有的濾除高頻信號的濾波效果，符合設(shè)計(jì)要求。5設(shè)計(jì)總結(jié)本次課程設(shè)計(jì)是數(shù)字濾波器的設(shè)計(jì)，是在上學(xué)期dsp原理即數(shù)字信號處理理論知識(shí)的基礎(chǔ)上，對dsp芯片的首次完整應(yīng)用。通過完整的實(shí)踐過程，對dsp芯片組成的基本系統(tǒng)的相關(guān)電路有了比較深的認(rèn)識(shí)，熟悉了dsp芯片的系統(tǒng)設(shè)計(jì)和應(yīng)用開發(fā)流程，并利用圖書館、網(wǎng)絡(luò)、詢問同學(xué)等方式查找資料和解決相關(guān)的難題。數(shù)字濾波器是dsp的典型應(yīng)用，學(xué)會(huì)了有助于觸類旁通，利于進(jìn)一步的學(xué)習(xí)研究，能做到理解其他基于dsp的系統(tǒng)的功能和工作原理。掌握了基于dsp的應(yīng)用開發(fā)，開闊了視野，增長了知識(shí)，是進(jìn)入現(xiàn)代數(shù)字信號處理領(lǐng)域重要技能，乃至大規(guī)模集成電路的開發(fā)也是會(huì)用到的

24、基礎(chǔ)，今后要予以重視并積極努力去學(xué)習(xí)。回想這幾天實(shí)踐過程,真的很有感觸。的確是這樣，從初學(xué)dsp匯編語言到今天可以簡單但是獨(dú)立的運(yùn)用；從不知道怎么解決一個(gè)簡單的語法錯(cuò)誤，到今天可以清楚地了解復(fù)雜語句之間的工作過程，進(jìn)步源于一點(diǎn)一滴的積累，這或許正是我們這次獨(dú)立地完成課程設(shè)計(jì)最大的收獲?；仡櫿麄€(gè)數(shù)字濾波器的設(shè)計(jì)過程，感慨很幼稚但也有屬于自己的深刻。由于dsp初次應(yīng)用的不甚了解，真的感覺自己走了不少的彎路。最主要的就是對實(shí)驗(yàn)思路的不了解。本次設(shè)計(jì)結(jié)合了ccs和matlab軟件實(shí)驗(yàn)平臺(tái)，真正了解了設(shè)計(jì)的過程，思路也變得相對簡易和明確。當(dāng)然，真正的動(dòng)手去做才發(fā)現(xiàn)要明白的東西絕不僅僅是書上的那點(diǎn)東西，很

25、多實(shí)踐上東西，如果真正的明白了，比理論上弄懂要深刻的多。比如在編寫程序的過程中，對卷積的編程必須基于對理論知識(shí)的深刻了解，同時(shí)實(shí)踐對理論也具有能動(dòng)的反作用，加深對卷積問題的理解。只有注意積累每次的體會(huì)和錯(cuò)誤，下次碰到該注意的地方才懂得小心，碰到出過錯(cuò)誤的地方才會(huì)格外注意不要再犯同樣的錯(cuò)誤，這種方法應(yīng)該適用做任何一件事。在設(shè)計(jì)中遇到一些困難和問題，在向老師請教和與同學(xué)的討論中，解決了問題，覺得很有收獲。當(dāng)你用心的去做一件事的時(shí)候，這件事就不會(huì)輕易的用不知對錯(cuò)的結(jié)果去應(yīng)付你，就像你不會(huì)去拿這樣的態(tài)度去對待他，做一件事和想一件事也不會(huì)是一樣，如果當(dāng)時(shí)我只是停留在第一天的空想，我就不會(huì)在編程和仿真的過

26、程中發(fā)現(xiàn)和改正那么多的錯(cuò)誤。實(shí)踐之所以高于理論，或許正是因?yàn)樗鼮槲覀兲峁┝烁喾稿e(cuò)誤和改正錯(cuò)誤的機(jī)會(huì)吧。最后衷心感謝老師的悉心指導(dǎo)，董老師的教導(dǎo)是我們前進(jìn)的方向。參考文獻(xiàn)1 范壽康 m北京：電子工業(yè)出版社，20062 萬永格 m北京：科學(xué)出版社，20073 程佩青 m北京：清華大學(xué)出版社，20014 高西全 m北京：電子工業(yè)出版社，20095 范壽康 m北京：電子工業(yè)出版社，20096 胡廣書 m北京：清華大學(xué)出版社，2005附錄fir數(shù)字濾波器程序：.title fir.asm .mmregs .global _c_int00li: .usect li , 41shuai: .usect s

27、huai, 41in: .usect in, 200 out: .usect out, 240 .datatable: .word -0*32768/10000, -12*32768/10000, -17*32768/10000, 18*32768/10000,20*32768/10000 .word 46*32768/10000, 36*32768/10000, -25*32768/10000, -101*32768/10000,-110*32768/10000 .word 0*32768/10000, 175*32768/10000, 259*32768/10000, 105*32768/

28、10000,-253*32768/10000 .word -552*32768/10000, -428*32768/10000, 312*32768/10000, 1483*32768/10000,2566*32768/10000 .word 3006*32768/10000, 2566*32768/10000, 1483*32768/10000, 312*32768/10000,-428*32768/10000 .word -522*32768/10000, -253*32768/10000, 105*32768/10000, 259*32768/10000,175*32768/10000

29、.word 0*32768/10000, -110*32768/10000, -101*32768/10000, -25*32768/10000,36*32768/10000 .word 46*32768/10000, 20*32768/10000, -8*32768/10000, -17*32768/10000,-12*32768/10000,-0*32768/10000 .text_c_int00: ssbx frct stm #li,ar5 ;ar5指向沖激響應(yīng) rpt #40 mvpd table,*ar5+ stm #shuai,ar4 ;ar4指向緩沖區(qū)間 rptz a,#40 s

30、tl a,*ar4+ stm #in,ar3 ;ar3指向輸入 stm #out,ar2 ;ar2指向輸出 stm #(li+40),ar5 ; stm #(shuai+40),ar4 stm -1,ar0 ld #shuai,dp mvdd *ar3+,*ar4 ;向緩沖去送數(shù) stm #240,brc rptb #table ;卷積算法 stm #41,bk rptz a,#40 mac *ar4+0%,*ar5+0%,a sth a,*ar2+ table: mvdd *ar3+,*ar4+0% .end-o fir.out-m fir.mapmemory page 0: eprom:

31、org=0e000h len=1000h vecs : org=0ff80h len=0080h page 1: saram: org=0060h len=0040h daram: org=0101h len=1380h sections .text : eprom page 0 .vectors : vecs page 0 .data : eprom page 0 .bss : saram page 1 li : align(128) daram page 1 shuai : align(128) daram page 1 in : daram page 1 out : daram page

32、 1 .title vectors.asm .ref _c_int00 .sect .vectors b _c_int00 .endfir數(shù)字濾波器程序： .titleiir.asm .mmregs .def _c_int_00x .usect x,3y .usect y,3b0 .usect b0,3a0 .usect a0,3 output .usect output,256input .usect intput,256 .datatable .word 0 .word 0 .word 0 .word 0 .word 10000*16384/10000 .word 19996*16384/

33、10000 .word 10000*16384/10000 .word 3753*16384/10000 .word -10541*16384/10000 .text_c_int_00: ssbx frct stm #x,ar1 rpt #1 mvpd #table,*ar1+ stm #y,ar1 rpt #1 mvpd #table+2,*ar1+ stm #b0,ar1 rpt #2 mvpd table+4,*ar1+ stm #a0,ar1 rpt #1 mvpd #table+7,*ar1+ stm #x+2,ar2 stm #a0+1,ar3 stm #y+1,ar4 stm #

34、b0+2,ar5 stm #3,bk stm #-1,ar0 stm #input,ar6 stm #output,ar7 stm 255,brciir rptb next-1 ld *ar6+,a stl a,*ar2 mpy *ar2+0%,*ar5+0%,a mac *ar2+0%,*ar5+0%,a mac *ar2,*ar5+0%,a mac *ar4+0%,*ar3+0%,a mac *ar4+0%,*ar3+0%,a mar *ar3+0% sth a,*ar4 sth a,*ar7+next: b next .endiir.obj-o iir.out-m iir.map-e _

35、c_int_00memory page 0:rom1 :org=0100h,len=1000h rom2 :org=2000h,len=0500h page 1: daram0 :org=1000h,len=0500h daram1 :org=1500h,len=0010h daram2 :org=1600h,len=0010h daram3 :org=1700h,len=0010h daram4 :org=1800h,len=0010h sections .text:rom1 page0 .data : daram0 page 1 x : daram1 page 1 y : daram2 page 1 b0 : daram3 page 1 x0 : daram4 page 1 .title vectors.asm .ref _c_int_00 .sect .vectors b _c_int_00 .endmatlab程序：%窗函數(shù)法設(shè)計(jì)fir低通濾波器 clear all ws=0.5*pi wp=0.3*pi wc=(ws+wp)/2%由過度帶決定n ow=ws-wpn=8*pi/owalfa=(n-1)/2; nm=0:n-1%理想低通濾波器 hd=sin(wc*(nm-alfa)/(pi*