小弟是初学者,最近在学模糊控制,这是我的程序
%Fuzzy Controller
clear all;
close all;
a=newfis('fuzzf');
f1=6/3;
a=addvar(a,'input','e',[-3*f1,3*f1]); %Parameter e
a=addmf(a,'input',1,'NB','zmf',[-3*f1,-1*f1]);
a=addmf(a,'input',1,'NM','trimf',[-3*f1,-2*f1,0]);
a=addmf(a,'input',1,'NS','trimf',[-3*f1,-1*f1,1*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,3*f1]);
a=addmf(a,'input',1,'PB','smf',[1*f1,3*f1]);
f2=6/3;
a=addvar(a,'input','ec',[-3*f2,3*f2]); %Parameter ec
a=addmf(a,'input',2,'NB','zmf',[-3*f2,-1*f2]);
a=addmf(a,'input',2,'NM','trimf',[-3*f2,-2*f2,0]);
a=addmf(a,'input',2,'NS','trimf',[-3*f2,-1*f2,1*f2]);
a=addmf(a,'input',2,'Z','trimf',[-2*f2,0,2*f2]);
a=addmf(a,'input',2,'PS','trimf',[-1*f2,1*f2,3*f2]);
a=addmf(a,'input',2,'PM','trimf',[0,2*f2,3*f2]);
a=addmf(a,'input',2,'PB','smf',[1*f2,3*f2]);
f3=7/3;
a=addvar(a,'output','u',[-3*f3,3*f3]); %Parameter u
a=addmf(a,'output',1,'NB','zmf',[-3*f3,-1*f3]);
a=addmf(a,'output',1,'NM','trimf',[-3*f3,-2*f3,0]);
a=addmf(a,'output',1,'NS','trimf',[-3*f3,-1*f3,1*f3]);
a=addmf(a,'output',1,'Z','trimf',[-2*f3,0,2*f3]);
a=addmf(a,'output',1,'PS','trimf',[-1*f3,1*f3,3*f3]);
a=addmf(a,'output',1,'PM','trimf',[0,2*f3,3*f3]);
a=addmf(a,'output',1,'PB','smf',[1*f3,3*f3]);
rulelist=[1 1 7 1 1; %Edit rule base
1 2 7 1 1;
1 3 6 1 1;
1 4 6 1 1;
1 5 5 1 1;
1 6 4 1 1;
1 7 4 1 1;
2 1 7 1 1;
2 2 7 1 1;
2 3 6 1 1;
2 4 5 1 1;
2 5 5 1 1;
2 6 4 1 1;
2 7 3 1 1;
3 1 6 1 1;
3 2 6 1 1;
3 3 6 1 1;
3 4 5 1 1;
3 5 4 1 1;
3 6 3 1 1;
3 7 3 1 1;
4 1 6 1 1;
4 2 6 1 1;
4 3 5 1 1;
4 4 4 1 1;
4 5 3 1 1;
4 6 2 1 1;
4 7 2 1 1;
5 1 5 1 1;
5 2 5 1 1;
5 3 4 1 1;
5 4 3 1 1;
5 5 3 1 1;
5 6 2 1 1;
5 7 2 1 1;
6 1 5 1 1;
6 2 4 1 1;
6 3 3 1 1;
6 4 2 1 1;
6 5 2 1 1;
6 6 2 1 1;
6 7 1 1 1;
7 1 4 1 1;
7 2 4 1 1;
7 3 2 1 1;
7 4 2 1 1;
7 5 2 1 1;
7 6 1 1 1;
7 7 1 1 1];
a=addrule(a,rulelist);
%showrule(a) % Show fuzzy rule base
a1=setfis(a,'DefuzzMethod','mom'); % Defuzzy
writefis(a1,'fuzzf'); % save to fuzzy file "fuzzf.fis" which can be
% simulated with fuzzy tool
a2=readfis('fuzzf');
disp('fuzzy controller table:e=[-6,+6],ec=[-6,+6]');
Ulist=zeros(13,13);
for i=1:13
for j=1:13
e(i)=-7+i;
ec(j)=-7+j;
Ulist(i,j)=evalfis([e(i),ec(j)],a2);
end
end
Ulist=ceil(Ulist)
figure(1);
plotfis(a2);
figure(2);
plotmf(a,'input',1);
figure(3);
plotmf(a,'input',2);
figure(4);
plotmf(a,'output',1);
运行后输出:fuzzy controller table:e=[-6,+6],ec=[-6,+6]
Ulist =
8 7 8 6 5 5 5 4 3 2 1 1 1
7 7 7 6 5 5 5 4 3 2 1 1 1
8 7 8 6 5 4 3 3 3 2 1 -2 -2
6 6 6 6 5 4 3 2 2 1 -1 -2 -2
5 5 5 5 5 4 3 2 1 -1 -2 -2 -2
5 5 5 4 4 2 2 1 -1 -1 -3 -3 -3
5 5 5 4 3 2 1 -1 -2 -3 -4 -4 -4
4 4 4 2 2 1 -1 -1 -2 -3 -4 -4 -4
3 3 3 2 1 -1 -2 -2 -2 -3 -4 -4 -4
3 3 2 1 -1 -1 -3 -3 -3 -3 -4 -5 -5
3 3 1 -1 -2 -3 -4 -4 -4 -4 -4 -6 -7
1 1 1 -1 -4 -4 -4 -4 -4 -5 -6 -6 -6
1 1 1 -1 -4 -4 -4 -4 -4 -5 -7 -6 -7
我想问一下,u的论域是[-7,+7],为什么会有8??
当我运行evalfis([e(1),ec(1)],a2)
输出
ans =
7.0000
这又是为什么呢?如何解决?
MATLAB的去模糊方法用的是加权平均还是最大隶属度?
有时候我单独求一个u值的时候,会出来小数,比如
evalfis([e(1),ec(2)],a2)
ans =
6.4400
是不是MATLAB自带的模糊程序有取整功能?
[本帖最后由 zhihuihaiwan 于 2009-4-14 19:16 编辑]
本帖子中包含更多资源
您需要 登录 才可以下载或查看,没有帐号?注册
x