模糊控制运算

一、例子介绍
输入量为两个，手长和脚长 。通过模糊逻辑控制得出变态程度的输出；
手长和脚长都是从70cm+n*10cm到130cm；
变态程度为1到7种；
模糊规则为总和140和150 为1；160和170为2。。。以此类推；
当出现手脚差值大于30cm变态程度转化为1；
二、代码实现
%模糊控制器设计
a=newfis(‘fuzzf’); %创建新的模糊推理系统

%输入1
f1=1;
a=addvar(a,‘input’,‘手长’,[0,100f1]);
%添加 e 的模糊语言变量
a=addmf(a,‘input’,1,‘70cm’,‘zmf’,[0,20f1]);
%添加 e 的模糊语言变量的隶属度函数（z型）
a=addmf(a,‘input’,1,‘80cm’,‘trimf’,[20f1,25f1,30f1]);
%隶属度函数为三角形
a=addmf(a,‘input’,1,‘90cm’,‘trimf’,[30f1,35f1,40f1]);
a=addmf(a,‘input’,1,‘100cm’,‘trimf’,[40f1,45,50f1]);
a=addmf(a,‘input’,1,‘110cm’,‘trimf’,[50f1,60f1,70f1]);
a=addmf(a,‘input’,1,‘120cm’,‘trimf’,[70,80f1,85f1]);
a=addmf(a,‘input’,1,‘130cm’,‘smf’,[85f1,100*f1]);

%输入2
f2=1;
a=addvar(a,‘input’,‘脚长’,[0,100f2]);
%添加 ec 的模糊语言变量
a=addmf(a,‘input’,2,‘70cm’,‘zmf’,[0f2,20f2]);
a=addmf(a,‘input’,2,‘80cm’,‘trimf’,[20f2,25f2,30f1]);
a=addmf(a,‘input’,2,‘90cm’,‘trimf’,[30f2,35f2,40f2]);
a=addmf(a,‘input’,2,‘100cm’,‘trimf’,[40f2,45f2,50f2]);
a=addmf(a,‘input’,2,‘110cm’,‘trimf’,[50f2,60f2,70f2]);
a=addmf(a,‘input’,2,‘120cm’,‘trimf’,[70,80f2,85f2]);
a=addmf(a,‘input’,2,‘130cm’,‘smf’,[85f2,100*f2]);

%输出
f3=1.5;
a=addvar(a,‘output’,‘变态程度’,[0,100f2]);
%添加 u 的模糊语言变量
a=addmf(a,‘output’,1,‘rabbish’,‘zmf’,[0f2,20f2]);
a=addmf(a,‘output’,1,‘just soso’,‘trimf’,[20f2,25f2,30f1]);
a=addmf(a,‘output’,1,‘a little’,‘trimf’,[30f2,35f2,40f2]);
a=addmf(a,‘output’,1,‘hai xing’,‘trimf’,[40f2,45f2,50f2]);
a=addmf(a,‘output’,1,‘good’,‘trimf’,[50f2,60f2,70f2]);
a=addmf(a,‘output’,1,‘very good’,‘trimf’,[70,80f2,85f2]);
a=addmf(a,‘output’,1,‘nb’,‘smf’,[85f2,100*f2]);

%规则库
rulelist=[1 1 1 1 1; %编辑模糊规则，后俩个数分别是规则权重和AND OR选项
1 2 1 1 1;
1 3 1 1 1;
1 4 1 1 1;
1 5 1 1 1;
1 6 1 1 1;
1 7 1 1 1;

           2 1 1 1 1;
           2 2 2 1 1;
           2 3 2 1 1;
           2 4 3 1 1;
           2 5 3 1 1;
           2 6 1 1 1;
           2 7 1 1 1;
           
           3 1 1 1 1;
           3 2 2 1 1;
           3 3 3 1 1;
           3 4 3 1 1;
           3 5 4 1 1;
           3 6 4 1 1;
           3 7 1 1 1;
           
           4 1 1 1 1;
           4 2 3 1 1;
           4 3 3 1 1;
           4 4 4 1 1;
           4 5 4 1 1;
           4 6 5 1 1;
           4 7 5 1 1;
           
           5 1 1 1 1;
           5 2 3 1 1;
           5 3 4 1 1;
           5 4 4 1 1;
           5 5 5 1 1;
           5 6 5 1 1;
           5 7 6 1 1;
           
           6 1 1 1 1;
           6 2 1 1 1;
           6 3 4 1 1;
           6 4 5 1 1;
           6 5 5 1 1;
           6 6 6 1 1;
           6 7 6 1 1;
           
           
           7 1 1 1 1;
           7 2 1 1 1;
           7 3 1 1 1;
           7 4 5 1 1;
           7 5 6 1 1;
           7 6 6 1 1;
           7 7 7 1 1;];

a=addrule(a,rulelist); %添加模糊规则函数
showrule(a) %显示模糊规则函数
a1=setfis(a,‘DefuzzMethod’,‘centroid’); %设置解模糊方法
writefis(a1,‘fuzzf’); %保存模糊系统
a2=readfis(‘fuzzf’); %从磁盘读出保存的模糊系统
disp(‘fuzzy Controller table:e=[0,100],ec=[0,100]’);%显示矩阵和数组内容

%推理
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);
三、运行结果