一级倒立摆的建模与控制分析

一级倒立摆的系统分析一、倒立摆系统的模型建立如图1-1所示为一级倒立摆的物理模型图1-1 一级倒立摆物理模型对于上图的物理模型我们做以下假设：M：小车质量m：摆杆质量b：小车摩擦系数l：摆杆转动轴心到杆质心的长度I：摆杆惯量F：加在小车上的力x：小车位置ɸ：摆杆与垂直向上方向的夹角θ：摆杆与垂直向下方向的夹角（考虑到摆杆初始位置为竖直向下）图1-2是系统中小车和摆杆的受力分析图。

其中，N和P为小车与摆杆相互作用力的水平和垂直方向的分量。

注意：实际倒立摆系统中的检测和执行装置的正负方向已经完全确定，因而矢量方向定义如图所示，图示方向为矢量正方向。

图1-2 小车及摆杆受力分析分析小车水平方向受力，可以得到以下方程：M ẍ＝Ｆ-ｂẋ-Ｎ (1-1)由摆杆水平方向的受力进行分析可以得到以下方程：N =md 2dt 2(x +l sin θ) (1-2)即： N =mẍ+mlθcos θ−mlθ2sin θ (1-3)将这个等式代入式（1-1）中，可以得到系统的第一个运动方程： (M +m )ẍ+bẋ+mlθcos θ−mlθ2sin θ=F (1-4)为推出系统的第二个运动方程，我们对摆杆垂直方向上的合力进行分析，可以得出以下方程： P −mg =md 2dt 2(l cos θ) (1-5)P −mg =− mlθsin θ−mlθ2cos θ (1-6) 利用力矩平衡方程可以有：−Pl sinθ−Nl cosθ=Iθ (1-7)注意：此方程中的力矩方向，由于θ=π+ɸ，cosɸ=−cosθ，sinɸ=−sinθ，所以等式前面含有负号。

合并两个方程，约去P和N可以得到第二个运动方程：(I+ml2)θ+mgl sinθ=−mlẍcosθ (1-8)设θ=π+ɸ，假设ɸ与1（单位是弧度）相比很小，即ɸ＜＜1，则可以进行近似处理：cosθ=−1，sinθ=−ɸ，（dθdt ）2=0。

