FIR滤波器程序的编写.docVIP

目录 1. 引言 2. FIR滤波器基本原理 2.1两种FIR滤波器的框图及原理 2.2用DSP实现的数字FIR滤波器的优点 3. FIR滤波器的设计方法 4. FIR滤波器的c54x实现 5. FIR滤波器的实例设计 5.1编写滤波器序列 5.2汇编源程序主程序的编写 5.3中断向量表的编写 5.4相关的连接命令文件 5.5程序在CCS C500C中编译 6. 设计心得 7. 参考文献  1. 引言 在数字信号处理中，滤波占有及其重要的地位。数字滤波是谱分析、 雷达信号处理、通信信号处理等应用中的基本处理算法，它能够满足滤 波器对幅度和相位的严格要求，解决了模拟滤波器所无法克服的电压漂 移、温度漂移和噪声等问题，同时有很高的可编程性和灵活性。数字滤 波是DSP最基本的应用领域，一个DSP芯片执行数字滤波算法的能力反 应了这种芯片的功能强弱，用 DSP实现的数字滤波器可以十分方便地改 变滤波器的特性。 2. FIR滤波器基本原理 2.1两种FIR滤波器的框图及原理 对于FIR数字滤波器来说，它的实现一般采用非递归的算法，图2.1 给出的是直接型和转置型FIR滤波器的机构图。 (b)转置型横向滤波器 图2.1直接型和转置型 FIR滤波器结构图 FIR 数字滤波器的差分方程表达式为： y( n)= En h(i)x( n-i) 其中N为FIR滤波器的阶数，FIR滤波器的单位脉冲响应h(n)是一个有 限长序列。 在数字信号处理应用中往往需要设计线性相位的滤波器， FIR 滤波器在 保证幅度特性满足技术要求的同时，很容易做到严格的线性相位特性。 为了是滤波器满足线性相位条件，要求其单位脉冲响应 h(n) 为实序列， 且满足偶对称或奇对称条件，即 h(n)=h(N-1-n) 或 h(n)=-h(N-1-n) 。这 样，当N为偶数时，偶对称线性相位FIR滤波器的差分方程表达式为： y(n)=刀N/2-1 h(i)(x( n-i)+x(N-1- n-i)) 由上可见， FIR 滤波器不断地对输入样本 x(n) 延时后，再作乘法累 加运算，将滤波结果y(n)输出，因此，FIR实际上是一种乘法累加运算。 而对于线性相位 FIR 而言，利用线性相位 FIR 滤波器系数的对称特性， 可以采用结构精简的 FIR 结构将乘法器数目减少一半，但加法器数目不 变，这时乘法操作数和加法操作次数是 1： 2。事实上，线性相位 FIR 滤 波器是用得最多的FIR滤波器。 2.2用DSP实现的数字FIR滤波器的优点 在数字滤波器中， FIR 滤波器的最主要的特点是没有反馈回路，故 不存在不稳定的问题；同时，可以在幅度特性随意设置的同时，保证精 度、严格的线性相位。稳定和线性相位特性是 FIR 滤波器的突出优点。 FIR 滤波器常表示为直接型和转置型两种结构，其功能是等效的。 在用专用DSP或 FPGA等硬件设计长阶数FIR滤波器时，常采用转置型结 构，因为这种结构便于设计成大量规则排列、 形式相同的乘法 /累加器阵 列，避免了直接型中的多输入加法器在长阶数时的实现困难。 在用通用定点DSP实现FIR滤波器时，更多的是采用直接型结构，对于 C54X来说，其单乘法器、双累加器结构更适合于对称 FIR实现；如果采 用转置结构，因为定点DSP仅乘法/累加器是32位或40位的，数据存储 器都是16位的，用16位存储器保证FIR滤波器每级乘法/累加器输出的 中间结果比较麻烦，程序效率也降低了。 3. FIR 滤波器的设计方法 FIR 滤波器的设计方法主要有窗函数法和频率采样法，其中，窗函 数法是最基本的方法。窗函数法的过程是：设期望的滤波器理想响应为 Hd(ejw)，需要寻找一个传递函数 H (ej「=工N-1h(n)e -j ®去逼近H(ej w), 其中最直接的方法就是将 H(ejJ的时域响应hd( n)用一个矩形窗FN (n) 进行截断，从而得到一个长度为 N的序列h(n)，即令 h(n)= h d(n) R n (n) 理想响应Hd(ej°)与其时域响应hd(n)也有Hd(ej°) = 2^hd(n) e-j ®的关系， 通常hd(n)是一个关于原点n = 0对称的无限长序列，截断后的h (n)还 应进行右移以形成一个因果的 FIR滤波器，即在n0时，h (n)= 0。为 了保证h (n)的线性相位，还应使h (n)关于中心点对称或反对称。很 明显，h (n)与期望的响应之间的误差随着窗长的增加而减小。 利用上述各种窗函数，DSP设计者可以利用Matlab工具很方便地设 计出逼近理想特性的FIR滤波器，然后将此FIR系数放入DSP程序中。 如下例： 设计一个FIR低通滤波器，其设计参数：滤波器阶数为 40,截至频 率 3 p =