信号处理综合设计性实验

课 程 设 计

设计内容：

1. 已知一个LTI系统的差分方程为：

y(n)?1.143y(n?1)?0.4128y(n?2)?0.0675x(n)?0.1349x(n?1)?0.0675x(n?2)

(1)初始条件y(-1)=1,y(-2)=2，输入x(n)=u(n)，计算系统的零输入响应； (2)当下面三个信号分别通过系统，分别计算系统的响应：

??7?x1(n)?cos(n)u(n) x2(n)?cos(n)u(n) x3(n)?cos(n)u(n)

10510(3)指出这是一个什么特性的系统。

2. 已知一个因果LTI离散系统的系统函数为：

0.03571?0.1428z?1?0.2143z?2?0.1428z?3?0.0357z?4 H(z)??1?2?3?41?1.035z?0.8264z?0.2605z?0.04033z（1）求系统的单位脉冲响应h(n)

(2)当信号通过系统:x(n)??1?cos(n)?cos(n)?u(n),计算系统的响应。

42(3)详细列出系统响应求解的步骤和原理。

??????思考题：

(1) 系统函数零极点对系统频响的影响？ (2) 离散系统响应的求解方法有哪些？

(3) 在MATLAB中用于离散系统求解的命令有哪些？各基于什么求解方法？

- 37 -

课 程 设 计

第4章 数字滤波器设计

通信与电子信息当中，在对信号作分析与处理时，常会用到有用信号叠加无用噪声的问题。这些噪声信号有的是与信号同时产生的，有的是在传输过程中混入的，在接收的信号中，必须消除或减弱噪声干扰，这是信号处理中十分重要的问题。根据有用信号与噪声的不同特性，消除或减弱噪声，提取有用信号的过程就称为滤波。滤波器的种类很多，实现方法也多种多样，本章利用Matlab来进行数字滤波器的设计。

数字滤波器是一离散时间系统，它对输入序列x(n)进行加工处理后，输出序列y(n)，并使y(n)的频谱与x(n)的频谱相比发生某种变化。

由DSP理论得知，无限长冲激响应（IIR）需要递归模型来实现，有限长冲激响应（FIR）滤波器可以采用递归的方式也可采用非递归的方式实现。本章把FIR与IIR滤波器分别用Matlab进行分析与设计。数字滤波器的结构参看《数字信号处理》一书。

数字滤波器的设计一般经过三个步骤： 1．给出所需滤波器的技术指标。

2．设计一个H（Z），使其逼近所需要的技术指标。 3．实现所设计的H（Z）。

4.1 IIR数字滤波器设计

设计IIR数字滤波器的任务就是寻求一个因果、物理可实现的系统函数H(z)，使它的频响H(ejω)满足所希望得到的低通频域指标，即通带衰减Ap、阻带衰减Ar、通带截频ωp、阻带截频ωr。而其它形式的滤波器由低通的变化得到。

采用间接法设计IIR数字滤波器就是按给定的指标，先设计一个模拟滤波器，进而通过模拟域与数字域的变换，求得物理可实现的数字滤波器。从模拟滤波器变换到数字滤波器常用的有：脉冲响应不变法和双线性变换法。IIR滤波器的设计过程如下 数字

?模拟频域指标?设计模拟滤波器H(S) ?设计数字滤波器H(z)

- 38 -

课 程 设 计

1. 模拟滤波器简介

模拟滤波器的设计方法已经发展得十分成熟，常用的高性能模拟低通滤波器有巴特沃斯型、切比雪夫型和椭圆型，而高通、带通、带阻滤波器则可以通过对低通进行频率变换来求得。必须指出，这三种滤波器都是非线性的相频特性。具体设计方法参见《数字信号处理》一书。

