数字滤波器的DSP实现实验.doc.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.PPTFIR滤波器的设计与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数字滤波器的原理。 熟悉IIR数字滤波器从设计到实现的具体过程及其应用。 学习定点DSP在算术运算时的定标方法。 了解Matlab在DSP中的应用及在DSP中直接调用由Matlab生成的数据文件的方法。 掌握IIR数字滤波的DSP实现方法。 6. 理解汇编语言编写dsp程序的方法 二、实验设备 计算机，CC4.1版软件，Matlab软件 三、实验原理 数字滤波就是将输入的信号序列按规定的算法处理之后，得到所希望的输出序列的过程。数字滤波的作用是滤除信号中某一部分频率分量。信号经过滤波处理，就相当于信号频谱与滤波器的频率响应相乘的结果。从时域来看，就是输入信号与滤波器的冲激响应作卷积和。数字滤波器在各种???域由广泛的应用，例如数字音响、音乐和语音合成、噪声消除、数据压缩、频率合成、谐波消除、过载检测、相关检测等。 无限冲激响应IIR数字滤波器一般具有如下差分方程的形式： （5.1） 其中x(n)是输入序列，y(n)是输出序列，bk和ak为滤波器系数，N是滤波器的阶数。 其系统函数可表示为： （5.2） IIR数字滤波器含有反馈支路（环路），一般呈现递归结构，即时刻n的输出不仅取决于n时刻的输入及以前的M个输入，还要取决于以前的N个输出。其单位脉冲响应h(n)是无限长的。 IIR数字滤波器是通过对采样数据进行数学运算来达到频域滤波的目的。其实现方法有两种：一种是频域法，即利用FFT算法对输入信号进行离散傅立叶变换，分析其频谱，然后根据所希望的频率特性进行滤波，再利用傅立叶反变换恢复出时域信号；另一种方法是时域法，这种方法是对离散抽样数据作差分数学运算来达到滤波目的。本实验采用的就是基于时域法来实现IIR数字滤波器。 IIR数字滤波器具有直接型、级联型、并联型和格型等多种结构形式。直接II型的IIR滤波器结构如图5.2所示。本实验就是采用标准型（即直接II型）结构实现的，其表达式如下所示。 （5.3） （5.4） 图5.2 直接II型的IIR滤波器结构 (5.3)式执行直接II型IIR滤波器的反馈环节，其结果存储于w(n)中。(5.4)式计算前馈环节，完成IIR滤波并将结果置于y(n)中。这两个式子就是本实验实现IIR滤波器的编程依据。使用LTD和MPY指令可以实现这两个式子中的乘加运算。 与较高阶的FIR滤波器相比，达到相同性能所需要的IIR滤波器阶数一般较低。因此，在IIR滤波器实现中所需要的MIPS较低。但是，IIR滤波器的定点编程不像FIR滤波器那样简单。IIR滤波器的直接形式实现对于系数量化极为敏感，大多数实际实现采用二阶级联来降低量化误差（每一个二阶基本节又采用直接型结构实现）。在用DSP来实现时，第一节的输出y1(k)可以存放在累加器ACC中作为第二节的输入，从而避免中间结果的截尾误差。一般情况下，应用缩放来防止溢出、应用级联结构、