DSP应用程序设计.ppt

《DSP应用程序设计.ppt》由会员分享，可在线阅读，更多相关《DSP应用程序设计.ppt（48页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第第6章章应用程序设计应用程序设计6.1 FIR滤波器的DSP实现在数字信号处理中，滤波占有极其重要的地位。数字滤波是语音处理、图像处理、模式识别、频谱分析等应用中的基本处理算法。6.1.4 FIR滤波器的滤波器的DSP实现实现FIR滤波器的输出表达式：滤波器的输出表达式：y(n)=b0 x(n)+b1x(n-1)+bn-1 x(n-N+1)基本算法基本算法:采用采用乘法累加运算乘法累加运算。即不。即不断地输入样本断地输入样本x(n)，经过经过z-1延时后延时后,再进行乘法再进行乘法-累加，最后输出滤波结累加，最后输出滤波结果果y(n)。(指令：指令：MAC/MACD)Z-1 算法算法 的实现

2、的实现常用：常用：线性缓冲区法线性缓冲区法和和循环缓冲区法循环缓冲区法线性缓冲区法：线性缓冲区法：顶部为低地址单元，存放最新样本顶部为低地址单元，存放最新样本底部为高地址单元，存放最老样本底部为高地址单元，存放最老样本指针指针ARx指向缓冲区底部指向缓冲区底部线性缓冲区线性缓冲区只能映射到只能映射到DARAM存储器存储器（1次读和次读和1次写）次写）数据的存放：输入数据：存放在数据存储器系数表：存放在程序存储器图6.1.8常用指令：DELAYMACD例：用线性缓冲区法实现例：用线性缓冲区法实现FIR滤波器滤波器数据存储器数据存储器数据存储器数据存储器暂存暂存暂存暂存y(n)y(n)x(n)x(

3、n)x(n-1)x(n-1)x(n-2)x(n-2)x(n-3)x(n-3)x(n-4)x(n-4)x(n-5)x(n-5)x(n-6)x(n-6)程序存储器程序存储器程序存储器程序存储器b b6 6b b5 5b b4 4b b3 3b b2 2b b1 1b b0 0 x x x xCOEFCOEFCOEFCOEF双操作数寻址指令：双操作数寻址指令：MACD*AR2-MACD*AR2-,COEFCOEF,A AAR2AR2AR2AR2x(n-6)x(n-6)x(n-5)x(n-5)程序清单:.title “FIR1.ASM”.mmregs .def start x .usect “x”，8

4、 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 .textstart:SSBX FRCT STM#x+7,AR2 STM#6,AR0 LD#x+1,DP PORTR PA1,x+1FIR1:RPTZ A,#6 MACD *AR2-,COEF,A STH A,*AR2 PORTW*AR2+,PA0 BD FIR1 PORTR PA1,*A

5、R2+0 .end 循环缓冲区法：数据存储器数据存储器数据存储器数据存储器x(n)x(n)x(n)x(n)x(n-7)x(n-7)x(n-7)x(n-7)x(n-1)x(n-1)x(n-1)x(n-1)x(n-2)x(n-2)x(n-2)x(n-2)x(n-3)x(n-3)x(n-3)x(n-3)x(n-4)x(n-4)x(n-4)x(n-4)x(n-5)x(n-5)x(n-5)x(n-5)x(n-6)x(n-6)x(n-6)x(n-6)数据存储器数据存储器数据存储器数据存储器x(n)x(n)x(n)x(n)x(n+1)x(n+1)x(n+1)x(n+1)x(n-1)x(n-1)x(n-1)x

6、(n-1)x(n-2)x(n-2)x(n-2)x(n-2)x(n-3)x(n-3)x(n-3)x(n-3)x(n-4)x(n-4)x(n-4)x(n-4)x(n-5)x(n-5)x(n-5)x(n-5)x(n-6)x(n-6)x(n-6)x(n-6)ARxARxARxARx数据存储器数据存储器数据存储器数据存储器x(n)x(n)x(n)x(n)x(n+1)x(n+1)x(n+1)x(n+1)x(n-1)x(n-1)x(n-1)x(n-1)x(n-2)x(n-2)x(n-2)x(n-2)x(n-3)x(n-3)x(n-3)x(n-3)x(n-4)x(n-4)x(n-4)x(n-4)x(n-5)x

