汕头大学电子系DSP实验报告6.docVIP

汕 头 大 学 实 验 报 告 学院: 系:专业:年级:成绩:实验目的FIR/IIR滤波器的原理及使用方法； 了解使用Matlab语言设计FIR/IIR 滤波器的方法； 了解DSP 对FIR/IIR 滤波器的设计及编程方法； 熟悉对FIR/IIR 滤波器的调试方法. 实验原理 FIR滤波器的原理： 在dsp中用循环缓冲区法来实现FIR算法： IIR滤波器的原理 传递函数H(z)为： 其对应差分方程： 编程时，开辟四个缓冲区，存放输入、输出变量和滤波器的系数 问题 （一）解释程序，并分析CCS显示的实验结果 FIR： .mmregs .global start .def start,_c_int00 KS .set 256 ；定义块重复次数，及输出的个数 N .set 16 ；定义卷积的阶次 COEF_FIR .sect COEF_FIR .word 62,188,86,-764,-1453,625,6202,11439 ；定义滤波器系数 .data INPUT .copy FIRIN.INC ;输入数据在数据区0x2400 OUTPUT .space 1024 ;输出数据在数据区0x2500 x_new .usect DATA1,N/2 x_old .usect DATA2,N/2 size .set N/2 .text _c_int00 b start nop nop start: SSBX FRCT ;设置FRCT（小数方式位） STM #x_new,AR2 ;AR2指向New缓冲区第一个单元 STM #x_old+(size-1),AR3 ;AR3指向Old缓冲区最后一个单元 STM #-1,AR0 STM #INPUT,AR4 ;模拟输入数据指针AR4初始化 STM #OUTPUT,AR5 ;滤波器输出数据指针AR5初始化 STM #KS-1,BRC ;;//块重复次数 RPTBD LOOP-1 STM #size,BK ;循环缓冲区块大小BK=size LD *AR4+,A STL A,*AR2 ;输入样本值 FIR_FILTER: ADD *AR2+0%,*AR3+0%,A ;AH=x(n)+x(n-15) RPTZ B,#size-1 ;;//重复执行下一条指令size次 firs *AR2+0%,*AR3+0%,COEF_FIR;B=B+AH*h(0),AH=x(n-1)+x(n-14) STH B,*AR5+ ;保存滤波输出数据到AP5所指向单元 MAR *+AR2(2)% ;修正AR2，指向New缓冲区最老的数据替代 MAR *AR3+% ;修正AR3，指向OLD缓冲区最老的数据替代 MVDD *AR2,*AR3+0% ;用New缓冲区最老的数据替代Old缓冲区中最老的数据 LOOP: EEND B EEND .end IIR： .title diir.asm .mmregs .global start .def start, _c_int00 N .set 6 INPUT .copy iirin.inc table .word 63, 0, -188, 0, 188, 0, -63 .word 0, 11675, 0, 9663, 0, 2948 .data BN .usect BN,N+1 AN .usect AN,N+1 INBUF .usect INBUF,256 OUTPUT .usect OUTPUT,256 .text .asg AR0, INDEX_P .asg AR2, XN_P .asg AR3, ACOEF_P .asg AR4, YN_P .asg AR5, BCOEF_P _c_int00 b start nop nop start: SSBX FRCT SSBX OVM SSBX SXM STM #BN+N,AR1 RPT #N MVPD #table,*AR1- STM #AN+N-1,AR1 RPT #N-1 MVPD #table+N+1,*