MATLAB环境下的模糊推理程序化方法

鲁棒性强、能够很好地克服传动系统中模型参数变化和非线性等不确定因素的优点，因此，模糊控制在实际控制系统中得到广泛应用。模糊控制器普遍采用的是离线进行模糊推理产生的模糊控制表，然后把控制表存储在单片机中进行在线查表控制，但模糊推理计算量大，这给模糊控制器的设计和调整带来了许多困难。此外，模糊控制存在的一个较大的缺点是稳态精度不高。为此，人们提出了许多改进的方法[１]，其中最直接的方法是增加模糊量化论域。但是，当量化论域中元素、控制规则的条数很多时，模糊推理计算控制表的运算量大大增加。本文在ＭＡＴＬＡＢ环境下，研究了模糊推理算法，研制了模糊推理程序。只要给出模糊控制器的输入和输出量的隶属函数矩阵、控制规则矩阵、模糊量的特征向量矩阵，运行该程序就可得到模糊控制表。该软件使用方便，为模糊控制的研究和应用提供了方便的条件。

１ 模糊推理概述

    在模糊控制中，关键是要求得模糊关系矩阵，应用最多的是离线进行模糊推理产生控制表，模糊推理是模糊控制器的核心。这里以两入一出、采用ＣＲＩ推理方法及重心法解模糊的模糊控制器为例，简要介绍模糊推理过程。

    设模糊集合 Ａｉ ＝（ａ１ｉ．．．ａｎｉ）∈Ｆ(Ｘ)，Ｂｉ ＝（ｂ１ｉ．．．ｂｍｉ）∈Ｆ(Ｙ)，Ｃ ＝（ｃ１ｉ．．．ｃｔｉ）∈Ｆ(Ｚ) (ｉ＝１,．．．,ｐ)，模糊规则为“ＩＦ Ａｉ ＡＮＤ Ｂｉ ＴＨＥＮ Ｃｉ”（ｉ＝１，．．．,ｐ）。其中ｐ为模糊变量的个数，Ｆ（Ｘ）、Ｆ（Ｙ）、Ｆ（Ｚ）分别为论域Ｘ、Ｙ、Ｚ上的模糊集。由此可求得总模糊关系：

    对某一模糊特征向量Ａ*、Ｂ*，进一步可求出输出的模糊量：

    最后，通过解模糊判决可求出输出的精确量:

    注：Ａｉ×Ｂｉ 是模糊向量的笛卡儿积，由公式（４）计算：

    其中,“Ｔ”表示转置，“ｏ”表示矩阵合成运算，使用最常见的 “∨—∧”算子计算，见公式：

    这里需要指出的是，把笛卡儿乘积看作一个ｎ×１与一个１×ｍ的矩阵合成。计算结果是ｎ×ｍ模糊阵，而（Ａｉ×Ｂｉ）Ｔ１表示把这个ｎ×ｍ模糊阵按行“拉直”成ｎｍ元模糊行向量，再转置成ｎｍ元模糊列向量；（Ａ*×Ｂ*）Ｔ２的含义是：把Ａ*×Ｂ*这个ｎ×ｍ模糊阵按行“拉直”成ｎｍ元模糊行向量。

２ 基于ＭＡＴＬＡＢ实现模糊推理程序设计

    在ＭＡＴＬＡＢ环境下，充分利用其矩阵处理能力[２]，可实现“∨—∧”、“∨—· ”及 “∧—”模糊算子运算、重心法（或最大隶属度法）解模糊。由于模糊推理过程涉及到合成、求模糊关系及解模糊运算等，采用自定义函数实现各个功能的模块化。图１为主函数程序框图；图２为求模糊关系程序框图；图３为解模糊程序框图；图４为合成程序框图。

    运行此程序，只要输入矩阵Ａ和Ｂ及相应矩阵Ｃ、控制矩阵Ｑ、特征向量矩阵Ａ*和Ｂ*及输出量化论域矩阵ＱＦ，就能得到模糊控制表。

３ 计算实例