7、(n-5)x(n-5)x(n-5)x(n+2)x(n+2)x(n+2)x(n+2)ARxARxARxARxARxARxARxARx循环缓冲区法循环缓冲区法从最老样本开始取数从最老样本开始取数一个一个y值求出后，值求出后，输入最新样本替代最输入最新样本替代最老样本老样本，其他数据位置不变，其他数据位置不变双操作数寻址：双操作数寻址：DARAM（2次读）次读）数据的存放：数据的存放：输入数据：存放在数据存储器输入数据：存放在数据存储器系数表：存放在数据存储器系数表：存放在数据存储器用用循环寻址循环寻址，BK（循环缓冲区长度（循环缓冲区长度N）要求：要求：缓冲区起始地址的缓冲区起始地址的k个最低有效

8、位必须为个最低有效位必须为0，且满足，且满足2kN。例如：例如：N=31，k的最小值为的最小值为5，则缓冲，则缓冲区的起始地址：区的起始地址：XXXX XXXX XXX0 0000B例：用循环缓冲区法实现例：用循环缓冲区法实现FIRFIR滤波器滤波器 .title“FIR2.ASM”.mmregs .def start .bss y,1xn .usect“xn”,7b0 .usect“b0”,7PA0 .set 0PA1 .set 1 .datatable:.word 1*32768/10 .word 2*32768/10 .word 3*32768/10 .word 4*32768/10 .

9、word 5*32768/10 .word 6*32768/10 .word 7*32768/10 .text start:SSBX FRCT STM#b0,AR1 RPT#6 MVPD table,*AR1+STM#xn+6,AR2 STM#b0+6,AR3 STM#7,BK STM#-1,AR0 LD#xn,DP PORTR PA1,xnFIR2:RPTZ A,#6 MAC*AR2+0%,*AR3+0%,A STH A,y PORTW y,PA0 BD FIR2 PORTR PA1,*AR2 (PORTR PA1,*AR2+0%).end y?链接命令文件：链接命令文件：vectors.o

10、bjfir.obj-o fir.out-m fir.map-estartMEMORYPAGE 0：EPROM：o=E000H l=1000H VECS：o=FF80H l=0080HPAGE 1：SPRAM：o=0060H l=0020H DARAM：o=0080H l=1380H SECTIONS.text：EPROM PAGE 0 .vectors：VECS PAGE 0.data：EPRAM PAGE 0.bss：SPRAM PAGE 1xn：align(8)DARAM PAGE 1b0：align(8)DARAM PAGE 1 6.1.5 FIR滤波器的设计实例滤波器的设计实例6.1.

11、5 FIR滤波器的设计实例滤波器的设计实例 设计一个设计一个FIRFIR低通滤波器，其设计低通滤波器，其设计参数：参数：滤波器阶数：滤波器阶数：N=40N=40；截止频率：截止频率：wp=0.35，ws=0.4。滤波器系数可由滤波器系数可由MATLABMATLAB中的中的fir2fir2函数函数产生，函数调用格式：产生，函数调用格式：1.1.由给定的设计参数确定滤波器的系数由给定的设计参数确定滤波器的系数 f=0 0.3535 0.4 1；m=1 1 0 0；b=fir2(39,f,m)滤波器的幅频、相频特性滤波器的幅频、相频特性,其格式：其格式：freqz(b,512512,10001000

12、)第第第第6 6 6 6章章章章 应用程序设计应用程序设计应用程序设计应用程序设计滤波器的频率特性滤波器的频率特性滤波器的频率特性滤波器的频率特性:1.1.1.1.由给定的设计参数确定滤波器的系数由给定的设计参数确定滤波器的系数由给定的设计参数确定滤波器的系数由给定的设计参数确定滤波器的系数 第第第第6 6 6 6章章章章 应用程序设计应用程序设计应用程序设计应用程序设计 采用循环缓冲区实现采用循环缓冲区实现采用循环缓冲区实现采用循环缓冲区实现FIRFIRFIRFIR滤波器的源程序如下：滤波器的源程序如下：滤波器的源程序如下：滤波器的源程序如下：2 2 2 2汇编源程序汇编源程序汇编源程序汇编

