1.创建新的模糊推理系统。
2.随后添加变量e和ec,并添加隶属度函数。
3.添加输出u的模糊语言变量并添加隶属函数。
4.添加模糊规则库。
5.利用a=addrule添加模糊规则函数,showwrule用于显示模糊规则函数。推理得出Ulist
代码
%模糊控制器设计
a=newfis('fuzzf'); %创建新的模糊推理系统
%输入1
f1=10;
a=addvar(a,'input','e',[-30*f1,30*f1]);
%添加 e 的模糊语言变量
a=addmf(a,'input',1,'NB','zmf',[-30*f1,-10*f1]);
%添加 e 的模糊语言变量的隶属度函数(z型)
a=addmf(a,'input',1,'NM','trimf',[-30*f1,-20*f1,0]);
%隶属度函数为三角形
a=addmf(a,'input',1,'NS','trimf',[-30*f1,-10*f1,10*f1]);
a=addmf(a,'input',1,'Z','trimf',[-20*f1,0,20*f1]);
a=addmf(a,'input',1,'PS','trimf',[-10*f1,10*f1,30*f1]);
a=addmf(a,'input',1,'PM','trimf',[0,20*f1,30*f1]);
a=addmf(a,'input',1,'PB','smf',[10*f1,30*f1]);
%输入2
f2=10;
a=addvar(a,'input','ec',[-30*f2,30*f2]);
%添加 ec 的模糊语言变量
a=addmf(a,'input',2,'NB','zmf',[-30*f2,-10*f2]);
a=addmf(a,'input',2,'NM','trimf',[-30*f2,-20*f2,0]);
a=addmf(a,'input',2,'NS','trimf',[-30*f2,-10*f2,10*f2]);
a=addmf(a,'input',2,'Z','trimf',[-20*f2,0,20*f2]);
a=addmf(a,'input',2,'PS','trimf',[-10*f2,10*f2,30*f2]);
a=addmf(a,'input',2,'PM','trimf',[0,20*f2,30*f2]);
a=addmf(a,'input',2,'PB','smf',[10*f2,30*f2]);
%输出
f3=3;
a=addvar(a,'output','u',[-30*f3,3*f3]);
%添加 u 的模糊语言变量
a=addmf(a,'output',1,'NB','zmf',[-30*f3,-10*f3]);
a=addmf(a,'output',1,'NM','trimf',[-30*f3,-20*f3,0]);
a=addmf(a,'output',1,'NS','trimf',[-30*f3,-10*f3,10*f3]);
a=addmf(a,'output',1,'Z','trimf',[-20*f3,0,20*f3]);
a=addmf(a,'output',1,'PS','trimf',[-10*f3,10*f3,30*f3]);
a=addmf(a,'output',1,'PM','trimf',[0,20*f3,30*f3]);
a=addmf(a,'output',1,'PB','smf',[10*f3,30*f3]);
%规则库
rulelist=[1 1 1 1 1; %编辑模糊规则,后俩个数分别是规则权重和AND OR选
2 1 3 1 1;
3 1 4 1 1;
4 1 2 1 1;
5 1 3 1 1;
6 1 4 1 1;
7 1 3 1 1;
1 2 1 1 1;
2 2 3 1 1;
3 2 4 1 1;
4 2 2 1 1;
5 2 3 1 1;
6 2 4 1 1;
7 2 3 1 1;
1 3 1 1 1;
2 3 3 1 1;
3 3 4 1 1;
4 3 2 1 1;
5 3 3 1 1;
6 3 4 1 1;
7 3 3 1 1;
1 4 1 1 1;
2 4 3 1 1;
3 4 4 1 1;
4 4 2 1 1;
5 4 3 1 1;
6 4 4 1 1;
7 4 3 1 1;
1 5 1 1 1;
2 5 3 1 1;
3 5 4 1 1;
4 5 2 1 1;
5 5 3 1 1;
6 5 4 1 1;
7 5 3 1 1;
1 6 1 1 1;
2 6 3 1 1;
3 6 4 1 1;
4 6 2 1 1;
5 6 3 1 1;
6 6 4 1 1;
7 6 3 1 1;
1 7 1 1 1;
2 7 3 1 1;
3 7 4 1 1;
4 7 2 1 1;
5 7 3 1 1;
6 7 4 1 1;
7 7 3 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)=-4+i;
ec(j)=-4+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);
结果截图