zoukankan      html  css  js  c++  java
  • 直线型一阶倒立摆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仿真文件:可点击

  • 相关阅读:
    PAT 甲级 1004 Counting Leaves
    C#列的一些操作
    EXT.NET 一些用法
    ext.net 所有版本
    鱼线拉力值
    CRM 总是弹出登录窗口
    台钓线组及各种绑法
    windows server 2012 r2 安装IIS失败
    [iOS]UITableView刷新
    [Linux]APUE读书笔记: 文件和目录
  • 原文地址:https://www.cnblogs.com/deciduousmap/p/12193724.html
Copyright © 2011-2022 走看看