模糊控制在工程中的应用.docx

1、模糊控制在工程中的应用第 8 章 模糊控制在工程中的应用8.1倒立摆系统的 T-S 模型模糊控制模糊控制在工业过程控制、 机器人控制、 运载工具控制及家电产品等领域有 着广泛的应用，本章重点介绍几个例子。本节课介绍基于 T-S 模糊模型的倒立摆控制。8.1.1倒立摆系统概述倒立摆系统是一个典型的非线性、 强耦合、 多变量和不稳定的高阶系统， 许 多抽象的控制理论概念都可以通过倒立摆实验直观的表现出来， 是控制理论教学 的理想实验设备和进行控制理论研究的典型实验平台， 也是新成果、 新方法的验 证平台，开发平台，一直受到教学和科研人员的广泛关注。因此，从其肇始之日 至今的半个世纪的发展历程中，

2、先后出现了形式各异的倒立摆， 大致可以分为以 下五大类 :直线倒立摆、平面倒立摆、斜轨道和圆轨道倒立摆以及并行倒立摆。(1)直线倒立摆 直线倒立摆是由可以沿直线导轨运动的小车以及一端铰接于小车之上的匀 质长杆组成的系统， 如图 1.1 所示。对于单级倒立摆和二级倒立摆系统的研究已 经历了很长的历程， 并且有很多控制成功的报道。 在此基础上， 三级倒立摆的研 究也取得了很大进展， 不仅在系统仿真方面， 而且在实物实验中， 都出现了控制 成功的范例。北京师范大学李洪兴教授分别于 2001 年 6 月和 2002 年 8 月完成 了四级倒立摆系统的仿真和实物实验， 是目前世界上控制成功的多级倒立摆系

3、统 中级数最多的Mntxnr胃线倒立摆(2)平面倒立摆如果小车在水平面内自由运动，即为二维倒立摆系统。图1.2是一种旋臂式 二维单级倒立摆的示意图：通过两个电机Ma和Mb分别控制后臂和前臂来控制摆 杆支点在水平面的自由运动，并进一步控制摆杆的平衡。其中一为4个测量 角度的位置传感器。还有一种小车式二维倒立摆：使用两个电机分别控制X轴和 丫轴的运动，使得摆杆支点在水平面内自由运动，并进一步控制摆杆的平衡。(3)斜轨道和圆轨道倒立摆如果小车运动轨迹不是水平的直线，而是在倾斜的轨道上或圆形的轨道上运 动，即为斜轨道或圆轨道的倒立摆系统。其中因斜轨道型二级倒立摆系统与实际 的控制问题模型相近，对其进行

4、的研究也比较广泛。斜轨道二级倒立摆如图 1.3所示，其轨道与水平方向成a的夹角。图1.4为圆轨道单级倒立摆的示意图。电 机带动旋臂旋转，将摆杆的支点限制于以旋臂长度为半径的圆形轨道上， 并控制 倒立摆的两级摆杆平衡斜轨道倒立擢(4)并行倒立摆对一些其他类型的倒立摆系统，也有人对其进行了研究。如图1.5所示的并行倒立摆系统。所谓并行倒立摆系统，就是在同一个小车上安装两根互相独立的单级摆杆，通过驱动小车来实现同时保持两根摆杆平衡的控制目标。(5)旋转式倒立摆旋转式倒立摆系统是不通过小车，直接利用电机转动进行控制的倒立摆系 统。图1.6是电机带动旋臂控制摆杆角度的倒立摆。旋转倒直摆8.1.2T-S模

5、糊模型模糊模型在形式上表现为一系列“ IF THEN ”的模糊规则的组合。按照模 糊规则后件不同的结构，我们可以将模糊模型分为三种，即Mamdani模糊模型、 模糊关系模型和T-S模糊模型。其中T-S模糊模型得到了众多学者的广泛关注， 也是本文研究的重点。Takagi和Sugeno于1985年提出了著名的T-S模糊模型，旨在开发从给定 的输入-输出数据集产生模糊规则的系统化方法。在 T-S模糊模型中，规则后 件是模型输入的函数：R : If x is Athen y = fj(x), i =1,2,.,K其中，xRm是输入(前件)变量；y R是输出(后件)变量；R表示第i条 规则；K是规则库中

