1、2020年7月28日,DSP原理及應(yīng)用,1,內(nèi)容提要 數(shù)字信號(hào)處理主要面向密集型的運(yùn)算，包括乘法-累加、數(shù)字濾波和快速傅里葉變換等。 本章結(jié)合數(shù)字信號(hào)處理和通信中最常見、最具有代表性的應(yīng)用，介紹通用數(shù)字信號(hào)處理算法的DSP實(shí)現(xiàn)方法，主要包括： 有限沖激響應(yīng)（FIR）數(shù)字濾波器 無限沖激響應(yīng)（IIR）數(shù)字濾波器 快速傅里葉變換（FFT） 正弦信號(hào)發(fā)生器。,第6章 應(yīng)用程序設(shè)計(jì),2020年7月28日,DSP原理及應(yīng)用,2,6.1 FIR濾波器的DSP實(shí)現(xiàn),在數(shù)字信號(hào)處理中，濾波占有極其重要的地位。數(shù)字濾波是語音處理、圖像處理、模式識(shí)別、頻譜分析等應(yīng)用中的基本處理算法。用DSP芯片實(shí)現(xiàn)數(shù)字濾波除了具

2、有穩(wěn)定性好、精確度高、不受環(huán)境影響等優(yōu)點(diǎn)外，還具有靈活性好等特點(diǎn)。 數(shù)字濾波器是DSP的基本應(yīng)用，本節(jié)主要討論FIR濾波器的基本結(jié)構(gòu)、設(shè)計(jì)方法和DSP實(shí)現(xiàn)方法。,2020年7月28日,DSP原理及應(yīng)用,3,6.1 FIR濾波器的DSP實(shí)現(xiàn),6.1.1 FIR濾波器的基本結(jié)構(gòu),數(shù)字濾波是將輸入的信號(hào)序列，按規(guī)定的算法進(jìn)行處理，從而得到所期望的輸出序列。,(6.1.1),x(n): 輸入序列，y(n): 輸出序列，bi : 濾波器系數(shù) N: 濾波器的階數(shù)。,FIR濾波器的差分方程為：,2020年7月28日,DSP原理及應(yīng)用,4,6.1.1 FIR濾波器的基本結(jié)構(gòu),對(duì)式（6.1.1）進(jìn)行z變換，可得

3、FIR濾波器的傳遞函數(shù)：,(6.1.2),FIR濾波器的結(jié)構(gòu) ：,2020年7月28日,DSP原理及應(yīng)用,5,6.1.1 FIR濾波器的基本結(jié)構(gòu),FIR濾波器的單位沖激響應(yīng)h(n)為有限長(zhǎng)序列。,偶對(duì)稱線性相位FIR濾波器的差分方程：,N偶數(shù),(6.1.3),若h(n)為實(shí)數(shù)，且滿足偶對(duì)稱或奇對(duì)稱的條件，則FIR濾波器具有線性相位特性。 偶對(duì)稱：h(n)= h(N-1-n)； 奇對(duì)稱：h(n)= -h(N-1-n)。,2020年7月28日,DSP原理及應(yīng)用,6,6.1.1 FIR濾波器的基本結(jié)構(gòu),在數(shù)字濾波器中，F(xiàn)IR濾波器具有如下幾個(gè)主要特點(diǎn)：, FIR濾波器無反饋回路，是一種無條件穩(wěn)定系統(tǒng)

