LQR的一倒立摆仿真研究

综

合

性

实

验

设

计

题目：直线一级倒立摆

班级：０８０２

姓名：郭长春

指导老师：张白莉学号：200807211065摘要

倒立摆是一个复杂的多变量强祸合不稳定非线性的系统，借助于这样一个系统可以有效的检测各种控制理论的好坏。本论文研究了倒立摆系统的建模方法，采用Lagrange方程对一级倒立摆系统进行建模。针对倒立摆的状态空间方程，采用MatLab软件对倒立摆的稳定性，能控性，能观性和可控度做了分析;采用最优控制理论对倒立摆系统进行稳摆控制并进行仿真，并利用该理论将倒立摆的数学模型加以改进，并对改进后的倒立。经过仿真分析确认改进后的倒立摆系统控制效果得到了明显的改善。

正文

1．倒立摆系统研究的背景

随着各种控制理论的发展，理论的正确性，它在实际应用中的可行性，都需要用一个典型的对象来验证，这种控制对象必须能够比较各种控制理论的优劣。倒立摆系统就是我们需要的理想的实验控制平台。倒立摆控制系统，是一个复杂的多变量强祸合不稳定非线性系统。对倒立摆系统进行研究，能有效的反映控制策略中的许多典型问题，，针对倒立摆而研究的控制方法在双足机器人的行走平衡方面、军工方面、火箭发射过程的垂直度控制、卫星在飞行中的姿态控制以及日常生产过程和生活领域中的应用都非常广泛。因此，倒立摆系统的研究，有着重要的工业背景和实际意义。

2．直线一级倒立摆系统的建模

直线一级倒立摆的工作原理可简述为：用一种强有力的控制方法使小车以一定的规律来回运动，从而使摆杆在垂平面内稳定下来。这样的系统就是倒立摆的控制系统如下图所示若小车不动，摆杆会由于重力作用倒下来若在水平方向给小车一个力，则摆杆受到一个力矩，这个力矩使摆杆朝与小车相反的方向运动，通过规律性的改变小车的受力方向使得摆杆在竖直方向左右摆动从而实现摆杆在竖直方向的动态平衡。在忽略了空气阻力和各种摩擦之后，可将直线一级倒立摆系统抽象成小车和匀质杆组成的系统。如下如所示。

图倒立摆系统

其中：

M 小车的质量

m 摆杆的质量

b 小车的摩擦系数

l 摆杆转动轴心到杆质心的长度

I 摆杆的惯量

F 加在小车上的力

X 小车的位置

φ摆杆与垂直向上方向的夹角

θ摆杆与垂直向下方向的夹角（考虑到摆杆初始位置为竖直向下）

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

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

分析小车水平方向所受的合力，可以得到以下方程：

-

=

M-

F

N

x b

x

