内容发布更新时间 : 2024/11/7 16:48:37星期一 下面是文章的全部内容请认真阅读。

2．分析直接计算线性卷积和利用FFT计算线性卷积的时间。

六、实验报告要求

1.简述实验原理及目的。

2.记录调试运行情况及所遇问题的解决方法。 3. 给出实验结果，并对结果作出分析。 4．简要回答思考题。

11

实验五 IIR数字滤波器的设计

一、实验目的

1．掌握双线性变换法及脉冲相应不变法设计IIR数字滤波器的具体设计方法及其原理，熟悉用双线性变换法及脉冲响应不变法设计低通、高通和带通IIR数字滤波器的计算机编程。

2．观察双线性变换及脉冲响应不变法设计的滤波器的频域特性，了解双线性变换法及脉冲响应不变法的特点。

3．熟悉Butterworth滤波器、Chebyshev滤波器和椭圆滤波器的频率特性。

二、实验原理

1．利用模拟滤波器设计IIR数字滤波器方法

（1）根据所给出的数字滤波器性能指标计算出相应的模拟滤波器的设计指标。 （2）根据得出的滤波器性能指标设计出相应的模拟滤波器的系统函数H(S)。

（3）根据得出的模拟滤波器的系统函数H(S)，经某种变换得到对该模拟滤波器相应的数字仿真系统——数字滤波器。

将模拟滤波器转换成数字滤波器的实质是，用一种从s平面到z平面的映射函数将Ha(s)转换成H(z)。对这种映射函数的要求是：(1) 因果稳定的模拟滤波器转换成数字滤波器，仍是因果稳定的。 (2)数字滤波器的频率响应模仿模拟滤波器的频响，s平面的虚轴映射z平面的单位圆，相应的频率之间成线性关系。脉冲响应不变法和双线性变换法都满足如上要求。

2．脉冲响应不变法

用数字滤波器的单位脉冲响应序列h(n)模仿模拟滤波器的冲激响应ha(t),让h(n)正好等于ha(t)的采样值，即h(n)=ha(nT)，其中T为采样间隔。

3．双线性变换法

s平面与z平面之间满足以下映射关系：

1?z?1s? ?11?zs平面的虚轴单值地映射于z平面的单位圆上，s平面的左半平面完全映射到z平面的单位圆内。双线性变换不存在混叠问题。

双线性变换时一种非线性变换??tg(?/2)，这种非线性引起的幅频特性畸变可通过预畸而得到校正。

以低通数字滤波器为例，将设计步骤归纳如下：

(1)确定数字滤波器的性能指标：通带临界频率fp、阻带临界频率fs；通带内的最大衰减Ap；阻带内的最小衰减As；

(2)确定相应的数字角频率，ωp=2πfp；ωs=2πfs； (3)计算经过预畸的相应模拟低通原型的频率，tg；

(4)根据Ωp和Ωs计算模拟低通原型滤波器的阶数N，并求得低通原型的传递函数Ha(s)； (5)用上面的双线性变换公式代入Ha(s)，求出所设计的传递函数H(z)； (6)分析滤波器特性，检查其指标是否满足要求。

三、主要实验仪器及材料

微型计算机、Matlab6.5教学版、TC编程环境。

12

四、实验内容

1．知识准备

在实验前复习数字信号处理理论课中有关滤波器设计的知识，认真阅读本实验的原理部分。

2．编制用脉冲响应不变法和用双线性变换法设计IIR数字滤波器的程序。采样周期、通带和阻带临界频率以及相应的衰减等参数在程序运行时输入；根据这些输入参数，计算阶数N、传递函数；输出分子分母系数；绘制幅频特性曲线，绘制点数为50点。

（1）利用脉冲响应不变法设计Butterworth、Chebyshev和椭圆数字低通滤波器，要求满足?p?0.2?，?s?0.6?，Ap?2dB，Aa?15dB。

（2）利用双线性变换法设计Butterworth、Chebyshev和椭圆数字滤波器，要求满足

?p?0.2?，?s?0.6?，Ap?2dB，Aa?15dB

（3）设计Butterworth、Chebyshev和椭圆高通数字滤波器，3dB数字截止频率为

?c?0.2?rad，阻带下边频?s?0.05?rad，阻带衰减As?48dB。

五、实验思考

1．双线性变换法中Ω和ω之间的关系是非线性的，在实验中你注意到这种非线性关系了吗？从那几种数字滤波器的幅频特性曲线中可以观察到这种非线性关系？