4、； FIR濾波器可以設(shè)計(jì)成具有線性相位特性。,2020年7月28日,DSP原理及應(yīng)用,7,6.1 FIR濾波器的DSP實(shí)現(xiàn),6.1.2 FIR濾波器的MATLAB設(shè)計(jì),MATLAB是一種功能強(qiáng)、效率高、便于進(jìn)行科學(xué)和工程計(jì)算的交互式軟件包，它集數(shù)值分析、矩陣運(yùn)算、信號(hào)處理和圖形顯示于一體，為用戶提供了方便、友好的界面環(huán)境。 MATLAB中的工具箱（Toolbox）包含了許多實(shí)用程序。它提供了多種FIR濾波器設(shè)計(jì)方法。 標(biāo)準(zhǔn)頻率響應(yīng)設(shè)計(jì)法fir1加窗線性相位FIR濾波器 任意頻率響應(yīng)設(shè)計(jì)法fir2,2020年7月28日,DSP原理及應(yīng)用,8,6.1 FIR濾波器的DSP實(shí)現(xiàn),6.1.3 FIR濾

5、波器的DSP實(shí)現(xiàn),FIR濾波器的輸出表達(dá)式：,y(n)=b0 x(n)+b1x(n-1)+ +bN-1x(n-N+1),（6.1.4）,bi為濾波器系數(shù)，x(n)為濾波器在n時(shí)刻的輸入，y(n)為n時(shí)刻的輸出。,基本算法: 采用乘法累加運(yùn)算。即不斷地輸入樣本x(n)， 經(jīng)過z-1延時(shí)后,再進(jìn)行乘法-累加，最后輸出濾 波結(jié)果y(n)。,2020年7月28日,DSP原理及應(yīng)用,9,6.1.3 FIR濾波器的DSP實(shí)現(xiàn),1. z-1算法的實(shí)現(xiàn),常用的方法:,線性緩沖區(qū)法 循環(huán)緩沖區(qū)法。,(1) 線性緩沖區(qū)法,又稱延遲線法。,特點(diǎn)：, 在數(shù)據(jù)存儲(chǔ)器中開辟一個(gè)N單元的緩沖區(qū)(滑窗)，用來存放最新的N個(gè)輸

6、入樣本；, 從最老樣本開始取數(shù)，每取一個(gè)數(shù)后，樣本向下移位；, 讀完最后一個(gè)樣本后，輸入最新樣本并存入緩沖區(qū)的頂部。,2020年7月28日,DSP原理及應(yīng)用,10,1. z-1算法的實(shí)現(xiàn),(1) 線性緩沖區(qū)法,緩沖區(qū):,頂部為低地址單元，存放最新樣本；,緩沖區(qū)頂部,最新樣本,底部為高地址單元，存放最老樣本；,緩沖區(qū)底部,最老樣本,指針ARx指向緩沖區(qū)底部。,ARx,2020年7月28日,DSP原理及應(yīng)用,11,(1) 線性緩沖區(qū)法,求y(n)的過程：,取數(shù)、移位和運(yùn)算:, 以ARx為指針，按x(n-7)x(n)的順序取數(shù)，每取一次數(shù)后，數(shù)據(jù)向下移一位，并完成一次乘法累加運(yùn)算；, 當(dāng)經(jīng)過8次取數(shù)

7、、移位和運(yùn)算后，得y(n)；, 求得y(n)后，輸入新樣本x(n+1)，存入緩沖區(qū)頂部單元；, 修改指針ARx，指向緩沖區(qū)的底部。,（N=8）,2020年7月28日,DSP原理及應(yīng)用,12,(1) 線性緩沖區(qū)法,求y(n)的過程：,算法：,ARx,x(n-7),ARx,x(n-6),y(n)=y7=b7x(n-7)+0,x(n-6),ARx,x(n-5),y(n)=y6=b6x(n-6)+y7,x(n-5),ARx,x(n-4),y(n)=y5=b5x(n-5)+y6,x(n-4),ARx,x(n-3),y(n)=y4=b4x(n-4)+y5,x(n-3),ARx,x(n-2),y(n)=y3