13、源程序 .title title title title “FIRFIRFIRFIR.ASMASMASMASM”.mmregsmmregsmmregsmmregs .bssbssbssbss y y y y，1 1 1 1K_FIR_BFFRK_FIR_BFFRK_FIR_BFFRK_FIR_BFFR .set 40set 40set 40set 40PA0 PA0 PA0 PA0 .set 0set 0set 0set 0PA1 PA1 PA1 PA1 .set 1set 1set 1set 1FIR_COFF_TABLE FIR_COFF_TABLE FIR_COFF_TABLE FIR_

14、COFF_TABLE .usectusectusectusect “FIR_COFFFIR_COFFFIR_COFFFIR_COFF”,40 40 40 40 ;定义数据存储空间定义数据存储空间定义数据存储空间定义数据存储空间D_DATA_BUFFER D_DATA_BUFFER D_DATA_BUFFER D_DATA_BUFFER .usectusectusectusect “FIR_BFRFIR_BFRFIR_BFRFIR_BFR”，40 40 40 40 ;定义数据存储空间定义数据存储空间定义数据存储空间定义数据存储空间 .datadatadatadataCOFF_FIR_STARTC

15、OFF_FIR_STARTCOFF_FIR_STARTCOFF_FIR_START：.word word word word-7 7 7 7*32768/10000,332768/10000,332768/10000,332768/10000,3*32768/1000032768/1000032768/1000032768/10000 ;b b b b0 0 0 0,b b b b1 1 1 1 .word 14word 14word 14word 14*32768/10000,1032768/10000,1032768/10000,1032768/10000,10*32768/1000032

16、768/1000032768/1000032768/10000 ;b b b b2 2 2 2,b b b b3 3 3 3 .word 3word 3word 3word 3*32768/10000,32768/10000,32768/10000,32768/10000,-7 7 7 7*32768/10000 32768/10000 32768/10000 32768/10000 ;b b b b38383838b b b b39393939第第第第6 6 6 6章章章章 应用程序设计应用程序设计应用程序设计应用程序设计 .texttexttexttext .def FIR_INITdef

17、 FIR_INITdef FIR_INITdef FIR_INIT .def FIR_TASKdef FIR_TASKdef FIR_TASKdef FIR_TASKFIR_INITFIR_INITFIR_INITFIR_INIT：SSBX FRCT SSBX FRCT SSBX FRCT SSBX FRCT ;设置小数乘法设置小数乘法设置小数乘法设置小数乘法 STM#FIR_COFFSTM#FIR_COFFSTM#FIR_COFFSTM#FIR_COFF，AR5 AR5 AR5 AR5 ;AR1AR1AR1AR1指向指向指向指向b b b b0 0 0 0单元单元单元单元 RPT#K_FIR

18、_BFFRRPT#K_FIR_BFFRRPT#K_FIR_BFFRRPT#K_FIR_BFFR-1 1 1 1 ;设置传输次数设置传输次数设置传输次数设置传输次数 MVPD#COFF_FIR_STARTMVPD#COFF_FIR_STARTMVPD#COFF_FIR_STARTMVPD#COFF_FIR_START,*AR5+*AR5+*AR5+*AR5+;系数系数系数系数b b b bi i i i传输至数据区传输至数据区传输至数据区传输至数据区 STM#D_DATA_BUFFERSTM#D_DATA_BUFFERSTM#D_DATA_BUFFERSTM#D_DATA_BUFFER，AR4

19、AR4 AR4 AR4;D_DATA_BUFFERD_DATA_BUFFERD_DATA_BUFFERD_DATA_BUFFER缓冲区清缓冲区清缓冲区清缓冲区清0 0 0 0 RPTZ A RPTZ A RPTZ A RPTZ A，#K_FIR_BFFR#K_FIR_BFFR#K_FIR_BFFR#K_FIR_BFFR-1 1 1 1 STL A STL A STL A STL A，*AR4+AR4+AR4+AR4+STM#STM#STM#STM#(D_DATA_BUFFER+K_FIR_BFFRD_DATA_BUFFER+K_FIR_BFFRD_DATA_BUFFER+K_FIR_BFFRD

