• 直线型一阶倒立摆2---建模


    三、直线型一阶倒立摆模型建立

        一级倒立摆系统是一个不稳定的系统,需要对其进行机理建模。 在研究过程中,应忽略空气摩擦、等,而后可将倒立摆系统进行抽象化,认为其由小车和匀质刚性杆两部分组成并对这两部分进行如图所示的受力分析:

     

    其中为小车的质量和摆杆的质量;b、Fx分别为小车的摩擦系数、施加在小车上的作用力和小车的位置[8];I和 分别为摆杆的惯量和摆杆转动轴心到质心的长度; 和  分别为摆杆与竖直向上方向和竖直向下方向的夹角;N P 分别为摆杆作用力的水平与竖直分量。

     

    小车水平方向的合力: 

                                                          Mfrac{mathrm{d^{2}x} }{mathrm{d} t^{2}}=F-Bfrac{mathrm{dx} }{mathrm{d} t}-N        (1)     

     摆杆水平方向的合力:  

                                                            N=mfrac{mathrm{d^{2} (x-lsin	heta )} }{mathrm{d} t^{2}}             (2)      

    摆杆水平方向运动方程:

                                                            (M+m)x^{''}+bx^{'}+ml	heta ^{''}cos	heta -ml	heta ^{'}sin	heta =F        (3)

    摆杆力矩平衡方程:

                                                               -PIsin	heta -Nlsin	heta =I	heta ^{''}           (4)

    摆杆在竖直方向的合力:

                                                                        P=mg+ml	heta ^{''}sin	heta +ml	heta ^{'}cos	heta         (5)

    可得到摆杆在竖直方向的运动方程:

                                                                         (1+ml^{2})	heta ^{''}+mglsin	heta =-mlx^{''}cos	heta         (6)

    摆杆竖直方向运动方程:

                                                                            I+Ml^{''}phi -mglphi =mlx             (7)

    将作用力 F 用 u代替,同时进行线性化,即得到:     

                                                                             (I+ml^{2})phi ^{''}-mglphi =mlx^{''}      (8)

                                                                              (M+m)x^{^{''}}+bx^{'}-mlphi ^{''}=u    (9)                  

    其中	heta =pi +phiphi为小角度。

    质量均匀分布的摆杆,对于式(8)有                       I=frac{1}{3}ml^{2}                                   (11)

    由式(8)、(11)得到:                                               (frac{4}{3}ml^{2})phi ^{''}-mglphi=mlx^{''}       (12)

    对质量均匀摆杆,取X=[egin{matrix} x &x^{'} &phi & phi ^{'} end{matrix}],u^{'}=x^{''}可得到线性一阶直线倒立摆状态空间描述:

                                                                                 egin{bmatrix} x^{'}\ x^{''}\ phi ^{'}\ phi ^{''} end{bmatrix}=egin{bmatrix} 0 & 1 & 0 & 0\ 0 & 0 & 0 & 0\ 0& 0& 0 &1 \ 0& 0 &frac{3g}{4l} & 0 end{bmatrix} egin{bmatrix} x\ x^{'}\ phi \ phi ^{'} end{bmatrix}+egin{bmatrix} 0\ 1\ 0\ frac{3}{4l} end{bmatrix}u^{'}         (13)

                                                                                 y=egin{bmatrix} x\ phi end{bmatrix}=egin{bmatrix} 1 & 0& 0& 0\ 0& 0& 1& 0 end{bmatrix}egin{bmatrix} x\ x^{'}\ phi \ phi ^{'} end{bmatrix}                         (14)

     

        对系统进行可控性分析,由控制矩阵Qc=[B  AB  A^{2}B  A^{3}B]=egin{bmatrix} 0 &1 &0 &0 \ 1&0 &0 &0 \ 0 & frac{3}{4l} &0 &frac{9g}{16l^{2}} \ frac{3}{4l} &0 & frac{9g}{16l^{2}} &0 end{bmatrix},用matlab计算可知,Qc的秩为4,系统可控。系统可进行状态变量的极点配置。

       对系统进行可观测性分析,由观测矩阵Qo=[C  CA  CA^{2}  CA^{3}],用matlab计算可知,Qo的秩为4,系统可观测。系统可设计观测器,并且观测器可控。

     

    其它博文链接:直线型一阶倒立摆1---概念篇

    有需要直线型一阶倒立摆的VREP仿真文件:可点击

  • 相关阅读:
    oracle PL/SQL(procedure language/SQL)程序设计之异常(exception)
    oracle PL/SQL(procedure language/SQL)程序设计之触发器(trigger)
    oracle PL/SQL(procedure language/SQL)程序设计之函数+过程+包
    oracle PL/SQL(procedure language/SQL)程序设计之游标cursors
    一个python拖库字段的小脚本
    perl6 拖库脚本
    python2 处理urllib/urllib2错误并打印源码
    perl6 修改文件并覆盖
    perl6文件操作
    python urllib2练习发送简单post
  • 原文地址:https://www.cnblogs.com/deciduousmap/p/12193724.html
Copyright © 2020-2023  润新知