(基于IIR的滤波器)选编.doc

《DSP技术与应用》课程设计 题目：基于DSP的IIR滤波器的设计 学院 电子与信息工程学院 专业 电子信息工程 班级 电工13本1 姓名 李 杭 学号 130916013 指导老师 刘明华老师 时间 2016年06月 评分 1.数字滤波器的优越性 DSP 芯片是一种特别适合数字信号处理运算的微处理器，主要用来实时、快速地实现各种数字信号处理算法。数字信号处理由于具有精度高、灵活性强等优点，已广泛应用于图像处理、数字通信、雷达等领域。对数字滤波器而言,从实现方法上有FIR滤波器和无限冲激响应(IIR)滤波器之分。由于FIR滤波器只有零点因此这一类系统不像IIR系统那样易取得比较好的通带与阻带衰减特性。但是IR系统 2、高效率（因为结构简单、系数小、乘法操作较少） 3、与模拟滤波器有对应关系 4、可以解析控制，强制系统在指定位置为零点 数字滤波器的实现方法 目前数字滤波器的主要实现方法有： 1．在通用的微型计算机上用软件实现。 2．用 DSP(Digital Signal Processing)处理器实现 3．用固定功能的专用信号处理器实现 4．用FPGA等可编程器件来开发数字滤波算法。 本课题选择在TM320C5509 DSP综合实验开发系统平台上进行开发。最后采用窗函数法在CCS环境下，利用MATLAB仿真得到的滤波器系数，编程实现IIR数字滤波器，得出滤波结果波形，并对结果进行分析与总结。 数字滤波器的介绍 数字滤波器是对数字信号实现滤波的线性时不变系统。数字滤波实质上是一种运算过程，实现对信号的运算处理。输入数字信号（数字序列）通过特定的运算转变为输出的数字序列，因此，数字滤波器本质上是一个完成特定运算的数字计算过程，也可以理解为是一台计算机。描述离散系统输出与输入关系的卷积和差分方程只是给数字信号滤波器提供运算规则，使其按照这个规则完成对输入数据的处理。时域离散系统的频域特性: , 数字滤波器根据其冲激响应函数的时域特性，可分为两种，即无限长冲激响应(IIR)数字滤波器和有限长冲激响应(FIR)数字滤波器。IIR 数字滤波器的特征是，具有无限持续时间冲激响应，需要用递归模型来实现，其差分方程为： 系统函数为： 设计IIR滤波器的任务就是寻求一个物理上可实现的系统函数H(z)，使其频率响应H(z)满足所希望得到的频域指标，即符合给定的通带截止频率、阻带截止频率、通带衰减系数和阻带衰减系数。 IIR滤波器的设计方法及原理 IIR数字滤波器是一种离散时间系统，其系统函数为： 2 IIR滤波器的设计过程 根据以上IIR数字滤波器设计方法，下面运用双线性变换法基于MATLAB设计一个IIR带通滤波器，其中带通的中心频率为ωp0=0.5π，通带截止频率ωp1=0.4π，ωp2=0.6π;通带最大衰减αp=3dB;阻带最小衰减αs=15dB;阻带截止频率ωs2=0.7π 设计步骤： (1)根据任务,确定性能指标:在设计带通滤波器之前,首先根据工程实际的需要确定滤波器的技术指标: 带通滤波器的阻带边界频率关于中心频率ωp0几何对称，因此ws1=wp0- (ws2-wp0)=0.3π 通带截止频率wc1=0.4π,wc2=0.6π;阻带截止频率wr1=0.3π，wr2=0.7π；阻带最小衰减αs=3dB和通带最大衰减αp=15dB; Ω=2/T*tan(w/2)对带通数字滤波器H(z)的数字边界频率预畸变，得到带通模拟滤波器H(s)的边界频率主要是通带截止频率ωp1,ωp2;阻带截止频率ωs1，ωs2的转换。 为了计算简便，对双线性变换法一般T=2s 通带截止频率wc1=(2/T)*tan(wp1/2)=tan(0.4π/2)=0.7265 wc2=(2/T)*tan(wp2/2)=tan(0.6π/2)=1.3764 阻带截止频率wr1=(2/T)*tan(ws1/2)=tan(0.3π/2)=0.5095 wr2=(2/T)*tan(ws2/2)=tan(0.7π/2)=1.9626 s=3dB和通带最大衰减αp=15dB; (3)运用低通到带通频率变换公式λ=(((Ω^2)-(Ω0^2))/(B*Ω))将模拟带通滤波器指标转换为模拟低通滤波器指标。 B=wc2-wc1=0.6499 normwr1=(((wr1^2)-(w0^2))/(B*wr1))=2.236 normwr2=(((wr2^2)-(w0^2))/(B*wr2))=2.236 normwc1=(((wc1^2)-(w0^2))/(B*wc1))=1 normwc2=(((wc2^2)-(w0^2))/(B*wc2))=1 得出，normwc=1，normwr=2.236 模拟低通滤波器指标：normwc=1，normwr=2.