数字滤波器的设计 MATLAB.docx

PAGE 2 PAGE 2 实验报告 题目： 数字滤波器的设计 课程名称： 成 绩： 目录 TOC \o 1-3 \h \z \u 32377 1 实验目的 2 270 2 实验环境 2 28323 3 实验原理 2 22520_WPSOffice_Level2 3.1 IIR数字滤波器设计 2 7756_WPSOffice_Level2 3.2 FIR滤波器的设计 6 7756_WPSOffice_Level1 4 实验内容 11 20532_WPSOffice_Level2 4.1实验步骤 11 23357_WPSOffice_Level2 4.2实验结果及分析 13 20532_WPSOffice_Level1 5.实验思考题 26 PAGE PAGE 25 1 实验目的 1、了解两种工程上最常用的变换方法：脉冲响应不变法和双线性变换法。 2、掌握双线性变换法设计 IIR 滤波器和熟悉用窗函数设计 FIR 数字滤波器的原理及方法。熟悉用双线性设计法设计低通、带通和高通IIR 数字滤波器的计算机程序。 3、熟悉用双线性变换法设计数字 Butterworth 和 Chebyshev 滤波器的全过程。熟悉线性相位 FIR 滤波器的幅频特性和相位特性。 4、观察用双线性变换法设计的滤波器的频域特性，并与脉冲响应不变法相比较，了解双线性变换法的特点。了解多项式乘积和多项式乘方运算的计算机编程方法，了解各种不同窗函数对滤波器性能的响应。 2 实验环境 1.一台计算机。含Matlab软件。 3 实验原理 3.1 IIR数字滤波器设计 从模拟滤波器设计IIR数字滤波器具有四种方法:微分差分变换法、脉冲响应不变法、双线性变换法、匹配z变换法;在工程上常用的是其中两种:脉冲响应不变法、双线性变换法。脉冲响应不变法需要经历如下基本步骤：由已知系统传输函数H(s)计算系统冲激响应；对h(t)进行等间隔采样得到h(n)=h(nT);由h(n)获得数字滤波器的系统响应H(z)。这种方法非常直观，其算法宗旨是保证所设计的IIR滤波器的脉冲响应和响应模拟滤波器的冲激响应在采样点上完全一致。而双线性变换法的设计准则是使数字滤波器的频率响应与参考模拟滤波器的频率响应相似。脉冲响应不变法一个重要的特点是频率坐标的变换是线性的(ω=ΩT)，其缺点是有频谱的周期延拓效应，存在频谱混淆的现象。为了克服脉冲响应不变法可能产生的频谱混淆，提出了双线性变换法，它依靠双线性变换式:?s=1+Z-1 Ω=tg(ω/2) ω=2arctgΩ (3-1) 由上面的关系式可知，当0→∞时，∞终止在折叠频率∞=π处，整个jΩ轴单值地对应于单位圆的一周。因此双线性变换法不同于脉冲响应不变法，不存在频谱混淆的问题。从式(3-1)还可以看出，两者的频率不是线性关系。这种非线性关系使得通带截止频率、过渡带的边缘频率的相对位置都发生了非线性畸变。这种频率的畸变可以通过预畸来教正。用双线性变换法设计数字滤波器时，一般总是先将数字滤波器的各临界频率经过式(3-1)的频率预畸，求得相应参考模拟滤波器的各临界频率，然后设计参考模拟滤波器的传递参数，最后通过双线性变换式求得数字滤波器的传递函数。这样通过双线性变换，正好将这些频率点映射到我们所需要的位置上。在用双线性变换法设计滤波器的过程中，我们也可以通过原型变换，直接求得归一化参考模拟滤波器原型参数，从而使得设计更加简化。下表是IIR低通、带通、高通滤波器设计双线性原型变换公式的总结。 表（3-1） 变换类型 变换关系式 参考模拟模型 备注 低通变换 s= Ω=tg(ω/2) ω=2πfT f:模拟频率 T：采样周期 cos ω1 高通变换 s= Ω=|ctg( 带通变换 s= Ω=| 在本实验中，我们只涉及Butterworth和Chebyshev两种滤波器的设计，相应的这两种参考模拟原形滤波器的设计公式见下表： 表（3-2） 类型 巴特沃斯（Butterworth） 切比雪夫（Chebyshev） 阶数 N=lg N≥ 传递 函数 H s H xk=- 有关 参数 3dB频率 Ωc=Ωs/ Ωc=Ωp/ 波纹系数 ε=10 备注 通带临界频率Ωp及衰减αp，阻带临界频率Ωs及衰减αs 综上所述，以巴特沃斯低通数字滤波器设计为例，可以将双线性变换法设计数字滤波器的步骤归纳如下: 1、确定数字滤波器的性能指标。这些指标包括:通带、阻带临界频率fp、fs；通带内的最大衰减αp 2、确定相应的