３.１ 常规模糊控制器

    下面以常规模糊控制器为例来验证此程序的正确性。在常规模糊控制器中，输入语言变量偏差Ｅ、偏差变化率ＣＥ、输出语言变量控制量Ｕ的赋值表及Ｅ、ＣＥ→Ｕ的控制规则表见文献[３]。文献[３]中的原控制表如表１所示，利用本文编制的模糊推理程序（采用“∨—∧”算子、最大隶属度解模糊）得到的控制表如表２所示。

    在表１、表２中最上面一行为输入量误差Ｅ的量化等级，最左一列为输入量误差变化ＣＥ的量化等级。表１与表２中的数据存在差异，这是因为表１中的带*的数据均是修正后的结果。这一点可以通过特征展开近似推理方法[４]进行验证。具体验证方法如下：

表1 E、CE→U 控制表

	-6	-5	-4	-3	-2	-1	-0	+0	1	2	3	4	5	6
-6	7	6*	7	7*	4	4	4	4	2*	1	0	0	0	0
-5	6*	6*	6*	6*	4	4	4	4	2*	2*	0	0	0	0
-4	7	6*	7	6*	4	4	4	4	2*	1	0	0	0	0
-3	6*	6*	6*	6*	5*	5*	5*	5*	2*	2*	0*	-2*	-2*	-2*
-2	7	6*	7	6*	4	4	1	1	0	0	-3*	-4	-4	-4
-1	7	6*	7	6*	4	4	1	1	0	-3*	-3*	-4	-4	-4
0	7	6*	7	6*	4	1*	0	0	-1*	-4	-6*	-7	-6*	-7
1	4	4	4	3*	1*	0	-1	-1	-4	-4	-6*	-7	-6*	-7
2	4	4	4	2*	0	0	-1	-1	-4	-4	-6*	-7	-6*	-7
3	2*	2*	2*	0*	0*	0*	-1*	-1*	-3*	-3*	-6*	-6*	-6*	-6*
4	0	0	0	-1*	-1	-3*	-4	-4	-4	-4	-6*	-7	-6*	-7
5	0	0	0	-1*	-1*	-2*	-4	-4	-4	-4	-6*	-6*	-6*	-6*
6	0	0	0	-1*	-1	1*	-4	-4	-4	-4	-6*	-7	-6*	-7

表2 E、CE→U 控制表

	-6	-5	-4	-3	-2	-1	-0	+0	1	2	3	4	5	6
-6	7	6.5	7	6.5	4	4	4	4	1.5	1	0	0	0	0
-5	6.5	6.5	6.5	6.5	4	4	4	4	1.5	1.5	0	0	0	0
-4	7	6.5	7	6.5	4	4	4	4	1.5	1	0	0	0	0
-3	6.5	6.5	6.5	6.5	4	4	3	3	1	1	-3	-3	-3	-3
-2	7	6.5	7	6.5	4	4	1	1	0	0	-4	-4	-4	-4
-1	7	6.5	7	6.5	4	4	1	1	0	0	-4	-4	-4	-4
0	7	6.5	7	6.5	4	4	0	0	-4	-4	-6.5	-7	-6.5	-7
1	4	4	4	4	0	0	-1	-1	-4	-4	-6.5	-7	-6.5	-7
2	4	4	4	4	0	0	-1	-1	-4	-4	-6.5	-7	-6.5	-7
3	3	3	3	3	-1	-1	-3	-3	-4	-4	-6.5	-6.5	-6.5	-6.5
4	0	0	0	0	-1	-1.5	-4	-4	-4	-4	-6.5	-7	-6.5	-7
5	0	0	0	0	-1.5	-1.5	-4	-4	-4	-4	-6.5	-6.5	-6.5	-6.5
6	0	0	0	0	-1	-1.5	-4	-4	-4	-4	-6.5	-7	-6.5	-7

