模糊控制课程设计报告书

模糊控制

学 院：电气工程学院 姓 名： 明 学 号：

任课教师： 刁 晨

班 级： 09级自动化3班

单倒置摆控制系统的状态空间设计

1 / 14

. . . .

一． 设计题目

1. 介绍

单倒置摆系统的原理图，如图1所示。设摆的长度为L、质量为m，用铰链安装在质量为M的小车上。小车有一台直流电动机拖动，在水平方向对小车施加控制力u,相对参考系产生位移z。若不给小车施加控制力，则倒置摆会向左或向右倾倒，因此，它是一个不稳定系统。控制的目的是，当倒置摆无论出现向左或向右倾倒时，通过控制直流电动机，使小车在水平方向运动，将倒置摆保持在垂直位置上。

m z l θ u M

2. 用途

倒立摆系统以其自身的不稳定性为系统的平衡提出了难题,也因此成为自动控制实验中验证控制算法优劣的极好的实验装置。单倒立摆的系统结构、数学模型以与系统的稳定性和可控性,对倒立摆进行了成功的控制,并在MATLAB 中获得了良好的仿真效果。倒立摆控制理论将在半导体与精密仪器加工、机器人技术、伺服控制领域、导弹拦截控制系统、航空器对接技术等方面具有广阔的开发利用前景。

3. 意义

倒立摆是一种典型的快速、多变量、非线性、绝对不稳定系统. 人们试图寻找同的控制方法以实现对倒立摆的控制,以便检验或说明该方法对严重非线性和绝对不稳定系统的控制能力 。同时，由于摩擦力的存在，该系统具有一定的不确定性。对这样一个复杂系统的研究在理论上将涉与系统控制中的许多关键问题：如非线性问题、鲁棒性问题、镇定问题、随动问题以与跟踪问题等都可以以它为例进行研究。

二． 被控对象的模型

为简化问题，工程上往往忽略一些次要因素。这里，忽略摆杆质量、执行电动机惯性以与摆轴、轮轴、轮与接触面之间的摩擦与风力。设小车瞬时位置为z，倒置摆出现的偏角为

2 / 14

. . . .

θ，则摆心瞬时位置为(z+lsinθ)。在控制力u的作用下，小车与摆均产生加速运动，根据牛顿第二定律，在水平直线运动方向的惯性力应与控制力u平衡，则有

d2zd2M2m2(zlsinθ)u dtdt即

•••••2(Mm)zmlθcosθ -mlθsinθu

d2[m2(zlsinθ)]lcosθmglsinθ dt即

•••••2 (1)

由于绕摆轴旋转运动的惯性力矩与重力矩平衡，因而有

zcoslcoslsincosgsin

•2 (2)

式(1)、式(2)两个方程都是非线性方程，需作线性化处理。由于控制的目的是保持倒置摆直立，因此，在施加合适u的条件下，可认为、均接近零，此时sin≈，cos≈1，且可忽略项，于是有

•2(Mm)zmlu

••••••••

(3)

(4)

zlg

联立求解式(3) 、式(4)，可得

mg1u (5) MM••(Mm)1gu (6)

MlMl消去中间变量，可得输入变量为u、输出变量为z的系统微分方程为

(Mm)g••1••g(4)zzuu (7)

MlMMl••z在此问题中，设上述两变量的论域为2x12和5radsx25rads,则设计步骤为

第1步。首先，对x1在其论域上建立三个隶属度函数，即如图 1所示的正值（P）、零（Z）和负值（N）。然后，对x2在其论域上亦建立3个隶属度函数，即图2所示的正值（P）、零（Z）和负值（N）。

3 / 14

. . . .

图2-3 输入x1的分区

图2-4输入x2的分区

第2步。为划分控制空间（输出），对uk在其论域上建立5个隶属度函数，

24uk24，如图3（注意，图上划分为7段，但此问题中只用了5段）。

4 / 14

. . . .

图2-5输出u的分区

第3步。用表1所示的3*3规则表的格式建立9条规则（即使我们可能不需要这么多）。本系统中为使倒立摆系统稳定，将用到和ddt。表中的输出即为控制作用u(t)。 X1 x2 P Z N P PB P Z Z P Z N N Z N NB 表1模糊控制规则表 第4步。我们可用表1中规则导出该控制问题的模型。并用图解法来推导模糊运算。假设初始条件为

