版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、广西大学实验报告纸序号学号姓名贡献排名实验报告分数1(组长):130212023712(组员):1302120238 2实验项目线性系统状态空间分析与运动解【实验时间】2015年11月12日【实验地点】课外(宿舍)【实验目的】1、学会利用MATLAB实现离散系统传递函数模型的生成2、学会利用MATLAB将连续系统离散化【实验设备与软件】1、MATLAB/Simulink数值分析软件2、计算机一台【实验原理】1、求矩阵特征值和特征向量命令格式V J=eig(A)Cv=eig(A)说明:V特征向量,J是Jordan型,cv是特征值列向量2、求运动的方法(1)利用Laplace逆变换-适合于连续/离
2、散线性系统采用ilaplace/iztrans对传递函数求逆,这种方法一般是零输入情况下求响应。(2)用连续(离散)状态转移矩阵表示系统解析解-适合于线性定常系统对连续定常系统有:假设初始时刻为零,LTI系统的解析解为。若u(t)是单位阶跃输入,则上述解可写成。进一步简化为:对离散线性定常系统有:(3)状态方程的数值分析方法-适合于连续线性系统和非线性系统采用直接数值积分很容易的处理各种定常/时变和线性/非线性系统。有很多数值积分方法,其中有一类预测-修正数值积分方法+自适应步长调整的算法比较有效。在MATLAB/Simulink中包含的多种有效的、适用于不同类型的ODE求解算法,典型的是Ru
3、nge-Ktuta算法,其通常使用如下的函数格式:t,x=ode45(odefun,ti,tf,x0,options)-采用四阶、五阶Runge-Ktuta算法t,x=ode23(odefun,ti,tf,x0,options)-采用二阶、三阶Runge-Ktuta算法说明:a.这两个函数是求解非刚性常微分方程的函数。 b.参数options为积分的误差设置,取值为相对误差reltol和绝对误差abstol;ti,tf求解的时间范围;x0是初值是初值向量;t,x是解。(4)利用CotrolToolBox的离散化求解函数-适合于TLI系统用step()/impulse()函数求取阶跃输入/冲激输
4、入时系统的状态响应:当系统G是连续的情况下:调用y,t,x=step/impulse(G )会自动对连续系统G选取采样时间范围和周期;调用y,t,x=step/impulse(G ,ti:Ts:tf)由用户自己定义对连续系统G的样时间范围和周期;当系统G是离散的情况下:调用y,t,x=step/impulse(G )会按离散系统G给出的采样周期计算;调用y,t,x=step/impulse(G ,ti:Ts:tf)是Ts必须与离散系统G的采样时间范围和周期一致。另外lsim()函数调用格式:y,x,t=lsim(G,u,ti,TS,tf,x0)零输入响应调用函数initial(),格式:y,x
5、,t=(G,x0)(5)利用simulink环境求取响应-适用于所有系统求取响应使用simulink求取线性或非线性系统的响应,调用格式如下:t,x,y=sim(XX.mdl,ti:Ts:tf,options,u)【实验内容】已知线性系统:已知线性系统1、利用Matlab求零状态下的阶跃响应(包括状态和输出),生成两幅图:第一幅绘制各状态响应曲线并标注;第二幅绘制输出响应曲线。状态响应曲线:A=-21 19 -20;19 -21 20;40 -40 -40;B=0;1;2;C=1 0 2;D=0; %输入状态空间模型各矩阵,若没有相应值,可赋空矩阵X0=0;0;0; % 输入初始状态sys=s
6、s(A,B,C,D); %构造传递函数y,x,t=step(sys); % 绘以时间为横坐标的状态响应曲线图plot(t,x);grid;title('状态响应曲线')输出响应程序:A=-21 19 -20;19 -21 20;40 -40 -40; B=0;1;2;C=1 0 2;D=0;X0=0;0;0num,den=ss2tf(A,B,C,D,1); sys=tf(num,den);step(sys)gridtitle('输出响应曲线')图一(状态响应曲线)图二(输出响应曲线)2、利用Matlab求零状态下的冲激响应(包括状态和输出),生成两幅图:第一幅绘
7、制各状态响应曲线并标注;第二幅绘制输出响应曲线。状态响应曲线程序:A=-21 19 -20;19 -21 20;40 -40 -40;B=0;1;2;C=1 0 2;D=; %输入状态空间模型各矩阵,若没有相应值,可赋空矩阵x0=0;0;0; % 输入初始状态sys=ss(A,B,C,D); %构造传递函数y,x,t= impulse(sys);plot(t,x);grid;title('状态响应曲线')输出响应曲线程序:A=-21 19 -20;19 -21 20;40 -40 -40;B=0;1;2;C=1 0 2;D=0;X0=0;0;0num,den=ss2tf(A,B
8、,C,D,1);sys=tf(num,den);impulse(sys);grid;title('')图三(状态响应曲线')图四(输出响应曲线)3、若控制输入为,且初始状态为,求系统的响应,要求a.在simulink只能够画出模型求响应,生成两幅图:第一幅绘制各状态响应曲线并标注;第二幅绘制输出响应曲线。程序如下:t=0:0.01:5;u=(1+exp(-t).*cos(5*t).*(t<3)+1*(t>=3);t=t'u=u'ut=t,u;t1,x,y=sim('shiyan5.mdl',t,ut);plot(t1,x)fi
9、gure(2);plot(t1,y)创建的模型图如下:图五(模型图)b.编写.m文件求响应,生成两幅图:第一幅绘制各状态响应曲线并标注;第二幅绘制输出响应曲线。状态响应曲线:t=0:0.02:5;u=(1+exp(-t).*cos(5*t).*(t<3)+1*(t>=3);t=t'u=u'A=-21 19 -20;19 -21 20;40 -40 -40;B=0;1;2;C=1 0 2;D=0; %输入状态空间模型各矩阵,若没有相应值,可赋空矩阵 X0=0.2;0.2;0.2; % 输入初始状态u=(t=0); %就是个条件判断,只有t=0的时候,u才为“1”sys
10、=ss(A,B,C,D); %构造传递函数plot(t,x);grid;title('状态响应曲线')输出响应曲线:plot(t,y);grid;title('输出响应曲线')图六(状态响应曲线)图七(输出响应曲线)4、以阶跃输入情况下的,分析各模块对响应有什么影响。图八(阶跃输入时)阶跃输入的图像到答稳定时间快,曲线平滑5、求系统的传递函数在MATLAB软件Command Window窗口中输入以下程序A=-21 19 -20;19 -21 20;40 -40 -40;B=0;1;2;C=1 0 2;D=0;num,den=ss2tf(A,B,C,D,1);p
11、rintsys(num,den)程序运行结果为图七6、若采用K增益负反馈,绘制闭环根轨迹图,并对根轨迹加以描述说明。A=-21 19 -20;19 -21 20;40 -40 -40;B=0;1;2;C=1 0 2;D=0;num,den=ss2tf(A,B,C,D,1);rlocus(num,den);gridtitle('K增益负反馈闭环根轨迹图')图九(K增益负反馈闭环根轨迹图)采用K增益负反馈,画出如图所示的根轨迹图。由图可知,共有3条根轨迹,第一条最终趋于原点;第二条收敛在2060之间;第三条最终趋于无穷远处。7、在Matlab中绘制Bode图和Nyquist图,并对
12、图给予说明。绘制Bode图:A=-21 19 -20;19 -21 20;40 -40 -40;B=0;1;2;C=1 0 2;D=0;sys=tf(num,den)bode(num,den)gridtitle('Bode图')汇出的波特图如图所示,由图可知,对复制响应分析可得,交越频率在转折频率之后,故复制的变化主要发生在低频段。对相频特性进行分析,可知此系统的相频特性角度均为负值,并且最后的相角是趋于-90度的。绘制Nyquist图:nyquist(sys)title('Nyquist图 ')图十(波特图)图十一(Nyquist图)画出的奈奎斯特图如上所示,根据此图可知,此系统是稳定的系统,由奈奎斯特
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- GB/T 23025-2024信息化和工业化融合管理体系生产设备运行管控信息模型分类与应用指南
- 电工电子技术(第3版) 课件 2.1 正弦交流电的基本特征
- 项目投标与招标管理规范制度
- 总承包公司永临结合做法选用图册
- 维修员工工作总结
- 心理健康教育直播课心得体会范文(30篇)
- 物业管理制度15篇
- 中国著名书法家简介
- 【培训课件】节约里程法
- 倾斜角与斜率课件
- 2024山东省化工行业职业技能大赛(化工总控工)试题库-下(判断、简答题)
- 2023年重庆轨道交通招聘考试试题及答案
- 历史人教部编版八年级(上册)22.抗日战争的胜利课件(25张)2024版新教材
- 19 父爱之舟 (教学设计)-2024-2025学年小学语文五年级上册统编版
- 2024年教资考试时政高频考点141条
- 15《搭船的鸟》(教学设计)2024-2025学年统编版语文三年级上册
- 个人简历模板(5套完整版)
- 2024至2030年中国传染病医院产业发展动态及未来前景展望报告
- 2024年保安员资格考试初级理论知识判断题库及答案共150题
- 内部追偿合同协议书
- 中国近现代史纲要智慧树知到答案2024年北京师范大学等跨校共建
评论
0/150
提交评论