DSP原理与应用 实验五.doc

实验IR滤波器设计的DSP实现 反馈通道和前向通道的差分方程 ： 反馈： w(n) = x(n) + a1w(n-1) + a2w(n-2) 前向： y(n) = b0w(n) + bw(n-1)+ b2w(n-2) 计算步骤：① 计算反馈通道；② 计算前向通道；③ 输出结果y(n)；④ 重复循环。 设计一个二阶标准型低通滤波器, H(z)= ， 输入信号:x(n)=[07FFFH,0,0,0,0],画出输出波形z(n)及其幅频谱。 存储单元分配： 1)XN 5个单元，存放输入数据。 2)x0暂存x(n),w(n)或y(n)。 x1暂存w(n-1)。x2暂存w(n-2)。 3)系数表：B2,B1,B0,A2,A1(共5个单元) 4) z5个单元，存放输出数据。 ;源程序 ch5p1.asm .mmregs .def aa .bss z,5 XN .usect XN,5 ;输入数据单元 x0 .usect x,1 ;暂存x(n)，w(n)，y(n)的单元 x1 .usect x,1 ;暂存w(n-1)的单元 x2 .usect x,1 ;暂存w(n-2)的单元 B2 .usect COEF,1 ;滤波系数表数据单元 B1 .usect COEF,1 B0 .usect COEF,1 A2 .usect COEF,1 A1 .usect COEF,1 .data tableXN: .word 07FFFH,0,0,0,0 ;输入数据 tablex: .word 0,0 ;w(n-1)和w(n-2)初始值 tableBA: .word 676*32768/100000 ;系数表/10 .word 1352*32768/100000 .word 676*32768/100000 .word 4142*32768/100000 .word -14142*32768/100000 .text aa: LD #x0,DP ;设置数据页：DP=1 SSBX FRCT ;小数运算 STM #XN,AR1 RPT #4 MVPD #tableXN,*AR1+ ; 将程序存储空间中标号为”tableXN”的地址里的数据存入AR1所指向的数据存储空间，且地址AR1加1。 STM #x1,AR1 RPT #1 MVPD #tablex,*AR1+ ; 将程序存储空间中标号为”tablex”的地址里的数据存入AR1所指向的数据存储空间，且地址AR1加1。 STM #B2,AR1 RPT #4 MVPD #tableBA,*AR1+ ; 将程序存储空间中标号为”tableBA”的地址里的数据存入AR1所指向的数据存储空间，且地址AR1加1。 STM #z,AR2 ;AR2指向z STM #5,AR4 ;AR4=5 (循环次数) STM #XN,AR3 ;AR3指向XN IIR1: LD #0,B ;B=0 LD *AR3+,B ;将一个输入数据装入B STL B,@x0 ;将B中的数据存入x0 LD @x0,16,A ;将x0中的数据(x(n))左移16位，装入A LD @x1,T ; T=(x1) MAC @A1,A ;A=A+A1*(x1) LD @x2,T ;T=(x2) MAC @A2,A ;A=A+A2*(x2) STH A,@x0 ;将A中的数据(w(n))存入x0 ;以上是反馈通道计算 MPY @B2,A ;A=B2*T [即A=B2*(x2)=B2*w(n-2)] LTD @x1 ;T=(x1)