FIR滤波器程序的编写.docVIP

PAGE 13 目录 TOC \o 1-3 \h \z \u 1.引言 1 2.FIR滤波器基本原理 1 2.1 两种FIR滤波器的框图及原理 1 2.2 用DSP实现的数字FIR滤波器的优点 2 3.FIR滤波器的设计方法 3 4.FIR滤波器的c54x实现 4 5. FIR滤波器的实例设计 7 5.1 编写滤波器序列 7 5.2 汇编源程序主程序的编写 8 5.3 中断向量表的编写 10 5.4 相关的连接命令文件 12 5.5 程序在CCS C5000中编译 12 6. 设计心得 13 7. 参考文献 14  1.引言 在数字信号处理中,滤波占有及其重要的地位。数字滤波是谱分析、雷达信号处理、通信信号处理等应用中的基本处理算法，它能够满足滤波器对幅度和相位的严格要求，解决了模拟滤波器所无法克服的电压漂移、温度漂移和噪声等问题，同时有很高的可编程性和灵活性。数字滤波是DSP最基本的应用领域，一个DSP芯片执行数字滤波算法的能力反应了这种芯片的功能强弱，用DSP实现的数字滤波器可以十分方便地改变滤波器的特性。 2.FIR滤波器基本原理 2.1 两种FIR滤波器的框图及原理 DDDh D D D h0 h1 hN-2- hN-1 多输入加法器 输出y（n） 输入x(n ) h0 h0 h1 hN-2- hN-1 D D D 输出y（n） 输入x(n ) （b）转置型横向滤波器 图2.1 直接型和转置型FIR滤波器结构图 FIR数字滤波器的差分方程表达式为： y(n)=∑Nh(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-1h(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(ejω)，需要寻找一个传递函数H (ejω)＝ΣN-1h(n)e-jω去逼近Hd(ejω)，其中最直接的方法就是将Hd(ejω)的时域响应hd(n)用一个矩形窗RN（n）进行截断，从而得到一个长度为N的序列h(n)，即令 h(n)= hd(n) RN（n） 理想响应Hd(ejω)与其时域响应hd(n)也有Hd(ejω)＝Σ∞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低通