最新DSP原理与应用.pptVIP

图7-6 的信号流图 * hhjkl 图7-7 二阶IIR滤波器的直接型Ⅱ的实现 由于直接型Ⅱ对于给定的传递函数具有最小可能的延迟数、加法器数和乘法器数，所以被称为标准型。 * hhjkl 7.4.2 高阶IIR滤波器的结构 高阶IIR滤波器的差分方程和系统函数分别为： * hhjkl 图7-8 高阶IIR滤波器的直接型Ⅱ实现（L=M+1） * hhjkl 图7-9 高阶IIR滤波器的串联型结构 图7-10 第k个二阶节 * hhjkl 对于 ，有 其中， 注意： * hhjkl 7.4.3 IIR滤波器的C语言实现 例7-11,采用二维数组编写的IIR滤波器C语言程序。 temp=xin; /*xin为IIR滤波器的输入 */ for(k=0;kN_IIR;k++) /*N_IIR为IIR滤波器二阶节的个数 { w[0][k]=temp-a[1][k]*w[1][k]-a[2][k]*w[2][k]; /*这里,temp为本二阶节的输入,也是上一级二阶节的输出*/ temp=b[0][k]*w[0][k]+b[1][k]*w[1][k]+b[2][k]*w[2][k]; /*这里,temp为本二阶节的输出,也是下一级二阶节的输入*/ w[2][k]=w[1][k]; w[1][k]=w[0][k]; } xoutput=temp;/*xoutput为IIR滤波器的输出 */ a[i][k]（i=1,2）和b[i][k] （i=0,1,2）分别为第k（k=0,1,K）个二阶IIR滤波器环节的系数。 w[i][k] (i=0,1,2)是第k个二阶IIR滤波器环节中对应于图7-10中(m= 0, 1,2)的中间信号。 * hhjkl 7.4.4 IIR滤波器的汇编语言实现 例7-12，采用指针编写的IIR滤波器汇编语言程序。 ;iir.asm ;函数原型: ;void iir(int *,unsigned int,int *,int *,unsigned int,int *); ;入口参数： ;参数0:AR0 - 输入信号缓冲区指针 ;参数1:T0 - 输入信号缓冲区的样本数 ;参数2:AR2 - 输出信号缓冲区指针 ;参数3:AR1 – IIR滤波器系数数组指针 ;参数4:T1 - 二阶IIR节的个数 ;参数5:AR3 – 延迟线指针 * hhjkl ;第k节IIR滤波器： ;w(n)=x(n)-a1k*w(n-1)-a2k*w(n-2) ;y(n)=b0k*w(n)+b1k*w(n-1)+b2k*w(n-2) ; ;存储器分配: ; 暂存单元[2* N_IIR] 系数 [5* N_IIR] ; AR3 - w1i AR7 - a1k ; w1j a2k ; : b2k ; w2i b0k ; w2j b1k ; : : ; 标定：Q14 格式 * hhjkl .global _iir .sect iir_code _iir pshm ST1_55 ; 保存ST1,ST2,ST3 pshm ST2_55 pshm ST3_55 psh T3 ; 保存T3 pshboth XAR7 ; 保存AR7 or #0x340,mmap(ST1_55) ;设置FRCT,SXMD,SATD bset SMUL ;置位SMUL sub #1,T0 ; 样本数 - 1 mov T0,BRC0 ; 设置外循环计数器 sub #1,T1,T0 ; 二阶节个数 - 1 mov T0,BRC1 ; 设置内循环计数器 mov T1,T0 ;设置循环缓冲区大小 sfts T0,#1 mov mmap(T0),BK03 ;BK03=2*二阶节个数 sfts T0,#1 add T1,T0 * hhj