倒立摆系统的状态空间极点配置控制设计

《倒立摆系统的状态空间极点配置控制设计》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、倒立摆系统的状态空间极点配置控制设计电子与信息工程学院自动化094丛长龙摘要：为实现多输入、多输出、高度非线不稳定的倒立摆系统平衡稳定控制，将倒立摆系统的非线性模型进行近似线性化处理，获得系统在平衡点附近的线性化模型。利用牛顿—欧拉方法建立直线型一级倒立摆系统的数学模型。在分析的基础上，基于状态反馈控制中极点配置法对直线型倒立摆系统设计控制器。由MATLAB仿真表明采用的控制策略是有效的，设计的控制器对直线型一级倒立摆系统的平衡稳定性效果好，提高了系统的干扰能力。关键词：倒立摆、极点配置、MATLAB仿真引言：倒立摆是进行控制理论研究的典型试验平台，由于倒立摆本身所具有的高阶次、不稳定、

2、非线性和强耦合性，许多现代控制理论的研究人员一直将他视为典型的研究对象，不断从中发掘出新的控制策略和控制方法。控制器的设计是倒立摆系统的核心内容，因为倒立摆是一个绝对不稳定的系统，为使其保持稳定并且可以承受一定的干扰，基于极点配置法给直线型一级倒立摆系统设计控制器1．数学模型的建立倒立摆系统其本身是自不稳定的系统，实验建模存在着一定的困难。在忽略掉一些次要的因素之后，倒立摆系统就是一典型的运动的刚体系统，可以在惯性坐标系中应用经典力学理论建立系统动力学方程。下面采用牛顿-欧拉方法建立直线型一级倒立摆系统的数学模型。1.1微分方程的数学模型在忽略了空气阻力和各种摩擦力之后，可将直线一级倒立

3、摆系统抽象成小车和匀质杆组成的系统，如图1所示：12倒立摆系统的状态空间极点配置控制设计电子与信息工程学院自动化094丛长龙图1：直线一级倒立摆模型设系统的相关参数定义如下：M：小车质量m：摆杆质量b：小车摩擦系数l：摆杆转动轴心到杆质心的长度I：摆杆质量F：加在小车上的力x：小车位置Φ：摆杆与垂直方向上方向的夹角θ：摆杆与垂直方向下方向的夹角（摆杆的初始位置为竖直向下）如下图2所示为小车和摆杆的受力分析图。其中，N和P为小车与摆杆相互作用力的水平和垂直方向的分量。12倒立摆系统的状态空间极点配置控制设计电子与信息工程学院自动化094丛长龙图2：小车和摆杆受力分析图应用牛顿方法来建立系统

4、的动力学方程过程如下：分析小车水平方向所受的合力，可以得到以下的方程：由摆杆水平方向的受力进行分析可以得到下面的等式：将此等式代入上述等式中，可以得到系统的第一个运动方程：为了推出系统的第二个运动方程，我们对摆杆垂直方向上的合力进行分析，可以得到下面的方程：力矩平衡方程如下：注意：此方程中力矩的方向，由于故等式前面有负号。合并这两个方程，约去P和N，得到第二个运动方程：设θ=π+φ，当摆杆与垂直向上方向之间的夹角φ与1（单位是弧度）相比很小时，即Φ<<1时，则可以进行如下近似处理：12倒立摆系统的状态空间极点配置控制设计电子与信息工程学院自动化094丛长龙线性化后得到该系统数学模型的微分

5、方程表达式：1.2状态空间数学模型控制系统的状态空间方程可写成如下形式：解代数方程可得如下解：整理后可得系统的状态空间方程：12倒立摆系统的状态空间极点配置控制设计电子与信息工程学院自动化094丛长龙对于质量均匀分布的摆杆，其转动惯量为：代入微分方程模型中得：化简后可得：设则有：12倒立摆系统的状态空间极点配置控制设计电子与信息工程学院自动化094丛长龙1.3实际系统模型实际系统参数如下：M：小车质量，0.5Kg;m：摆杆质量，0.2Kg;b：小车摩擦系数，0.1N/m/sec；l：摆杆转动轴心到杆质心的长度，0.3m;I：摆杆质量，0.006Kg·m·m；T：采样时间，0.005s。将

6、上述系统参数代入可得系统实际模型。摆杆角度和小车位移的传递函数：摆杆角度和小车加速度之间的传递函数：摆杆角度和小车所受外界作用力的传递函数：12倒立摆系统的状态空间极点配置控制设计电子与信息工程学院自动化094丛长龙以外界作用力作为输入的系统状态方程：以小车加速度作为输入的系统状态方程：2.状态空间极点配置12倒立摆系统的状态空间极点配置控制设计电子与信息工程学院自动化094丛长龙经典控制理论的研究对象主要是单输入单输出的系统，控制器设计时一般需要有关被控对象的较精确模型，现代控制理论主要是依据现代数学工具，将经典控制理论的概念扩展到多输入多输出系统。极点配置法通过设计状态反馈控制器将多

7、变量系统的闭环系统极点配置在期望的位置上，从而使系统满足瞬态和稳态性能指标。设计要求：用极点配置方法设计控制器，使得在小车上施加0.1N的阶跃信号时，闭环系统的响应指标为：（1）要求系统调整时间小于3s（2）稳态时摆杆与垂直方向的夹角变化小于0.1弧度状态方程为：选择控制信号：可解得：直接利用MATLAB极点配置函数[K，PREC，MESSAGE]=PLACE(A，B，P)来计算。选取调整时间ts=2.0s，阻尼比为ξ=0.5，可得