信号与系统 分析实验报告
实验项目名称：离散线性时不变系统分析；
连续时间系统分析
所属课程名称：信号与系统实验教程 实 验 类 型 ：验证型 指 导 教 师 ：
实 验 日 期 ：2013.06.04 班 级 ： 学 号 ： 姓 名 ：
离散线性时不变系统分析
一、实验目的
1. 掌握离散线性时不变系统的单位序列响应、单位阶跃响应和任意激励下响应的MATLAB求解方法。
2. 掌握离散线性时不变系统的频域分析方法； 3. 掌握离散线性时不变系统的复频域分析方法；
4. 掌握离散线性时不变系统的零极点分布与系统特性的关系。
二、实验原理及方法
1.
离散线性时不变系统的时域分析
描述一个N阶线性时不变离散时间系统的数学模型是线性常系统差分方程，N阶线性时不变离散系统的差分方程一般形式为
（2.1） NMk?0?aky(n?k)??bix(n?i)i?0也可用系统函数来表示
Y(z)H(z)??X(z)?bziM?i （2.2）
?azkk?0i?0N?kb(z)b0?b1z?1?b2z?2???bMz?M??a(z)1?a1z?1?a2z?2???aNz?N系统函数H(z)反映了系统响应和激励间的关系。一旦上式中ak，bi的数据确定了，系统的性质也就确定了。特别注意a0必须进行归一化处理，即a0?1。
对于复杂信号激励下的线性系统，可以将激励信号在时域中分解为单位序列或单位阶跃序列的线性叠加，把这些单元激励信号分别加于系统求其响应，然后把这些响应叠加，即可得到复杂信号作用于系统的零状态响应。因此，求解系统的单位序列响应和单位阶跃响应尤为重要。由图2-1可以看出一个离散LSI系统响应与激励的关系。
x(n)h(n)y(n)?x(n)*h(n)
X(z)H(z)Y(z)?X(z)H(z)图2-1 离散LSI系统响应与激励的关系
(1) 单位序列响应（单位响应）
单位响应h(n)是指离散线性时不变系统在单位序列?(n)激励下的零状态响应，因此
h(n)满足线性常系数差分方程(2.1)及零初始状态，即
?ah(n?k)??b?(n?i)kik?0i?0NM， h(?1)?h(?2)???0 (2.3)
按照定义，它也可表示为
h(n)?h(n)??(n) (2.4) 对于离散线性时不变系统，若其输入信号为x(n)，单位响应为h(n)，则其零状态响应
yzs(n)为
yzs(n)?x(n)*h(n) (2.5)
可见，h(n)能够刻画和表征系统的固有特性，与何种激励无关。一旦知道了系统的单位响应h(n)，就可求得系统对任何输入信号x(n)所产生的零状态响应y(n)。
zsMATLAB提供了专门用于求离散系统冲激响应的函数impz()，其调用格式有 [h,n]=impz(b,a)
求解离散系统的单位响应，其中b?[b,b,b?,Mb，]a?[1,a1,a2,?,aN]，012,n?[0,1,2,?]?；
[h,n]=impz(b,a,N)
求解离散系统的单位响应，采样点数由N确定，n?[0,1,2,?,N-1]?；
impz(b,a) ：在当前窗口，用stem(n,h)绘出图形。 （2）单位阶跃响应
单位阶跃响应s(n)是指离散离散时不变系统在单位阶跃序列u(n)激励下的零状态响应，它可以表示为
ns(n)?u(n)?h(n)?m????h(m) (2.6)
上式表明，离散线性时不变系统的单位阶跃响应是单位响应的累加和，系统的单位阶跃响应和系统的单位响应之间有着确定的关系，因此，单位阶跃响应也能完全刻画和表征一个线性时不变系统。
MATLAB提供了专门用于求离散系统单位阶跃响应的函数stepz( )，其调用格式有 [s,n]=stepz(b,a) ：求解离散系统的单位阶跃响应，其中，b?[b,b,b?,,b]012Ma?[1,a1,a2,?,aN]，n?[0,1,2,?]?；
[s,n]=stepz(b,a,N) ：求解离散系统的单位阶跃响应，采样点数由N确定，n?[0,1,2,?,N-1]?；
stepz(b,a) ：在当前窗口，用stem(n,s)绘出图形。 （3）任意激励下的零状态响应
已经知道，离散线性时不变系统可用常系数线性差分方程(2.1)式来描述，Matlab提供的函数dlsim( )能对上述差分方程描述的离散线性时不变系统的响应进行仿真，该函数不仅能绘制指定时间范围内的系统响应波形图，而且还能求出系统响应的数值解。其调用格式有
dlsim(b,a, x) ：求解输入序列为x的零状态响应
需要特别强调的是，Matlab总是把由分子和分母多项式表示的任何系统都当作是因果系统。所以，利用impz (b,a)，stepz(b,a)，dlsim(b,a,x)函数求得的响应总是因果信号。
同时，卷积和也是线性时不变系统求解零状态响应的重要工具之一。假设系统的输入信号为x(n)，单位响应为h(n)，则系统的零状态响应y(n)可由（2.5）式求解。Matlab提
zs供了专门用于求离散系统卷积和的函数conv( )，其调用格式有
y=conv(x,h) ：求解序列x，h的卷积和，若序列x的长度为n1，序列h的长度为n2，卷积和y的长度为n1+n2-1。这一点需要特别注意，否则，作图时容易造成横纵坐标长度不匹配。
（4）带初始状态的任意激励下的全响应
任意激励下的离散线性时不变系统的全响应为零输入响应和零状态响应之和，表示为
y(n)?yzi(n)?yzs(n) （2.7）
Matlab提供了用于求离散系统全响应的函数filter( )，其调用格式有 y=filter( b,a,x) ：求解零状态响应；
y=filter( b,a,x,zi) ：求解初始条件为zi的系统的全响应，zi向量的长度为max(length(a),length(b))-1，返回值为系统的全响应。
z = filtic(b,a,y,x)：将初始状态转换为初始条件，其中x?[x(?1),x(?2),x(?3),?,x(?m)]，
y?[y(?1),y(?2),y(?3),?,y(?n)]；
z = filtic(b,a,):将初始状态转换为初始条件x?0,y?[y(?1),y(?2),y(?3),?,y(?n)] 2 离散线性时不变系统的复频域（Z域）分析
（1）利用Z变换解差分方程
在前面图2-1中表示了离散系统的响应与激励的关系，由图可知，系统的响应既可以用时域的方法求解，也可以用Z域的方法求解。当已知系统输入序列的Z变换X(z)，系统函数H(z)时，系统响应序列的Z变换可由Y(z)?X(z)H(z)求出。Matlab提供了用于求序列Z变换和Z反变换的函数，其调用格式有
X=ztrans(x)：求无限长序列x的Z变换，返回Z变换的表达式，注意这里x，X都是符号表达式；
x=iztrans(X)：求X（z）的Z反变换x(n)，返回Z反变换的表达式，注意这里x，X都是符号表达式；
[r,p,c]=residuez(b,a)：把b(z)/a(z)展开成部分分式；
[b,a]=residuez(r,p,c):根据部分分式的r、p、c数组，返回有理多项式。 （2）系统的零极点分布与系统因果性和稳定性的关系
因果系统的单位响应h(n)一定满足当n?0时，h(n)=0，那么其系统函数H(z)的收敛域一定包含?点，即?点不是极点，极点分布在某个圆的圆内，收敛域是圆外区域。
，对照z变换定义，系统稳定要求收敛域包含单位圆。系统稳定要求 ?n????|h(n)|??如果系统因果且稳定，收敛域包含?点和单位圆，那么收敛域可表示为：
r?|z|??,0?r?1 (2.8)
MATLAB提供了用于求系统零极点的函数，其调用格式有 roots()：利用多项式求根函数来确定系统函数的零极点位置；
roots(a)：求极点位置，a为系统函数H(z)分母多项式所构成的系数向量； roots(b)：求零点位置，b为系统函数H(z)分子多项式所构成的系数向量；
zplane(b,a)：绘制由行向量b和a构成的系统函数的零极点分布图；
zplane(z,p)：绘制由列向量z确定的零点、列向量p确定的极点构成的零极点分布图。
（3）系统的零极点分布与系统频率响应的关系 将式(2.2)因式分解，得到
(2.9) M?1H(z)?AmN?(1?dkz?1)k?1?(1?cmz?1)式中，A?ba，c是H(z)的零点，d是其极点。A参数影响频率响应的幅度大小，
mk00影响系统特性的是零点c和极点d的分布。下面采用几何方法研究系统零极点分布对系统mk频率特性的影响。
将式(2.9)的分子、分母同乘以zN?M，得到：