由摆杆水平方向的受力进行分析可以得到下面等式：

)sin (22

θl x dt

d m N += 即：

θθθθsin cos 2 ml ml x

m N -+= 把这个等式代入式(3-1)中，就得到系统的第一个运动方程：

F ml ml x b x

m M =-+++θθθθsin cos )(2 为了推出系统的第二个运动方程，我们对摆杆垂直方向上的合力进行分析，可以得到下面方程：

)cos (22

θl dt

d m mg P =- θθθθ

cos sin 2 ml ml mg P --=- 力矩平衡方程如下：

θ

θθ I Nl Pl =--cos sin 注意：此方程中力矩的方向，由于θφθφφπθsin sin ,cos cos ,-=-=+=，故等式前面有负号。

合并这两个方程，约去P 和N ，得到第二个运动方程：

θθθcos sin )(2x

ml mgl ml I -=++ 设φπθ+=（φ是摆杆与垂直向上方向之间的夹角），假设φ与1（单位是弧度）相比很小，即1<<φ，则可以进行近似处理：0)(

,sin ,1cos 2=-=-=dt

d θφθθ。用u 来代表被控对象的输入力F ，线性化后两个运动方程如下： 2(+)()I ml mgl mlx M m x bx ml u

????-=?++-=? 对式(3-9)进行拉普拉斯变换，得到

2111011(+)()J m L mgl m Lr m m r r m L u

????-=?++-=? 3.传递函数模型

推导传递函数时假设初始条件为0。

由于输出为角度φ，求解方程组的第一个方程，可以得到：

21121()()[]()J m L g R s s m L s

+=-Φ

或 2

122111()()()s m Ls R s J m L s m gL

Φ=+- （18-10） 如果令x

v =，则有： 122111()()()s m L V s J m L s m gL

Φ=+- 把上式代入方程组的第二个方程，得到：

222111101211()()()()()()J m L g J m L g m m s s ml s s U s m L

s m L s ????+++-++Φ-Φ=???????? 整理后得到传递函数：

212432110111()()()()m L s s q J m L m m m gL m gL U s s s s s q q q

Φ=+++-- （18-11） 其中 2201111()()()q m m J m L m L ??=++-??

4.状态空间方程模型

设系统状态空间方程为：

Du CX y Bu AX X

+=+=

方程组 对φ ,x

解代数方程，得到解如下： 222211111222

1001101011010111011222101011010110110()()()()()()()()()x x J m L m gL J m L x x u J m m m m L J m m m m L J m m m m L m L m gL m m m L x u J m m m m L J m m m m L J m m m m L ???

??=??-++?=++?++++++?=??-+=++?++++++? 整理后得到系统状态空间方程：

222211111222101011010110101110112210101101011001000()0

0()()()00010()00()()(x x J m L m gL J m L x x J m m m m L J m m m m L J m m m m L m L

m gL m m m L J m m m m L J m m m m L J m ????????????-++????????++++++????=+??????????????-+??????++++??2101)u m m m L ????????????????++??

u x x x y ??????+?????

?????????????=??????=0001000001φφφ （18-13） 5.系统的物理参数

代入倒立摆系统的参数。

小车的质量：M=1kg

摆杆的质量：m=0.1kg

小车的摩擦系数：b=0.1N/(m ·1-s )

摆杆转动轴心到杆质心的长度：l =0.5m

重力加速度：g=10/2

s

摆杆惯量：I=0.003kg ·2m

采样频率：0.005s.

将上述数据带入（18-10）、（18-11）、（18-13） 摆杆角度与小车位移之间的传递函数：5

.0028..005.0)()(22

-=s s s x s φ (18-14) 摆杆角度与小车所受外力之间的传递函数：

)()(s u s φ=s

s s s s 78253.16078.19099821.078253.12342

--+ (18-15) 系统的空间状态表达式为：

u x x x x ?????

???????+????????????????????????--=??????

????????766784452.1098939929.000434628975.191766784.0010000883392.00989399.000010φφφφ （18-15） u x x x y ??????+??????????????????=??????=0001000001φφφ 6.经典控制器设计——PID 控制

PID 控制器因其结构简单，容易调节，且不需要对系统建立精确的模型，在控制上应用较广

PID 的控制器传递函数为：

PID

PID I P D I P D PID den num S K S K S K S K K S K S G =++=++=2)( （18-17） close all

clear all

%设定系统参数

M=1;

m=0.1;

b=0.1;

l=0.5;

I=0.003;

g=9.8;

q=(M+m)*(I+m*l^2)-(m*l)^2

%输入倒立摆的传递函数

num=[m*l/q 0 0];

den=[l b*(I+m*l^2)/q -(M+m)*m*g*l/q -b*m*g*l/q 0];

%计算并显示系统未校正的极点

[r,p,k]=residue(num,den);

disp('the poles are:');

p

%设定PID 参数

kp=40;

ki=1;

kd=10;

numPID=[kd kp ki];

denPID=[1 0];

%计算经过PID 校正后的闭环传递函数

numc=conv(num,denPID);

denc=polyadd(conv(den,denPID),conv(num,numPID));

%j 计算并显示校正后系统的极点

[Rc,Pc,Kc]=residue(num,denc);

disp('the poles are:');

Pc

%纺真原系统在干扰脉冲下的输出曲线图

t=0:0.005:5;

yc=impulse(numc,denc,t);

plot(t,yc);

grid on

title('mpulse Response');

xlabel('Time-sec');

ylabel('Amplitude');

文件中用到两个多项式之和的函数polyadd ，他不是MATLAB 工具因此必须把他考到polyadd.m 文件中并把该文件保存到MATLAB 安装路径下的work 目录下, polyadd.m:

function[poly]=polyadd(poly1,poly2)

if length(poly1)

short=poly1;

long=poly2;

else

short=poly2;

long=poly1;

end

mz=length(long)-length(short);

if mz>0

poly=[zeros(1,mz),short]+long;

else

poly=long+short;

end

运行系统得到闭环极点为:

p =

-6.2264

6.1194

-0.0909

系统位于s平面右半部的极点可见是不稳定的.矫正后闭环极点为: Pc =

-32.3411

-3.1918

-0.0002

-0.0002

-0.0002

矫正后系统稳定.矫正后的单位脉冲响应如下图.

参考文献：

[1}孙灵芳，孔辉，刘长国等.倒立摆系统及研究现状.机床与液压.2008,36(07):306-309.

[2]罗晶，陈平.一阶倒立摆的PID控制.实验室研究与探索.2005, 24(11) :26^28.

[3]杨平，徐春梅，贺茂康等.直线二级倒立摆的PID实时控制.上海电力学院学报.2008,24(03):236-247.

[4〕吴昊，秦志，朱学峰等.直线柔性连接两级倒立摆控制器设计.控制工程.2003,10(06):497 ^-499.

[5〕杨平，徐春梅，王欢等.直线型一级倒立摆状态反馈控制设计及实现.上海电力学院学报，.2007,23(01):21-25.

[6〕袁性忠，姜新建.基于滑模变结构的倒立摆系统稳定控制.控制理论与应用.2004,21(05):720-723.

[7]姚利娜，王宏，周靖林等.倒立摆系统的变结构控制方案.控制理论与应用.2004,21(05):724-727.

[8]周文彬，朱学峰，陈华艳。基于留优遗传算法的倒立摆最优LQR控制.控制工程.2004, 11 (S,) :160^-163.

[10] Hiroshi Wakuya,Yoichiro Sankai，Yuzuru Morita.Seasonal variation of neuro-controller

training period

一阶倒立摆控制系统

一阶直线倒立摆系统 姓名： 班级： 学号：

目录 摘要 (3) 第一部分单阶倒立摆系统建模 (4) （一）对象模型 (4) （二）电动机、驱动器及机械传动装置的模型 (6) 第二部分单阶倒立摆系统分析 (7) 第三部分单阶倒立摆系统控制 (11) （一）内环控制器的设计 (11) （二）外环控制器的设计 (14) 第四部分单阶倒立摆系统仿真结果 (16) 系统的simulink仿真 (16)

摘要： 该问题源自对于娱乐型”独轮自行车机器人”的控制，实验中对该系统进行系统仿真，通过对该实物模型的理论分析与实物仿真实验研究，有助于实现对独轮自行车机器人的有效控制。 控制理论中把此问题归结为“一阶直线倒立摆控制问题”。另外，诸如机器人行走过程中的平衡控制、火箭发射中的垂直度控制、卫星飞行中的姿态控制、海上钻井平台的稳定控制、飞机安全着陆控制等均涉及到倒立摆的控制问题。 实验中通过检测小车位置与摆杆的摆动角，来适当控制驱动电动机拖动力的大小，控制器由一台工业控制计算机（IPC）完成。实验将借助于“Simulink封装技术——子系统”，在模型验证的基础上，采用双闭环PID控制方案，实现倒立摆位置伺服控制的数字仿真实验。实验过程涉及对系统的建模、对系统的分析以及对系统的控制等步骤，最终得出实验结果。仿真实验结果不仅证明了PID方案对系统平衡控制的有效性，同时也展示了它们的控制品质和特性。 第一部分单阶倒立摆系统建模

（一） 对象模型 由于此问题为”单一刚性铰链、两自由度动力学问题”，因此，依据经典力学的牛顿定律即可满足要求。 如图1.1所示，设小车的质量为0m ，倒立摆均匀杆的质量为m ，摆长为2l ，摆的偏角为θ，小车的位移为x ，作用在小车上的水平方向上的力为F ，1O 为摆杆的质心。 图1.1 一阶倒立摆的物理模型 根据刚体绕定轴转动的动力学微分方程，转动惯量与角加速度乘积等于作用于刚体主动力对该轴力矩的代数和，则 1）摆杆绕其重心的转动方程为 sin cos y x l F J F l θθθ=-&& （1-1） 2）摆杆重心的水平运动可描述为 2 2(sin )x d F m x l dt θ=+ （1-2） 3）摆杆重心在垂直方向上的运动可描述为 2 2(cos )y d F mg m l dt θ-= （1-3） 4）小车水平方向运动可描述为 202x d x F F m dt -= （1-4）

倒立摆系统的建模及Matlab仿真资料

第1 页共11 页 倒立摆系统的建模及Matlab仿真 1.系统的物理模型 考虑如图(1)所示的倒立摆系统。图中，倒立摆安装在一个小车上。这里仅考虑倒立摆在图面内运动的二维问题。 图(1)倒立摆系统 假定倒立摆系统的参数如下。 摆杆的质量：m=0.1g l=1m小车的质量：摆杆的长度：2重力加速度：g=9.8m/M=1kg s摆杆的质量在摆杆的中心。 设计一个控制系统，使得当给定任意初始条件(由干扰引起)时，最大超调量?≤10%，调节时间ts ≤4s ，通过小车的水平运动使倒立摆保持在垂直位置。 2.系统的数学模型 2.1建立倒置摆的运动方程并将其线性化。 为简化问题，在数学模型中首先假设：1)摆杆为刚体；2）忽略摆杆与支点之间的摩擦；3）忽略小车与接触面间的摩擦。 ?）,在u设小车瞬时位置为z,摆心瞬时位置为（作用下，小车及摆均产生加速远 动，sin?lz根据牛顿第二定律，在水平直线远动方向的惯性力应与u平衡，于是有 22dzd?)?sinu?M?m(zl22dtdt???2????z(M?mml?)cos?mlusin? 即：??①

绕摆轴转动的惯性力矩与重力矩平衡，因而有． 第2 页共11 页 2??d??? sin??lcosm(z?lsinmgl)??2dt?????22???????即： nis?l?ocgcosincoszs?ls??② 以上两个方程都是非线性方程，为求得解析解，需作线性化处理。由于控制的目的是保持倒立摆直?2?????且可忽略则，立，在试驾合适的外力条件下，假定θ很小，接近于零时合理的，1sincos??,项。于是有 ???M?zm?u?ml??)(③ ????g?z?l??④联立求解可得1mg?u?z????MM 1)?m(M????u??MlMl 列写系统的状态空间表达式。2.2??T xx,x,x,，选取系统变量则 xx,x,xx?,42134123xx??211mgux???x?32MM x?x?431)(M?mu?x?x? 34MlMl 即00100????z??1mg??????000?z?????d MM??Bu?Ax?xux????????00001???dt????1gm?(M)????000??????? MlMl??????Cx?0?y?xx1001代入数据计算得到：0100????000?1??????T0D,?0??1BA?,?001,C100??1000??00011?? 11 页3 页共第 3.设计控制器3.1判断系统的能控性和稳定性 1100????0011????23BBAABAB?Q?故被控对象完全可控， rank()=4,Q kk??11?0?10??011?10???22???11?。出现大于零的特征值，故被，，0 解得特征值为 0由特征方程0??11I?A?)(控对象不稳定3.2确定希望的极点， 另一对为远极点，认为系统性能主要由主导，选其中一对为主导极点和希望的极点n=4ss21极点决定，远极点只有微小影响。根据二阶系统的关系式，先确定主导极点???42??1????10.?e??t1.67?有，闭环可得；取误差带，于是取，则6.?059?0.02.?0? pns??n2????1?js??=-10.8j,远极点选择使它和原点的距离大于主导极点与原点 距离主导极点为?n,21s??15倍，取的54,33.3采用状态反馈方法使系统稳定并配置极点 ??kkkk?k；状态反馈系统的状态方程，馈状态反的控制规律为为kxu??3102?，其

(完整版)一级倒立摆系统分析

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

杆相互作用力的水平和垂直方向的分量。注意：实际倒立摆系统中的检测和执行装置的正负方向已经完全确定，因而矢量方向定义如图所示，图示方向为矢量正方向。 图1-2 小车及摆杆受力分析 分析小车水平方向受力，可以得到以下方程： M x?＝Ｆ-ｂx?-Ｎ (1-1) 由摆杆水平方向的受力进行分析可以得到以下方程： N =m d 2dt (x +l sin θ) (1-2) 即： N =mx?+mlθcos θ?mlθ2sin θ (1-3) 将这个等式代入式（1-1）中，可以得到系统的第一个运动方程： (M +m )x?+bx?+mlθcos θ?mlθ2sin θ=F (1-4) 为推出系统的第二个运动方程，我们对摆杆垂直方向上的合力进行分析，可以得出以下方程： P ?mg =m d 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θ=?mlx?cosθ (1-8) 设θ=π+?，假设?与1（单位是弧度）相比很小，即?＜＜1，则 可以进行近似处理：cosθ=?1，sinθ=??，（dθ dt ） 2 =0。用u来 代表被控对象的输入力F，线性化后的两个运动方程如下： {(I+ml2)??mgl?=mlx? (M+m)x?+bx??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 ?g s ]Φ(s) (1-11) 或改写为：Φ(s) X(s)=mls2 (I+ml2)s2?mgl (1-12) 如果令v=x?,则有：Φ(s) V（s）=ml (I+ml2)s2?mgl (1-13) 如果将上式代入方程组的第二个方程，可以得到： (M+m)[(I+ml2) ml ?g s ]Φ(s)s2+b[(I+ml2) ml +g s ]Φ(s)s?mlΦ(s)s2= U(s) (1-14) 整理后可得传递函数： Φ(s) U(s)= ml q s2 s4+b(I+ml 2) q s3?(M+m)mgl q s2?bmgl q s (1-15)

最优化方法课程设计实验报告_倒立摆

倒立摆控制系统控制器设计实验报告

成员：陈乾睿 2220150423 郑文 2220150493 学院：自动化 倒立摆控制系统控制器设计实验 一、实验目的和要求 1、目的 （1）通过本设计实验，加强对经典控制方法（LQR控制器、PID控制器）和智能控制方法（神经网络、模糊控制、遗传算法等）在实际控制系统中的应用研究。（2）提高学生有关控制系统控制器的程序设计、仿真和实际运行能力. （3）熟悉MATLAB语言以及在控制系统设计中的应用。 2、要求 （1）完成倒立摆控制系统的开环系统仿真、控制器的设计与仿真以及实际运行结果 （2）认真理解设计内容，独立完成实验报告，实验报告要求：设计题目，设计的具体内容及实验运行结果，实验结果分析、个人收获和不足，参考资料。程序

清单文件。 二、实验内容 倒立摆控制系统是一个典型的非线性系统，其执行机构具有很多非线性，包括：死区、电机和带轮的传动非线性等。 本设计实验的主要内容是设计一个稳定的控制系统，其核心是设计控制器，并在MATLAB/SIMULINK环境下进行仿真实验，并在倒立摆控制实验平台上实际验证。 算法要求：使用LQR以外的其它控制算法。 三、倒立摆系统介绍 倒立摆是机器人技术、控制理论、计算机控制等多个领域、多种技术的有机结合，其被控系统本身又是一个绝对不稳定、高阶次、多变量、强耦合的非线性系统，可以作为一个典型的控制对象对其进行研究。倒立摆系统作为控制理论研究中的一种比较理想的实验手段，为自动控制理论的教学、实验和科研构建一个良好的实验平台，以用来检验某种控制理论或方法的典型方案，促进了控制系统新理论、新思想的发展。由于控制理论的广泛应用，由此系统研究产生的方法和技术将在半导体及精密仪器加工、机器人控制技术、人工智能、导弹拦截控制系统、航空对接控制技术、火箭发射中的垂直度控制、卫星飞行中的姿态控制和一般工业应用等方面具有广阔的应用开发前景。 倒立摆的形式和结构各异，但所有的倒立摆都具有以下的特性：非线性，不确定性，耦合性，开环不稳定性，约束限制。 经过相关论文和文献的查询，我们决定采用模糊控制的方法进行倒立摆的控制。

一级倒立摆地Simulink仿真

单级倒立摆稳定控制 直线一级倒立摆系统在忽略了空气阻力及各种摩擦之后，可抽象成小车和匀质摆杆组成的系统，如图1所示。 图1 直线一级倒立摆系统 图2 控制系统结构 假设小车质量M =0.5kg ，匀质摆杆质量m=0.2kg ，摆杆长度2l =0.6m ，x (t )为小车的水平位移，θ为摆杆的角位移，2 /8.9s m g =。控制的目标是通过外力u (t)使得摆直立向上（即0)(=t θ）。该系统的非线性模型为： u ml x m M ml mgl x ml ml J +=++=++22)sin ()()cos (sin )cos ()(θθθθθθθ ，其中231ml J =。 解： 一、 非线性模型线性化及建立状态空间模型 因为在工作点附近（0,0==θ θ ）对系统进行线性化，所以 可以做如下线性化处理：32 sin ,cos 13!2!θθθθθ≈-≈-

当θ很小时，由cos θ、sin θ的幂级数展开式可知，忽略高次项后， 可得cos θ≈1，sin θ≈θ，θ’^2≈0； 因此模型线性化后如下： （J+ml^2）θ’’+mlx ’’=mgl θ (a) ml θ’’+(M+m) x ’’=u (b) 其中23 1ml J = 取系统的状态变量为,,,,4321θθ ====x x x x x x 输出T x y ][θ=包括小车位移和摆杆的角位移. 即X=????????????4321x x x x =????? ???????''θθx x Y=??????θx =??????31x x 由线性化后运动方程组得 X1’=x ’=x2 x2’=x ’’=m m M mg 3)(43-+-x3+m m M 3)(44-+u X3’ =θ’=x4 x4’=θ’’=ml l m M g m M 3)(4)(3-++x3+ml l m M 3)(43-+-u 故空间状态方程如下： X ’=????????????'4'3'2'1x x x x =????????????????? ?-++-+-03)(4)(300100003)(4300 0010ml l m M g m M m m M mg ????????????4321x x x x + ???????? ??????????-+--+ml l m M m m M 3)(4303)(440 u

倒立摆实验报告

倒立摆实验报告 机自82 组员：李宗泽 李航 刘凯 付荣

倒立摆与自动控制原理实验 一．实验目的: 1.运用经典控制理论控制直线一级倒立摆,包括实际系统模型的建立、根轨迹分析和控制器设计、频率响应分析、PID 控制分析等内容. 2.运用现代控制理论中的线性最优控制LQR 方法实验控制倒立摆 3.学习运用模糊控制理论控制倒立摆系统 4.学习MATLAB工具软件在控制工程中的应用 5.掌握对实际系统进行建模的方法，熟悉利用MATLAB 对系统模型进行仿真，利用学习的控制理论对系统进行控制器的设计，并对系统进行实际控制实验，对实验结果进行观察和分析，非常直观的感受控制器的控制作用。 二. 实验设备 计算机及等相关软件 固高倒立摆系统的软件 固高一级直线倒立摆系统,包括运动卡和倒立摆实物 倒立摆相关安装工具 三．倒立摆系统介绍 倒立摆是机器人技术、控制理论、计算机控制等多个领域、多种

技术的有机结合，其被控系统本身又是一个绝对不稳定、高阶次、多变量、强耦合的非线性系统，可以作为一个典型的控制对象对其进行研究。倒立摆系统作为控制理论研究中的一种比较理想的实验手段，为自动控制理论的教学、实验和科研构建一个良好的实验平台，以用来检验某种控制理论或方法的典型方案，促进了控制系统新理论、新思想的发展。由于控制理论的广泛应用，由此系统研究产生的方法和技术将在半导体及精密仪器加工、机器人控制技术、人工智能、导弹拦截控制系统、航空对接控制技术、火箭发射中的垂直度控制、卫星飞行中的姿态控制和一般工业应用等方面具有广阔的利用开发前景。 倒立摆已经由原来的直线一级倒立摆扩展出很多种类，典型的有直线倒立摆环形倒立摆，平面倒立摆和复合倒立摆等，本次实验采用的是直线一级倒立摆。 倒立摆的形式和结构各异，但所有的倒立摆都具有以下的特性: 1) 非线性2) 不确定性3) 耦合性4) 开环不稳定性5) 约束限制 倒立摆控制器的设计是倒立摆系统的核心内容，因为倒立摆是一个绝对不稳定的系统，为使其保持稳定并且可以承受一定的干扰，需要给系统设计控制器，本小组采用的控制方法有：PID 控制、双PID 控制、LQR控制、模糊PID控制、纯模糊控制 四．直线一级倒立摆的物理模型： 系统建模可以分为两种：机理建模和实验建模。实验建模就是通过在研究对象上加上一系列的研究者事先确定的输入信号，激励

