




下载本文档
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、求解系统的状态方程一、实验设备PC计算机, MATLAB 软件,控制理论实验台二、实验目的(1) 掌握状态转移矩阵的概念。学会用 MATLAB 求解状态转移矩阵(2) 学习系统齐次、非齐次状态方程求解的方法, 计算矩阵指数,求状态响应;(3) 通过编程、 上机调试,掌握求解系统状态方程的方法,学会绘制输出响应和状态响应曲线;(4) 掌握利用MATLAB 导出连续状态空间模型的离散化模型的方法。三、实验原理及相关基础(1) 参考教材P99101 “3.8利用MATLAB求解系统的状态方程”(2) MATLAB 现代控制理论仿真实验基础(3) 控制理论实验台使用指导四、实验内容(1)求下列系统矩阵
2、 A 对应的状态转移矩阵a)(b)代码:syms lambdaA=lambda 0 0;0 lambda 0;0 0 lambda;syms t;f=expm(A*t)c)代码:syms t;syms lambda;A=lambda 0 0 0;0 lambda 1 0;0 0 lambda 1;0 0 0 lambda;f=expm(A*t)(2) 已知系统a) 用MATLAB 求状态方程的解析解。 选择时间向量 t,绘制系统的状态响应曲线 观察并记录这些曲线。(1)代码:A=0 1; -2 -3;B=3;0;C=1 1;D=0;u=1;syms t;f=expm(A*t);% 状态转移矩阵
3、x0=0;s1=f*B*u;s2=int(s1,t,0,t)% 状态方程解析解状态曲线:( 2 ) A=0 1;-2 -3;syms t;f=expm(A*t);X0=1;0;t=0:0.5:10;for i=1:length(t);g(i)=double(subs(f(1),t(i);endplot(t,g)(3)状态转移矩阵syms lambdaA=lambda 0 0;0 lambda 0;0 0 lambda;syms tf=expm(A*t)精品资料b) 计算系统在初始状态作用下状态响应和输出响应的数值解 (用函数initial( ),绘制系统的状态响应曲线和输出响应曲线。 观察并记
4、录这些响应曲线, 然后将这 一状态响应曲线与 a)中状态响应曲线进行比较。代码:A=0 1; -2 -3;B=3;0;C=1 1;D=0;G=ss(A,B,C,D);t=0:0.5:10;x0=1;0y0,t,x0=initial(G,x0,t);plot(t,x0,'-',t,y0,'-')。记c) 根据 b) 中所得的状态响应的数值解,绘制系统的状态轨迹(用命令 plot(x(:,1), x(:,2)录系统状态转移的过程,结合 a)和 b)中的状态响应曲线分析这一过程。代码:A=0 1; -2 -3;B=3;0;C=1 1;D=0;t=0:0.01:10;x
5、0=1;0;G=ss(A,B,C,D)y,t,x=initial(G,x0,t);plot(x(:,1),x(:,2)2) 令初始状态为零,输入为 u(t)=1(t).a) 用MATLAB 求状态方程的解析解。 选择时间向量 t ,绘制系统的状态响应曲线。 观察并记录这些曲线。 代码: A=0 1; -2 -3; B=3;0; C=1 1;D=0;G=ss(A,B,C,D);y,t,x=step(G);plot(t,x)b) 计算系统在初始状态作用下状态响应和输出响应的数值解 , 绘制系统的状态响应曲线和输出响应曲线。观察并记录这些响应曲线,然后将这一状态响应曲线与 a).中状态响应曲线进行比
6、较。代码:A=0 1; -2 -3;B=3;0;C=1 1;D=0;G=ss(A,B,C,D);G=ss(A,B,C,D);t=0:0.5:10;x0=1;-1;y0,t,x0=initial(G,x0,t);plot(t,x0,'-',t,y0,'-')c) 根据 b) 中所得的状态响应的数值解,绘制系统的状态轨迹。记录系统状态转 移的过程,结合 a) 和 b) 中的状态响应曲线分析这一过程。代码:A=0 1; -2 -3;B=3;0;C=1 1;D=0;t=0:0.5:10;G=ss(A,B,C,D);x0=0 0;y0,t,x0=initial(G,x0,
7、t);plot(t,x0,'-',t,y0,'-')绘制系统的状态响应曲线、 输出响应曲线和状态轨迹。 观察和分析这些响应曲线和状态轨迹是否是( 1)和( 2)中的响应曲线和状态轨迹的叠加代码:A=0 1; -2 -3;B=3;0;C=1 1;D=0;t=0:0.01:10;x0=1 -1;G=ss(A,B,C,D);y,t,x=initial(G,x0,t);plot(t,x)4) 令初始状态为零,输入为 u(t)=3sin(5t) 。计算状态响应和输出响应的数值解 lsim( ) ,并绘制系统的状态响应曲线、输出响应曲线和状态轨迹。(用函数代码:A=0 1;
8、 -2 -3;B=3;0;C=1 1;D=0;t=0:0.01:10;u=3*sin(5*t);G=ss(A,B,C,D); y,t,x=lsim(G,u,t); plot(t,x)3)已知系统1) 当输入为 u(t)= (t)时,用函数 initial( ) 和 impulse( ) 求解系统的状态响应和输 出响应的 数值解,并绘制系统的状态响应曲线、输出响应曲线和状态轨迹。状态响应:A=0,1,0;0,0,1;-6,-11,-6;B=0;0;1;C=6,0,0;D=0;t=0:0.5:10;G=ss(A,B,C,D);x0=1,0,-1;y,t,x=initial(G,x0,t); u=o
9、nes(size(t);plot(t,x,t,y)y,t,x=initial(G,x0,t);精品资料输出响应:A=0,1,0;0,0,1;-6,-11,-6;B=0;0;1;C=6,0,0;D=0;t=0:0.01:10;u=ones(size(t);G=ss(A,B,C,D);y,t,x=lsim(G,u,t);plot(t,x)2) 当输入为 u(t)=1(t) 时,用函数 initial( )和 step( )求解系统的状态响应和输出 响应的数值解,并绘制系统的状态响应曲线、输出响应曲线和状态轨迹。状态响应A=0,1,0;0,0,1;-6,-11,-6;B=0;0;1;C=6,0,0;
10、D=0;t=0:0.5:10;G=ss(A,B,C,D);x0=1,0,-1;u=step(G);plot(t,x,t,y)输出响应:A=0,1,0;0,0,1;-6,-11,-6;B=0;0;1;C=6,0,0;D=0;G=ss(A,B,C,D);y,t,x=step(G);plot(t,x)3) 当输入为 u(t)= t 时,用函数 initial( )和 lsim( )求解系统的状态响应和输出响应的数值解,并绘制系统的状态响应曲线、输出响应曲线和状态轨迹状态响应A=0,1,0;0,0,1;-6,-11,-6;B=0;0;1;C=6,0,0;D=0;t=0:0.5:10;G=ss(A,B,
11、C,D);x0=1,0,-1;u=t;plot(t,x,t,y)输出响应A=0,1,0;0,0,1;-6,-11,-6;B=0;0;1;C=6,0,0;D=0;t=0:0.01:10;u=t;G=ss(A,B,C,D);y,t,x=lsim(G,u,t);plot(t,x)4)当输入为时,用函数 initial( ) 和 lsim( ) 求解系统的状态响应和输出响应的数值解,并绘制系统的状态响应曲线、输出响应曲线和状态轨迹 . 状态响应A=0,1,0;0,0,1;-6,-11,-6;B=0;0;1;C=6,0,0;D=0;t=0:0.5:10;G=ss(A,B,C,D);x0=1,0,-1;u
12、=sin(t);plot(t,x,t,y)输出响应A=0,1,0;0,0,1;-6,-11,-6;B=0;0;1;C=6,0,0;D=0;t=0:0.01:10;u=sin(t);G=ss(A,B,C,D);y,t,x=lsim(G,u,t);plot(t,x)y,t,x=initial(G,x0,t);精品资料(4) 已知一个连续系统的状态方程是若取采样周期秒 0.05 T1)试求相应的离散化状态空间模型;代码:syms T;A=0 1;-25 -4;B=0;1;Gz,Hz=c2d(A,B,T) exp(- 2*T - 21(1/2)*T*i)/2 + exp(- 2*T + 21(1/2)
13、*T*i)/2 + (21(1/2)*exp(- 2*T -21(1/2)*T*i)*i)/21 - (21(1/2)*exp(- 2*T + 21(1/2)*T*i)*i)/21,(21(1/2)*exp(- 2*T - 21(1/2)*T*i)*i)/42 - (21(1/2)*exp(- 2*T + 21(1/2)*T*i)*i)/42 - (21(1/2)*exp(- 2*T -21(1/2)*T*i)*25*i)/42 + (21(1/2)*exp(- 2*T + 21(1/2)*T*i)*25*i)/42, exp(- 2*T -21(1/2)*T*i)/2 + exp(- 2*T
14、 + 21(1/2)*T*i)/2 - (21(1/2)*exp(- 2*T - 21(1/2)*T*i)*i)/21+ (21(1/2)*exp(- 2*T + 21(1/2)*T*i)*i)/211/25 - exp(- 2*T + 21(1/2)*T*i)/50 - (21(1/2)*exp(- 2*T - 21(1/2)*T*i)*i)/525 +(21(1/2)*exp(- 2*T + 21(1/2)*T*i)*i)/525 - exp(- 2*T - 21(1/2)*T*i)/50(21(1/2)*exp(- 2*T - 21(1/2)*T*i)*i)/42 - (21(1/2)*exp(- 2*T + 21(1/2)*T*i)*i)/422) 分析不同采样周期下,离散化状态空间模型的结果A=0 1;-25 -4;B=0;1;Gz,Hz=c2d(A,B,0.05)Gz =0.9709 0.0448-1.1212 0.7915Hz =0.0012 0.0448五、实验心
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 氢液化器可靠性分析及其安全控制策略研究
- 复杂环境下汽轮发电机绝缘监测系统设计与研究
- 基于大语言模型的智能合约漏洞检测研究
- 秦皇岛市人民医院机械清创技术考核
- 2025年版旅游团队租车合同协议范本发布平台
- 基于慷慨教师的知识蒸馏及人脸交换轻量化研究
- 2025年电子信息产品维修合同条件
- 2025年仓储管理合同协议书
- 2025云南考察调(流)动县中医医院中层干部7人考试模拟试题及答案解析
- 2025年焊工试题库及答案
- 人教版初中物理八年级上册《质量》说课(附教学反思、板书)课件
- 工程欠款协议书范文
- 2025年中国鞋撬市场调查研究报告
- 事故隐患内部报告奖励制度
- 2024年广东省中考满分作文《当好自己故事的主角》2
- DB37-T 5001-2021 住宅工程外窗水密性现场检测技术规程
- 冬季心血管疾病的预防
- 心衰患者出入量管理研究进展
- 安全部经理竞聘汇报
- 《物料摆放规范》课件
- 《智能建造技术与装备》 课件 第二章 BIM技术与应用
评论
0/150
提交评论