20、_DATA_BUFFER+K_FIR_BFFR-1 1 1 1)，AR4AR4AR4AR4 STM#STM#STM#STM#(FIR_COFF_TABLE+K_FIR_BFFRFIR_COFF_TABLE+K_FIR_BFFRFIR_COFF_TABLE+K_FIR_BFFRFIR_COFF_TABLE+K_FIR_BFFR-1 1 1 1)，AR5AR5AR5AR5第第第第6 6 6 6章章章章 应用程序设计应用程序设计应用程序设计应用程序设计 STM#STM#STM#STM#-1 1 1 1，AR0AR0AR0AR0 ;设置双操作数减量设置双操作数减量设置双操作数减量设置双操作数减量 LD

21、#D_DATA_BUFFERLD#D_DATA_BUFFERLD#D_DATA_BUFFERLD#D_DATA_BUFFER，DPDPDPDP ;设置页指针设置页指针设置页指针设置页指针 PORTR PA1PORTR PA1PORTR PA1PORTR PA1，D_DATA_BUFFERD_DATA_BUFFERD_DATA_BUFFERD_DATA_BUFFER ;输入输入输入输入x x x x(n n n n)FIR_TASKFIR_TASKFIR_TASKFIR_TASK：STM#K_FIR_BFFRSTM#K_FIR_BFFRSTM#K_FIR_BFFRSTM#K_FIR_BFFR，B

22、KBKBKBK RPTZ A RPTZ A RPTZ A RPTZ A，#K_FIR_BFFR#K_FIR_BFFR#K_FIR_BFFR#K_FIR_BFFR-1 1 1 1 ;重复操作重复操作重复操作重复操作 MAC MAC MAC MAC *AR4+0%AR4+0%AR4+0%AR4+0%，*AR5+0%AR5+0%AR5+0%AR5+0%，A A A A ;双操作数乘法双操作数乘法双操作数乘法双操作数乘法-累加累加累加累加 STH ASTH ASTH ASTH A，yyyy ;暂存暂存暂存暂存y y y y(n n n n)PORTW y PORTW y PORTW y PORTW y

23、，PA0PA0PA0PA0 ;输出输出输出输出y y y y(n n n n)BD FIR_TASK BD FIR_TASK BD FIR_TASK BD FIR_TASK ;循环循环循环循环 PORTR PA1PORTR PA1PORTR PA1PORTR PA1，*AR4+0%AR4+0%AR4+0%AR4+0%;输入最新样本，并修正输入最新样本，并修正输入最新样本，并修正输入最新样本，并修正AR3AR3AR3AR3 .endendendend第第第第6 6 6 6章章章章 应用程序设计应用程序设计应用程序设计应用程序设计 在在在在C5402DSPC5402DSPC5402DSPC5402

24、DSP硬件系统中，用户可使用的存储器硬件系统中，用户可使用的存储器硬件系统中，用户可使用的存储器硬件系统中，用户可使用的存储器资源：资源：资源：资源：3 3 3 3汇编源程序的链接命令文件汇编源程序的链接命令文件汇编源程序的链接命令文件汇编源程序的链接命令文件 程序存储空间：程序存储空间：程序存储空间：程序存储空间：片内片内片内片内0080H0080H0080H0080H 3FFFH3FFFH3FFFH3FFFH，片外片外片外片外48000H48000H48000H48000H 4FFFFH4FFFFH4FFFFH4FFFFH；数据存储空间：数据存储空间：数据存储空间：数据存储空间：片内片内片

