试验二：数字滤波器的DSP实现试验所给试验例程：IIR_FIR文件夹.DOC

实验二：数字滤波器的DSP实现实验 所给实验例程： IIR_FIR文件夹 Demo_IIR文件夹 IIR4_asm文件夹 或 IIR4_C文件夹 采用宏替换的方法实现的4阶IIR滤波器（采用直接型） IIR_asm_m文件夹：信号及滤波器系数产生、滤波器设计效果演示的matlab程序 IIR4_CAS_asm文件夹 对上面直接型的改进（采用级联型） Cca文件夹 采用C语言调用4阶IIR汇编程序 IIR_CAS_asm_m文件夹：信号及滤波器系数产生的matlab程序 阶次N的确定.ppt 滤波器设计时的采用matlab确定滤波器阶数的方法介绍 Demo_FIR文件夹 FIR15_C文件夹 采用C语言实现15阶FIR滤波器（截止频率为350Hz），应按照实验1的方法通过探点引人信号 Matlab文件夹：信号及滤波器系数产生、滤波器设计效果演示的matlab程序 Lab2_FIR.PPT FIR滤波器的设计与DSP实现的介绍 exp10a_IIR 对应“实验2a IIR滤波器的DSP实现”，为4阶IIR滤波器实现的汇编语言程序 exp11a_FIR_line 对应“实验2b FIR滤波器的DSP实现”，为11阶FIR滤波器实现的汇编语言程序（采用线性缓冲区实现） 注：重点做蓝色字体对应的实验。即做IIR的DSP实现实验时先做IIR4_asm（或IIR4_C），再做exp10a_IIR；做FIR的DSP实现实验时先做FIR15_C，再做exp11a_FIR_line。实验的原理、步骤及结果观察方法参考下面的实验指导。 实验2a IIR滤波器的DSP实现 一、实验目的四、实验内容与步骤IIR数字滤波器的内容。 2. 运行所给的例程Demo_IIR\IIR4_asm，对滤波结果进行观察。理解并读懂例程： （1）如何观察滤波结果？为什么要观察数据存储器0x200与0x300处？如果要使输出波形在数据存储器0x1000处开始，如何修改？ （2）如果要对512点输入信号滤波，如果修改？ （3）如果要编写5阶512点IIR程序，如果修改？ （4）如果本实验要对f1=5KHz、f2=30KHz、fs=200KHz混叠信号滤波，如果修改？ （5）所给例程采用宏调用实现，若采用循环方法或函数调用方法实现，该如何实现？ 2． 给定IIR数字低通滤波器的主要技术指标：通带截止频率1kHz，通带最大衰减=4dB，阻带截止频率1.5kHz，阻带最小衰减15dB，采样频率为10kHz。要求用双线性变换法设计IIR数字低通滤波器。 3. 打开Matlab，在命令窗口下输入FDATool，将出现滤波器设计窗口，如图5.3所示。按照图5.3所示对滤波器的参数进行设置，得到IIR滤波器的系数。 图5.3 滤波器设计主界面 点击图标，观察该滤波器的极零图bd =[ 0.0052 0.0207 0.0311 0.0207 0.0052 ] ad =[ 1.0000 -2.3358 2.2608 -1.0229 0.1807 ] 然后将滤波器的各个系数用213进行定标，并将其取整后转换成十六进制数。这可以在matlab的命令窗口下输入“BD=dec2hex(fix(bd.*8192))”获得，其结果如下： BD =[ 2A 0A9 0FE 0A9 2A ] 即b0=02Ah，b1＝0A9h，b2=0FEh，b3=0A9h，b4=02Ah。 同理，也将滤波器系数的分母系数取整后转换成十六进制数。 也可以用Matlab语言通过编程来设计IIR低通滤波器，见参考程序：配套实验10“exp10a_IIR /matlab”。 图5.4 极零图 图5.5 滤波器系数导出 注意：若为负数，则应该先不管符号，定标并取整后转换成二进制数，再求反加1，变成补码；正数补码即为其本身。例如， -a1=2.3358 2.3358×213＝4ABEh -a2=-2.2608 2.2608×213＝4858h 其补码为-2.2608＝0B7A8h -a3=1.0229 1.0229×213＝20BBh -a4=-0.1807 0.1807×213＝5C8h 其补码为-0.1807＝0FA38h 4. 利用Matlab生成一个由频率分别是f1=500Hz和f2=2000Hz，幅度分别是5和2.5的两正弦波叠加构成的混叠信号，并将其离散化（采样频率fs=10kHz），再取128点该离散序列，存入数据文件testdata.dat中（参见下面