哈工大一阶倒立摆

哈尔滨工业大学 控制科学与工程系 控制系统设计课程设计报告

姓名：院（系）： 专业：自动化班号： 任务起至日期： 2014 年9 月9 日至 2014 年9 月20 日 课程设计题目：直线一级倒立摆控制器设计 已知技术参数和设计要求： 本课程设计的被控对象采用固高公司的直线一级倒立摆系统GIP-100-L。 系统内部各相关参数为: M小车质量0.5kg; m摆杆质量0.2kg; b小车摩擦系数0.1N/m/sec; l摆杆转动轴心到杆质心的长度0.3m; I摆杆惯量0.006kg*m*m; T采样时间0.005秒。 设计要求： 1.推导出系统的传递函数和状态空间方程。用Matlab进行阶跃输入仿真，验证系统的稳定性。 2.设计PID控制器，使得当在小车上施加0.1N的脉冲信号时，闭环系统的响应指标为： （1）稳定时间小于5秒； （2）稳态时摆杆与垂直方向的夹角变化小于0.1弧度。 3.设计状态空间极点配置控制器，使得当在小车上施加0.2m的阶跃信号时，闭环系统的响应指标为： （1）摆杆角度错误！未找到引用源。和小车位移x的稳定时间小于3秒 （2）x的上升时间小于1秒 （3）错误！未找到引用源。的超调量小于20度（0.35弧度） （4）稳态误差小于2%。 工作量： 1.建立直线一级倒立摆的线性化数学模型； 2.倒立摆系统的PID控制器设计、Matlab仿真及实物调试； 3.倒立摆系统的极点配置控制器设计、Matlab仿真及实物调试。