25、内片内0080H0080H0080H0080H 3FFFFH3FFFFH3FFFFH3FFFFH，片外片外片外片外4000H4000H4000H4000H 7FFFH7FFFH7FFFH7FFFH。第第第第6 6 6 6章章章章 应用程序设计应用程序设计应用程序设计应用程序设计FIRFIRFIRFIR滤波器源程序的链接命令文件如下：滤波器源程序的链接命令文件如下：滤波器源程序的链接命令文件如下：滤波器源程序的链接命令文件如下：/*SOLUTION FILE FOR FIR SOLUTION FILE FOR FIR SOLUTION FILE FOR FIR SOLUTION FILE FOR

26、 FIR.CMD CMD CMD CMD*/vectorsvectorsvectorsvectors.objobjobjobj firfirfirfir.objobjobjobj -o firo firo firo fir.outoutoutout -m firm firm firm fir.map map map map -e fir_inite fir_inite fir_inite fir_init MEMORY MEMORY MEMORY MEMORY PAGE 0 PAGE 0 PAGE 0 PAGE 0：EPROMEPROMEPROMEPROM：org=0E000H org=0E0

27、00H org=0E000H org=0E000H lenlenlenlen=1000H=1000H=1000H=1000H VECS VECS VECS VECS：org=0FF80H org=0FF80H org=0FF80H org=0FF80H lenlenlenlen=0080H=0080H=0080H=0080H PAGE 1 PAGE 1 PAGE 1 PAGE 1：SPRAMSPRAMSPRAMSPRAM：org=0060H org=0060H org=0060H org=0060H lenlenlenlen=0020H=0020H=0020H=0020H DARAM DARA

28、M DARAM DARAM：org=0080H org=0080H org=0080H org=0080H lenlenlenlen=1380H=1380H=1380H=1380H 第第第第6 6 6 6章章章章 应用程序设计应用程序设计应用程序设计应用程序设计FIRFIRFIRFIR滤波器源程序的链接命令文件如下：滤波器源程序的链接命令文件如下：滤波器源程序的链接命令文件如下：滤波器源程序的链接命令文件如下：SECTIONSSECTIONSSECTIONSSECTIONS .texttexttexttext：EPROM PAGE 0 EPROM PAGE 0 EPROM PAGE 0 EPR

29、OM PAGE 0 .vectorsvectorsvectorsvectors：VECS PAGE 0 VECS PAGE 0 VECS PAGE 0 VECS PAGE 0 .datadatadatadata：EPRAM PAGE 0 EPRAM PAGE 0 EPRAM PAGE 0 EPRAM PAGE 0 .bssbssbssbss：SPRAM PAGE1 SPRAM PAGE1 SPRAM PAGE1 SPRAM PAGE1 FIR_BFR FIR_BFR FIR_BFR FIR_BFR：alignalignalignalign(128128128128)DARAM PAGE 1 D

30、ARAM PAGE 1 DARAM PAGE 1 DARAM PAGE 1 FIR_COFF FIR_COFF FIR_COFF FIR_COFF：alignalignalignalign(128128128128)DARAM PAGE 1 DARAM PAGE 1 DARAM PAGE 1 DARAM PAGE 1 (3)3)3)3)系数对称系数对称系数对称系数对称FIRFIRFIRFIR滤波器的实现滤波器的实现滤波器的实现滤波器的实现 系数对称系数对称系数对称系数对称FIRFIRFIRFIR滤波器具有线性相位的特性，在数字信号处理中滤波器具有线性相位的特性，在数字信号处理中滤波器具有线性相

31、位的特性，在数字信号处理中滤波器具有线性相位的特性，在数字信号处理中应用十分广泛。应用十分广泛。应用十分广泛。应用十分广泛。设滤波器设滤波器设滤波器设滤波器N=8N=8N=8N=8，若系数若系数若系数若系数b b b bn n n n=b b b bN N N N-1-1-1-1-n n n n，则为对称则为对称则为对称则为对称FIRFIRFIRFIR滤波器。其输出滤波器。其输出滤波器。其输出滤波器。其输出方程：方程：方程：方程：y y y(n n n)=)=)=b b b0 00 x x x(n n n)+)+)+b b b1 11x x x(n n n-1)+-1)+-1)+b b b2

