- 1、本文档共4页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
信号与系统实验连续LTI系统的时域分析
实验二 连续LTI系统的时域分析
一. 实验目的
1. 加深对线性时不变系统中零状态响应概念的理解,掌握其求解方法;
2. 掌握求解给定连续系统的冲激响应和阶跃响应的方法。
二. 实验原理
1.连续系统零状态响应的数值解
线性时不变 (LTI) 连续时间系统用常系数线性微分方程进行描述,系统的零状态响应就是在系统初始状态为零条件下微分方程的解。MATLAB控制系统工具箱提供了一个lism函数来求解连续时间系统的零状态响应,其调用格式为
y = lism(sys,f,t)
其中t表示计算系统响应的时间抽样点向量,f是系统输入信号向量,sys是LTI系统模型,用来表示微分方程、状态方程。在求解微分方程时,微分方程的LTI系统模型sys要借助MATLAB中的tf函数来获得,其调用格式为
sys = tf(b,a)
其中a、b分别为微分方程左端和右端各项的系统向量。
例如系统方程
该方程左边、右边的系数向量分别为,。
例1:描述某线性时不变系统的方程为
试求:当时,系统的零状态响应。
解:实现所要求运算的m文件如下,
a = [1 4 4]; %将y(t)各阶导数的系数放在向量a中
b = [1 3]; %将f(t)各阶导数的系数放在向量b中
sys = tf(b, a); %求系统模型sys
td = 0.01; %定义时间间隔
t = 0 : td : 10; %定义时间向量
f = exp(-t); %将f(t)表示出来
y = lsim(sys, f, t); %求系统的零状态响应
plot(t, y); %绘出零状态响应的波形
xlabel(t(sec)); %给出x坐标的标签
ylabel(y(t)); %给出y坐标的标签
grid on %在图上显示方格
程序运行结果见图1。
图1. 零状态响应波形
2.连续系统的冲激响应和阶跃响应
连续系统的冲激响应、阶跃响应分别是输入信号为和所对应的零状态响应。MATLAB控制系统工具箱专门提供了两个函数求解连续系统的冲激响应和阶跃响应。
冲激响应:y = impulse(sys, t) ;
阶跃响应:y = step(sys, t);
其中t表示计算系统响应的时间抽样点向量,sys是LTI系统模型。
例2:描述某线性时不变系统的方程为
试求系统的冲激响应h(t)和阶跃响应g(t)。
解:实现所要求运算的m文件如下,
b=[20 100];? %输入微分方程右边的系数行向量;
a=[1 2 100]; %输入微分方程左边的系统行向量;
sys=tf(b,a);
t=0:0.02:4; %定义时间向量
figure(1)
impulse(sys,t); %计算冲激响应并绘制波形
figure(2)
step(sys,t); %计算阶跃响应并绘制波形
运算结果如图2(a)、(b)所示。
图2(a) 冲激响应波形
图2(b) 阶跃响应波形
三. 实验内容
1.已知系统的微分方程为
若,计算系统的零状态响应、冲激响应和阶跃响应,并画出波形图。
2.已知系统的微分方程为
分别求系统的零状态响应和阶跃响应,并在一个窗口中画出两个图形,观察结果。如果输入信号改为,那么结果又将如何?
四. 实验报告
1. 报告组成部分:实验目的、实验原理、实验内容、实验结果及分析。
2. 在实验结果中列出全部m文件。
3. 简述心得体会。
文档评论(0)