以表１的第１０行第７列为例，即Ｅ*＝（０,０,０,０,０,０,１,０,０,０,０,０,０,０）,ＣＥ*＝（０,０,０,０,０,０,０,０,０,１,０,０,０）为其特征向量,由特征展开近似推理方法[４]的式:

    可计算出α３＝０．３,α４＝１,β５＝０．７,β６＝０．７,β７＝０．１,(α３∧β５）Ｕ４＝０．３Ｕ４,（α３∧β６）Ｕ３＝０．３Ｕ３,（α３∧β７）Ｕ３＝０．１Ｕ３,(α４∧β５）Ｕ３＝０．７Ｕ３,（α４∧β６）Ｕ２＝０．７Ｕ２,（α４∧β７）Ｕ２＝０．１Ｕ２,由其余５０条规则得到的αｉ或βｊ总有一个为０，故这一对模糊输入得到的模糊输出为：

    Ｕ＝０．３Ｕ４∪０．７Ｕ３∪０．３Ｕ３∪０．１Ｕ３∪０．７Ｕ２∪０．１Ｕ２
    ＝(０,０．１４,０．４９,０．７,０．４９,０．５６,０．７,０．３,０．１５,０,０,０,０,０,０)

    经过解模糊得到精确输出为：

    ｕ＝((－４)＋(－１))／２＝－２．５

３.２ 增加模糊量化论域的模糊控制器

    以增加模糊量化论域的模糊控制器为例进行模糊推理，计算出控制表。取误差Ｅ、误差变化ＣＥ和控制量Ｕ的量化论域均为：

    {－１０，－９，－８，－７，－６，－５，－４，－３，－２，－１，０，１，２，３，４，５，６，７，８，９，１０}

    取误差Ｅ、误差变化ＣＥ和控制量Ｕ的模糊论域均为：

    {ＮＶＢ，ＮＢ，ＮＭ，ＮＳ，ＮＶＳ，ＺＯ，ＰＶＳ，ＰＳ，ＰＭ，ＰＢ，ＰＶＢ}

    误差Ｅ、误差变化ＣＥ和控制量Ｕ的隶属函数均如图５所示。

    借鉴常规模糊控制器设计经验，得到语言控制规则表,如表３所示。表中共有１２１条控制规则，其中一些规则可以合并。但利用计算机进行推理计算，这些规则就没必要合并了。

表3 E、CE→U 模糊控制规则表

E    U CE	NVB	NB	NM	NS	NVS	ZO	PVS	PS	PM	PB	PVB
NVB	PVB	PVB	PB	PM	PM	PS	PS	PVS	ZO	ZO	ZO
NB	PVB	PB	PB	PM	PM	PS	PS	PVS	ZO	ZO	ZO
NM	PB	PB	PM	PM	PS	PS	PVS	PVS	ZO	ZO	ZO
NS	PB	PM	PM	PS	PS	PVS	PVS	ZO	NVS	NVS	NVS
NVS	PM	PM	PS	PS	PVS	PVS	ZO	NVS	NVS	NS	NS
ZO	PM	PS	PS	PVS	PVS	ZO	NVS	NVS	NS	NS	NS
PVS	PS	PS	PVS	PVS	ZO	NVS	NVS	NS	NS	NM	NM
PS	PS	PVS	PVS	ZO	NVS	NVS	NS	NS	NM	NM	NM
PM	PVS	PVS	ZO	NVS	NVS	NS	NS	NM	NM	NB	NB
PB	PVS	ZO	NVS	NVS	NS	NS	NM	NM	NB	NB	NVB
PVB	ZO	NVS	NVS	NS	NS	NM	NM	NB	NB	NVB	NVB

    利用研制的模糊推理计算程序，采用“∨—？”得到的Ｅ，ＣＥ→Ｕ控制表，如表４所示。在表中最上面一行为误差Ｅ的量化等级，最左一列为误差变化ＣＥ的量化等级。

表4 控制表

u	-10	-9	-8	-7	-6	-5	-4	-3	-2	-1	0	1	2	3	4	5	6	7	8	9	10
-10	9.7	9.7	9.7	9.7	8	7	6	6	6	5	4	4	4	3	2	1	0	0	0	0	0
-9	9.7	8.7	8.7	8.7	8	7	6	6	6	5	4	4	4	3	2	1	0	0	0	0	0
-8	9.7	8.7	8	8	8	7	6	6	6	5	4	4	4	3	2	1	0	0	0	0	0
-7	8.7	8.7	8	7	7	7	6	5	5	5	4	3	3	3	2	1	0	0	0	0	0
-6	8	8	8	7	6	6	6	5	4	4	4	3	2	2	2	1	0	0	0	0	0
-5	8	7	7	7	6	5	5	5	4	3	3	3	2	1	1	0	-1	-1	-1	-1	-1
-4	8	7	6	6	6	5	4	4	4	3	2	2	2	1	0	-1	-2	-2	-2	-2	-2
-3	7	7	6	5	5	5	4	3	3	3	2	1	1	0	-1	-1	-2	-3	-3	-3	-3
-2	6	6	6	5	4	4	4	3	2	2	2	1	0	-1	-2	-2	-2	-3	-4	-4	-4
-1	6	5	5	5	4	3	3	3	2	1	1	0	-1	-1	-2	-3	-3	-3	-4	-4	-4
0	6	5	4	4	4	3	2	2	2	1	0	-1	-2	-2	-2	-3	-4	-4	-4	-4	-4
1	5	5	4	3	3	3	2	1	1	0	-1	-1	-2	-3	-3	-3	-4	-5	-5	-5	-5
2	4	4	4	3	3	2	2	1	0	-1	-1	-1	-1	-3	-4	-4	-4	-5	-6	-6	-6
3	4	3	3	3	2	1	1	0	-1	-1	-2	-3	-3	-3	-4	-5	-5	-5	-6	-6	-6
4	4	3	2	2	2	1	0	-1	-2	-2	-2	-3	-4	-4	-4	-5	-6	-6	-6	-6	-6
5	3	3	2	1	1	0	-1	-1	-2	-3	-3	-3	-4	-5	-5	-5	-6	-7	-7	-7	-7
6	2	2	2	1	0	-1	-2	-2	-2	-3	-4	-4	-4	-5	-6	-6	-6	-7	-8	-8	-8
7	2	1	1	0	-1	-1	-2	-3	-3	-3	-4	-5	-5	-5	-6	-7	-7	-7	-8	-8.7	-8.7
8	2	1	0	-1	-2	-2	-2	-3	-4	-4	-4	-5	-6	-6	-6	-7	-8	-8	-8	-8.7	-9.7
9	1	0	-1	-1	-2	-3	-3	-3	-4	-5	-5	-5	-6	-7	-7	-7	-8	-8.7	-8.7	-8.7	-9.7
10	0	-1	-2	-2	-2	-3	-4	-4	-4	-5	-6	-6	-6	-7	-8	-8	-8	-8.7	-9.7	-9.7	-9.7

４ 仿真结果

    这里以对象

    为例，分别使用表２、表４中的数据构成模糊控制器，响应曲线如图６中曲线１、２所示，其中采样周期 Ｔ＝０．１ｓ。从图６中可以看出，增加模糊量化论域的模糊控制器的稳态精度明显提高。

    此程序方便地实现了模糊推理的功能，免去了手工计算的麻烦。而且，不论变量的量化等级为多大，只要给定模糊控制器的输入和输出量的隶属函数矩阵、控制规则矩阵、模糊量的特征向量矩阵，运行该软件就可快速地获得相应的模糊控制表。另外，该软件可实现“∨—∧”、“∨—· ”、 “—∧”模糊算子运算及重心法、最大隶属度原则模糊判决方法。因此，本文为模糊控制的研究和应用提供了方便的条件。

编辑： 引用地址：http://www.eeworld.com.cn/designarticles/sensor/200605/1421.html

本网站转载的所有的文章、图片、音频视频文件等资料的版权归版权所有人所有，本站采用的非本站原创文章及图片等内容无法一一联系确认版权者。如果本网所选内容的文章作者及编辑认为其作品不宜公开自由传播，或不应无偿使用，请及时通过电子邮件或电话通知我们，以迅速采取适当措施，避免给双方造成不必要的经济损失。