%模糊控制器设计
clear all
close all
clc;
a=newfis('fuzzyf'); %建立一个模糊系统，名为a
f1=1;
%定义误差输入的隶属度函数
a=addvar(a,'input','e',[-3*f1,3*f1]);
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=1;
%定义误差输入的隶属度函数
a=addvar(a,'input','ec',[-3*f2,3*f2]);
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=1.5;
%定义输出量的隶属度函数
a=addvar(a,'output','u',[-3*f3,3*f3]);
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; % e ec u 规则使用的权重 前两个变量之间与为1，0为或
1 2 7 1 1;
1 3 6 1 1;
1 4 6 1 1;
1 5 5 1 1;
1 6 5 1 1;
1 7 4 1 1;
2 1 7 1 1;
2 2 7 1 1;
2 3 6 1 1;
2 4 6 1 1;
2 5 5 1 1;
2 6 4 1 1;
2 7 4 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 2 1 1;
4 1 6 1 1;
4 2 5 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 4 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 3 1 1;
7 3 3 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); %生成规则库
a1=setfis(a,'DefuzzMethod','mom');%反模糊化
writefis(a1,'fuzzf');
a2=readfis('fuzzf');
ruleview(a2);
disp('......................................................');
disp(' fuzzy controller table:e=[-3,+3],ec=[-3,+3] ');
disp('......................................................');
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
figure(1)
plotfis(a2); %画模糊控制推理图
figure(2)
plotmf(a,'input',1);
figure(3)
plotmf(a,'input',2);
figure(4)
plotmf(a,'output',1);
- 1
- 2
- 3
- 4
前往页