6、的规则数；A是第i条规则的前件模糊集合，通过隶属函数 来定义：(x)：Rm 0,1其前件命题“Xis A ”通常表示成对于x单独成分定义的不变模糊集合简单命题的 逻辑组合，通常为下面的组合形式：Ri : If x1 is A|i and x is A； and and Xm is 从，Then % = f (xx；, Xm) ； i = 1, 2, , K (8-1)其中A为模糊子集，隶属函数可以取三角形、梯形或者高斯型； f(X1,X；，Xm)是后件的精确函数，通常是输入变量 Xj的多项式，也可以是任意函数。当f(X1,X；, Xm)为一阶多项式且带有常数项，即：Yi = aT x bi (

7、8-2)其中，ai是参数变量，bi是标量补偿。我们称这种模糊模型为仿射 T-S模糊模 型。当bi =0,i =1,2/ ,K时，结论函数成为一种特殊形式，这时模型称为齐次 T-S模型(或线性T-S模型)：If xis Ai then 比=aTx, i =1,2,.,K这种模型与仿射T-S模型相比，逼近非线性系统能力是有限的。当a： =0,i =1,2,，K时，模型结论部分是一常数，所得到的模型称为零阶 T-S模糊模型，也称为单点T-S模糊模型：If x is Athen % = b, i =1,2,.,K这个模型也可以看成是语言模糊模型结论模糊集简化成单一值的特殊形式。 一般意义上的T-S模糊

8、模型就是指仿射T-S模糊模型。由于T-S模糊模型的后件为线性函数的形式，因此避免了繁琐的去模糊化过程，系统的输出可以表示为：W1 Zt =卩必 亠q亠Z?二 qx - q2y D加权平均z WtZW2Z2=3也图8-1 T-S模糊模型推理过程图8-1为一阶T-S模糊模型的模糊推理过程。由于每条规则都有一个精确输 出，通过加权平均可以得到整体输出，从而避免 Mamdani模型所需耗时的去模 糊过程。有时，一个简单的T-S模糊模型可以产生复杂的行为，下面是一个两输入系统的例子例2两输入单输出T-S模糊模型如果 X 小，and Y 小，则 z=-x+y+1;如果 X 小, and Y 大，贝U z=

9、-y+3;如果 X 大，and Y 小,则 z=-x+3;如果 X 大，and Y 大，贝U z=x+y+2;图8-2 (a)给出了输入X和丫的隶属函数，图8-2 (b)是所产生的输入一 输出曲面。曲面是复杂的，但仍然可以看出，曲面由几个平面构成，每个平面由 一条模糊规则的输出方程来描述。LET-小大rr-3-10X135*111小大-10.50-510.5-335rp-?B m 3 am r Dnjfnrm Rp 3 rae图8-2例2中T-S模糊模型：(a)隶属函数；(b)输入一输出曲面如果T-S模糊系统的输出作为它自身的一个输入，则得到所谓的动态 T-S模糊系统。具体地讲，一个动态 T-

10、S模糊系统采用的模糊规则，最常用的形式为非线性自回归NARX模型。在系统输入输出变量的阶次己知时，可将式 转述为Ri : If y(k) is A； and y(k1) is A； and and y(km 1) is 您，and u(k) is B1 and u(k1) is B2 and and u(k-n 1) is B；m n(8-6)Then yi (k 1) a；y(k - p 1)亠二 b；u(k - q 1) C ;p d q =1模糊系统的输出为各子系统输出的加权平均，即K wi yi (k 1)y(k 1)= k . (8-7)Z wii 二其中模糊化采用单点模糊化，清晰化

11、采用加权平均法。W，是第i条规则的适应度，运算采用求积法，即m nwi A；y(kp 1)【B；u(kq 1). (8-8)pd qd8.1.3单级倒立摆系统的数学模型我们实验中采用的是固高科技有限公司生产的 L1IP型直线单级倒立摆。单 级倒立摆的结构如图3所示。其中M代表金属小车，m代表摆杆。摆杆与小车 自由连接。小车由电机通过皮带驱动并沿着直线轨道运动。 控制目的是通过电机 驱动小车，使摆杆动态平衡在倒立状态，这一过程即平衡控制；在此前提下，将 小车控制到指定位置，即位置控制。图中，r :小车的位移，单位(m)；9 :摆杆与垂直位置方向夹角，单位(rad);M :小车的质量，单位(kg)

12、;m :摆杆的质量，单位(kg)；l :摆杆的中心到转轴的长度，单位(m)；J :摆杆对重心的转动惯量，单位(kg m2)；u :电机对小车施加的作用力，单位(N);F :小车所受的等效摩擦力，单位(N);卩：小车所受的等效摩擦系数，单位(kg/s)；f :摆杆所受的摩擦阻力矩系数，单位(kg m2/s)；规定图示中r的方向为位移的正方向，顺时针方向为摆角的正方向。则该单级倒立摆动力学非线性方程组为i广 “ m Omlrcos 日+ (ml +J)T=mglsi n fB .- 2 . (&9)(M +m)r+mlB cos=u+ml日 sinT r当选取的状态变量为xr ； X2 - v ；

