资源描述：
《倒立摆系统的建模及matlab仿真》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、倒立摆系统的建模及MatIab仿真1.系统的物理模型考虑如图（1）所示的倒立摆系统。图中，倒立摆安装在一个小车上。这里仅考虑倒立摆在图面内运动的二维问题。uv//////////////////图（1）倒立摆系统假定倒立摆系统的参数如下。摆杆的质量：m=0.lg摆杆的长度：/=1nH、午.的质量：M=lkg重力加速度:g=9.8m/?摆杆的质量在摆杆的中心。设计一个控制系统，使得当给定任意初始条件（由干扰引起）时，最大超调量S<10%,调节时间ts^4s,通过小车的水平运动使倒立摆保持在垂直位置。2.系统的数学
2、模型2.1建立倒置摆的运动方程并将其线性化。为简化问题，在数学模型中首先假设：1）摆杆为刚体；2）忽略摆杆与支点之间的摩擦；3）忽略小车与接触面间的摩擦。设小车瞬时位置为z，摆心瞬时位置为（z+/sin~），在u作用下，小车及摆均产生加速远动，根据牛顿第二定律，在水平直线远动方向的惯性力应与u平衡，于是有M-++/sin0）=udtdt~即：,2Z7?-^T(Z+/sind)/cos汐二mg/sin沒(M+m)z+ml3cos3-mid2sin6-u绕摆轴转动的惯性力矩与重力矩平衡，因而冇即：zcos^+/^c
3、os23-Id2sin^cos^=gsin沒②以上两个方程都是非线性方程，为求得解析解，需作线性化处理。由于控制的I；的是保持倒立摆直立，在试驾合适的外力条件下，假定0很小，接近于零时合理的，Kijsin^«^cos^«l,且可忽略少沒项。于是有99(M+m)z+mid=uz+lO=gO联立求解可得MlMl2.2列写系统的状态空间表达式选取系统变量&，'，〜々，*=则X2X3=XA.(M+m)1xd=x,u4Ml•Ml"0100’01•00哪01•ax=—z—Mx+Mdte00010•e00(M+m)g01Ml
4、Mlu=Ax+Buy-x-[iooo]x=Cx代入数据计算得到:010000-10000100II0,B=[010-l]r,C=[l000],Z)=01.设计控制器3.1判断系统的能控性和稳定性Qk=[bABA2BA3b]010一110-101110•110rank(Qk)=4,故被控对象完全可控由特征方程
5、2/-A
6、=/l2a2-11)=0解得特征值为0,0，土a/TT。出现大于零的特征值，故被控对象不稳定3.2确定希望的极点希望的极点n=4,选其中一对为主导极点&和*2,另一对为远极点，认为系统性能主要由
7、主导极点决定，远极点只有微小影响。根据二阶系统的关系式，先确定主导极点主导极点为<7<0.10.59,于是取《=0.6;取误差带A=0.02有〔=,贝ljty”=1.67，闭环1±0.8j,远极点选择使它和原点的距离大于主异极点与原点距离的5倍，取^^=-153.3采用状态反馈方法使系统稳定并配置极点状态反馈系统的状态方程为状态反馈的控制规律为w=-h，k=[A:okxk2k3]；x={A-BK)x^Bv,其特征多项式为-(A-BK)=/I4+(k}-+(k{)-k2-l1)A2-10^,2-10^()⑤希望
8、特征多项式为u+15)2(2+1-0.8j)(又+1+0.8y)=24+3223+286.64#+499.22+369比较以上两式系数，解得状态反馈矩阵[-36.9-49.92-334.54-81.92]1.设计全维观测器0000-100-1，rank(2P=4,故被控对象完全可观4.1判断系统的能观性10Qs=[cATC（，）2c（，）3c]=2004.2确定观测器的反馈増益全维观测器的动态方程为x=(A-GC)x+Bv+GCx；其特征多项式为
9、2/-(A-GC)
10、=/14+g()乂3+(&-11)A2+(-
11、1lg0-g2)2+(-l1§丨-g3)⑦取观测器的希望极点为：-45,-45,-3+3j,-3-3j；则希望特征多项式为(A+15)2U+1-0.8+1+0.87)=/+9623+2583A2+13770A+34650⑧比较以上两式系数，解得观测器反馈矩阵G=[962594-14826-64984]"2.降维状态观测器的设计5.1建立倒置摆三维子系统动态方程设小车位移Z由输出传感器测量，因而无需估计，可以设计降维(三维)状态观测器，通过重新排列被控系统变量的次序，把需由降维状态观测器估计的状态变量与输出传感器
12、测得的状态变量分离开。将Z作为第四个状态变量，则被控系统的状态方程和输出方程变换为參z1e0曇d十-1z0_0-1000010011001000u•ze•ezd~dt〉，=[0001]離Zeez简记为:II161
13、办2-I0-10001，Ai2=[O0O]7',b{=[iO-l]70110"x2=z=yfA2]=[l0O]r，A2i=O，b2=Q，/,=1被控系统的n-q维子系统动态方
