文档介绍：实验六离散线性时不变系统分析
一、实验目的
1. 掌握离散LSI系统的单位序列响应、单位阶跃响应和任意激励下响应的MATLAB求解方法。
2. 掌握离散LSI系统的频域分析方法;
3. 掌握离散LSI系统的复频域分析方法;
4. 掌握离散LSI系统的零极点分布与系统特性的关系。
二、实验原理及方法
离散LSI系统的时域分析
描述一个N阶线性时不变离散时间系统的数学模型是线性常系统差分方程,N阶LSI离散系统的差分方程一般形式为
()
也可用系统函数来表示
()
系统函数反映了系统响应和激励间的关系。一旦上式中,的数据确定了,系统的性质也就确定了。特别注意必须进行归一化处理,即。
对于复杂信号激励下的线性系统,可以将激励信号在时域中分解为单位序列或单位阶跃序列的线性叠加,把这些单元激励信号分别加于系统求其响应,然后把这些响应叠加,即可得到复杂信号作用于系统的零状态响应。因此,求解系统的单位序列响应和单位阶跃响应尤为重要。由图6-1可以看出一个离散LSI系统响应与激励的关系。
图6-1 离散LSI系统响应与激励的关系
(1) 单位序列响应(单位响应)
单位响应是指离散LSI系统在单位序列激励下的零状态响应,因此满足线性常系数差分方程()及零初始状态,即
, ()
按照定义,它也可表示为
()
对于离散LSI系统,若其输入信号为,单位响应为,则其零状态响应为
()
可见,能够刻画和表征系统的固有特性,与何种激励无关。一旦知道了系统的单位响应,就可求得系统对任何输入信号所产生的零状态响应。
MATLAB提供了专门用于求离散时间系统冲激响应的函数impz(),其调用格式有
[h,n]=impz(b,a)
求解离散系统的单位脉冲响应,其中,
,;
[h,n]=impz(b,a,N)
求解离散系统的单位响应,采样点数由N确定,;
impz(b,a) :在当前绘图窗口,用stem(n,h)绘出图形。
(2)单位阶跃响应
单位阶跃响应是指离散LTI系统在单位阶跃序列激励下的零状态响应,它可以表示为
()
上式表明,离散LSI系统的单位阶跃响应是单位脉冲响应的累加和,系统的单位阶跃响应和系统的单位响应之间有着确定的关系,因此,单位阶跃响应也能完全刻画和表征一个LSI系统。
MATLAB提供了专门用于求离散系统单位阶跃响应的函数stepz( ),其调用格式有
[s,n]=stepz(b,a) :求解离散系统的单位阶跃响应,其中,,;
[s,n]=stepz(b,a,N) :求解离散系统的单位阶跃响应,采样点数由N确定,;
stepz(b,a) :在当前窗口,用stem(n,s)绘出图形。
任意激励下的零状态响应
已经知道,离散LSI系统可用常系数线性差分方程()式来描述,Matlab提供的函数dlsim( )能对上述差分方程描述的离散LSI系统的响应进行仿真,该函数不仅能绘制指定时间范围内的系统响应波形图,而且还能求出系统响应的数值解。其调用格式有
dlsim(b,a, x) :求解输入序列为x的零状态响应
需要特别强调的是,Matlab总是把由分子和分母多项式表示的任何系统都当作是因果系统。所以,利用impz (b,a),stepz(b,a),dlsim(b,a,x)函数求得的响应总是因果信号。
同时,卷积和也是LSI系统求解零状态响应的重要工具之一。假设系统的输入信号为,单位响应为,则系统的零状态响应可由()式求解。Matlab提供了专门用于求离散系统卷积和的函数conv( ),其调用格式有
y=conv(x,h) :求解序列x,h的卷积和,若序列x的长度为n1,序列h的长度为n2,卷积和y的长度为n1+n2-1。这一点需要特别注意,否则,作图时容易造成横纵坐标长度不匹配。
例如已知输入序列,系统的单位响应,求系统的零状态响应。参考程序如下:
n1=0:19;
x=.^n1;
subplot(231)
stem(n1,x,'.');
title('x(n)');
n2=0:9;
h=ones(1,10);
subplot(232)
stem(n2,h,'.');
title('h(n)');
y=conv(x,h); %求卷积和
m1=n1(1)+n2(1); %卷积和的起始位置
m2=n1(length(n1))+n2(length(n2)); %卷积和的结束位置
n3=m1:m2; % 卷积和的长度,这是需要大家特别注意的地方。
subplot(233)
stem(n3,y,'.');
title('y(n)');
运行结果如图6-2所示。
图