2．能否利用公式完成脉冲响应不变法的数字滤波器设计？为什么？

六、实验报告要求

1.简述实验原理及目的。

2.按照实验步骤及要求，比较各种情况下的滤波性能。 3.总结实验所得主要结论。 4.简要回答思考题。

13

实验六 FIR数字滤波器的设计

一、实验目的

1．熟悉FIR滤波器的设计基本方法

2．掌握用窗函数设计FIR数字滤波器的原理与方法，熟悉相应的计算机高级语言编程。 3．熟悉线性相位FIR滤波器的幅频特性和相位特性。 4．了解各种不同窗函数对滤波器性能的影响。

二、实验原理与方法

FIR滤波器的设计问题在于寻求一系统函数H(z)，使其频率响应H(ej?)逼近滤波器要求的理想频率响应Hd(ej?)，其对应的单位脉冲响应hd(n)。

1．用窗函数设计FIR滤波器的基本方法

设计思想：从时域从发，设计h(n)逼近理想hd(n)。设理想滤波器Hd(ej?)的单位脉冲响应为hd(n)。以低通线性相位FIR数字滤波器为例。

Hd(e)?hd(n)?12?j?n????h????d(n)ejn?

?Hd(ej?)ejn?d?

hd(n)一般是无限长的，且是非因果的，不能直接作为FIR滤波器的单位脉冲响应。要想得

到一个因果的有限长的滤波器h(n)，最直接的方法是截断h(n)?hd(n)w(n)，即截取为有限长因果序列，并用合适的窗函数进行加权作为FIR滤波器的单位脉冲响应。按照线性相位滤波器的要求，h(n)必须是偶对称的。对称中心必须等于滤波器的延时常数，即

?h(n)?hd(n)w(n) ?a?(N?1)/2?用矩形窗设计的FIR低通滤波器，所设计滤波器的幅度函数在通带和阻带都呈现出振荡

现象，且最大波纹大约为幅度的9%，这个现象称为吉布斯（Gibbs）效应。为了消除吉布斯效应，一般采用其他类型的窗函数。

2．典型的窗函数

（1）矩形窗(Rectangle Window)

w(n)?RN(n)

其频率响应和幅度响应分别为：

sin(N?/2)?j?j?W(e)?esin(?/2)N?12，WR(?)?sin(N?/2)

sin(?/2) 14

（2）三角形窗(Bartlett Window)

N?1?2n，0?n??2 w(n)??N?12nN?1?2?，?n?NN?12?2N?122?sin(N?/4)??j?j?e其频率响应为：W(e)???N?sin(?/2)?（3）汉宁(Hanning)窗，又称升余弦窗

12n?w(n)?[1?cos()]RN(n)

2N?1?j2?2?W(e)?{0.5WR(?)?0.25[WR(??)?WR(??)]}eN?1N?1j?N?1?2其频率响应和幅度响应分别为：

W(?)?0.5WR(?)?0.25[WR(??2?2?)?WR(??)] N?1N?1（4）汉明(Hamming)窗，又称改进的升余弦窗

w(n)?[0.54?0.46cos(2n?)]RN(n) N?12?2?)?WR(??)]N?1N?1

其幅度响应为：

W(?)?0.54WR(?)?0.23[WR(??（5）布莱克曼(Blankman)窗，又称二阶升余弦窗

w(n)?[0.42?0.5cos(其幅度响应为：

2n?4n?)?0.08cos()]RN(n) N?1N?1W(?)?0.42WR(?)?0.25[WR(??2?2?)?WR(??)] N?1N?14?4??0.04[WR(??)?WR(??)]N?1N?1

（6）凯泽(Kaiser)窗

I0(?1?[1?2n/(N?1)]2)w(n)?，0?n?N?1I0(?)

其中：β是一个可选参数，用来选择主瓣宽度和旁瓣衰减之间的交换关系，一般说来，β越大，过渡带越宽，阻带越小衰减也越大。I0(·)是第一类修正零阶贝塞尔函数。 若阻带最小衰减表示为As??20log10?s，β的确定可采用下述经验公式：

As?21?0?(As?21)0.4?0.07886(As?21)21?As?50 ?0.5842?0.1102(As?8.7)As?50?若滤波器通带和阻带波纹相等即δp=δs时，滤波器节数可通过下式确定：

15