哈尔滨工业大学 (1) 控制系统设计课程设计报告 (1) 一．实验设备简介 (3) 二．直线一阶倒立摆数学模型的推导 (6) 2.1概述 (6) 2.2数学模型的建立 (7) 2.3一阶倒立摆的状态空间模型： (9) 2.4实际参数代入： (10) 三．定量、定性分析系统的性能 (11) 3.1 对系统的稳定性进行分析 (11) 3.2 对系统的稳定性进行分析： (12) 四. 实际系统的传递函数与状态方程 (13) 五. 系统阶跃响应分析 (14) 六．一阶倒立摆PID控制器设计 (15) 6.1 PID控制分析 (15) 6.2 PID控制参数设定及MATLAB仿真 (17) 6.3 PID控制实验 (18) 七．状态空间极点配置控制器设计 (19) 7.1 状态空间分析 (20) 7.2 极点配置及MA TLAB仿真 (21) 7.3 利用爱克曼公式计算 (21) 八．课程设计心得与体会 (22) 一．实验设备简介 倒立摆控制系统：Inverted Pendulum System （IPS） 倒立摆控制系统是一个复杂的、不稳定的、非线性系统，是进行控制理论教学及开展各种控制实验的理想实验平台。对倒立摆系统的研究能有效的反映控制中的许多典型问题：如非线性问题、鲁棒性问题、镇定问题、随动问题以及跟踪问题等。通过对倒立摆的控制，用来检验新的控制方法是否有较强的处理非线性和不稳定性问题的能力。同时，其控制方法在军工、航天、机器人和一般工业过程领域中都有着广泛的用途，如机器人行走过程中的平衡控制、火箭发射中的垂直度控制和卫星飞行中的姿态控制等。 倒立摆是进行控制理论研究的典型实验平台。倒立摆是机器人技术、控制理论、计算机控制等多个领域、多种技术的有机结合，其被控系统本身又是一个绝对不稳定、高阶次、多变量、强耦合的非线性系统，可以作为一个典型的控制对象对其进行研究。最初研究开始于二十世纪50 年代，麻省理工学院（MIT）的控制论专家根据火箭发射助推器原理设计出一级倒立摆实验设备。近年来，新的控制方法不断出现，人们试图通过倒立摆这样一个典型的控制对象，检验新的控制方法是否有较强的处理多变量、非线性和绝对不稳定系统的能力，从而从中找出最优秀的控制方法。

