【正文】 AR6+,A 。加載輸入數(shù)據(jù) fir_filter: STL A,*AR40% 。*AR40% 是否更正確？ RPTZ A,(CIR_BFFR_LENGTH1) MAC *AR4+0%, *AR5+0%, A 。 STH A,*AR7+ 。 fir_filter_loop NOP NOP NOP .end b) 建 立 匯 編源文件 （續(xù)） b) 建立匯編源文件 （續(xù)） 。RPT CIR_BFFR_LENGTH1 。將 FIR系數(shù)從程序 區(qū) 。MVPD FIR_COFF_START, *AR5+ 。將系 數(shù)從程序區(qū)移到數(shù)據(jù)區(qū) STM K_FIR_INDEX, AR0 。循環(huán)緩沖區(qū)步 長(zhǎng)是 1→AR0 STM d_data_buffer, AR4 。循環(huán)緩沖區(qū)首地址裝入 AR4 RPTZ A, CIR_BFFR_LENGTH 。重復(fù) 16次 ,清零循環(huán)緩沖區(qū) STL A, *AR4+ 。清零循環(huán)緩沖區(qū) STM (d_data_buffer+CIR_BFFR_LENGTH1), AR4 STM coff_data, AR5 。指 向?yàn)V波器系數(shù)緩 沖 區(qū) 頂 部 。 。下面的子程序的功 能：使 用 MAC指令執(zhí)行 FIR濾波 。累加器 A(濾波器輸出 )=A+h(i)*x(ni), i=0,1,2,3.....15。 。 fir_task: 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 *AR6+,A 。加載輸入數(shù)據(jù) fir_filter: STL A,*AR40% 。*AR40% 是否更正確？ RPTZ A,(CIR_BFFR_LENGTH1) MAC *AR4+0%, *AR5+0%, A 。 STH A,*AR7+ 。 fir_filter_loop .end 59 指 向 循 環(huán)緩沖 區(qū) 底 部加載最新值 10()????Nkkb x n kb) 建 立匯編源文件 （續(xù)） 。 。下面的子程序的功能： 使用 MAC指令執(zhí)行 FIR濾波 。累加器 A(濾波器輸出 )=A+h(i)*x(ni), i=0,1,2,3.....15。 。 fir_task: STM K_FRAME_SIZE1, BRC 。塊 重復(fù) 256次 ,濾波 256輸入數(shù)據(jù) RPTBD fir_filter_loop1 STM CIR_BFFR_LENGTH, BK 。設(shè)置 FIR循環(huán)緩沖區(qū)大 小 16 LD *AR6+,A 。 從輸 入數(shù)據(jù)緩沖 區(qū)準(zhǔn)備加載數(shù)據(jù)到循環(huán)緩沖區(qū) fir_filter: STL A,*AR4+% 。*AR40%使 x循環(huán)緩沖順序 與 FIR圖 同 , 這里相反 RPTZ A, (CIR_BFFR_LENGTH1) MAC *AR4+0%, *AR5+0%, A 。濾波器輸 出 = A + h(i)*x(ni), STH A, *AR7+ 。 保存結(jié)果到輸出數(shù)據(jù)緩沖區(qū) fir_filter_loop .end 60 10()????Nkkb x n k加 載數(shù) 據(jù)到循環(huán)緩沖 區(qū)替代最老數(shù)據(jù) FIR濾波器在 C54x上的實(shí)現(xiàn)舉例 b) 建立匯編源程序及 中 斷向量 表文件 （續(xù)） .title .ref fir_init .sect vectors RESET: BD fir_init 。 復(fù) 位時(shí) 跳轉(zhuǎn)到 fir濾波程序 fir_init開始 NOP NOP .space 124*16 。沒有使用其它中斷 , 保留空間 .end 61 b) 建立匯編源程序及鏈接命令文件 （續(xù)） MEMORY { PAGE 0: EPROG: =0x1400, len=0x7C00 VECT: =0xFF80, len=0x80 PAGE 1: USERREGS: =0x0060, len=0x1C /*沒有用 */ IDATA: =0x0080, len=0x50 IDATA1: =0x0180, len=0x1280 EDATA: =0x1400, len=0x8000 } SECTIONS { .text : EPROG PAGE 0 .fir_prog : EPROG PAGE 0 vectors : VECT PAGE 0 fir_coff : EPROG PAGE 0 fir_bfr : IDATA1 PAGE 1 fir_vars : IDATA1 PAGE 1 STACK : EDATA PAGE 1 } 62 FIR濾波器在 C54x上的實(shí)現(xiàn)舉例 目標(biāo)：設(shè)計(jì)一個(gè) 15階 FIR低通濾 波 器，截止頻率 ② 用 CCS編程實(shí)現(xiàn) FIR濾 波； a) 打開 CCS，新建立一工程文件 。 b) 建立 匯 編源文件 ,中 斷向量 表 和 鏈接命令文件 ，并添 加到 。 c) 在 project菜單下選擇 rebuild all選項(xiàng)。 d) 在 file 菜單下，選擇 load program 選項(xiàng)，將生成的 文件裝載到 DSP中 。 63 1） 產(chǎn)生驗(yàn)證數(shù)據(jù)文 件 (輸入數(shù)據(jù) ) clear all x(1:20)=0 %行向量 20元素皆為 0 x(2:2:1000,:)=400 %1000 20矩陣 ,奇數(shù)行為 0,偶數(shù)行皆為 400 x=x? % 20 1000矩 陣 ,奇 數(shù)列為 0,偶 數(shù)列皆 為 400 x=x(:) %列向 量 20200元素 ,間 隔 20元素 0與 400切換 fid1=fopen(??,?wt?) %?t?以文本方式打開可寫入文 ?w?) fprintf(fid1,?1651 1 0 0 0\n?) %向文件 ??寫文 件頭 a1=fprintf(fid1,39。%d\n39。,x) fclose(fid1) FIR濾波器在 C54x上的實(shí)現(xiàn)舉例 ③ 用 MATLAB產(chǎn)生驗(yàn)證數(shù)據(jù)文件， 用 CCS觀察輸入輸出波形。 64 %向文件 ??寫數(shù)據(jù) 1） 產(chǎn)生驗(yàn)證數(shù)據(jù)文 件 (輸入數(shù)據(jù) ) clear all x(1:20)=0 %行向量 20元素皆為 0 x(2:2:1000,:)=400 %1000 20矩陣 ,奇數(shù)行為 0,偶數(shù)行皆為 400 x=x? % 20 1000矩 陣 ,奇 數(shù)列為 0,偶 數(shù)列皆 為 400 x=x(:) %列向 量 20200元素 ,間 隔 20元素 0與 400切換 fid1=fopen(??,?wt?) %以文本方式打開可寫入文件 fprintf(fid1,39。1651 1 0 0 0\n39。) %輸出文件頭 a1=fprintf(fid1,39。%d\n39。,x) fclose(fid1) FIR濾波器在 C54x上的實(shí)現(xiàn)舉例 ③ 用 MATLAB產(chǎn)生驗(yàn)證數(shù)據(jù)文件， 用 CCS觀察輸入輸出波形。 65 CCS支持的 .dat文件的格 式 文件頭 為 ： 定數(shù) 數(shù)據(jù)格式 起 始地址 頁(yè) 類型 數(shù) 據(jù)塊大小 1651 1 3000 0 C8 固 定標(biāo)識(shí) ：又稱 MagicNumber，是固定數(shù) 1651 數(shù)據(jù)格式 ： 1十六進(jìn)制 2十進(jìn)制 3十進(jìn)制長(zhǎng)整型 4 十進(jìn) 制浮點(diǎn) 型 起始地 址 ： 為 裝載數(shù)據(jù)在 存儲(chǔ) 器 的 基 地址 0x3000 頁(yè)類型 ： 0程 序 存儲(chǔ) 器 1 數(shù) 據(jù) 存儲(chǔ) 器 長(zhǎng)度 ： 裝入數(shù)據(jù)的長(zhǎng) 度 0xC8(十進(jìn)制數(shù) 200) FIR濾波器在 C54x上的實(shí)現(xiàn)舉例 ③ 用 MATLAB產(chǎn)生驗(yàn)證數(shù)據(jù)文件， 用 CCS觀察輸入輸出波形。 2)選 擇 FileDataLoad… 打開之前 Matlab 生成 的文件 3)將 Address設(shè)為 input，Length為 C8（ 200d），Page為 data 4)運(yùn)行程序 , 點(diǎn)擊 Run按鈕。 66 FIR濾波器在 C54x上的實(shí)現(xiàn)舉例 5） 可以在 view菜單下選擇 graph/time frequence，設(shè) 置好相應(yīng)的參數(shù)，來觀測(cè)輸入和輸出數(shù)據(jù)的波形。 67 FIR濾波器在 C54x上的實(shí)現(xiàn)舉例 5)可 以在 view菜單下選擇 graph/time frequence，設(shè)置好相應(yīng)的參數(shù)，來觀測(cè)輸入和輸出數(shù)據(jù)的波形。 68 FIR濾波器在 C54x上的實(shí)現(xiàn)舉例 5)可 以在 view菜單下選擇 graph/time frequence，設(shè)置好相應(yīng)的參數(shù)，來觀測(cè)輸入和輸出數(shù)據(jù)的波形。 69 第 9章 C54x綜合應(yīng)用系統(tǒng)設(shè)計(jì) 5)可 以在 view菜單下選擇 graph/time frequence，設(shè)置好相應(yīng)的參數(shù)，來觀測(cè)輸入和輸出數(shù)據(jù)的波形。 70 第 9章 C54x綜合應(yīng)用系統(tǒng)設(shè)計(jì) IIR濾波器的 DSP實(shí)現(xiàn) IIR濾波器的結(jié)構(gòu)簡(jiǎn)單，運(yùn)算量小，具有經(jīng)濟(jì)、靈活、高效等特點(diǎn) 。 IIR濾波器的基本原理和設(shè)計(jì)方法 1. IIR濾波器的特點(diǎn) IIR濾波器差分方程的一般形式為： 01( ) ( ) ( )MNkkkky n b x n k a y n k??? ? ? ???71 IIR濾波器的 DSP實(shí)現(xiàn) IIR濾波器具有無限長(zhǎng)的單位脈沖響應(yīng)，在結(jié)構(gòu)上存在反饋回路，即是遞歸型的。 當(dāng) M=N時(shí)，在零初始條件下，對(duì)上式進(jìn)行 z變換，得到傳遞函數(shù)： 101111. . . (). . . NN NNkNNkNkb z b z b z zH z Cz a z a z p?? ? ???? ? ? ?01( ) ( ) ( )MNkkkky n b x n k a y n k??? ? ? ???山東大學(xué)生物醫(yī)學(xué)工程劉忠國(guó) 72 IIR濾波器的 DSP實(shí)現(xiàn) 對(duì) 于 因果 IIR濾波器，系統(tǒng)穩(wěn)定的條件如下： ?若 ，當(dāng) 時(shí)， ，系統(tǒng)穩(wěn)定； | | 1kp ? n ?? ( ) 0hn ?IIR濾波器結(jié)構(gòu)具有多種形式： (1)直接 I型 (也稱直接型 ) (2)直接 II型 (也稱標(biāo)準(zhǔn)型 ) (3)級(jí)聯(lián) 型 (多個(gè)二階基本 節(jié)的級(jí)聯(lián)