32、22x x x(n n n-2)+-2)+-2)+b b b3 33x x x(n n n-3)-3)-3)+b b b3 33x x x(n n n-4)+-4)+-4)+b b b2 22x x x(n n n-5)+-5)+-5)+b b b1 11x x x(n n n-6)+-6)+-6)+b b b0 00 x x x(n n n-7)-7)-7)=b b0 0 x x(n n)+)+x x(n n-7)+-7)+b b1 1 x x(n n-1)+-1)+x x(n n-6)-6)+b b2 2 x x(n n-2)+-2)+x x(n n-5)+-5)+b b3 3 x x(n

33、 n-3)+-3)+x x(n n-4)-4)需要：需要：需要：需要：4 4 4 4次乘法、次乘法、次乘法、次乘法、7 7 7 7次加法次加法次加法次加法对称对称对称对称对称对称FIRFIRFIRFIRFIRFIR滤波器的实现方法：滤波器的实现方法：滤波器的实现方法：滤波器的实现方法：滤波器的实现方法：滤波器的实现方法：两个两个两个两个两个两个N/2N/2N/2N/2N/2N/2长度的循环缓冲区：长度的循环缓冲区：长度的循环缓冲区：长度的循环缓冲区：长度的循环缓冲区：长度的循环缓冲区：NewNewNewNewNewNew和和和和和和OldOldOldOldOldOld ；注意：注意：注意：注意

34、：注意：注意：循环缓冲区的首地址循环缓冲区的首地址循环缓冲区的首地址循环缓冲区的首地址循环缓冲区的首地址循环缓冲区的首地址NewNewOldOldx(n)x(n)x(n-3)x(n-3)x(n-2)x(n-2)x(n-1)x(n-1)x(n-4)x(n-4)x(n-5)x(n-5)x(n-6)x(n-6)x(n-7)x(n-7)设置循环缓冲区指针：设置循环缓冲区指针：设置循环缓冲区指针：设置循环缓冲区指针：设置循环缓冲区指针：设置循环缓冲区指针：AR1AR1AR1AR1AR1AR1指向指向指向指向指向指向NewNewNewNewNewNew区中的最新数据，区中的最新数据，区中的最新数据，区中的

35、最新数据，区中的最新数据，区中的最新数据，AR2AR2AR2AR2AR2AR2指向指向指向指向指向指向OldOldOldOldOldOld区中的最老数据；区中的最老数据；区中的最老数据；区中的最老数据；区中的最老数据；区中的最老数据；AR1AR1AR1AR1AR2AR2AR2AR2 在程序存储器中设置系数表；在程序存储器中设置系数表；在程序存储器中设置系数表；在程序存储器中设置系数表；在程序存储器中设置系数表；在程序存储器中设置系数表；系数表系数表b0b0b1b1b2b2b3b3COEFCOEFCOEFCOEFx(n)x(n)x(n-7)x(n-7)AR1AR1AR1AR1AR2AR2AR2A

36、R2系数对称系数对称系数对称系数对称系数对称系数对称FIRFIRFIRFIRFIRFIR滤波器指令：滤波器指令：滤波器指令：滤波器指令：滤波器指令：滤波器指令：格式：格式：FIRS FIRS XmemXmem，YmemYmem，PmadPmad 功能：功能：PmadPmadPmadPmadPARPARPARPAR；当当当当(RC)(RC)(RC)(RC)0 0 0 0，则，则，则，则B+AHB+AHB+AHB+AH (Pmem)(Pmem)(Pmem)(Pmem)B B B B，(XmemXmemXmemXmem)+(Ymem)+(Ymem)+(Ymem)+(Ymem)16161616A A

37、A A，PAR+1PAR+1PAR+1PAR+1PARPARPARPAR，RC-1 RC-1 RC-1 RC-1 RCRCRCRC 其中，其中，其中，其中，其中，其中，PmemPmemPmemPmemPmemPmem是通过是通过是通过是通过是通过是通过PARPARPARPARPARPAR寻址。寻址。寻址。寻址。寻址。寻址。程序清单：程序清单：程序清单：程序清单：.mmregsmmregs .def start .def start .bssbss y y,1 1x_new:.x_new:.usectusect“DATA1”“DATA1”,4 4x_old:.x_old:.usectusect“

