- 1、本文档共11页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
MATLAB软件对LTI连续系统时域进行分析仿真
目录
前言 1
项目概况 1
正文 2
3.1设计目的与意义 2
3.1.1目的 2
3.1.2意义 2
3.2目标与总体方案 2
3.2.1目标 2
3.2.2总体方案 2
3.3方法与内容 2
3.3.1设计原理 2
3.3.2 内容 2
3.4 结论 10
致谢 10
参考文献 10
前言
MATLAB是由美国athworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中,为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案,并在很大程度上摆脱了传统非交互式程序设计语言(如C、Fortran)的编辑模式,代表了当今国际科学计算软件的先进水平。《信号与线性系统分析》是电子信息类专业本科生的一门必修专业基础课。通过本门课的学习,使学生能够熟练掌握信号的表示、运算,了解系统的性质及表示;掌握连续系统的时域分析、连续系统的时域响应的求法;熟练掌握卷积、卷积和;了解信号的分解、掌握周期信号的傅立叶级数及频谱、熟练掌握傅立叶变换的性质、线性时不变系统的频域分析、取样定理;掌握拉普拉斯变换、熟练掌握拉普拉斯变换的性质、复频域分析;掌握离散信号的z变换、熟练掌握z变换的性质;在学生完成本课程的学习后,能够进行线性系统的分析工作,并为后续《数字信号处理》、《通信原理》、《自动控制原理》、《随机信号分析》等课程的学习奠定理论基础。通过本,,MATLAB软件,其中为单位阶跃函数。初始状态为0.
方程取拉普拉斯变换:
经整理得传递函数:
然后,利用MATLAB设计程序并用Simulink进行系统仿真。
3.3.2 内容
根据系统传递函数构建如图3-1所示的冲激响应仿真模型
图3-1冲激响应仿真模型
启动仿真就可以在示波器窗口中看到仿真结果的变化曲线,所求冲激响应如图3-2所示:
图3-2冲激响应
根据系统传递函数构建如图3-3所示的阶跃响应仿真模型
图3-3阶跃响应仿真模型
启动仿真就可以在示波器窗口中看到仿真结果的变化曲线,所求阶跃响应如图3-4所示:
图3-4所求阶跃响应
其中,在Simulink仿真系统建立的时候,双击Transfer Fcn模块,弹开其参数设置对话框,在分子、分母栏中填写所需的系数如图3 -5
图3-5 Transfer Fcn模块参数设置对话框
求解单位冲激响应及阶跃响应,其程序为:
clear all
a=[1 7 12];b=[2];
subplot(2,1,1)
impulse(b,a)
title('系统的冲激响应波形h(t)')
grid on
subplot(2,1,2)
step(b,a)
title('系统的阶跃响应波形g(t)')
grid on
图3-6单位冲激响应及阶跃响应波形
给定一个激励信号,建立仿真模型,仿真分析系统的响应。求解所设计的系统的零输入响应,零状态响应。,, ,零状态响应函数为:,构造如图3-7
图3-7零状态响应SIMULINK模型
修改MATLAB Fcn模块设置函数为exp(-u)如图3-8
图3-8 Fcn模块设置
双击Transfer Fcn弹出其参数设置对话框,在分子分母栏中填写所需的系数如图3-9
图3-9 Transfer Fcn参数设置对话框
③启动仿真,即可在示波器窗口中看到仿真结果的变化曲线如图3-10
图3-10零状态响应波形
零输入响应程序为:
clear all
r=[1 7 12]; %特征多项式系数矩阵
p=roots(r); %计算特征根
v=rot90(vander(p)); %形成范德蒙矩阵
y0=[1;1]; %初始条件矩阵
c=v\y0; %计算系数矩阵
t=linspace(0,10,101);
for k=1:length(p)
y_ji(k,:)=exp(p(k)*t); %生成解空间基信号
end
yt=c.'*y_ji; %计算零输入响应yzi
plot(t,yt,'b');
grid on
图3-11零输入响应波形
求解零状态响应,其程序为:
clear all
Num=2; %分子系数;
Den=[1 7 12]; %分母系数;
S_tf=tf(Num,Den); %利用分子分母系数形成LTI系统传输函数H(s);
t=0:0.01:;5
et=exp(-t); %激励为指数信号
figure(1);
lsim(S_tf,et,t); %计算系统在指数信号作用下的响应
grid on
零状态响应波形,如图3-12
图3-12零状态响应
频率响应程序为:
clear all
文档评论(0)