【正文】 ***************** 。 circular address, so when save out_wave 。 The following codes are used to initalize McBSP0 ! 。 choose XCR10 stm 40h,xcr10 stm 5,spsa0 。 The following codes are initalized AC01, ALL AC01 setup same ! 。 A=05,B=40(0x28) idle 1 。 DA 0 dB) ld 3,a idle 1 。 send 8th regs 0 ld 0h,a stm 0,38h stm 2021h,39h 。 TIMER start, each about 100HZ .... orm 8,imr 。 each AC01 has 25 dots to D/A ! mvdm wave_ptr,ar7 。******************************************************************** 。 The following codes are serived for vc540239。************************************************************************ 。 0 a, then repeat 34 times mac *ar2+0%,*ar3+0%,a 。 save FIR result to wave buffer ! stl a,temp mvmd ar7,out_wave_ptr 。 The following data is used by make wave ! when using, must copy to 。 Filter Length 38 。C034 .word 20 。C026 .word 157 。C018 .word 9162 。C010 .word 585 。C002 .word 7 。C004 .word 7 。C012 .word 701 。C020 .word 10230 。C028 .word 773 。C036 .word 7 。************************************************************************ 。 output buffer 0x1000 (show from 0x1020) con_fir: 29 popm ar3 popm ar2 popm ar0 popm al popm ah popm ar7 popm bk popm st0 popm st1 ret 。 save FIR result in temp ld temp,2,a 。 new AD data a, remove two LSBs stm 1,ar0 stm N,bk 。 set xf=0 show_con: ld 5,a 。 has new AD data popm st0 popm st1 rete 。 not equal ,jump ! stm 0,mcbsp_init 。 transmit ! 。 stop TIMER ! stm 0ffffh,prd 。 send 7th regs 0 ld 3,a idle 1 。 send 3 ld 00405h,a idle 1 。 flp=10M/2/40/A= 25k (all pass) 。 enable all int ! 。 choose RCR10 stm 40h,39h stm 3,spsa0 。 if is_new_data ==1 , has new AD data by slave ac01 ! 24 stm 1000h,ar0 stm 1800h,ar5 stm 400h,bk ld *ar0,t ld 0,b 。 save out wave buffer39。 stm 97ffh,58h 。/* 。 faction multiply stm 10h,26h 。 interrupte vector table ! 。 路漫漫其修遠(yuǎn)兮，吾將上下而求索。濾波后經(jīng) D/A 輸出模擬量，結(jié)果證實(shí)該程序可實(shí)現(xiàn)對(duì)采樣率為 25kHz 的模擬信號(hào)進(jìn)行實(shí)時(shí)濾波處理。由于串口發(fā)送中斷將每 （ 25KHz ）產(chǎn)生一次，所以我們將一個(gè)周期的方波信號(hào)分 25 次送出，這樣經(jīng) D/A 變化后便可得到 1KHz 的方波。注意 FIR 濾波系數(shù)存 ；放在程序存貯區(qū)， filter_coff 為系數(shù)起始地址 MAR *ar2% ； 修改 Buffer1 指針 15 MAR *+ar3(2)% ； 修改 Buffer2 指針 圖 4 輸入序列的存儲(chǔ)形式 DSP 實(shí)現(xiàn) 利用 DES5402PP 實(shí)驗(yàn)箱的模擬信號(hào)輸出通道產(chǎn)生一個(gè) 1KHz 的方波 ,然后利用信號(hào)輸入通道對(duì)產(chǎn)生的方波進(jìn)行低通濾波 ,得到一個(gè) 1KHz 的正弦信號(hào)，并使用 CCS 的圖形顯示工具顯示輸入和輸出的波形。長(zhǎng)度為 N 的線性相位 FIR 濾波器的輸出表達(dá)式為： ?)(ny ? ?? ?????12/0 )]}1([)(){(NK kNnxknxkh 32 利用 FIRS 指令 ,需要將輸入數(shù)據(jù)緩沖分成兩個(gè) ,循環(huán)緩沖區(qū)大小寄存器的值設(shè)為N/2。 AR0=1 STM N,BK 。 13 使用帶 MAC 指令的循環(huán)尋址實(shí)現(xiàn) FIR 濾波器 為了能正確使用循環(huán)尋址，必須先初始化 BK，塊長(zhǎng)為 N。在一些特殊的場(chǎng)合，要求的信號(hào)處理速度極高，用通用 DSP 芯片很難實(shí)現(xiàn)。該設(shè)備采用廣泛使用的 TMS320VC5402 為核心，配備包括 AD/DA、 鍵盤(pán)、液晶等 I|O 在內(nèi)的完整的DSP 系統(tǒng) ]4[ 。 FIR 濾波器的 MATLAB 設(shè)計(jì) MATLAB 是 英文 Matrix Laboratory(矩陣實(shí)驗(yàn)室 )的縮寫(xiě)，它是由美國(guó)MathWorks 公司推出的用于數(shù)值計(jì)算和圖形處理的數(shù)學(xué)計(jì)算環(huán)境， 一套用于科學(xué)計(jì)算的可視化高性能語(yǔ)言與軟件環(huán)境。目前常用的 FIR 濾波器的 設(shè)計(jì)方法主要有三種，分別是窗函數(shù)法、頻率采樣法和切比雪夫等波紋逼近優(yōu)化設(shè)計(jì)法。運(yùn)用線性相位響應(yīng)指標(biāo)進(jìn)行濾波器設(shè)計(jì)具有如下優(yōu)點(diǎn)： ① 只包含實(shí)數(shù)算法，不涉及復(fù)數(shù)運(yùn)算； ② 不存在延遲失真，只有固定數(shù)量的延遲； ③ 長(zhǎng)度為 N 的濾波器（階數(shù)為 N1），計(jì)算量為 N/2 數(shù)量級(jí)。因此，指標(biāo)的形6 式一般在頻域中給出幅度和相位響應(yīng)。低通有限沖激響應(yīng)濾波器（低通 FIR 濾波器）有其獨(dú)特的優(yōu)點(diǎn)，由于 FIR 系統(tǒng)只有零點(diǎn)，因此，系統(tǒng)總是穩(wěn)定的，而且容易實(shí)現(xiàn)線性相位和允許實(shí)現(xiàn)多通道濾波器。本文首先介紹了 FIR 數(shù)字濾波器的原理，接下給出了在 MATLAB 環(huán)境下設(shè)計(jì)濾波器的方法，同時(shí)對(duì)本論文中的 FIR 數(shù)字濾波器給以設(shè)計(jì)，最后用 DSP 來(lái)實(shí)現(xiàn)。 本文的內(nèi)容提要 本文第一部分是引言，第二部分介紹濾波器的基本原理和設(shè)計(jì)方法 ；第三部分是 FIR 濾波器的 DSP 實(shí)現(xiàn)；最后，第四部分是結(jié)論，對(duì)本文作了總結(jié)。它提供對(duì)幅度響應(yīng)函數(shù)的要求，一般應(yīng)用于 FIR 濾波器的設(shè)計(jì)。之后，利用數(shù)字濾波器的設(shè)計(jì)方法，設(shè)計(jì)出一個(gè)實(shí)際濾波器模型來(lái)逼近給定的目標(biāo)。 2 頻率抽樣法 設(shè)計(jì)的基本思想是：從頻域出發(fā)，把給定的理想的頻率響應(yīng) )( ?jd eH 以等間隔抽樣，所得到的 h(k)作逆離散傅氏變換，從而求得 h(n)，并用與之相應(yīng)的頻率響應(yīng) )( ?jeH 去逼近理想的頻率響應(yīng) )( ?jd eH 。在設(shè)計(jì)數(shù)字濾波器時(shí)，通常采用 MATLAB 來(lái)進(jìn)行輔助設(shè)計(jì)和仿真。 CCS 專門(mén)為 TI 的DSP 設(shè)計(jì)，擴(kuò)充了基本的開(kāi)發(fā)工具，增加了調(diào)試和實(shí)時(shí)分析功能，主要有以下特點(diǎn)： 1 項(xiàng)目管理和編輯工具，用來(lái)對(duì)項(xiàng)目進(jìn)行創(chuàng)建、配置、生成和運(yùn)行，支持通過(guò) makefile 和命令行來(lái)生成可執(zhí)行代碼； 2 代碼生成工具包括 C 編譯器、匯編器和鏈接器等工具，另外 CCS 還提供了可視化的鏈接工具 Visual Linker，幫助用戶來(lái)調(diào)整代碼和數(shù)據(jù)內(nèi)存，實(shí)現(xiàn)程序內(nèi)存的最優(yōu)化； 3 CCS 提供了豐富的調(diào)試工具，支持多種模擬器和仿真版，支持在程序運(yùn)行中探測(cè)存儲(chǔ)器特定地址的值，用戶可以方便地在 PC 和目標(biāo)版上交換數(shù)據(jù)； 4 提供了豐富的函數(shù)庫(kù)支持，如 CSL、 DSPLIB 以及 IMGLIB 等； 5 支持實(shí)時(shí)內(nèi)核 DSP/BIOS，極大的便利了程序的開(kāi)發(fā)，使開(kāi)發(fā)者能夠把精12 力集中在對(duì)算法的研究上 ]5[ 。 在上述幾種方法中，第 ② 種方法的缺點(diǎn)是速度較慢，一般可用于 DSP 算法的模擬。 圖 3 FIR 濾波器存儲(chǔ)器里的數(shù)據(jù)存儲(chǔ)方式 FIR 濾波器存儲(chǔ)器里的數(shù)據(jù)存儲(chǔ)方式中濾波系數(shù)指針初始化時(shí)指向 h(N1) ，經(jīng)過(guò)一次 FIR 濾波計(jì)算后，在循環(huán)尋址的作用下，仍然指向 h(N1) 。重復(fù) MAC 指令 N 次 ,先將 A 清零 14 MAC *FIR_DATA_P+0%,*FIR_COFF_P+0%,A 。每次進(jìn)行濾波之前 ,應(yīng)先將緩沖區(qū) 1 頂部的數(shù)據(jù)移到緩沖區(qū) 2 的底部 ,新來(lái)的一個(gè)樣本存儲(chǔ)到緩沖區(qū) 1 中時(shí) ,并對(duì)緩沖區(qū) 1 指針 AR2 加 1(使用循環(huán)尋址 )。 AC01 提供一個(gè) 14bit 的 D/A 和一個(gè) 14bit 的 A/D 通道。完成串口設(shè)置后，還需要修改中斷向量表以便正確響應(yīng)串口 0 的接收和發(fā)送中斷請(qǐng)求。在實(shí)際應(yīng)用中，只需按要求修改濾波器參數(shù)，并對(duì)程序作較少的改動(dòng)，即可實(shí)現(xiàn)不同截止頻率的低通濾波器，實(shí)用性較強(qiáng)。bootloader... .ref _c_int00 .mmregs gpiocr .set 3ch gpiosr .set 3dh dxr10 .set 23h drr10 .set 21h spsa0 .set 38h spcr10 .set 39h spcr20 .set 39h xcr10 .set 39h pcr0 .set 39h dxr11 .set 43h spsa1 .set 48h spcr11 .set 49h spcr21 .set 49h xcr11 .set