38、DATA2”“DATA2”,4 4size .set 4 size .set 4 PA0 .set 0PA0 .set 0PA1 .set 1PA1 .set 1 .data .dataCOEF:.word 1*32768/10COEF:.word 1*32768/10 .word 2*32768/10 .word 2*32768/10 .word 3*32768/10 .word 3*32768/10 .word 4*32768/10 .word 4*32768/10y(n)数据存储器数据存储器数据存储器数据存储器x(n)x(n-3)x(n-2)x(n-1)x(n-4)x(n-5)x(n-6

39、)x(n-7)y y y ySPRAMSPRAMSPRAMSPRAMx_newx_newx_newx_newDATA1DATA1DATA1DATA1x_oldx_oldx_oldx_oldDATA2DATA2DATA2DATA2程序存储器程序存储器程序存储器程序存储器b0b1b2b3COEFCOEFCOEFCOEF系数表系数表系数表系数表 .text.text.text.text.text.textstart:start:start:start:start:start:LD#x_new,DPLD#x_new,DPLD#x_new,DPLD#x_new,DPLD#x_new,DPLD#x_new

40、,DP SSBX FRCT SSBX FRCT SSBX FRCT SSBX FRCT SSBX FRCT SSBX FRCT STM#x_new,AR1 STM#x_new,AR1 STM#x_new,AR1 STM#x_new,AR1 STM#x_new,AR1 STM#x_new,AR1 STM#x_old+STM#x_old+STM#x_old+STM#x_old+STM#x_old+STM#x_old+(sizesizesizesizesizesize-1 1 1 1 1 1),AR2,AR2,AR2,AR2,AR2,AR2 STM#size,BK STM#size,BK STM#s

41、ize,BK STM#size,BK STM#size,BK STM#size,BK STM#-1,AR0 STM#-1,AR0 STM#-1,AR0 STM#-1,AR0 STM#-1,AR0 STM#-1,AR0 PORTR PA1,x_new PORTR PA1,x_new PORTR PA1,x_new PORTR PA1,x_new PORTR PA1,x_new PORTR PA1,x_newFIR3:ADD *AR1+0%,*AR2+0%,AFIR3:ADD *AR1+0%,*AR2+0%,AFIR3:ADD *AR1+0%,*AR2+0%,AFIR3:ADD *AR1+0%,*

42、AR2+0%,AFIR3:ADD *AR1+0%,*AR2+0%,AFIR3:ADD *AR1+0%,*AR2+0%,A RPTZ B,#(size-1)RPTZ B,#(size-1)RPTZ B,#(size-1)RPTZ B,#(size-1)RPTZ B,#(size-1)RPTZ B,#(size-1)FIRS *AR1+0%,*AR2+0%,COEF FIRS *AR1+0%,*AR2+0%,COEF FIRS *AR1+0%,*AR2+0%,COEF FIRS *AR1+0%,*AR2+0%,COEF FIRS *AR1+0%,*AR2+0%,COEF FIRS *AR1+0%,

43、*AR2+0%,COEF STH B,*(y)STH B,*(y)STH B,*(y)STH B,*(y)STH B,*(y)STH B,*(y)PORTW*(y),PA0 PORTW*(y),PA0 PORTW*(y),PA0 PORTW*(y),PA0 PORTW*(y),PA0 PORTW*(y),PA0 MAR *+AR1(2)%MAR *+AR1(2)%MAR *+AR1(2)%MAR *+AR1(2)%MAR *+AR1(2)%MAR *+AR1(2)%MAR *AR2+%MAR *AR2+%MAR *AR2+%MAR *AR2+%MAR *AR2+%MAR *AR2+%MVDD

44、*AR1,*AR2+0%MVDD *AR1,*AR2+0%MVDD *AR1,*AR2+0%MVDD *AR1,*AR2+0%MVDD *AR1,*AR2+0%MVDD *AR1,*AR2+0%BD FIR3 BD FIR3 BD FIR3 BD FIR3 BD FIR3 BD FIR3 PORTR PA1,*AR1 PORTR PA1,*AR1 PORTR PA1,*AR1 PORTR PA1,*AR1 PORTR PA1,*AR1 PORTR PA1,*AR1 .end .end .end .end .end .endHIGHPASS.set 0 BANDPASS.set 1LOWPASS