8、=b3x(n-3)+y4,x(n-2),ARx,x(n-1),y(n)=y2=b2x(n-2)+y3,x(n-1),ARx,x(n),y(n)=y1=b1x(n-1)+y2,x(n),y(n)=b0 x(n)+y1,PORTR,x(n+1),ARx,結(jié)果：, y(n),2020年7月28日,DSP原理及應(yīng)用,13,(1) 線性緩沖區(qū)法,求y(n+1)的過程：,算法：,結(jié)果：, y(n),取數(shù)順序：,x(n-6)x(n+1),x(n-5),x(n-4),x(n-3),x(n-2),x(n-1),x(n),x(n+1),最新樣本：,x(n+2),x(n+2), y(n+1), y(n+2), y(

9、n+3), y(n+4), y(n+5), y(n+6), y(n+7),ARx,2020年7月28日,DSP原理及應(yīng)用,14,(1) 線性緩沖區(qū)法,Z-1的運(yùn)算是通過執(zhí)行存儲(chǔ)器延時(shí)指令來實(shí)現(xiàn)的。 即將數(shù)據(jù)存儲(chǔ)器中的數(shù)據(jù)向較高地址單元移位來進(jìn)行 延時(shí)。,其指令：,DELAY Smem ；(Smem) Smem+1,DELAY *AR3- ；AR3指向源地址,將延時(shí)指令與其他指令結(jié)合使用，可在同樣的 機(jī)器周期內(nèi)完成這些操作。例如：,LD + DELAY LTD MAC + DELAY MACD,2020年7月28日,DSP原理及應(yīng)用,15,(1) 線性緩沖區(qū)法,注意：用線性緩沖區(qū)實(shí)現(xiàn)z-1運(yùn)算時(shí)

10、，緩沖區(qū)的數(shù)據(jù) 需要移動(dòng)，這樣在一個(gè)機(jī)器周期內(nèi)需要一次 讀和一次寫操作。因此，線性緩沖區(qū)只能定 位在DARAM中。,優(yōu)點(diǎn)：,在存儲(chǔ)器中新老數(shù)據(jù)的位置直觀明了。,2020年7月28日,DSP原理及應(yīng)用,16,1. z-1算法的實(shí)現(xiàn),(2)循環(huán)緩沖區(qū)法,特點(diǎn)：, 在數(shù)據(jù)存儲(chǔ)器中開辟一個(gè)N個(gè)單元的緩沖區(qū) (滑窗)，用來存放最新的N個(gè)輸入樣本； 從最新樣本開始取數(shù)； 讀完最老樣本后，輸入最新樣本來代替最老 樣本，而其他數(shù)據(jù)位置不變； 用BK寄存器對(duì)緩沖區(qū)進(jìn)行間接尋址，使緩沖 區(qū)地址首尾相鄰。,2020年7月28日,DSP原理及應(yīng)用,17,(2)循環(huán)緩沖區(qū)法,緩沖區(qū)：,頂層為低地址單元，存放最新樣本；,

11、x(n),底層為高地址單元，存放最老樣本；,x(n-7),x(n-1),x(n-2),x(n-3),x(n-4),x(n-5),x(n-6),ARx指向最新樣本單元。,ARx,計(jì)算過程：, 以ARx為指針，按順序取數(shù)，并修改指針；,x(n),ARx,x(n-1),ARx,x(n-2),ARx,x(n-3),ARx,x(n-4),ARx,x(n-5),ARx,x(n-6),ARx,x(n-7),ARx, 每取1次數(shù)后，完成1次乘法累加計(jì)算；,y(n)=y0=b0 x(n)+0,y(n)=y1=b1x(n-1)+y0,y(n)=y2=b2x(n-2)+y1,y(n)=y3=b3x(n-3)+y2,

