[计算机软件及应用]第6章 应用程序设计531.pptVIP

2003.2.16 DSP原理及应用 第6章 应用程序设计 2. FIR滤波器的实现 (2) 用循环缓冲区实现FIR滤波器 设N=7，FIR滤波器的算法： y(n)=b0x(n)+b1x(n-1)+…+b5x(n-5) +b6x(n-6) DARAM y(n) x(n) x(n-1) x(n-2) x(n-3) x(n-4) … x(n-5) x(n-6) DARAM b0 b1 b2 b3 b4 b5 b6 y b0 xn 第6章 应用程序设计 (2) 用循环缓冲区实现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 .word 3*32768/10 .word 4*32768/10 .word 5*32768/10 .word 6*32768/10 .word 7*32768/10 DARAM y(n-1) x(n) x(n-1) x(n-2) x(n-3) x(n-4) … b0 b1 b2 b3 b4 … x(n-5) x(n-6) b5 b6 ;源文件标题 ;定义MMR寄存器符号名 ;定义模块 ;给y保留1个空间 y ;给xn段保留7个空间 xn ;给b0段保留7个空间 b0 ;PA0赋值为0 ;PA1赋值为1 ;从ROM的table定义数据 ;定义0.1 ;定义0.2 ;定义0.3 ;定义0.4 ;定义0.5 ;定义0.6 ;定义0.7 ROM table 0.1 0.2 0.3 0.4 0.5 0.6 0.7 第6章 应用程序设计 .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,@xn FIR2: RPTZ A,#6 MAC *AR2+0%,*AR3+0%,A STH A,@y PORTW @y,PA0 BD FIR2 PORTR PA1,*AR2+0% .end DARAM y(n-1) x(n) x(n-1) x(n-2) x(n-3) x(n-4) … b0 b1 b2 b3 b4 … x(n-5) x(n-6) b5 b6 y xn b0 … A ;设置小数乘法 ;AR1指向b0 AR1→ ;设置传输次数 ROM table 0.1 0.2 0.3 0.4 0.5 0.6 0.7 ;系数传输至数据区 0.1 0.1 AR1→ 0.2 0.2 AR1→ 0.3 0.3 AR1→ 0.4 0.4 AR1→ 0.5 0.5 AR1→ 0.6 0.6 AR1→ 0.7 0.7 AR1→ ;AR2指向x(n-6)单元 AR2→ ;AR3指向b6单元 AR3→ ;设置缓冲区长度 ;设置双操作数增量 ;设置页指针 ;输入x(n) x(n) ;A清0,设置迭代次数 00 0000 0000 ;双操作数乘法累加 x(n-6) 0.7 0.7x(n-6)+0 AR2→ AR3→ x(n-5) 0.6 0.6x(n-5)+A AR2→ AR3→ x(n-4) 0.5 0.5x(n-4)+A AR2→ AR3→ x(n-3) 0.4 0.4x(n-3)+A AR2→ AR3→ x(n-2) 0.3 0.3x(n-2)+A AR2→ AR3→ x(n-1) 0.2 0.2x(n-1)+A AR2→ AR3→ x(n) 0.1 0.1x(n)+A AR2→ AR3→ ;存储y(n) y(n) ;输出y(n) ;循环 ;输入最新x(n+1)，修正AR2 x(n+1) AR2→ 第6章 应用程序设计 链接命令文件： FIR2.obj vectors.obj -o FIR2.out -m FIR2.map -e star