（1）巴特沃斯低通滤波器的特点是：通、阻带均为单调下降。这种单调下降的特性使得系统的误差分布不均匀。在设计中，如果在通带满足指标，阻带指标就过于好了。或者说，在阻带满足指标，通带指标也过于好了。这样不利于以最小的阶数来满足设计指标。 （2）切比雪夫模拟低通滤波器 分成Ι型和ΙΙ型：Ι型是通带等波纹，阻带单调下降； ΙΙ型滤波器是通带单调下降，阻带等波纹。切比雪夫的等波纹特性使得可以用较小的阶数设计出符合指标的滤波器。

（3）椭圆滤波器特点是通、阻带均为等波纹。椭圆滤波器的误差均匀分布在通、阻带，比较上面几种滤波器，在同样衰耗指标下，椭圆滤波器所要求的阶数最小。在同样阶数情况下，椭圆滤波器的通带到阻带的变化最陡峭、性能最好。 2．模拟滤波器映射到数字滤波器

从模拟滤波器变换到数字滤波器常用的有：冲激响应不变法和双线性变换法。 （1）冲激响应不变法

冲激响应不变法的基本思路为：设模拟滤波器的单位冲激响应为ha（t），数字滤波器的单位脉冲响应h(n)，是对ha（t）的均匀取样，即

?? h(n)= h(nT) ???? H(z) H(S)?????h(t)??

冲激响应不变法由于保持了模拟滤波器的冲激响应在取样时刻的取样值，所以具有较好的时间特性。但由于从模拟域（s域）变换到数字域（z域）的映射关系不是一一对应的关系，因而频谱混叠难以避免, 所以在实际应用中有它的局限性。 3． Matlab滤波器设计函数简介

拉氏逆变换抽样z变换 - 39 -

课 程 设 计

使用MATLAB信号处理工具箱中提供的函数可以很容易地实现IIR滤波器设计。常用的IIR数字滤波器设计函数有： （1）IIR滤波器阶数（order）选择

buttord - 巴特沃斯（Butterworth）滤波器阶数选择 cheb1ord - 切比雪夫（Chebyshev）I型滤波器阶数选择 cheb2ord - 切比雪夫（Chebyshev）II型滤波器阶数选择 ellipord - 椭圆（Elliptic）滤波器阶数选择 buttord用来确定巴特沃斯滤波器的阶数：

[N,Wc] = buttord (Wp,Wr,Ap,Ar,options)

Ap,Ar为通带最大衰减和阻带最小衰减，以dB为单位； Wp,Wr为归一化通带截频和阻带截频。

Options=’low’, ’high’, ’bandpass’, ’stop’, ’s’, 分别对应低通、高通、带通、带阻、模拟滤波器，默认情况下，为低通或带通。对于带通和带阻滤波器，Wp=[Wp1,Wp2]，Wp1和 Wp2为通带的上、下截频，Wr=[Wr1,Wr2]，Wr1和Wr2分别为阻带的上、下截频。若为模拟滤波器定阶，令Options= ’s’, 则Wp和Wr以（归一或非归一）给出，因此可大于1。

输出参数N为滤波器的阶数，Wc为归一化3dB截频。若Ap=3dB，则Wc= Wp。 其他滤波器的定阶函数类似buttord的格式：

切比雪夫I型滤波器[N,Wc] = cheb1ord(Wp, Wr, Ap, Ar) 切比雪夫II型滤波器[N,Wc] = cheb2ord(Wp, Wr, Ap, Ar) 椭圆滤波器[N,Wc] = ellipord(Wp, Wr, Ap, Ar)

（2）完整的滤波器设计函数，根据滤波器阶数和通带截频设计出实际数字滤波器的H(z)

（模拟域到数字域的变换只采用双线性变换法）： butter - 巴特沃斯（Butterworth）滤波器设计 cheby1 -切比雪夫（Chebyshev）I型滤波器设计 cheby2 - -切比雪夫（Chebyshev）II型滤波器设计 ellip - 椭圆（Elliptic）滤波器设计

maxflat - 通用的巴特沃斯（Butterworth）低通滤波器设计

- 40 -