鲁棒控制系统设计

鲁棒控制设计报告

学院

专业

报告人

目 目

1 绪论 (2)

1.1控制系统设计背景 (2)

1.2本文主要工作分配 (3)

2 一级倒立摆模型建立 (4)

2.1一级倒立摆的工作原理 (4)

2.2一级倒立摆的数学模型 (4)

3 H∞鲁棒控制器设计 (6)

3.1基于Riccati方程的H∞控制 (7)

3.2基于LMI的H∞控制 (7)

4 一级倒立摆系统的仿真 (9)

4.1一级倒立摆控制系统设计 (9)

4.2闭环控制系统仿真及分析 (10)

5 结论 (13)

1 绪论

1.1控制系统设计背景

一级倒立摆系统是一个典型非线性多变量不稳定系统，在研究火箭箭身的姿态稳定控制、机器人多自由度运动稳定设计、直升机飞行控制等多种领域中得到了广泛的应用，因此以倒立摆作为被控对象进行控制方法的研究具有重要的现实意义。为解决一级倒立摆系统的非线性、强耦合、多变量、自然不稳定问题，本文利用H∞鲁棒控制实现对一级倒立摆的控制。

Mg

图1.1 一级倒立摆系统结构图

本文采用的直线一级倒立摆的基本系统如图1.1所示，它是由沿直线导轨运动的小车以及一端固定于小车上的材质均匀的摆杆组成，它是一个不稳定的系统，当倒立摆出出现偏角θ后，如果不给小车施加控制力，倒立摆会倾倒。所以本文采用H∞鲁棒控制方法的目的是通过调节水平力F的大小控制小车的运动，使倒立摆处于竖立的垂直位置。控制指标为：倒立摆系统的从初始状态调节到小车停留在零点、并使摆杆的摆角为0的稳定状态。

1.2本文主要工作分配

第一章：对一级倒立摆系统的特点、结构以及控制要求进行阐述。

第二章：根据一级倒立摆的结构，利用机理建模法建立被控对象的精确数学模型，并在系统平衡点处进行线性化，得到系统简化的状态方程。

第三章：首先H∞鲁棒控制的基本原理，然后分别利用Riccati方程和LMI 方法设计H∞状态反馈控制器。

第四章：首先使用MATLAB计算基于Riccati方程的H∞状态反馈控制器和基于LMI的H∞状态反馈控制器，然后进行闭环控制系统的仿真并控制系统的性能分析。

第五章：对本次设计进行总结。

2 一级倒立摆模型建立

2.1一级倒立摆的工作原理

如图1.1所示，倒立摆装置主要由摆杆、小车以及导轨组成。导轨的一端装有用来测量小车位移的电位计，摆杆与小车的连接处安装测量摆角的装置，小车可以沿着有界轨道直线移动，同时摆杆可以在垂直平面内自由运动。直流电动机通过传送带拖动小车运动，从而使倒立摆稳定在竖立的垂直位置。为简化系统分析，在实际模型建立过程中，忽略空气流动的阻力以及各种摩擦力，这样可以将倒立摆抽象为由小车和均匀材质的刚性摆杆组成的系统。小车质量为M ，摆杆质量为m ，小车位置x ，作用在小车上力大小为F ，摆杆的长度为2L l =，均匀材质的摆杆质心是摆杆的中心。

2.2一级倒立摆的数学模型

被控对象的数学模型是过程中的输入量和输出量之间的函数关系，常用的有机理建模法和实验建模两种方法。本文采用的是机理建模的方法，根据过程的内在机理，利用相关的平衡方程，获得所需要的数学模型。

对摆杆进行受力分析，转动惯量与加速度的乘积等于刚体主动力对该轴力矩的代数和，则摆杆绕其重心的转动方程为：

sin cos y x J F l F l θθθ=- (2.1)

摆杆在水平方向上受到的合力为：

2

2(sin )x d F m x l dt

θ=+ (2.2) 摆杆在垂直方向上受到的合力为：

2

2(sin )y d F mg m x l dt

θ-=+ (2.3) 小车在水平方向上受到合力：

22y d x F F M dt

-= (2.4)

将等式(2.2)(2.3)分别带入等式(2.1)和(2.4)中：

()2

cos sin J ml mlx mgl θθθ++= (2.5) ()()

2cos sin F M m x ml θθθθ=++⋅-⋅ (2.6)

整理得到系统精确模型为： ()()()()()()()2222222222222222sin sin cos cos lgsin cos sin cos cos J ml F ml J ml m l g x J ml M m m l M m m mlF m l J ml M m m l θθθθθθθθθθθθ⎧+++-⎪=++-⎪⎨+--⎪=⎪++-⎩

(2.7) 式中，转动惯量2

3

ml J =。 由等式(2.7)得知，一级直线倒立摆系统的动力学模型为非线性微分方程，因此选择工作点00θ=、00x =对系统进行线性化，即可近似认为0θ≈、sin θθ≈、cos 1θ≈，得到进一步的简化模型： ()()()()()()2222222g M m m l ml F J M m mMl J M m mMl J ml m gl x F J M m mMl J M m mMl θθθ+⎧=-⎪++++⎪⎨+⎪=-+⎪++++⎩

(2.8) 以摆杆与竖直向上方向的偏角θ，小车的位移x 、摆杆摆角变化θ和小车的速度x 作为四个状态变量，考虑控制输入干扰ω，将(2.8)转化为状态方程的形式：

12x Ax B B u ω=++ (2.9)

式中，u F =，x x x θθ⎡⎤⎢=⎥

⎢⎥⎢⎥⎢⎥⎣⎦，1200100001000000k A k ⎡⎤⎢⎥⎢⎢⎢=⎥⎥⎥⎣⎦，21100B b b ⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦=，42300B k k ⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦

=，其中1b 、2b 为不大于1的正数，()()12

g M m m l k J M m mMl +=++，()22

22m gl k J M m mMl =-++，()32ml k J M m mMl =-++，()()2

42J ml k J M m mMl +=++。