12、y(n)=y4=b4x(n-4)+y3,y(n)=y5=b5x(n-5)+y4,y(n)=y6=b6x(n-6)+y5,y(n)=y7=b7x(n-7)+y6,算 法：,2020年7月28日,DSP原理及應(yīng)用,18,(2)循環(huán)緩沖區(qū)法,算 法：,計(jì)算過程：, 每取1次數(shù)后，完成1次乘法累加計(jì)算；, 求得y(n)后，輸入新樣本替代最老樣本；,x(n+1）,x(n+1), 修改指針ARx，指向最新樣本單元。,ARx, 求y(n):,取數(shù)順序：x(n)x(n-7),最新樣本：x(n+1),ARx：指向x(n+1)單元, 求y(n+1):,取數(shù)順序：x(n+1)x(n-6),最新樣本：x(n+2),A

13、Rx：指向x(n+2)單元, 求y(n+2):,取數(shù)順序：x(n+2)x(n-5),最新樣本：x(n+3),ARx：指向x(n+3)單元,2020年7月28日,DSP原理及應(yīng)用,19,(2)循環(huán)緩沖區(qū)法,循環(huán)緩沖區(qū)的優(yōu)點(diǎn)：, 緩沖區(qū)數(shù)據(jù)不需要移動(dòng)； 可以使用SARAM存儲(chǔ)器。,實(shí)現(xiàn)N個(gè)循環(huán)緩沖區(qū)單元首尾相鄰，可用BK寄存器按模間接尋址來實(shí)現(xiàn)。,常用指令：, *ARx%,; 增（減）量、按模修正ARx addr=ARx，ARx=circ(ARx1), *ARx0%,;增（減） AR0、按模修正ARx addr=ARx，ARx=circ(ARxAR0), *+ARx(1K)%,;加(1K)、按模修

14、正ARx addr=circ(ARx+1K)，ARx=circ(ARx+1K),2020年7月28日,DSP原理及應(yīng)用,20,(2)循環(huán)緩沖區(qū)法,circ是根據(jù)BK寄存器中的緩沖區(qū)長(zhǎng)度，對(duì)(ARx+1)、(ARx-1)、(ARx+AR0)、(ARx-AR0)和(ARx+1k)的值進(jìn)行取模，使指針ARx指向緩沖區(qū)，實(shí)現(xiàn)循環(huán)緩沖區(qū)首尾相鄰。,if 0index+step BK index=index+step else if index+stepBK index=index+step-BK else if index+step 0 index=index+step+BK,循環(huán)尋址 的算法：,ind

15、ex:存放在輔助寄存器中的地址指針； step:步長(zhǎng)，可正可負(fù)。,2020年7月28日,DSP原理及應(yīng)用,21,(2)循環(huán)緩沖區(qū)法, 用BK規(guī)定循環(huán)緩沖區(qū)的長(zhǎng)度N； 緩沖區(qū)起始地址的k個(gè)最低有效位必須為 0，且滿足2kN 且|step|BK 。,要求：,例如：若N=31，k的最小值為5，則緩沖區(qū)的起始地址：XXXX XXXX XXX0 0000B,若N=32，k的最小值為6，緩沖區(qū)的起始地址：XXXX XXXX XX00 0000B,2020年7月28日,DSP原理及應(yīng)用,22,(2)循環(huán)緩沖區(qū)法,例如：(BK)=N=8，(AR1)=0060h，用*AR1+%間接尋址。,第一次尋址后，AR1指

16、向0061h單元； 第二次尋址后，AR1指向0062h單元； 第八次尋址后，AR1指向0068h單元； 將BK按8取模，AR1回到0060h單元。,2020年7月28日,DSP原理及應(yīng)用,23,6.1.4 FIR濾波器的DSP實(shí)現(xiàn),2. FIR濾波器的實(shí)現(xiàn),C54x提供的MAC指令和循環(huán)尋址方式，可使FIR數(shù)字濾波器在單周期內(nèi)完成每個(gè)樣值的乘法 -累加計(jì)算。而每個(gè)樣值的乘法-累加計(jì)算,可采用RPTZ和MAC指令結(jié)合循環(huán)尋址方式來實(shí)現(xiàn)。,為了實(shí)現(xiàn)對(duì)應(yīng)項(xiàng)乘積運(yùn)算，輸入的樣值x(n)和濾波系數(shù)bi必須合理的存放，并正確初始化存儲(chǔ)塊和塊指針。樣值x(n)和濾波系數(shù)bi的存放可用線性緩沖區(qū)或循環(huán)緩沖區(qū)實(shí)

17、現(xiàn)。,2020年7月28日,DSP原理及應(yīng)用,24,2. FIR濾波器的實(shí)現(xiàn),(1) 用線性緩沖區(qū)實(shí)現(xiàn)FIR濾波器,設(shè)N=7，F(xiàn)IR濾波器的算法：,y(n)=b0 x(n)+b1x(n-1)+b5x(n-5) +b6x(n-6),x,b,雙操作數(shù)尋址指令：MACD *AR1-,b,A,功能：A=A+(AR1)(b), AR1-1AR1,(AR1)(AR1+1),AR1,x(n-6),b6,x(n-6),b6x(n-6),+A,b6x(n-6)+A,AR1,x(n-5),2020年7月28日,DSP原理及應(yīng)用,25,(1) 用線性緩沖區(qū)實(shí)現(xiàn)FIR濾波器,程序清單： .title “FIR1.AS

18、M” .mmregs .def start x .usect “x”，8 PA0 .set 0 PA1 .set 1 .data COEF： .word 1*32768/10 .word 2*32768/10 .word -4*32768/10 .word 3*32768/10 .word -4*32768/10 .word 2*32768/10 .word 1*32768/10,；自定義數(shù)據(jù)空間,x,暫存y(n),x(n-1),x(n-2),x(n-3),x(n-4),x(n-5),x(n-6),定義系數(shù)bi,COEF,b6,；定義b6=0.1,；定義b5=0.2 ；定義b4=-0.4 ；定

19、義b3=0.3 ；定義b2=-0.4 ；定義b1=0.2 ；定義b0=0.1,b5,b4,b3,b2,b1,b0,x(n),2020年7月28日,DSP原理及應(yīng)用,26,(1) 用線性緩沖區(qū)實(shí)現(xiàn)FIR濾波器,程序清單： .text start: SSBX FRCT STM #x+7,AR2 STM #6,AR0 LD #x+1,DP PORTR PA1,x+1 FIR1: RPTZ A,#6 MACD *AR2-,COEF,A STH A,*AR2 PORTW *AR2+,PA0 BD FIR1 PORTR PA1,*AR2+0% .end,;設(shè)置小數(shù)乘法,;設(shè)置AR2,AR2,;設(shè)置復(fù)位值A(chǔ)

20、R0=6,;設(shè)置緩沖區(qū)首地址,;輸入x(n),x(n),;A清0,設(shè)置迭代次數(shù),00 0000 0000,;7次乘法累加和移位,b6x(n-6),AR2,x(n-5),A+b5x(n-5),AR2,x(n-4),A+b4x(n-4),AR2,x(n-3),A+b3x(n-3),AR2,x(n-2),A+b2x(n-2),AR2,x(n-1),A+b1x(n-1),x(n),A+b0 x(n),AR2,AR2,y(n-1),;暫存y(n),y(n),;輸出y(n),修改AR2,AR2,;循環(huán),;輸入最新數(shù)據(jù), 修改AR2=AR2+AR0,x(n+1),AR2,2020年7月28日,DSP原理及應(yīng)

21、用,27,2. FIR濾波器的實(shí)現(xiàn),(2) 用循環(huán)緩沖區(qū)實(shí)現(xiàn)FIR濾波器,設(shè)N=7，F(xiàn)IR濾波器的算法：,y(n)=b0 x(n)+b1x(n-1)+b5x(n-5) +b6x(n-6),y,b0,xn,2020年7月28日,DSP原理及應(yīng)用,28,(2) 用循環(huán)緩沖區(qū)實(shí)現(xiàn)FIR濾波器,程序清單： .title “FIR2.ASM” .mmregs .def start .bss y,1 xn .usect “xn”,7 b0 .usect “b0”,7 PA0 .set 0 PA1 .set 1 .data table: .word 1*32768/10 .word 2*32768/10 .

22、word 3*32768/10 .word 4*32768/10 .word 5*32768/10 .word 6*32768/10 .word 7*32768/10,;源文件標(biāo)題,;定義MMR寄存器符號(hào)名,;定義模塊,;給y保留1個(gè)空間,y,;給xn段保留7個(gè)空間,xn,;給b0段保留7個(gè)空間,b0,;PA0賦值為0,;PA1賦值為1,;從ROM的table定義數(shù)據(jù),;定義0.1,;定義0.2,;定義0.3,;定義0.4,;定義0.5,;定義0.6,;定義0.7,0.1,0.2,0.3,0.4,0.5,0.6,0.7,2020年7月28日,DSP原理及應(yīng)用,29,.text start: S

23、SBX FRCT STM #b0,AR1 RPT #6 MVPD table,*AR1+ STM #xn,AR2 STM #b0,AR3 STM #7,BK STM #1,AR0 LD #xn,DP PORTR PA1,xn FIR2: RPTZ A,#6 MAC *AR2+0%,*AR3+0%,A STH A,y PORTW y,PA0 BD FIR2 PORTR PA1,*AR2(-1)% .end,;設(shè)置小數(shù)乘法,;AR1指向b0,;設(shè)置傳輸次數(shù),;系數(shù)傳輸至數(shù)據(jù)區(qū),0.1,0.2,0.3,0.4,0.5,0.6,0.7,;AR2指向x(n)單元,;AR3指向b0單元,;設(shè)置緩沖區(qū)長(zhǎng)度,

24、;設(shè)置雙操作數(shù)增量,;設(shè)置頁指針,;輸入x(n),;A清0,設(shè)置迭代次數(shù),;雙操作數(shù)乘法累加,;存儲(chǔ)y(n),;輸出y(n),;循環(huán),;輸入最新x(n+1)，修正AR2,2020年7月28日,DSP原理及應(yīng)用,30,鏈接命令文件： FIR2.obj vectors.obj -o FIR2.out -m FIR2.map -e start MEMORY PAGE0: EPROM:org=0E000h, len=1000h VECS: org=0FF80h, len=0080h PAGE1: SPRAM:org=0060h, len=0020h DARAM:org=0080h, len=1380h

25、 ,;選定的目標(biāo)文件 ;生成FIR2的輸出文件 ;生成FIR2的存儲(chǔ)器映像文件 ;定義源程序的入口地址 ;定義目標(biāo)存儲(chǔ)器空間 ;第0頁：程序存儲(chǔ)器 ;EPROM的起始地址：E000h 長(zhǎng)度：4K ; VECS的起始地址：FF80h 長(zhǎng)度：0080h ;第1頁：數(shù)據(jù)存儲(chǔ)器 ; SPRAM的起始地址：0060h 長(zhǎng)度：0020h ; DARAM的起始地址：0080h 長(zhǎng)度：1380h,存儲(chǔ)空間,0E000,4k,0EFFF,0FF80,80h,0FFEF,0060,20h,007F,0080,1380h,13EF,2020年7月28日,DSP原理及應(yīng)用,31,鏈接命令文件：,SECTIONS .t

26、ext:EPROM PAGE 0 .data:EPROM PAGE 0 .bss :SPRAM PAGE 1 xn:align(8)DARAM PAGE 1 b0:align(8)DARAM PAGE 1 .vections:VECS PAGE 0 ,;在存儲(chǔ)器中定義輸出段的位置 ;text段定位在程序存儲(chǔ)器 即源程序位于程序存儲(chǔ)器 ;系數(shù)區(qū)定義在程序存儲(chǔ)器 ;bss段定義在暫存器 ; 從xn起8個(gè)單元定義在DARAM ;從b0起8個(gè)單元定義在DARAM ;vections定義在VECS區(qū),2020年7月28日,DSP原理及應(yīng)用,32,2. FIR濾波器的實(shí)現(xiàn),(3) 系數(shù)對(duì)稱FIR濾波器的實(shí)

27、現(xiàn),系數(shù)對(duì)稱FIR濾波器具有線性相位的特性，在數(shù)字信號(hào)處理中應(yīng)用十分廣泛，常用于相位失真要求較高的場(chǎng)合。,設(shè)濾波器N=8，若系數(shù)bn=bN-1-n，則為對(duì)稱FIR濾波器。其輸出方程：,y(n) = b0 x(n)+b1x(n-1)+b2x(n-2)+b3x(n-3) +b3x(n-4)+b2x(n-5)+b1x(n-6)+b0 x(n-7),= b0 x(n)+x(n-7) + b1 x(n-1)+x(n-6) +b2 x(n-2)+x(n-5) +b3 x(n-3)+x(n-4) ,需要： 4次乘法 7次加法,2020年7月28日,DSP原理及應(yīng)用,33,(3) 系數(shù)對(duì)稱FIR濾波器的實(shí)現(xiàn),

28、對(duì)稱FIR濾波器的實(shí)現(xiàn)方法：, 在RAM中開辟兩個(gè)N/2長(zhǎng)度的循環(huán)緩沖區(qū)New和Old ，分別存放N/2個(gè)新數(shù)據(jù)和老數(shù)據(jù)；,x(n),x(n-3),x(n-2),x(n-1),x(n-4),x(n-5),x(n-6),x(n-7), 設(shè)置循環(huán)緩沖區(qū)指針： AR1指向New區(qū)中的最新數(shù)據(jù)， AR2指向Old區(qū)中的最老數(shù)據(jù)；,AR1,AR2, 在程序存儲(chǔ)器中設(shè)置系數(shù)表；,b0,b1,b2,b3,COEF, 進(jìn)行(AR1)+(AR2)AH加法運(yùn)算，并修改數(shù)據(jù)指針， AR1-1AR1，AR2-1AR2；,x(n),x(n),x(n-7),x(n-7),x(n) + x(n-7),x(n) + x(n-

29、7),AR1,AR2,2020年7月28日,DSP原理及應(yīng)用,34,對(duì)稱FIR濾波器的實(shí)現(xiàn)方法：, 累加器B清0，完成塊操作，重復(fù)執(zhí)行4次；,乘法累加：(AH)bi +BB；,修改系數(shù)指針：PAR+1PAR；,加法運(yùn)算：(AR1)+(AR2)AH；,修改數(shù)據(jù)指針：AR1-1AR1，,AR2-1AR2；, 保存和輸出結(jié)果;, 修正數(shù)據(jù)指針，,AR1指向New區(qū)的最老數(shù)據(jù)；,AR2指向Old區(qū)的最老數(shù)據(jù)。, 用New區(qū)的最老數(shù)據(jù)替代Old區(qū)的最老數(shù)據(jù)，輸入新數(shù)據(jù)替代New區(qū)的最老數(shù)據(jù) ；, 重復(fù)執(zhí)行 。,2020年7月28日,DSP原理及應(yīng)用,35,(3) 系數(shù)對(duì)稱FIR濾波器的實(shí)現(xiàn),系數(shù)對(duì)稱FIR濾波器指令：,格式： FIRS Xmem，Ymem，Pmem,功能： PmemPAR； 當(dāng)(RC)0，則B+AH(Pmem)B， (Xmem)+ (Ymem)16A， PAR+1PAR，(RC)-1 RC,其中， Pmem是通過PAR尋址。,2020年7月28日,DSP原理及應(yīng)用,36,(3) 系數(shù)對(duì)稱FIR濾波器的實(shí)現(xiàn),程序清單