【文章內(nèi)容簡介】 .word 4*32768/10。 .text x(n4) x(n5) x(n6) x(n7) AR4 x(n) x(n3) x(n2) x(n1) AR3 x_new x_old 47 FIR濾波器的 DSP實現(xiàn) begin: LD x_new, DP SSBX FRCT STM x_new, AR3 。AR3指 向新緩沖區(qū)第一個單元 STM x_old+(size1), AR4 。AR4指 向老緩沖區(qū)最后 1個單元 STM size, BK 。循 環(huán)緩沖區(qū)長度 STM 1, AR0 PORTR PA1,＃ x_new 。輸 入 x[n] FIR: ADD *AR3＋ 0%,*AR4＋ 0%, A 。AH=x[n]+x[n7] (第一次 ) RPTZ B, (size1) 。讓 B=0, 下條 指令執(zhí) 行 size(4)次 FIRS *AR3+0%,*AR4＋ 0%, COEF 。B=B+AH * h0, 。AH=x[n1]+x[n6], (RC)1→RC, (PAR ＋ 1)→PAR STH B, @y 。保存結(jié)果 x(n4) x(n5) x(n6) x(n7) AR4 x(n) x(n3) x(n2) x(n1) AR3 x_new x_old 1 2 3 4 5 AR3 AR4 ；以上程序執(zhí)行完畢后， AR3和 AR4指向位置如圖 48 + PORTW @y， PA0 ； 輸出結(jié)果 MAR * +AR3(2)% ； 修正 AR3， 指向 新緩沖區(qū)最老數(shù)據(jù) MAR * AR4+% ；修正 AR4， 指向 老緩沖區(qū)最老數(shù)據(jù) MVDD *AR3, *AR4+0% 。新 緩沖區(qū)向老緩沖區(qū)傳送一個數(shù) BD FIR PORTR PA1， * AR3 。輸 入新數(shù)據(jù)至新緩沖區(qū) .end x(n4) x(n5) x(n6) x(n7) AR4 x(n) x(n3) x(n2) x(n1) AR3 新 舊 ；前頁程序執(zhí)行完畢后，AR3和 AR4指向位置如圖 AR3 AR4 AR4 49 FIR濾波器的 DSP實現(xiàn) + ( +1 數(shù)字濾波器的 DSP實現(xiàn) FIR濾波器在 C54x上的 simulator實現(xiàn) 舉 例 例 1目標 : 設(shè)計一個 15階 FIR低通濾 波 器 , 截止頻率 (指數(shù)字頻率，對應(yīng)采樣頻率一半的數(shù)字頻率為 1)。 解題思路步驟： ① 借助 MATLAB求出濾波器系數(shù)； ② 用 CCS編程實現(xiàn) FIR濾 波； ③ 用 MATLAB產(chǎn)生驗證數(shù)據(jù)文件 (輸入波 形 )， 用 CCS觀察輸入輸出波形。 51 ① 用 MATLAB編程得出濾波器系數(shù) ； f = [0 1]。 %書中 f = [0 1] 對應(yīng)截止頻率約 m = [1 1 0 0]。 b=fir2(15, f, m) freqz(b,1,1000) 。傳遞函數(shù)分母為 1,ωc= 20log10(Hc)=6db 數(shù)字濾波器的 DSP實現(xiàn) FIR濾波器在 C54x上的實 現(xiàn) 舉 例 例 1目標 : 設(shè)計一個 15階 FIR低通濾 波 器 , 截止頻率 . 或者 用 Matlab 中的 Filter Design and Analysis Tool (FDATool)設(shè)計濾波器 52 ① 用 MATLAB編程得出濾波器系數(shù) ； f = [0 1]。 %書中 f = [0 1] 對應(yīng)截止頻率約 m = [1 1 0 0]。 b=fir2(15,f,m) freqz(b,1,1000) 數(shù)字濾波器的 DSP實現(xiàn) FIR濾波器在 C54x上的實 現(xiàn) 舉 例 目標：設(shè)計一個 15階 FIR低通濾 波 器，截止頻率 或者 用 Matlab 中的 Filter Design and Analysis Tool (FDATool)設(shè)計濾波器 b = 0 0 . 1 0 . 2 0 . 3 0 . 4 0 . 5 0 . 6 0 . 7 0 . 8 0 . 9 1 1 5 0 0 1 0 0 0 5 0 00N o r m a l i z e d Fr 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 00N o r m a l i z e d Fr e q u e n c y ( ? ? r a d / s a m p l e )Magnitude (dB)53 0 0 . 1 0 . 2 0 . 3 0 . 4 0 . 5 0 . 6 1 0 0 0 5 0 00N o r m a l i z e d Fr 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 1 050N o r m a l i z e d Fr e q u e n c y ( ? ? r a d / s a m p l e )Magnitude (dB) FIR濾波器在 C54x上的實現(xiàn)舉例 目標：設(shè)計一個 15階 FIR低通濾 波 器，截止頻率 ② 用 CCS編程實現(xiàn) FIR濾 波； a) 打開 CCS，新建立一工程文件 。 b) 建立 匯 編源文件 ,中 斷向量 表 和 鏈接命令文件 ，并添 加到 。 c) 在 project菜單下選擇 rebuild all選項。 d) 在 file 菜單下，選擇 load program 選項，將生成的 文件裝載到 DSP中 。 54 FIR濾波器在 C54x上的實現(xiàn)舉例 b) 建 立 匯 編源文件 , .mmregs CIR_BFFR_LENGTH .set 16 。濾波器系數(shù)個 數(shù) ,循 環(huán)緩沖區(qū)長度 K_FIR_INDEX .set 1 。循環(huán)緩沖區(qū)步長 K_FRAME_SIZE .set 256 。輸入數(shù)據(jù)緩 沖 區(qū)大小 FIR_COFF_START .sect fir_coff 。 16個 FIR濾 波器系數(shù) .word 9*32768/10000 .word 5*32768/10000 .word 65*32768/10000 .word 224*32768/10000 .word 509*32768/10000 .word 879*32768/10000 .word 1235*32768/10000 .word 1455*32768/10000 .word 1455*32768/10000 .word 1235*32768/10000 .word 879*32768/10000 .word 509*32768/10000 .word 224*32768/10000 .word 65*32768/10000 .word 5*32768/10000 .word 9*32768/10000 FIR_COFF_END 55 15階 FIR低通濾波器 .word 879*32768/10000 .word 509*32768/10000 .word 224*32768/10000 .word 65*32768/10000 .word 5*32768/10000 .word 9*32768/10000 FIR_COFF_END FIR濾波器在 C54x上的實現(xiàn)舉例 用 循環(huán)緩沖區(qū)實現(xiàn) FIR濾波器 FIR濾波 器輸入數(shù)據(jù) 循 環(huán)緩沖區(qū) 存儲器圖 數(shù) 據(jù)存儲器 濾波器系數(shù)緩沖區(qū) input 輸 入數(shù)據(jù)緩沖區(qū)（ 256個） AR6 output 輸 出 數(shù)據(jù)緩沖區(qū)（ 256個） AR7 56 ← AR7 ← AR7 58 FIR濾波器在 C54x上的實現(xiàn)舉例 coff_data .usect “fir_bfr”,20 。 濾 波器系 數(shù)緩沖區(qū) d_data_buffer .usect “fir_bfr” ,40 。 濾波運算輸入循環(huán)緩沖區(qū) input .usect fir_bfr,256 。輸入數(shù)據(jù)緩沖區(qū) output .usect fir_bfr,256 。輸出數(shù)據(jù)緩沖區(qū) .def fir_init 。初始化 FIR濾波 器程序地址標號 .def fir_task 。執(zhí)行 FIR濾波程序地 址 全 局 標號 。 。下 面程 序段的功能是 為 循環(huán)緩 沖 區(qū) 數(shù) 據(jù)和系數(shù)初始 化 。 .sect .fir_prog fir_init: STM input, AR6 。輸入數(shù)據(jù)緩沖區(qū)首址裝入 AR6 STM output, AR7 。輸出數(shù)據(jù)緩沖區(qū)首址裝入 AR7 STM coff_data, AR5 。系數(shù)緩沖 區(qū)首 址裝入 AR5 RPT CIR_BFFR_LENGTH1 。 161 MVPD FIR_COFF_START,*AR5+ 。將 FIR系數(shù)從程序區(qū)移到數(shù)據(jù)區(qū) STM K_FIR_INDEX,AR0 STM d_data_buffer,AR4 。循環(huán)緩沖區(qū)首地址裝入 AR4 RPTZ A,CIR_BFFR_LENGTH STL A,*AR4+ 。復(fù)位緩沖區(qū) STM (d_data_buffer+CIR_BFFR_LENGTH1),AR4 STM coff_data,AR5 。 。下面的子程序的功能為使用 MAC指令執(zhí)行 FIR濾波 。累加器 A(濾波器輸出 )=h(i)*x(ni),i=0,1,2,3.....15。因為 。 fir_task: NOP NOP NOP STM K_FRAME_SIZE1,BRC 。重復(fù) 256次，濾波256個輸入數(shù)據(jù) RPTBD fir_filter_loop1 STM CIR_BFFR_LENGTH,BK 。設(shè)置 FIR循環(huán)緩沖區(qū)大小 LD *