状态空间课设设计.docx
- 文档编号:24134200
- 上传时间:2023-05-24
- 格式:DOCX
- 页数:18
- 大小:329.53KB
状态空间课设设计.docx
《状态空间课设设计.docx》由会员分享,可在线阅读,更多相关《状态空间课设设计.docx(18页珍藏版)》请在冰豆网上搜索。
状态空间课设设计
西安石油大学
课程设计
电子工程学院电气专业
班
题目控制系统的状态空间法设计8
学生
指导老师
二○一○年十二月
《自动控制理论I》
《自动控制理论》课程设计任务书
题目
状态空间法设计
学生姓名
学号
2008
专业班级
设
计
内
容
与
要
求
一.设计内容:
已知一系统的状态方程为
⑴试判别系统的可控性并求状态反馈增益矩阵k,使得系统的闭环特征值(即闭环极点)为
=-10、
==-1+j
、
=-1-j
。
⑵试判别系统的可观性对系统设计状态观测器,使得系统的闭环极点为
=-4、
=-3+j、
=-3-j。
二.设计要求:
(1)编程绘制原系统的单位阶跃响应曲线,并计算出原系统的动态性能指标;
(2)编程判断该系统的可控性;
(3)编程设计校正方案(得到相应的状态反馈参数或状态观测器参数);
(4)利用SIMULINK绘制校正前、后系统状态模拟图;
(5)编程绘制校正后系统的单位阶跃响应曲线,并计算出校正后系统的动态性能指标;
(6)整理设计结果,得出设计结论并提交设计报告。
起止时间
2010年12月24日至2010年12月30日
指导老师签名
年月日
系(教研室)主任签名
年月日
学生签名
2010年1月1日
1.1任务书..........................................................2
1.2设计思想及内容..................................................4
2.21运用MATLAB编程.................................................5
2.2编制的程序......................................................5
2.22在SIMULINK中绘制状态图.......................................11
3.1结论........................................................13
4.1设计总结.....................................................13
5.1参考文献.....................................................14
2.1设计内容及思想:
1)内容:
已知一系统的状态方程为
⑴试判别系统的可控性并求状态反馈增益矩阵k,使得系统的闭环特征值(即闭环极点)为
=-10、
==-1+j
、
=-1-j
。
⑵试判别系统的可观性对系统设计状态观测器,使得系统的闭环极点为
=-4、
=-3+j、
=-3-j。
2)思想:
编程思想:
利用状态反馈任意极点配置闭环极点的充要条件是被控系统可控的;
利用输出至状态微分的反馈任意配置闭环极点的充要条件是被控系统可观测。
仿真思想:
设被控对象动态方程为
x’=Ax+Bu,y=Cx
构造一个动态方程与上式相同且能用计算机实现的模拟被控系统
x’1=Ax1+Bu,y1=Cx
x1,y1分别为模拟系统的状态向量和输出向量,是被控对象状态向量和输出向量的估值。
当模拟系统与被控对象的初始状态向量相同时,在同一输入作用下,有x1=x,可用x1作为状态反馈所需要的信息。
但是,被控对象的初始状态可能很不相同,模拟系统中积分器初始条件的设置只能估计,因而两个系统的初始状态总有差异,即使两个系统的A,B,C阵完全一样,也必定存在估计状态与被控对象实际状态的误差(x1-x),难以实现所需要的反馈。
但是,(x1-x)必定导致(y1-y)的存在,而被控系统的输出量总是可以用传感器测量的,于是根据这一般反馈控制原理,将(y1-y)负反馈至x’1处,控制(y1-y)尽快逼近于零,从而使(x1-x)尽快逼近于零,便可以利用x1来形成状态反馈。
H为观测器输出反馈阵,它把(y1-y)负反馈至x’1处,是为配置观测器极点,提高其动态性能。
2.2编制的程序:
2.2.1运用MATLAB编程:
1.可控性校正前:
A=[120;3-11;020];B=[0;0;1];
C=[-111];D=[0];%设定系数矩阵
p=[-10-1+i*sqrt(3)-1-i*sqrt(3)];%极点配置
CAM=ctrb(A,B);N=size(A);n=N
(1);
[num,den]=ss2tf(A,B,C,D);
G=tf(num,den);
step(G);
ifdet(CAM)~=0
rcam=rank(CAM);
ifrcam==n
disp('Systemiscontrolled')
elseifrcam disp('Systemisnotcontrolled') end elseifdet(CAM)==0 disp('Systemisnotcontrolled')%判定系统可控性 end ifrcam==n K=place(A,B,p)%极点配置,计算估计器增益矩阵K end 运行结果: Systemiscontrolled 由图可知系统是发散的 K= 77.500033.000012.0000 2.可观性校正前: A=[120;3-11;020];B=[1;0;0]; C=[-111];D=[0];%设定系数矩阵 p=[-4-3+i-3-i];%极点配置 CAM=ctrb(A,B);N=size(A);n=N (1); [num,den]=ss2tf(A,B,C,D); G=tf(num,den); step(G); ob=obsv(A,C) roam=rank(ob); ifroam==n disp('Systemisobservable') elseifroam~=n disp('Systemisnotobservable') end%判定系统可观测性 K=acker(A',C',p);H=K' AHC=A-H*C B%极点配置,计算观测器输出反馈矩阵H 运行结果: ob= -111 2-11 -17-1 Systemisobservable H= 27 24 13 AHC= 28-25-27 27-25-23 13-11-13 B= 1 0 0 3.可控性校正后: G1=ss(A-B*K,B,C,D) step(G1) 运行结果: a= x1x2x3 x1-11-10-12 x2-9-13-11 x3-12-10-12 b= u1 x10 x20 x31 c= x1x2x3 y1-111 d= u1 y10 Continuous-timemodel. >> 4.可观性校正后: G1=ss(AHC,B,C,D) step(G1) 运行结果: a= x1x2x3 x128-25-27 x227-25-23 x313-11-13 b= u1 x11 x20 x30 c= x1x2x3 y1-111 d= u1 y10 Continuous-timemodel 可控性校正后的系统阶跃响应: 根据上图可得: 上升时间 =0.594;峰值时间 =-0.306;调节时间 =4.14; =0.812;超调量 %=22.5%。 可观性校正后的系统阶跃响, 根据上图可得: 上升时间 =0.518;峰值时间 =-0.208;调节时间 =1.85; =0.548;超调量 %=3.88% (二)2.22在SIMULINK绘制状态图: 可控性校正前: 可控性校正前的阶跃响应曲线: 可控性校正后: 可控性校正后阶跃响应曲线 可观校正前 可观校正前的阶跃响应曲线 可观校正后 可观校正后阶跃响应曲线 3.1结论: 1.该状态方程所描述的系统的阶跃响应是发散的。 2.系统通过校正后,其响应是趋于稳定的,动态性能为: 上升时间 =0.594;峰值时间 =-0.306;调节时间 =4.14; =0.812;超调量 %=22.5% 3.此系统经过状态反馈后实现了题中要求的极点配置,同时系统经过状态反馈后可以按输入的意志达到期望的状态和输出响应。 4.1设计总结: 经过此次状态空间课程设计,我发现了Matlab在程序设计方面的方便之处,更体会到了Simulink软件在仿真方面的强大功能,调试方便,可以节省打来那个的编程时间。 更重要的是我学会了用Simulink仿真实际问题,它使我更加方便的观看到了电路的功能。 5.1参考文献: 薛朝妹,霍爱清,《自动控制理论》课程设计指导书,西安: 西安石油大学电子工程学院,2007年 瞿亮,基于的控制系统计算机仿真,北京: 清华大学出版社,2006年
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 状态 空间 设计