倒立摆系统的建模与matlab仿真

倒立摆系统的建模与matlab仿真

ID:23295933

大小:286.50 KB

页数:12页

时间:2018-11-06

倒立摆系统的建模与matlab仿真_第1页
倒立摆系统的建模与matlab仿真_第2页
倒立摆系统的建模与matlab仿真_第3页
倒立摆系统的建模与matlab仿真_第4页
倒立摆系统的建模与matlab仿真_第5页
资源描述:

《倒立摆系统的建模与matlab仿真》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库

1、word资料下载可编辑倒立摆系统的建模及Matlab仿真1.系统的物理模型考虑如图(1)所示的倒立摆系统。图中,倒立摆安装在一个小车上。这里仅考虑倒立摆在图面内运动的二维问题。图(1)倒立摆系统假定倒立摆系统的参数如下。摆杆的质量:m=0.1g摆杆的长度:l=1m小车的质量:M=1kg重力加速度:g=9.8m/摆杆的质量在摆杆的中心。设计一个控制系统,使得当给定任意初始条件(由干扰引起)时,最大超调量d≤10%,调节时间ts≤4s,通过小车的水平运动使倒立摆保持在垂直位置。2.系统的数学模型2.1建立倒置摆的运动方程并将其线性化。为简化问题,在数

2、学模型中首先假设:1)摆杆为刚体;2)忽略摆杆与支点之间的摩擦;3)忽略小车与接触面间的摩擦。设小车瞬时位置为z,摆心瞬时位置为(),在u作用下,小车及摆均产生加速远动,根据牛顿第二定律,在水平直线远动方向的惯性力应与u平衡,于是有专业技术资料word资料下载可编辑即:①绕摆轴转动的惯性力矩与重力矩平衡,因而有即:②以上两个方程都是非线性方程,为求得解析解,需作线性化处理。由于控制的目的是保持倒立摆直立,在试驾合适的外力条件下,假定θ很小,接近于零时合理的,则,且可忽略项。于是有③④联立求解可得2.2列写系统的状态空间表达式。选取系统变量,则即代

3、入数据计算得到:专业技术资料word资料下载可编辑3.设计控制器3.1判断系统的能控性和稳定性,rank()=4,故被控对象完全可控由特征方程解得特征值为0,0,。出现大于零的特征值,故被控对象不稳定3.2确定希望的极点希望的极点n=4,选其中一对为主导极点和,另一对为远极点,认为系统性能主要由主导极点决定,远极点只有微小影响。根据二阶系统的关系式,先确定主导极点可得,于是取;取误差带,则,闭环主导极点为=-10.8j,远极点选择使它和原点的距离大于主导极点与原点距离的5倍,取3.3采用状态反馈方法使系统稳定并配置极点状态反馈的控制规律为,;状态

4、反馈系统的状态方程为,其特征多项式为⑤希望特征多项式为⑥专业技术资料word资料下载可编辑比较以上两式系数,解得状态反馈矩阵4.设计全维观测器4.1判断系统的能观性,rank()=4,故被控对象完全可观4.2确定观测器的反馈增益全维观测器的动态方程为;其特征多项式为⑦取观测器的希望极点为:-45,-45,-3+3j,-3-3j;则希望特征多项式为⑧比较以上两式系数,解得观测器反馈矩阵5.降维状态观测器的设计5.1建立倒置摆三维子系统动态方程设小车位移z由输出传感器测量,因而无需估计,可以设计降维(三维)状态观测器,通过重新排列被控系统变量的次序,

5、把需由降维状态观测器估计的状态变量与输出传感器测得的状态变量分离开。将z作为第四个状态变量,则被控系统的状态方程和输出方程变换为专业技术资料word资料下载可编辑⑨简记为:式中,,,,,=0,,被控系统的n-q维子系统动态方程的一般形式为,式中,为子系统输出量。故倒置摆三维子系统动态方程为5.2.判断子系统的可观测性A1=[0-10;001;0110];C1=[100];专业技术资料word资料下载可编辑Qg1=obsv(A1,C1);r=rank(Qg1)运行Matlab程序;结果为r=3,故该子系统可观测降维状态观测器动态方程的一般形式为式中

6、h=。考虑被控对象参数,单倒置摆降维观测器动态方程的一般形式为5.3确定三维状态观测器的反馈矩阵h三维状态观测器的特征多项式为设希望的观测器闭环极点为-45,-3+3j,-3-3j,则希望特征多项式为比较以上两式系数,解得h=故所求三维状态观测器的动态方程为6.Matlab仿真分析6.1源程序通过Matlab对用全维状态观测器实现状态反馈的倒置摆系统进行仿真分析,下面是文件名为Inversion_pendulum_system.m的源程序%倒立摆系统建模分析专业技术资料word资料下载可编辑%a)判断系统能控性和能观性clearall;clcA=

7、[0100;00-10;0001;00110];B=[0;1;0;-1];C=[1000];D=0;Uc=ctrb(A,B);rc=rank(Uc);n=size(A);ifrc==ndisp('Thesystemiscontrolled.')elseifrc

8、b)判断系统稳定性P=poly(A),v=roots(P)Re=real(v);if(length(find(Re>0))~=0)dis

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。