计算智能--模糊控制器
文章目录
- 一、模糊控制器概述
- 二、模糊控制器的设计步骤
- 三、模糊控制器设计举例及matlab代码实现
- 四、运行结果及源代码
一、模糊控制器概述
在确定性控制系统中，根据控制器输出的个数，可分为单变量控制系统和多变量控制系统。在模糊控制系统中也可类似地划分为单变量模糊控制和多变量模糊控制。
二、模糊控制器的设计步骤
1、确定模糊控制器的结构
2、定义输入输出模糊集、隶属函数
3、建立模糊控制规则
4、模糊推理
5、反模糊化
三、模糊控制器设计举例及matlab代码实现
1、确定模糊控制器的结构（两个输入e和ec，一个输出u），即多变量模糊控制器。
2、定义输入输出模糊集、隶属函数
输入e的模糊集定义及隶属函数（e的输入范围为-9到9，有NB、NM、NS、Z、PS、PM、PB七个隶属函数）。
输入eu的模糊集定义及隶属函数（eu的输入范围为-9到9，有NB、NM、NS、Z、PS、PM、PB七个隶属函数）。
输出u的模糊集定义及隶属函数（u有NB、NM、NS、Z、PS、PM、PB七个隶属函数）。
3、建立模糊控制规则
对输入（e、eu）、输出（u）的隶属函数建立模糊规则。
模糊规则解释（本文中定义了36条模糊规则）：
2.If (e is NB) and (ec is NM) then (u is NB) (1)
如果输入e属于隶属函数NB，输入ec属于隶属函数NM，则输出u属于隶属函数NB。
10.If (e is NM) and (ec is PB) then (u is PS) (1)
如果输入e属于隶属函数NM，输入ec属于隶属函数PB，则输出u属于隶属函数PS。
4、模糊推理
通过模糊推理算法和模糊控制规则进行推理，得到最终的控制量。
5、反模糊化
将模糊推理结果转化为精确值。常用的反模糊化方法有三种，最大隶属度法、重心法、加权平均法，本例采用的是重心法。
四、运行结果及源代码
输入e和输入ec，经过模糊推理计算得到输出u：
输入e的7个隶属度函数图像：
输入ec的7个隶属度函数图像：
输出u的7个隶属度函数图像：
经过模糊推理得到的模糊矩阵：
源代码：
%模糊控制器设计
a=newfis('fuzzf'); %创建新的模糊推理系统
%输入1
f1=1;
a=addvar(a,'input','e',[-9*f1,9*f1]);
%添加 e 的模糊语言变量
a=addmf(a,'input',1,'NB','zmf',[-2*f1,3*f1]);
%添加 e 的模糊语言变量的隶属度函数（z型）
a=addmf(a,'input',1,'NM','trimf',[-4*f1,-1*f1,5*f1]);
%隶属度函数为三角形
a=addmf(a,'input',1,'NS','trimf',[-5*f1,3*f1,4*f1]);
a=addmf(a,'input',1,'Z','trimf',[-2*f1,0,2*f1]);
a=addmf(a,'input',1,'PS','trimf',[-1*f1,1*f1,3*f1]);
a=addmf(a,'input',1,'PM','trimf',[0,2*f1,6*f1]);
a=addmf(a,'input',1,'PB','smf',[6*f1,9*f1]);
%输入2
f2=1;
a=addvar(a,'input','ec',[-9*f2,9*f2]);
%添加 ec 的模糊语言变量
a=addmf(a,'input',2,'NB','zmf',[-5*f2,4*f2]);
a=addmf(a,'input',2,'NM','trimf',[-3*f2,3*f2,6*f2]);
a=addmf(a,'input',2,'NS','trimf',[-2*f2,4*f2,5*f2]);
a=addmf(a,'input',2,'Z','trimf',[-5*f2,0,3*f2]);
a=addmf(a,'input',2,'PS','trimf',[-3*f2,1*f2,4*f2]);
a=addmf(a,'input',2,'PM','trimf',[0,2*f2,7*f2]);
a=addmf(a,'input',2,'PB','smf',[2*f2,8*f2]);
%输出
f3=1.5;
a=addvar(a,'output','u',[-9*f3,9*f3]);
%添加 u 的模糊语言变量
a=addmf(a,'output',1,'NB','zmf',[-6*f3,3*f3]);
a=addmf(a,'output',1,'NM','trimf',[-4*f3,0,1*f3]);
a=addmf(a,'output',1,'NS','trimf',[-2*f3,1*f3,6*f3]);
a=addmf(a,'output',1,'Z','trimf',[-3*f3,2*f3,3*f3]);
a=addmf(a,'output',1,'PS','trimf',[-1*f3,4*f3,5*f3]);
a=addmf(a,'output',1,'PM','trimf',[0,3*f3,6*f3]);
a=addmf(a,'output',1,'PB','smf',[5*f3,7*f3]);
%规则库
rulelist=[1 1 1 1 1; %编辑模糊规则，后俩个数分别是规则权重和AND OR选项
1 2 1 1 1;
1 3 2 1 1;
1 4 2 1 1;
1 6 3 1 1;
1 7 7 1 1;
2 1 2 1 1;
2 2 6 1 1;
2 3 3 1 1;
2 7 5 1 1;
3 1 1 1 1;
3 2 3 1 1;
3 3 4 1 1;
3 4 2 1 1;
3 6 4 1 1;
3 7 5 1 1;
4 1 1 1 1;
4 4 2 1 1;
4 5 3 1 1;
4 6 3 1 1;
4 7 5 1 1;
5 1 2 1 1;
5 2 4 1 1;
5 5 3 1 1;
5 6 6 1 1;
5 7 5 1 1;
6 1 1 1 1;
6 2 3 1 1;
6 3 2 1 1;
6 6 4 1 1;
6 7 5 1 1;
7 1 1 1 1;
7 2 2 1 1;
7 5 3 1 1;
7 6 4 1 1;
7 7 5 1 1;
];
a=addrule(a,rulelist); %添加模糊规则函数
showrule(a) %显示模糊规则函数
a1=setfis(a,'DefuzzMethod','centroid'); %设置解模糊方法
writefis(a1,'fuzzf'); %保存模糊系统
a2=readfis('fuzzf'); %从磁盘读出保存的模糊系统
disp('fuzzy Controller table:e=[-3,+3],ec=[-3,+3]');%显示矩阵和数组内容
%推理
Ulist=zeros(7,7); %全零矩阵
for i=1:7
for j=1:7
e(i)=-3+i;
ec(j)=-2+j;
Ulist(i,j)=evalfis([e(i),ec(j)],a2); %完成模糊推理计算
end
end
% Ulist=ceil(Ulist) %朝正无穷方向取整
Ulist %朝正无穷方向取整
%画出模糊系统
figure(1); plotfis(a2);
figure(2);plotmf(a,'input',1);
figure(3);plotmf(a,'input',2);
figure(4);plotmf(a,'output',1);