45、.set 0.globalstart,fir.mmregsCOFF_FIR_START .sect coff_fir.if LOWPASS.includelowpasslowpass.inc.elseif BANDPASS.includebandpassbandpass.inc.else.includehighpasshighpass.inc.endif滤波器设计实验滤波器设计实验：K_FIR_BFFR.set32d_data_buffer.usect fir_bfr,64FIR_DP.usect fir_vars,0d_filin.usect fir_vars,1d_filout.usect

46、 fir_vars,100houtput.usect fir_vars,1.asgAR4,FIR_DATA_P.asgAR6,INBUF_P.asgAR7,OUTBUF_P.asgAR3,OUTBUF.sectfir_prog start:LD#FIR_DP,DPSTM#d_data_buffer,FIR_DATA_PRPTZ A,#K_FIR_BFFR-1STLA,*FIR_DATA_P+STM#d_filin,INBUF_PSTM#d_filout,OUTBUF_P STM#output,OUTBUFSTM#100h,BKfir_loop:NOP ;Add Breakpoint&porbe

47、 pointLD*INBUF_P,ACALL firSTHA,*OUTBUF_P+%STHA,*OUTBUFmain_end:bfir_loop fir:;SSBX SXM SSBX FRCT STM#d_data_buffer,FIR_DATA_PSTL A,*FIR_DATA_PSTM#(d_data_buffer+K_FIR_BFFR-1),FIR_DATA_P fir_task:RPTZ A,#K_FIR_BFFR-1 MACD *FIR_DATA_P-,COFF_FIR_START,ARET.endFir.cmdMEMORYPAGE 0:PROG:o=100h,l=2000h PAG

48、E 1:DATA1:o=2600h,l=100hDATA2:o=2100h,l=100hDATA3:o=2200h,l=100hDATA4:o=2300h,l=100hDATA5:o=2400h,l=100hDATA6:o=2500h,l=100hSECTIONScoff_fir:PROG PAGE 0fir_prog:PROGPAGE 0fir_vars:DATA1 PAGE 1 fir_coff:DATA2PAGE 1 fir_bfr :DATA3PAGE 1FFTFFT变换的变换的信号流程图信号流程图(N=8,N=8,共需共需3 3级蝶形运算级蝶形运算）：）：x(0)x(4)x(2)x(

49、6)x(1)x(5)x(3)x(7)WN0 x(7)x(6)x(5)x(4)x(3)x(2)x(1)x(0)-1 1-1 1-1 1WN0WN0WN0-1 1-1 1WN0WN2WN0WN2-1 1-1 1-1 1WN0WN1WN2WN3-1 1-1 1-1 1-1 1快速傅立叶变换快速傅立叶变换FFTFFT6.3 FFT的的DSP实现实现程序存储器的分配：图程序存储器的分配：图6.3.3数据存储器的分配：图数据存储器的分配：图6.3.4d_input:输入自然顺序输入自然顺序（x0r,x0i,x1r,x1i,x2r,x2i,)fft_data：比特反转输入序列：比特反转输入序列（x0r,x0

50、i,x4r,x4i,x2r,x2i,x6r,x6i,.)fft_out：FFT结果结果实现实现FFT算法的步骤：算法的步骤：实现输入数据的比特反转实现输入数据的比特反转位反转寻址（间接寻址）位反转寻址（间接寻址）实现实现N点复数点复数FFT功率谱的计算功率谱的计算输出输出FFT结果结果汇编语言程序：汇编语言程序：rfft_task：主调用程序：主调用程序bit_rev子程序：实现位码倒置子程序：实现位码倒置fft：完成：完成N点点FFT运算运算power：计算功率谱：计算功率谱bit_rev子程序子程序(N=8)bit_rev:STM#d_input,AR3STM#fft_data,AR7MV