用u来代表被控对象的输入力F，线性化后的两个运动方程如下：{(I+ml2)ɸ−mglɸ=mlẍ(M+m)ẍ+bẋ−mlɸ=u(1-9)假设初始条件为0，则对式（1-9）进行拉普拉斯变换，可以得到：{(I+ml2)Φ(s)s2−mglΦ(s)=mlX(s)s2(M+m)X(s)s2+bX(s)s−mlΦ(s)s2=U(s) (1-10) 由于输出为角度ɸ，求解方程组的第一个方程，可以得到：X(s)=[(I+ml2)ml −gs2]Φ(s) (1-11)或改写为：Φ(s)X(s)=mls2(I+ml2)s2−mgl(1-12)如果令v=ẍ,则有：Φ(s)V（s）=ml(I+ml2)s2−mgl(1-13)如果将上式代入方程组的第二个方程，可以得到：(M+m)[(I+ml2)ml −gs]Φ(s)s2+b[(I+ml2)ml+gs2]Φ(s)s−mlΦ(s)s2=U(s) (1-14) 整理后可得传递函数：Φ(s) U(s)=mlqs2s4+b(I+ml2)qs3−(M+m)mglqs2−bmglqs(1-15)其中q=[(M+m)(I+ml2)−(ml)2]假设系统状态空间方程为：X=AX+Buy=CX+Du (1-16) 方程组对ẍ,ɸ解代数方程，可以得到解如下：{ẋ=ẋẍ=−(I+ml2)bI(M+m)+Mml2ẋ+m2gl2I(M+m)+Mml2ɸ+(I+ml2)I(M+m)+Mml2uɸ=ɸɸ=−mlbI(M+m)+Mml2ẋ+mgl(M+m)I(M+m)+Mml2ɸ+mlI(M+m)+Mml2u(1-17)整理后可以得到系统状态空间方程：[ẋẍɸɸ]=[01000−(I+ml2)bI(M+m)+Mml2m2gl2I(M+m)+Mml200010−mlbI(M+m)+Mml2mgl(M+m)I(M+m)+Mml20][xẋɸɸ]+[(I+ml2)I(M+m)+Mml2mlI(M+m)+Mml2]uy=[xɸ]=[10000010][xẋɸɸ]+[0]u(1-18)由（1-9）的第一个方程为：(I+ml2)ɸ−mgl ɸ=mlẍ对于质量均匀分布的摆杆可以有：I=13ml2于是可以得到：(13ml2+ml2)ɸ−mgl ɸ=mlẍ化简可以得到：ɸ=3g4l ɸ+34lẍ(1-19)设X={x, ẋ, ɸ , ɸ}，u=ẍ则有：[ẋẍɸɸ]=[010000000001003g4l0][xẋɸɸ]+[134l]uy=[xɸ]=[10000010][xẋɸɸ]+[0]u(1-20)以上公式推理是根据牛顿力学的微分方程验证的。

直线一级倒立摆的建模及控制分析摘要：本文利用牛顿—欧拉方法，建立了直线型一级倒立摆系统的数学模型。

在分析的基础上, 采用状态反馈控制中极点配置法设计了用于直线型一级倒立摆系统的控制器。

此外，用MATLAB 仿真绘制了相应的曲线并做了分析。

一、问题描述倒立摆控制系统是机器人技术、控制理论、计算机控制等多个领域和多种技术的有机结合，其被控系统本身是一个绝对不稳定、高阶次、多变量、强耦合的非线性系统，是控制理论研究中较为理想的实验对象。

它为控制理论的教学、实验和科研构建了一个良好的实验平台，促进了控制系统新理论、新思想的发展。

倒立摆系统可以采用多种理论和方法来实现其稳定控制，如PID，自适应、状态反馈、智能控制等方法都己经在倒立摆控制系统上得到实现。

由于直线一级倒立摆的力学模型较简单，又是研究其他倒立摆的基础，所以本文利用所学的矩阵论知识对此倒立摆进行建模和控制分析。

二、方法简述本文利用牛顿—欧拉方法，建立了直线型一级倒立摆系统的数学模型。

在分析的基础上, 采用状态反馈控制中极点配置法设计了用于直线型一级倒立摆系统的控制器。

此外，用MATLAB 仿真绘制了相应的曲线并做了分析。

三、模型的建立及分析3.1 微分方程的推导在忽略了空气阻力，各种摩擦之后，可将直线一级倒立摆系统抽象成小车和匀质杆组成的系统，如下图1所示。

图1 直线一级倒立摆系统假设 M 为小车质量；m 为摆杆质量；b 为小车摩擦系数；l 为摆杆转动轴心到杆质心的长度；I 为摆杆惯量；F 为加在小车上的力；x 为小车位置；φ为摆杆与垂直向上方向的夹角；θ为摆杆与垂直向下方向的夹角。

图2是系统中小车和摆杆的受力分析图。

其中，N 和P 为小车与摆杆相互作用力的水平和垂直方向的分量。

值得注意的是: 在实际倒立摆系统中检测和执行装置的正负方向已确定, 因而矢量方向定义如图2所示, 图示方向为矢量正向。

(a) (b)图2 小车和摆杆的受力分析图分析小车水平方向所受的合力，可以得到以下方程：N x b F x M --= (1)由摆杆水平方向的受力进行分析可以得到下面等式：θθθθs i n c o s 2ml ml x m N -+= (2) 把这个等式代入上式中，就得到系统的第一个运动方程：()F ml ml x b x m M =-+++θθθθsin cos 2 (3)为了推出系统的第二个运动方程，我们对摆杆垂直方向上的合力进行分析，可以得到下面方程：θθθθc o s s i n 2 ml ml mg P --=- (4) 力矩平衡方程如下：θθθI Nl Pl =--cos sin (5)合并这(4)、(5)两个方程，约去P 和N ，得到第二个运动方程：()θθθc o s s i n 2x ml mgl ml I -=++ (6) 假设φ与1(单位是弧度)相比很小，即φ《1，则可以进行近似处理：0d d s i n 1c o s 2=⎪⎭⎫ ⎝⎛-=-=t θφθθ，， (7) 用u 来代表被控对象的输入力F ，线性化后两个运动方程如下：()()⎩⎨⎧=-++=-+u ml x b x m M xml mgl ml I φφφ 2 (8) 3.2 状态空间方程方程组(8)对φ，x 解代数方程，整理后的系统状态空间方程为： ()()()()()()()()u Mm l m M I m l Mm l m M I m lI x x Mm l m M I m M m gl Mm l m M I m lbMm l m M I gl m Mm l m M I b m l I x x ⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡++++++⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡+++++-+++++-=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡222222222200001000000010φφφφ u x x x y ⎥⎦⎤⎢⎣⎡+⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡=0001000001φφφ 对于质量均匀分布的摆杆有：3/2ml I =，于是可得：()x ml mgl ml ml =-+φφ223/ 化简得：xll g 4343+=φφ设}{x u x x X ==1，，，，φφ ，则有：14301004300100000000010u l x x l g x x⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡+⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡φφφφ10001000001u x x x y ⎥⎦⎤⎢⎣⎡+⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡=φφφ 3.3 实际系统模型实际系统模型参数: M =1.096 Kg ；m =0.109 Kg ；b =0.1 N/m/s ； l =0.25 m ；I =0.0034 kg ·m ·m ；采样频率 T =0.005 s 。

一阶倒立摆模糊控制实验报告本次实验旨在研究一阶倒立摆系统的模糊控制方法，通过对系统进行建模、设计控制器并进行仿真，最终评估控制效果。

实验过程主要包括系统建模、控制器设计、模糊控制器参数调节和性能评价四个步骤。

首先，我们对一阶倒立摆系统进行建模。

一阶倒立摆系统是一种具有非线性特性的控制系统，主要由电机、倒立摆、支撑杆等组成。

我们需要建立数学模型描述系统的动力学特性，包括倒立角度、倒立角速度、杆角度等状态变量，并考虑控制输入电压对系统的影响。

接着，我们设计模糊控制器。

模糊控制是一种基于模糊逻辑的控制方法，适用于非线性系统和模糊系统。

我们根据系统模型，设计模糊控制器的模糊规则、隶属函数等参数，以实现系统的稳定控制。

在设计过程中，我们需要考虑系统的性能指标，如超调量、稳态误差等。

第三步是模糊控制器参数调节。

通过仿真实验，我们可以对模糊控制器的参数进行调节，以使系统的性能达到最佳状态。

调节参数的过程需要考虑系统的稳定性、鲁棒性和响应速度，以达到控制效果的要求。

最后，我们对模糊控制系统进行性能评价。

通过对系统的响应曲线、稳定性、控制精度等指标进行分析，评价模糊控制器的控制效果。

我们可以比较模糊控制系统和传统控制系统的性能，探讨模糊控制在一阶倒立摆系统中的优势和局限性。

总的来说，本次实验通过研究一阶倒立摆系统的模糊控制方法，探讨了模糊控制在非线性系统中的应用。

通过实验，我们对模糊控制的基本原理和设计方法有了更深入的理解，同时也对一阶倒立摆系统的控制特性有了更清晰的认识。

希望通过实验的研究，能够为控制系统的设计和应用提供一定的参考和借鉴。

研究生《现代控制理论及其应用》课程小论文一级倒立摆的建模与控制分析学院：机械工程学院班级：机研131姓名：尹润丰学号： 2013212020162014年6月2日目录1. 问题描述及状态空间表达式建立............................ - 1 -1.1问题描述................................................................. - 1 -1.2状态空间表达式的建立..................................................... - 1 -1.2.1直线一级倒立摆的数学模型........................................... - 1 -1.2.2 直线一级倒立摆系统的状态方程 ...................................... - 5 -2.应用MATLAB分析系统性能.................................. - 6 -2.1直线一级倒立摆闭环系统稳定性分析......................................... - 6 -2.2 系统可控性分析.......................................................... - 7 -2.3 系统可观测性分析........................................................ - 8 -3. 应用matlab进行综合设计................................. - 8 -3.1状态反馈原理............................................................. - 8 -3.2全维状态反馈观测器和simulink仿真........................................ - 9 -4.应用Matlab进行系统最优控制设计......................... - 11 -5.总结.................................................... - 13 -1.问题描述及状态空间表达式建立1.1问题描述倒立摆是机器人技术、控制理论、计算机控制等多个领域、多种技术的有机结合，其被控系统本身又是一个绝对不稳定、高阶次、多变量、强耦合的非线性系统，可以作为一个典型的控制对象对其进行研究。