13、 X3二r ； X4 - v， 丫二rjT为输出向量。(8-9)式可化为以下一阶非线性方程组，Xi = 3X2 =沧2 2 2 2 2 2得到倒立摆系统的线性状态方程:- mlfX4 cosx2 m l g sin x2 cosx2(J + ml )x3+(J+ml )mlx4 sinx2+(J+ml )u(8-12)其中，相应矩阵如下所示:0001;B =00t32t33t34b30t42t43t44 _1 Ib4 一A =02t 20 1；C0|t_0 12 2 2a = (M m)(J ml ) - m l ;2 2 2/a ;t32 二一m 丨 g /a ; t33 二(J ml )A

14、 / a ; t34 = mlf t42 = (M m)mgl /a ; t43 二ml / a; t44 - -(M m) f /a ;2b3 =(J ml )/ a; b4 二-ml / a。本论文中所用倒立摆的实际参数为 M = 1.096kg , m = 0.109kg，卩二0.100kg/s，l = 0.250m , J= 0.0034kg m2, f = 0.00218kg m2/s。8.1.4倒立摆系统的T-S模糊控制器设计T-S模糊控制器的设计包括控制器的结构设计和参数设计。控制器的结构设计包括：1） 控制器的输入变量和输出变量；2） 模糊规则的数目；3） 前件语言变量的隶属函

15、数的类型;控制器的参数包括：1） 前件参数；2） 后件参数。1.输入变量和输出变量的确定我们使用T-S模糊控制器控制实际倒立摆系统T-S模糊控制器与倒立摆组成的闭环控制系统如图 8-4所示。倒立摆系统有4个状态变量：小车位移x，摆杆角度二，小车速度x，摆杆角速度二，即图8-4中x , angle, xdot, angledot,这些状态变量为T-S模糊控制器的输入。输出变量为小车所受的控制力T-S LogicController图8-4倒立摆的闭环控制系统2.模糊规则数目的确定在T-S模糊控制器中，将这4个状态变量的语言变量分别采用“正”、“负”个 语言值，从而组合出24 = 16条模糊控制规

16、则。每条规则采用 T-S型模糊规则。3.前件隶属函数形状的的确定我们可以采用高斯型函数作为前件隶属函数的类型。于是，T-S模糊控制器的结构已完全确定。4.T-S模糊控制器参数的确定若前件语言变量的隶属度函数采用高斯型函数， 可用中心和方差2个参数对其描述，则每个语言变量对应 4个参数，模糊控制器总共有 4个语言变量，因 此隶属函数对应16个前件参数。再加上每一个规则后件有 5个参数，对应16条规则，共有16x5=80个后件参数，需要优化设计的控制器参数总共有个 96。可以使用遗传算法，神经网络的算法优化这些参数。参数优化的过程见图8-5.图8-5模糊控制器参数的优化过程5.生成T-S模糊推理系

17、统通过上述方法生成16条规则的T-S模糊模型。图8-6表示生成的模糊控制器的输入变量：摆杆角度，摆杆角速度，小车位移，小车速度的隶属函数图8-6 T-S模糊模型输入变量的隶属函数16条规则的T-S模糊模型为:1.If x=in 1mf1 and xdot=in2mf1 and angle=in3mf1 and angledot=in4mf1,then y1=InputVector * -25.28 -21.07 69.55 11.81 0.3678 T2.If x=i n1mf1 a nd xdot=i n2mf1 a nd an gle=i n3mf1 a nd an gledot=i n4

18、 mf2,then y2=InputVector * -38.26 -20.69 69 11.21 -2.511 T3.If x=i n1mf1 a nd xdot=i n2mf1 a nd an gle=i n3mf2 a nd an gledot=in 4mf1,then y3=InputVector * -1.7 -22.37 54.66 11.77 4.813 T15.If x=i n1mf2 a nd xdot=in 2mf2 a nd an gle=i n3mf2 a nd an gledot=in 4mf1,then y15=InputVector * -2.86 -12.29 91.71 10.87 -10.81 T16.If x=i n1mf2 a nd xdot=in 2mf2 a nd an gle=i n3mf2 a nd an gledot=in 4mf2,then y16=InputVector * -15 -15.57 97.42 12.03 -10.24 T其中 InputVector= x, xdot, angle, angledot= x, x 8.1.5倒立摆系统的仿真控制