基于lmi的二级倒立摆的建模与仿真

《基于lmi的二级倒立摆的建模与仿真》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、基于LMI的二级倒立摆系统的H。鲁棒控制摘要倒立摆系统为典型的快速、多变量、非线性、绝对不稳定系统，且存在不确定因素。针对二级倒立摆系统中所受摩檫的不确定性，采用LMT方法，建立丫二级倒立摆模型，设计了H鲁棒控制器，给出了控制器的求解方法。仿真实验结果证明了该控制方法的有效性和可行性，并且具有很好的鲁棒稳定性和响应速度快的优越性，对高阶次不稳定系统具有很好的控制效果。关键词：二级倒立摆；线性矩阵不等式(LMI)鲁棒控制0引言现代控制工程所面临的问题极其复杂。实际的工程控制系统中，总是存在一定的不确定性。倒立摆即是一个包含不确定性的系统，也是控制理论的一个理想实验平台，对倒立摆系统

2、的研宄具有重要的理论和实际意义。本文采用线性矩阵不等式(LMI)方法，设计了二级倒立摆系统的鲁棒H状态反馈控制器，有效地克服了用求解两个联立的里卡辿方程获得H控制器时求解OO过程不容易收敛的困难，并且可降低控制器参数的数量级，使其在实控上易于实现。根据文献［1］中对LMI的处理方法，对二级倒立摆系统进行了仿真研宄，结果表明，这样的控制方法可使二级倒立摆系统具有很好的魯棒稳定性。1二级倒立摆系统建模1.1倒立摆系统结构图1是二级倒立摆的系统结构图，它由三部分组成：计算机、电气部分和机械部分。计算机部分有A/D、D/A转换模块，运动控制卡和PC机；电气部分主要有：光电编码器、直流功率

3、放大器、伺服电机和保护电路；机械部分有摆杆、轨道、运动小车和皮带轮等。图1二级倒立摆系统结构图1.2倒立摆系统特性分析倒立摆系统是典型的机械电子系统，具有如下特性：(1)欠冗余性。一般的倒立摆控制系统采用单电机驱动，无冗余结构。采用欠冗余的设计方法主要是在不失系统可靠性的前提下节约经济成本或节约有效的空间。(2)仿射非线性系统。倒立摆控制系统是一种典型的仿射非线性系统，可以用微分几何的方法进行分析。(3)不确定性。主要是指建立系统的数学模型时的参数误差、量测噪声以及机械传动过程中的非线性因素所导致的难以量化的部分。(4)耦合特性。倒立摆和小车之间，以及多级倒立摆的上下杆之间都是强

4、親合的。(5)开环不稳定性。倒立摆系统有两个平衡状态：竖直向下和竖直向上。竖直向下的状态是系统稳定的平衡点(考虑摩擦力的影响)，而竖直向上的状态是系统不稳定的平衡点，开环时微小的扰动都会使系统离开竖直向上的状态而进入竖直向下的状态。针对倒立摆系统的以上特性，在建模时，为了简单起见，一般忽略系统中些次要的、难以建模的一些因素，例如空气阻力、伺服电机由于安装而产生的静摩擦力、系统连接处的松弛程度、摆杆连接处质量分布的不均匀、传动皮带的弹性、传动齿轮的间隙等，将小车抽象为质点，摆杆抽象为匀质刚体，摆杆绕轴转动，这样可以通过力学原理建立系统较为精确的数学模型。为了研宄倒立摆系统的按制方法

5、，建立一个比较精确的倒立摆系统的线性模型是必不可少的。目前，对倒立摆系统建模一般采用两种方法：•十•顿力学分析方法和Lagrange方程法。1.3倒立摆系统数学建模为了简化二级倒立摆系统的数学模型，忽略空气流动作用在摆杆上的力矩干扰和上下摆转轴处的摩擦力矩，仅考虑小车与导轨的摩擦之后，可将倒立摆系统抽象成小车匀质杆和质量块组成，模型如图2所示。图中字母的意义和相关参数值如表1所示。图2二级倒立摆结构图表1二级倒立摆参数参数符号含义参数值M小车质量1.32kg上摆杆质量0.04kgm2下摆杆质量0.132kgm3质量块质量0.208kg下摆摆杆与垂直向上方向的夹角上摆摆杆与垂直向上

6、方向的夹角1,下摆摆杆转动中心到摆杆质心的距离0.09mb上摆摆杆转动中心到摆杆质心的距离0.27mF电机对小车的驱动力X小车位移颀flife1.3.1利用拉格朗日(Lagrange)方程推导动力学方程拉格朗日(Lagrange)方程为：(1)f,.为系统沿广义ddLdL_dt4机-1L为拉格朗H算子，q为系统的广义坐标，i-1,2,3...,n,坐标方向上的外力。L(q，q)=T(q，q)—V(q，q)(2)其中，T为系统的动能，V为系统的势能利用Lagrange方程，可以推导二级倒立摆的动力学方程，设系统的三个广义坐标分别是X，e2.第一步计算系统的动能：T=TM^Tm^Tm

7、2^Tm3(3)其中&，&，rm、，7；、分别表示小车，摆杆1,摆杆2,质量块的动能。而T,”'=Tm'+Tm'，其中7；是摆杆1质心平动动能，7：,是摆杆1绕质心转动动能。Tm、=H其中7；;是摆杆2质心平动动能，<2是摆杆2绕质心转动动能。2-2小车动能：1•2Tm=~Mx(4)摆杆1动能：=去'%4•卜^>m0)+J(Z,-Z^cos6,)]

8、=+去(5)乂因为=-J}w}2=-22m,(6)(7)1•2•參2•2Tnh=K,-m./.x^cos摆杆2动能:=r2d{x-