x101 和 x204rads

然后，我们在上例中取离散步长0k3，并用矩阵差分方程式导出模型的四部循环式。模型的每步循环式都会引出两个输入变量的隶属度函数，规则表产生控制作用u(k)的隶属度函数。我们将用重心法对控制作用的隶属度函数进行精确化，用递归差分方程解得新的x1和x2值为开始，并作为下一步递归差分方程式的输入条件。

分别为x1和x2的初始条件。从模糊规则表（表1）有 If(x1=P)and(x2=Z),then(u=P) If(x1=P)and(x2=N),then(u=Z) If(x1=Z)and(x2=Z),then(u=Z) If(x1=Z)and(x2=N),then(u=N)

表示了控制变量u的截尾模糊结果的并。利用重心法精确化计算后的控制值为u=-2。

在已知u=-2控制下，系统的状态变为

x11x10x203 x21x10x20u01

依次类推，可以计算出下一步的控制输出u(1)。模糊控制器能够满足倒立摆的运动控制。

三、模糊控制器的建立

3.1在MTALAB中的fuzzy 控制器的建立与封装

在命令窗口中输入：fuzzy然后回车可得出如下图所示：

5 / 14

. . . .

图3-1 模糊控制器设置界面

然后对其各个变量进行设置其步骤如下图3-2：

对输入变量X1进行设置如下图3-3所示：

6 / 14

. . . .

变量X2的设置如下图3-4所示：

输出量的设置图3-5所示：

7 / 14

. . . .

模糊规则控制表的设置如下图3-6所示：

8 / 14

. . . .

设置出来的效果图如图3-7(a),(b),(c)所示：

9 / 14

. . . .

（a）

10 / 14

. . . .

（b）

11 / 14

. . . .

（c)

3.2 最终在MATLAB中的搭建出来的框图如下：

图3-8 单级倒立摆在MTALAB中simulink仿真的框架图

主要的状态空间模块的参数设置如下：

12 / 14

. . . .

四、 仿真结果以与分析

通过（fuzzy)模糊控制模块，可以和包含模糊控制器的fis文件联系起来，还可以随时改变输入输出论域，隶属度函数以与模糊规则。仿真结果如下图：图4-1和图4-2。

图4-1

分析如下：从图4-1仿真图中可以看出，仿真时间大概在1秒左右趋于平衡，但是图中曲线最终稳定在-2.3左右，而不是在0附近稳定，猜想曲线应该最终稳定于0附近。

13 / 14

. . . .

图 4-2

分析如下：图4-1从图4-2仿真图中可以看出，仿真时间大概也在1秒左右趋于平衡，图中

曲线最终稳定在0.3左右，接近于0附近稳定，基本实现了仿真预期效果。

五、 结语

通过以上仿真实践的过程可以看出，对于倒立摆一类的包含非线性以与不稳定因素的控制系统模型，采用模糊控制的方法具有很强的适应性以与理想的控制效果，在此基础上，选用恰当的隶属度函数也是实现优化控制的重要条件。 在实际应用中，由于模型往往具有更多的复杂因素以与一些不可预测的干扰因素，所以有时需要简历更复杂的模糊控制方案以与更加精确的隶属度函数，有时也需要采用神经网络的算法结构来解决这样的实际问题，这些方法互相配合从而形成了较为完备的智能控制方案系统。

参考文献:

[1]力群.智能控制理论与应用 [2]胡寿松.自动控制原理简明教程

[3]永华. 新型PID 控制与其应用[M] . :机械工业, 2003.

[4]红兵,生权,王瑜. 基于RBF 网络二级倒立摆系统PID 控制[J ] . 微计算机信息,2007 (6) : 72 - 75.

[5]豹. 现代控制理论[M].:机械工业,1992.

[4]桂明，等．应用Mat]ab建模与仿真[M]．．科学，2000． [5]闻新，等模糊逻辑工具葙的分析与应用[M]，科学．2000． [6]志涌. 精通MATLAB [M ]. (第三版). :航空航天大学, 2000.

14 / 14