【正文】 ]+Cx[k1] 這是個(gè)二階差分方程，其單位沖擊響應(yīng)即為 sinkwT。 在本實(shí)驗(yàn)中將介紹另一種獲得正弦 /余弦信號(hào)的方法，即利用數(shù)字振蕩器用疊代方法產(chǎn)生正弦信號(hào)。相對(duì)于 IIR 濾波器， FIR 濾波器有著易于實(shí)現(xiàn)和系統(tǒng)絕對(duì)穩(wěn)定的優(yōu)勢(shì)，因此得到廣泛的應(yīng)用；對(duì)于時(shí)變系統(tǒng)濾波器 的研究則導(dǎo)致了以卡爾曼濾波為代表的自適應(yīng)濾波理論的產(chǎn)生。 數(shù)字信號(hào)處理（ DSP）是利用計(jì)算機(jī)或?qū)Ｓ锰幚碓O(shè)備，以數(shù)值計(jì)算的方法、對(duì)信號(hào)進(jìn)行采集、濾波、增強(qiáng)、壓縮、估值和識(shí)別等加工處理，借以達(dá)到提取信息和便于應(yīng) 用的目的，其應(yīng)用范圍涉及幾乎所有的工程技術(shù)領(lǐng)域。 if(index==1) index=L1。 n++。 out[i]=0。 int index=0。 /* Output buffer */ main() { unsigned int L=101。 恢復(fù) ST1, ST2 和 ST3 popm ST2_55 popm ST1_55 ret .end /* include define Fs 64000 /* 采樣頻率 */ define T 1/Fs /* 采樣時(shí)間 */ define f1 500 /* 正弦信號(hào) 1 頻率 */ define f2 2500 /* 正弦信號(hào) 2 頻率 */ define f3 3300 /* 正弦信號(hào) 3 頻率 */ define PI define w1 (2*PI*f1*T) /* 正弦信號(hào) 1 數(shù)字頻率 =2*pi*f1/Fs */ define w2 (2*PI*f2*T) /* 正弦信號(hào) 2 數(shù)字頻率 =2*pi*f2/Fs */ define w3 (2*PI*f3*T) /* 正弦信號(hào) 3 數(shù)字頻 率 =2*pi*f3/Fs */ define a1 /* 正弦信號(hào) 1 幅度 */ define a2 /* 正弦信號(hào) 2 幅度 */ define a3 /* 正弦信號(hào) 3 幅度 */ extern int fir(int *,int *,unsigned int,int )。設(shè)置外部循環(huán)次數(shù)為 L1 mpym *AR0+,*AR2+,AC0 。AR0 和 AR2 為循環(huán)緩沖區(qū)指針 mov T1,AR0 。置位 SMUL mov mmap(AR0),BSA01 。 參數(shù) 3: T1 輸入信號(hào)當(dāng)前值在循環(huán)緩沖區(qū)的序數(shù) 。 int fir(int *,int *, unsigned int， int) 。 2) 打開“ Debug” 選單，選擇 “ Halt” 選項(xiàng)或按 Shift+F5終止實(shí)驗(yàn)程序。 Display Type：選擇顯示形式 Graph Title：圖表名稱 Start Address：起始地址 Page：顯示數(shù)據(jù)模式 Acquisition Buffer Size：獲得緩沖器的大小 Index Increment：索引增量 Display Data Size：顯示數(shù)據(jù)的大小 DSP Data Type： DSP 數(shù)據(jù)類型 Sampling Rate：采樣頻率 按如下步驟打開圖形窗口顯示輸入信號(hào) 波形 a) Display Type: Singel Time， 如要看頻譜則選為 FFT Magnitude， Graph Title: d_filin Start Address:d_data_buffer Acquistition Buffer Size: 32 Display Data Size:32 DSP Data Type:16bit signed integer Autoscale:on 點(diǎn)擊下面的 OK 完成設(shè)置。將出現(xiàn)Break/Probe/Profile Points 的對(duì)話框，單擊 Probe Point 框中的選項(xiàng)之后，在 Connect 中選中 “FILE IN ： c:\ti\myprojects\fir\” 。以上步驟是為了完成在程序中加入斷點(diǎn)。 3） 將文件夾 lowpass拷貝到新建的工程文件夾內(nèi)。觀察濾波器輸入信號(hào)、輸出信號(hào)波形。 3. 向工程 high 中加載文件 、 、 。 （ 2） 為采用匯編語(yǔ)言編寫的 fir 濾波模塊。 Wn = 對(duì)應(yīng)的非歸一化頻率為 Fs/2，其中 Fs 為采樣頻率。 本設(shè)計(jì)是基于 FIR算法的 dsp 高通濾波器的設(shè)計(jì)，采用 DSP芯 片。在設(shè)計(jì) FIR 濾波器時(shí)，必須根據(jù)這些參數(shù)來確定 FIR 濾波器的系數(shù)及階數(shù)。 FIR 指的是單一采樣響應(yīng)的期間是有限長(zhǎng)的， IIR 指的是單一采樣響應(yīng)的期間是無限長(zhǎng)的。 二、實(shí)驗(yàn)原理 數(shù)字濾波器一直以來就是數(shù)字信號(hào)處理器最廣為人知的應(yīng)用。主要應(yīng)用是實(shí)時(shí)快速的實(shí)現(xiàn)各種數(shù)字信號(hào)處理算法，如卷積及各種變換等。目前的 DSP 芯片以其強(qiáng)大的數(shù)據(jù)處理功能在通信和其他信號(hào)處理領(lǐng)域得到廣泛注意并已成為開發(fā)應(yīng)用的熱點(diǎn)技術(shù)。 2． 掌握 FIR 數(shù)字濾波器的 DSP 實(shí)現(xiàn)方法。濾波器分成兩種形態(tài)：有限長(zhǎng)度脈沖響應(yīng)（ FIR）濾波器。濾波器特性可以由很多的參數(shù)來決定，這些參數(shù)包括帶通漣波、帶阻漣波、帶通截止頻率及帶阻截止頻率。這一過程本實(shí)驗(yàn)不詳述，本實(shí)驗(yàn)假定該過程已經(jīng)完成，并得到了濾波器的階數(shù)及系數(shù)。 Wn 為歸一化邊界頻率， 0 Wn ， Wn 處增益為 6dB。執(zhí)行該指令得到的系數(shù)為： b=[6,3,0,4,8,13,19,25,31,38,44,49,52,53,51,45,34,20,0,24,53,85,120,155,189,221,247,266,274,270,251,216,163,91,0,110,238,381,538,705,879,1055,1229,1397,1554,1695,1817,1915,1988,2033,30713,2033,1988,1915,1817,1695,1554,1397,1229,1055,879,705,538,381,238,110,0,91,163,216,251,270,274,266, 247,221,189,155, 120,85,53,24, 0,20,34,45,51,53,52, 49, 44, 38, 3, 25, 19, 13, 8, 4, 0, 3, 6] 采用函數(shù) plot(b)可以繪制濾波器的沖擊響應(yīng)圖： 0 20 40 60 80 100 120 0 . 200 . 20 . 40 . 60 . 811 . 2 采用函數(shù) freqz(b)可以繪制濾波器的頻率響應(yīng)圖： 0 0 . 1 0 . 2 0 . 3 0 . 4 0 . 5 0 . 6 0 . 7 0 . 8 0 . 9 1 1 0 0 0 0 5 0 0 005000N o r m a l i z e d F r e q u e n c y ( ? ? r a d / s a m p l e )Phase (degrees)0 0 . 1 0 . 2 0 . 3 0 . 4 0 . 5 0 . 6 0 . 7 0 . 8 0 . 9 1 1 0 0 5 0050N o r m a l i z e d F r e q u e n c y ( ? ? r a d / s a m p l e )Magnitude (dB) 2．實(shí)驗(yàn)源程序（清單） （ 1） 為主程序，用于設(shè)置濾波器、測(cè)試信號(hào)參數(shù)，產(chǎn)生測(cè)試信號(hào)，調(diào)用 fir 濾波程序。 high。 6．全速運(yùn)行程序，數(shù)秒鐘后暫停。 2） 將文件添加到工程中。 6） 打開 文件在 “fir_loop:” 程序段的 “NOP ” 上點(diǎn)擊右鍵選中 Toggle Breakpoint，重復(fù)單擊右鍵選中 Toggle Probe Point, 在 “main_end: ” 程序段內(nèi)的 “b fir_loop” 行 上單擊右鍵選中 Toggle Breakpoint 。就會(huì)在屏幕上再?gòu)棾?“File Input” 的對(duì)話框，雙擊C:\ti\myprojects\fir\lowpass\，單擊 Add Probe Point。 9） 選擇 VIEW 選單中 Graph 副選框內(nèi)的 Time/Frequency 打開對(duì)話框，修改以下內(nèi)容。 1) 打開“ Debug” 選單，選擇 “Animate” 選項(xiàng)或按 F12 運(yùn)行程序，觀察運(yùn)行結(jié)果。 該程序用于實(shí)現(xiàn) FIR 濾波器 ,可被 C 語(yǔ)言程序調(diào)用 。 參數(shù) 2: T0 FIR 濾波器的階數(shù) L 。設(shè)置 FRCT,SXMD,SATD bset SMUL 。設(shè)置循環(huán)緩沖區(qū)大小 or 0x5,mmap(ST2_55) 。T0=L2 mov T0,CSR 。 用 Q15 格式存放結(jié)果 popm ST3_55 。 /* input buffer */ int out[101]。 unsigned int n=0。i++) { in[i]=0。 signal += a3*cos((float)w3*n)。 index。數(shù)字信號(hào)處理（ DSP） 包括兩重含義：數(shù)字信號(hào)處理技術(shù)（ Digital Signal Processing）和數(shù)字信號(hào)處理器（ Digital Signal Processor）。其中，線性時(shí)不變的數(shù)字濾波器是最基本的類