倒立摆系统的建模及Matlab仿真

倒立摆系统的建模及Matlab 仿真 1.系统的物理模型 考虑如图(1)面内运动的二维问题。 图(1)倒立摆系统 假定倒立摆系统的参数如下。 摆杆的质量：m=0.1g 摆杆的长度：l =1m 小车的质量： M=1kg 重力加速度：g=9.8m/2s 摆杆的质量在摆杆的中心。 设计一个控制系统，使得当给定任意初始条件(由干扰引起)时，最大超调量δ ≤10%，调节时 间ts ≤4s ，通过小车的水平运动使倒立摆保持在垂直位置。 2.系统的数学模型 2.1建立倒置摆的运动方程并将其线性化。 为简化问题，在数学模型中首先假设：1)摆杆为刚体；2）忽略摆杆与支点之间的摩擦；3）忽略小车与接触面间的摩擦。 设小车瞬时位置为z,摆心瞬时位置为（θsin l z +）,在u 作用下，小车及摆均产生加速远动，根据牛顿第二定律，在水平直线远动方向的惯性力应与u 平衡，于是有 u l z dt d m dt z d M =++)sin (22 22θ 即： u ml ml z m M =-++θθθθsin cos )(2&&&&& ① 绕摆轴转动的惯性力矩与重力矩平衡，因而有

θθθsin cos )sin (22mgl l l z dt d m =??? ????+ 即： θθθθθθθsin cos sin cos cos 22g l l z =-+&&&&& ② 以上两个方程都是非线性方程，为求得解析解，需作线性化处理。由于控制的目的是保持倒立摆直 立，在试驾合适的外力条件下，假定θ很小，接近于零时合理的，则1cos ,sin ≈≈θθθ，且可忽略θ θ2&项。于是有 u ml z m M =++θ&&&& )( ③ θθg l z =+&&&& ④ 联立求解可得 u Ml Ml m M u M M mg z 1)(1 -+=+- =θθθ&&&& 2.2列写系统的状态空间表达式。 选取系统变量4321,,,x x x x ， []T x x x x x 4321,,,=则 u Ml x Ml m M x x x u M x M mg x x x 1 )(134433221-+= =+-==&&&& 即 []Cx x x y Bu Ax u Ml M x Ml g m M M mg z z dt d x ===+=?????? ? ???????-+?????????? ??? ? +- =???? ????????=000110100)(0 010 0000000 1 1θθ&&& 代入数据计算得到： [][]0,0001,1010,01100 1000010000 1 0==-=? ? ??? ? ??? ???-=D C B A T

倒立摆实验报告

目录 一、倒立摆系统介绍 (2) 1.1倒立摆系统简介 (2) 1.2 倒立摆组成及其原理 (2) 1.3 倒立摆特性 (3) 二、一级倒立摆 (3) 2.1一级倒立摆建模 (3) 2.2 一级倒立摆控制方法 (11) 2.2.1 单输入—单输出控制方法 (11) 超前滞后控制方法 2.2.2 单输入—多输出控制方法 (22) 双PID控制方法 2.2.3 多输入—多输出控制方法 (30) 极点配置法 二次线性最优控制法 三、二级倒立摆 (36) 3.1二级倒立摆建模 (36) 3.2 二级倒立摆控制方法 (46) 3.2.1 二次线性最优控制法 (46) 3.2.2 基于融合技术的模糊控制法 (48) 四、总结 (60) 五、参考文献 (63)

一、倒立摆系统介绍 1.1倒立摆系统简介 倒立摆是机器人技术、控制理论、计算机控制等多个领域、多种技术的有机结合，其被控系统本身又是一个绝对不稳定、高阶次、多变量、强耦合的非线性系统，可以作为一个典型的控制对象对其进行研究。最初研究开始于二十世纪50 年代，麻省理工学院（MIT）的控制论专家根据火箭发射助推器原理设计出一级倒立摆实验设备。近年来，新的控制方法不断出现，人们试图通过倒立摆这样一个典型的控制对象，检验新的控制方法是否有较强的处理多变量、非线性和绝对不稳定系统的能力，从而从中找出最优秀的控制方法。倒立摆系统作为控制理论研究中的一种比较理想的实验手段，为自动控制理论的教学、实验和科研构建一个良好的实验平台，以用来检验某种控制理论或方法的典型方案，促进了控制系统新理论、新思想的发展。由于控制理论的广泛应用，由此系统研究产生的方法和技术将在半导体及精密仪器加工、机器人控制技术、人工智能、导弹拦截控制系统、航空对接控制技术、火箭发射中的垂直度控制、卫星飞行中的姿态控制和一般工业应用等方面具有广阔的利用开发前景。平面倒立摆可以比较真实模拟火箭的飞行控制和步行机器人的稳定控制等方面的研究。 1.2倒立摆组成及其原理 倒立摆的组成包括计算机、运动控制卡、伺服系统、倒立摆本体和光电码盘、反馈测量元件等几大部分，组成一个闭环系统。对于直线型倒立摆，可以根据伺服电机自带的码盘反馈通过换算获得小车的位移，小车的速度信号可以通过差分法得到；各个摆杆的角度由光电码盘测得并直接反馈到控制卡，速度信号可以通过差分方法得到。计算机从运动控制卡中实时读取数据，确定控制策略（电机的输出力矩），并发送给运动控制卡。运动控制卡经过DSP 内部的控制算法实现该控制决策，产生相应的控制量，使电机转动，带动小车运动，保持摆杆平衡。

倒立摆仿真报告

计算机控制系统课题报告 1.倒立摆基本背景： 倒立摆，Inverted Pendulum ，是典型的多变量、高阶次，非线性、强耦合、自然不稳定系统。倒立摆系统的稳定控制是控制理论中的典型问题，在倒立摆的控制过程中能有效反映控制理论中的许多关键问题，如非线性问题、鲁棒性问题、随动问题、镇定、跟踪问题等。因此倒立摆系统作为控制理论教学与科研中典型的物理模型，常被用来检验新的控制理论和算法的正确性及其在实际应用中的有效性。从 20 世纪 60 年代开始，各国的专家学者对倒立摆系统进行了不懈的研究和探索。 倒立摆系统按摆杆数量的不同，可分为一级，二级，三级倒立摆等，多级摆的摆杆之间属于自由连接（即无电动机或其他驱动设备）。由中国的大连理工大学李洪兴教授领导的“模糊系统与模糊信息研究中心”暨复杂系统智能控制实验室采用变论域自适应模糊控制成功地实现了四级倒立摆。因此，中国是世界上第一个成功完成四级倒立摆实验的国家。 倒立摆的控制问题就是使摆杆尽快地达到一个平衡位置，并且使之没有大的振荡和过大的角度和速度。当摆杆到达期望的位置后，系统能克服随机扰动而保持稳定的位置。 2.倒立摆模型分析 倒立摆系统的输入为小车的位移（即位置）和摆杆的倾斜角度期望值，计算机在每一个采样周期中采集来自传感器的小车与摆杆的实际位置信号，与期望值进行比较后，通过控制算法得到控制量，再经数模转换驱动直流电机实现倒立摆的实时控制。直流电机通过皮带带动小车在固定的轨道上运动，摆杆的一端安装在小车上，能以此点为轴心使摆杆能在垂直的平面上自由地摆动。作用力F平行

于铁轨的方向作用于小车，使杆绕小车上的轴在竖直平面内旋转，小车沿着水平铁轨运动。当没有作用力时，摆杆处于垂直的稳定的平衡位置（竖直向下）。为了使杆子摆动或者达到竖直向上的稳定，需要给小车一个控制力，使其在轨道上被往前或朝后拉动。 我们的分析对象是一阶倒立摆。很多国内实验都说可以合理的假设空气阻力为0，但查阅了更多的文献和真正仿真做出模型并在网络上开源的一些实验后，我认为这是不正确的。空气阻力或许可以忽略，但是对于运动过程中的所有阻碍都忽略那就太为理想。也就是说，我们需要自己假设一个阻碍模型，即收到的所有阻力等效成一个包含速度，位姿等的广义函数。当然，我们的时间精力和所学知识都还有限，却也不想太过简单。我选取了一个阻力和速度成正比的函数关系，来在以后的建模和仿真过程中来模拟倒立摆所受的一切阻碍。 3.1 倒立摆物理建模：基于经典牛顿力学 受力分析如上图。 那我们在本实验中定义如下变量： M 小车质量 m 摆杆质量 b 小车摩擦系数 l 摆杆转动轴心到杆质心的长度（0.3 m）

2021年倒立摆实验报告(根轨迹)

*欧阳光明*创编 2021.03.07

I 摆杆惯量0.0034 kg*m*m g 重力加速度9.8 kg.m/s （2）直线一级倒立摆根轨迹校正控制原理 基于根轨迹法校正的基本思想是：假设系统的动态性能指标可由靠近虚轴的一对共轭闭环主导极点来表征，因此，可把对系统提出的时域性能指标的要求转化为一对期望闭环主导极点。确定这对闭环主导极点的位置后，首先根据绘制根轨迹的相角条件判断一下它们是否位于校正前系统的根轨迹上。如果这对闭环主导极点正好落在校正前系统的根轨迹上，则无需校正，只需调整系统的根轨迹增益即可；否则，可在系统中串联一个超前校正装置。 常见的校正器有超前校正、滞后校正以及超前滞后校正等。 2. 实验方法 （1）直线倒立摆建模、仿真与分析 利用牛顿-欧拉方法建立直线一级倒立摆系统的数学模型；依照根轨迹设计的步骤得到系统的控制器，利用MA TLAB Simulink中的工具进行仿真分析。 （3）直线一级倒立摆根轨迹校正控制 利用MATLAB Simulink来实现根轨迹校正控制参数设定和仿真，并利用该参数来设定只限一级倒立摆的根轨迹校正控制器值，分析和仿真倒立摆的运行情况。 3. 实验装置 直线单级倒立摆控制系统硬件结构框图如图1所示，包括计算机、I/O设备、伺服系统、倒立摆本体和光电码盘反馈测量元件等几大部分，组成了一个闭环系统。 图1 一级倒立摆实验硬件结构图 对于倒立摆本体而言，可以根据光电码盘的反馈通过换算获得小车的位移，小车的速度信号可以通过差分法得到。摆杆的角度由光电码盘检测并直接反馈到I/O设备，速度信号可以通过差分法得到。计算机从I/O设备中实时读取数据，确定控制策略（实际上是电

自动化实验-倒立摆实验-附仿真结果图复习课程

一、直线一级倒立摆的仿真 (一)直线一级倒立摆的数学建模 对于倒立摆系统，由于其本身是自不稳定的系统，实验建模存在一定的困难。但是忽略掉一些次要的因素后，倒立摆系统就是一个典型的运动的刚体系统，可以在惯性坐标系内应用经典力学理论建立系统的动力学方程。下面我们采用其中的牛顿－欧拉方法和拉格朗日方法分别建立直线型一级倒立摆系统的数学模型。 图2 直线一级倒立摆模型 φ摆杆与垂直向上方向的夹角； θ摆杆与垂直向下方向的夹角（考虑到摆杆初始位置为竖直向下）。 图3 小车及摆杆受力分析 分析小车水平方向所受的合力，可以得到以下方程： 由摆杆水平方向的受力进行分析可以得到下面等式：

把这个等式代入式1中，就得到系统的第一个运动方程： 为了推出系统的第二个运动方程，我们对摆杆垂直方向上的合力进行分析，可以得到下面方程： 力矩平衡方程如下： 注意：此方程中力矩的方向，由于θ=π+φ,cosφ= ?cosθ,sinφ= ?sinθ，故等式前面有负号。 合并这两个方程，约去P 和N，得到第二个运动方程： 设θ=π+φ（φ是摆杆与垂直向上方向之间的夹角），假设φ与1（单位是弧度）相比很小，即φ<<1，则可以进行近似处理： 。 用u 来代表被控对象的输入力F，线性化后两个运动方程如下： 对式9进行拉普拉斯变换，得到 注意：推导传递函数时假设初始条件为0。 由于输出为角度φ，求解方程组的第一个方程，可以得到：

或 如果令v = x，则有： 把上式代入方程组的第二个方程，得到： 整理后得到传递函数： 其中 设系统状态空间方程为： 方程组对解代数方程，得到解如下： 整理后得到系统状态空间方程：

一级倒立摆控制系统设计

基于双闭环PID控制的一阶倒立摆控制系统设计 一、设计目的 倒立摆是一个非线性、不稳定系统,经常作为研究比较不同控制方法的典型例子。设计一个倒立摆的控制系统，使倒立摆这样一个不稳定的被控对象通过引入适当的控制策略使之成为一个能够满足各种性能指标的稳定系统。 二、设计要求 倒立摆的设计要求是使摆杆尽快地达到一个平衡位置，并且使之没有大的振荡和过大的角度和速度。当摆杆到达期望的位置后，系统能克服随机扰动而保持稳定的位置。实验参数自己选定，但要合理符合实际情况，控制方式为双PID控制，并利用MATLAB进行仿真，并用simulink对相应的模块进行仿真。 三、设计原理 倒立摆控制系统的工作原理是：由轴角编码器测得小车的位置和摆杆相对垂直方向的角度，作为系统的两个输出量被反馈至控制计算机。计算机根据一定的控制算法，计算出空置量，并转化为相应的电压信号提供给驱动电路，以驱动直流力矩电机的运动，从而通过牵引机构带动小车的移动来控制摆杆和保持平衡。 四、设计步骤 首先画出一阶倒立摆控制系统的原理方框图 一阶倒立摆控制系统示意图如图所示： 分析工作原理，可以得出一阶倒立摆系统原理方框图：

一阶倒立摆控制系统动态结构图 下面的工作是根据结构框图，分析和解决各个环节的传递函数！ 1.一阶倒立摆建模 在忽略了空气流动阻力，以及各种摩擦之后，可将倒立摆系统抽象成小车和匀质杆组成的系统，如下图所示，其中： M ：小车质量 m ：为摆杆质量 J ：为摆杆惯量 F ：加在小车上的力 x ：小车位置 θ：摆杆与垂直向上方向的夹角 l ：摆杆转动轴心到杆质心的长度 根据牛顿运动定律以及刚体运动规律，可知： （1） 摆杆绕其重心的转动方程为 （2） 摆杆重心的运动方程为 得 sin cos ..........(1)y x J F l F l θθθ=- 2 22 2(sin ) (2) (cos ) (3) x y d F m x l d t d F mg m l d t θθ=+=-

倒立摆实验报告根轨迹

专业实验报告

（2）直线一级倒立摆根轨迹校正控制原理 基于根轨迹法校正的基本思想是：假设系统的动态性能指标可由靠近虚轴的一对共轭闭环主导极点来表征，因此，可把对系统提出的时域性能指标的要求转化为一对期望闭环主导极点。确定这对闭环主导极点的位置后，首先根据绘制根轨迹的相角条件判断一下它们是否位于校正前系统的根轨迹上。如果这对闭环主导极点正好落在校正前系统的根轨迹上，则无需校正，只需调整系统的根轨迹增益即可；否则，可在系统中串联一个超前校正装置。 常见的校正器有超前校正、滞后校正以及超前滞后校正等。 2. 实验方法 （1）直线倒立摆建模、仿真与分析 利用牛顿-欧拉方法建立直线一级倒立摆系统的数学模型；依照根轨迹设计的步骤得到系统的控制器，利用MATLAB Simulink中的工具进行仿真分析。 （3）直线一级倒立摆根轨迹校正控制 利用MATLAB Simulink来实现根轨迹校正控制参数设定和仿真，并利用该参数来设定只限一级倒立摆的根轨迹校正控制器值，分析和仿真倒立摆的运行情况。 3. 实验装置 直线单级倒立摆控制系统硬件结构框图如图1所示，包括计算机、I/O设备、伺服系统、倒立摆本体和光电码盘反馈测量元件等几大部分，组成了一个闭环系统。 图1 一级倒立摆实验硬件结构图 对于倒立摆本体而言，可以根据光电码盘的反馈通过换算获得小车的位移，小车的速度信号可以通过差分法得到。摆杆的角度由光电码盘检测并直接反馈到I/O设备，速度信号可以通过差分法得到。计算机从I/O设备中实时读取数据，确定控制策略（实际上是电机的输出力矩），并发送给I/O设备，I/O设备产生相应的控制量，交与伺服驱动器处理，然后使电机转动，带动小车运动，保持摆杆平衡。

一阶倒立摆控制系统设计

课程设计说明书 课程名称：控制系统课程设计设计题目：一阶倒立摆控制器设计院系：信息与电气工程学院班级： 设计者： 学号： 指导教师： 设计时间：2013年2月25日到2013年3月8号

课程设计（论文）任务书 指导教师签字：系（教研室）主任签字： 2013年3月5日

目录 一、建立一阶倒立摆数学模型 (4) 1. 一阶倒立摆的微分方程模型 (4) 2. 一阶倒立摆的传递函数模型 (6) 3. 一阶倒立摆的状态空间模型 (7) 二、一阶倒立摆matlab仿真 (9) 三、倒立摆系统的PID控制算法设计 (13) 四、倒立摆系统的最优控制算法设计 (23) 五、总结............................................................................................... 错误！未定义书签。 六、参考文献 (29)

一、建立一阶倒立摆数学模型 首先建立一阶倒立摆的物理模型。在忽略空气阻力和各种摩擦之后，可将直线一级倒立摆系统抽象成小车和匀质杆组成的系统，如图1所示。 系统内部各相关参数定义如下： M 小车质量 m 摆杆质量 b 小车摩擦系数 l 摆杆转动轴心到杆质心的长度 I 摆杆惯量 F 加在小车上的力 x 小车位置 φ摆杆与垂直向上方向的夹角 θ摆杆与垂直向下方向的夹角（考虑到摆杆初始位置为竖直向下）1.一阶倒立摆的微分方程模型 对一阶倒立摆系统中的小车和摆杆进行受力分析，其中，N和 P为小车与摆杆相互作用力的水平和垂直方向的分量。

图1-2 小车及摆杆受力图 分析小车水平方向所受的合力，可以得到以下方程： （1-1）由摆杆水平方向的受力进行分析可以得到下面等式： （1-2）即： （1-3） 把这个等式代入式(1-1)中，就得到系统的第一个运动方程： （1-4） 为了推出系统的第二个运动方程，我们对摆杆垂直方向上的合力进行分析，可以得到下面方程： （1-5） 即： （1-6）力矩平衡方程如下： （1-7） 由于所以等式前面有负号。

基于matlab的倒立摆仿真设计

基于matlab的倒立摆的仿真与设计姓名:贾永伟专业:测控技术与仪器学号:1123105950 年级:2011级 摘要：倒立摆系统是一个典型的快速、多变量、非线性、不稳定系统，对倒 立摆的控制研究无论在理论上和方法上都有深远的意义。 本论文以实验室原有的直线一级倒立摆实验装置为平台，重点研究其PID控制方法，设计出相应的PID控制器，并将控制过程在MATLAB上加以仿真。 关键词：一级倒立摆，PID，MATLAB仿真 一、倒立摆模型的研究意义 倒立摆系统的研究能有效的反映控制中的许多典型问题：如非线性问题、鲁棒性问题、镇定问题、随动问题以及跟踪问题等。通过对倒立摆的控制，用来检验新的控制方法是否有较强的处理非线性和不稳定性问题的能力。同时，其控制方法在军工、航天、机器人和一般工业过程领域中都有着广泛的用途，如机器人行走过程中的平衡控制、火箭发射中的垂直度控制都有重要意义 倒立摆控制系统是一个复杂的、不稳定的、非线性系统，是进行控制理论教学及开展各种控制实验的理想实卫星飞行中的姿态控制等。故其研究意义广泛。 二、倒立摆模型的数学建模 质量为m的小球固结于长度为L的细杆（可忽略杆的质量）上，细杆又和质量为M的小车铰接相连。由经验知：通过控制施加在小车上的力F（包括大小和方向）能够使细杆处于θ＝0的稳定倒立状态。在忽略其他零件的质量以及各种摩擦和阻尼的条件下，推导小车倒立摆系统的数学模型 分析过程如下： 如图所示，设细杆摆沿顺时针方向转动为正方向，水平向右方向为水平方向上的正方向。当细杆摆顺时针往右运动时水平方向施加的力应该为水平向右。 现对小车和细杆摆分别进行隔离受力分析：

单级倒立摆实验报告

单级倒立摆实验报告 1. 单级倒立摆系统的建模 单级倒立摆系统的建模可采用受力分析或Lagrange 方程建立得到。这里采用受力分析方法建模。如图所示： 根据牛顿第二定律： (cos )0Mx m x L u θθ++-= (2-1) cos sin 0mLx I mLg θθθ--= (2-2) 以摆杆偏角θ、角速度θ 、小车的位移x 和 小车速度x 为状态变量，即令： () T X x x θθ= (2-3) 同时假设倒立摆摆杆的垂直倾斜角度θ与1 （单位为rad ）相比很小，即1θ 。 则可以近似处理：cos θ≈1，sin 0θ≈，并 忽略高阶小量，则可得： 2222 ()()m L g I x u I m M mML I m M mML θ=+++++ （2-4） 22 ()()()mL m M g mL u I m M mML I m M mML θ θ+=-+++++ （2-5） 摆杆系统的状态方程为： 1222 2122 344122 ()()()()()x x m L g I x x u I m M mML I m M mML x x mL m M g mL x x u I m M mML I m M mML =???=+?++++? =??+=-+?++++? （2-6） 写成向量的形式为： X AX Bu y CX Du ?=+? =+? (2-7) 其中

0100000A 00010 00a b ?? ? ?= ? ? ??, 00c B d ?? ? ?= ? ??? ,10000010C ??= ???,00D ??= ??? （2-8） 参数a 、b 、c 、d 分别为： 222()m L g b I m M mML = ++ （2-9） 2 ()()mL m M g a I m M mML +=- ++ （2-10） 2 ()I c I m M mML = ++ （2-11） 2 ()mL d I m M mML =++ （2-12） 选择摆杆的倾斜角度θ和小车的水平位移x 作为系统的输出，则输出方程为： y CX = (2-13) 根据金棒-2型倒立摆系统实验平台的参数，m=0.2kg ，M=0.6kg ，L=0.158m ，I=0.001654kg.m 2 ，g=10N/kg.同时，这里建模时候使用的ｕ是以力作为输入信号的，实际上采用的是以电压作为输入信号，通过电机作了一定的转化，这里我们约定：先暂时以力作为输入信号，最后再统一处理。则有，a=2.3121，b=-58.5337，c=0.3830，d=7.3167。 因此，010000 2.31210A 00010058.53370?? ? ?= ? ?-??,00.383007.3167B ?? ? ?= ? ??? 2. 全状态反馈设计 2.1. 检验系统可控性 可控性矩阵纯ctrB=105 *0 00.00020.005300.00020.00530.148200.0001-0.00310.09370.0001-0.00310.0937 2.5164-????--? ?????-?? 显然rank(ctrB)=4,系统可控. 2.2. 反馈设计 要求：稳定调节时间3s n t s π ξω= <，摆角5θ< ，(5/90100) 5.56p σ

倒立摆系统建模及MATLAB仿真

倒立摆系统的建模及MATLAB仿真 通过建立倒立摆系统的数学模型,应用状态反馈控制配置系统极点设计倒立摆系统的控制器,实现其状态反馈,从而使倒立摆系统稳定工作。之后通过MA TLAB 软件中Simulink工具对倒立摆的运动进行计算机仿真,仿真结果表明,所设计方法可使系统稳定工作并具有良好的动静态性能。 倒立摆系统是1个经典的快速、多变量、非线性、绝对不稳定系统,是用来检验某种控制理论或方法的典型方案。倒立摆控制理论产生的方法和技术在半导体及精密仪器加工、机器人技术、导弹拦截控制系统和航空器对接控制技术等方面具有广阔的开发利用前景。因此研究倒立摆系统具有重要的实践意义,一直受到国内外学者的广泛关注。本文就一级倒立摆系统进行分析和研究,建立倒立摆系统的数学模型,采用状态反馈极点配置的方法设计控制器,并应用MA TLAB 软件进行仿真。 1 一级倒立摆系统的建模 1. 1 系统的物理模型 如图1 所示,在惯性参考系下,设小车的质量为M ,摆杆的质量为m ,摆杆长度为l ,在某一瞬间时刻摆角(即摆杆与竖直线的夹角)为θ,作用在小车上的水平控制力为f 。这样,整个倒立摆系统就受到重力,水平控制力和摩擦力的3 外力的共同作用。 图1 一级倒立摆物理模型 1. 2 系统的数学模型 在系统数学模型中,本文首先假设: (1) 摆杆为刚体。 (2)忽略摆杆与支点之间的摩擦。 (3)忽略小车与导轨之间的摩擦。 然后根据牛顿第二运动定律,求得系统的运动方程为:

方程(1) , (2) 是非线性方程,由于控制的目的是保持倒立摆直立,在施加合适的外力条件下,假定θ很小,接近于零是合理的。则sinθ≈θ,co sθ≈1 。在以上假设条件下,对方程线性化处理后,得倒立摆系统的数学模型: 1. 3 系统的状态方程 以摆角θ,角速度θ',小车的位移x ,速度x'为状态变量,输出为y 。即令: 则一级倒立摆系统的状态方程为: 2 控制器设计及MATLAB 仿真 2. 1 极点配置状态反馈的基本原理