离散线性时变系统分析.doc

《离散线性时变系统分析.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、实验六离散线性时不变系统分析一、实验目的1.掌握离散LSI系统的单位序列响应、单位阶跃响应和任意激励下响应的MATLAB求解方法。2.掌握离散LSI系统的频域分析方法；3.掌握离散LSI系统的复频域分析方法；4.掌握离散LSI系统的零极点分布与系统特性的关系。二、实验原理及方法1.离散LSI系统的时域分析描述一个N阶线性时不变离散时间系统的数学模型是线性常系统差分方程，N阶LSI离散系统的差分方程一般形式为（6.1）也可用系统函数来表示（6.2）系统函数反映了系统响应和激励间的关系。一旦上式中，的数据确定了，系统的性质也就确定了。特别注意必须进行归一化处理，即。对于复杂信号激

2、励下的线性系统，可以将激励信号在时域中分解为单位序列或单位阶跃序列的线性叠加，把这些单元激励信号分别加于系统求其响应，然后把这些响应叠加，即可得到复杂信号作用于系统的零状态响应。因此，求解系统的单位序列响应和单位阶跃响应尤为重要。由图6-1可以看出一个离散LSI系统响应与激励的关系。图6-1离散LSI系统响应与激励的关系(1)单位序列响应（单位响应）单位响应是指离散LSI系统在单位序列激励下的零状态响应，因此满足线性常系数差分方程(6.1)及零初始状态，即，(6.3)按照定义，它也可表示为(6.4)对于离散LSI系统，若其输入信号为，单位响应为，则其零状态响应为(6.5)可见

3、，能够刻画和表征系统的固有特性，与何种激励无关。一旦知道了系统的单位响应，就可求得系统对任何输入信号所产生的零状态响应。MATLAB提供了专门用于求连续系统冲激响应的函数impz()，其调用格式有[h,n]=impz(b,a)求解离散系统的单位响应，其中，，；[h,n]=impz(b,a,N)求解离散系统的单位响应，采样点数由N确定，；impz(b,a)：在当前窗口，用stem(n,h)绘出图形。（2）单位阶跃响应单位阶跃响应是指离散LTI系统在单位阶跃序列激励下的零状态响应，它可以表示为(6.6)上式表明，离散LSI系统的单位阶跃响应是单位响应的累加和，系统的单位阶跃响应和

4、系统的单位响应之间有着确定的关系，因此，单位阶跃响应也能完全刻画和表征一个LSI系统。MATLAB提供了专门用于求离散系统单位阶跃响应的函数stepz()，其调用格式有[s,n]=stepz(b,a)：求解离散系统的单位阶跃响应，其中，，；[s,n]=stepz(b,a,N)：求解离散系统的单位阶跃响应，采样点数由N确定，；stepz(b,a)：在当前窗口，用stem(n,s)绘出图形。(3)任意激励下的零状态响应已经知道，离散LSI系统可用常系数线性差分方程(6.1)式来描述，Matlab提供的函数dlsim()能对上述差分方程描述的离散LSI系统的响应进行仿真，该函数不仅

5、能绘制指定时间范围内的系统响应波形图，而且还能求出系统响应的数值解。其调用格式有dlsim(b,a,x)：求解输入序列为x的零状态响应需要特别强调的是，Matlab总是把由分子和分母多项式表示的任何系统都当作是因果系统。所以，利用impz(b,a)，stepz(b,a)，dlsim(b,a,x)函数求得的响应总是因果信号。同时，卷积和也是LSI系统求解零状态响应的重要工具之一。假设系统的输入信号为，单位响应为，则系统的零状态响应可由（6.5）式求解。Matlab提供了专门用于求离散系统卷积和的函数conv()，其调用格式有y=conv(x,h)：求解序列x，h的卷积和，若序列

6、x的长度为n1，序列h的长度为n2，卷积和y的长度为n1+n2-1。这一点需要特别注意，否则，作图时容易造成横纵坐标长度不匹配。（4）带初始状态的任意激励下的全响应任意激励下的离散LSI系统的全响应为零输入响应和零状态响应之和，表示为（6.7）在理论学习的过程中，同学们对低阶差分方程的求解已颇为头痛，高阶差分方程直接求解几乎不可能。Matlab提供了用于求离散系统全响应的函数filter()，其调用格式有y=filter(b,a,x)：求解零状态响应；y=filter(b,a,x,zi)：求解初始条件为zi的系统的全响应，zi向量的长度为max(length(a),lengt

7、h(b))-1，返回值为系统的全响应。z=filtic(b,a,y,x)：将初始状态转换为初始条件，其中，；z=filtic(b,a,y)：将初始状态转换为初始条件，其中，。2离散LSI系统的复频域（Z域）分析（1）利用Z变换解差分方程在前面图6-1中表示了离散系统的响应与激励的关系，由图可知，系统的响应既可以用时域的方法求解，也可以用Z域的方法求解。当已知系统输入序列的Z变换，系统函数时，系统响应序列的Z变换可由求出。Matlab提供了用于求序列Z变换和Z反变换的函数，其调用格式有X=ztrans(x