DSP实验一有限冲击响应滤波器.docVIP

实验一 有限冲击响应滤波器（FIR）算法实验 一、实验目的 （1）握用窗函数法设计FIR数字滤波器的原理和方法； （2）熟悉线性相位FIR数字滤波器特性； （3）了解各种窗函数对滤波器特性的影响。 二、实验设备 计算机，DSP实验箱， ICETEK-USB5100仿真器。 三、实验原理 （1）有限冲击响应数字滤波器的基础理论； （2）模拟滤波器原理（巴特沃斯滤波器、切比雪夫滤波器、椭圆滤波器、贝塞尔滤波器）； （3）数字滤波器系数的确定方法。 （4）源程序及注释 .include 240x.h .global start,FIR .data .include lowpass.inc ;低通滤波器系数（32个） N .set 32 ;32点FIR滤波 DPTRI .set 8000h ;输入波形地址 DPTRO .set 0a000h ;输出波形地址 i .set 1024 ;输入波形点数 .sect .vectors b start .bss buff,N ;系数缓存 .bss input,1 ;当前输入数据 .bss output,1 ;当前输出数据 .text start: ldp #0e0H splk #6fH,WDCR splk #5555H,WDKEY splk #0aaaaH,WDKEY ;关闭看门狗中断 splk #81feH,SCSR1 ;设置DSP运行频率40m ldp #buff splk #0,input splk #0,output mar *,ar0 lar ar0,#buff zac ;acc低16位清0 rpt #N-1 ;循环32次 sacl *+ ;系数缓存清零 lar ar2,#DPTRI ;ar2- 输入存储区 lar ar3,#DPTRO ;ar3- 输出缓冲区 lar ar7,#i ;ar7 处理波形点数 mar *,ar2 ;指定当前辅助存储器为ar2 done: lacl *+ ;读入当前输入值 sacl input sacl buff ;存入系数缓存 call FIR ;调用滤波子程序完成滤波处理 mar *,ar3 sach output,1 sach *+,1,7 ;滤波后的结果右移1位存到输出缓冲区 banz done,*-,2 ;点数-1，如未处理完循环 b $ ;子程序完成32点FIR滤波 FIR: mar *,ar1 lar ar1,#buff+N-1 mpyk #0 ;乘法器清零 zac ;累加器请零 rpt #N-1 ;计算32个累加和 macd coef,*- ;系数相乘、累加运算 apac ret .end 四、实验内容 1、编译并下载程序 2、打开观察窗口 选择菜单“View”、“Graph”、“Time/Frequency…”进行如下设置： 选择菜单“View”、“Graph”、“Time/Frequency…”进行如下设置： 在弹出的图形窗口中单击鼠标右键，选择“Clear Display”。 设置断点 在标号“done”语句设置软件断点（Toggle breakpoint）。 选择菜单“File”、“Data”、“Load…”；选中F:\ FIR目录中的文件sine.dat， 单击“打开”按钮；在“Load File into Memory”窗口中的“Address”项中输入DPTRI，在“Length”项中输入1024，单击“OK”并等待读入数据。 运行并观察结果 选择“Debug”菜单的“Animate”项，或按F12键运行程序； 观察“Input”、“Output”窗口中时域图形；观察滤波效果； 鼠标右键单击“Input”和“Output”窗口，选择“Properties…”项，设置“Display Type”为“FFT Magitude”，再单击“OK”按钮结束设置； 观察“Input”、“Output”窗口中频域图形；理解滤波效果。 停止程序运行并退出。 五、实验结果 输入波形为一个低频率的正弦波与一个高频的正弦波叠加而成。 实验现象时域图： 频域图： 通过观察频域和时域图，得知：输入波形中的低频波形通过了滤波器，而高频部分则被滤除。 实验三 数据存取实验 一、实验目的 （1）了解TMS320LF2407的内部存储器空间的分配及指令寻址方式； （2）学习用CodeComposer修改、填充DSP内存单元的方法； （3）学习操作TMS320LF240x内存空间的指令。 二、实验设备 计算机，DSP实验箱，ICETEK-USB5100仿真器。 三、实验内容 将第0页数据存储区60H~6FH中的数据转存到第4页数据存