一个fuzzybp训练的程序，值得初学者入门借鉴.txt

一个fuzzybp训练的程序

clc
clear
U=[1 2 3;2 1 3;3 1 2];%输入矩阵
D=[0.9 0.8 0.7;0.6 0.5 0.4;0.3 0.2 0.1];%期望输出矩阵
E=[0.1 0.3 0.56 0.25 0.31 0.4 0.25 0.25 0.24 0.12 0.13 0.14 0.36 0.56 0.456 0.364 0.58 0.546 0.456 0.43 0.156 0.36 0.23 0.25 0.24 0.26 0.28;
0.1 0.2 0.47 0.4 0.1 0.6 0.5 0.5 0.5 0.8 0.3 0.45 0.13 0.13 0.025 0.16 0.15 0.56 0.16 0.16 0.56 0.8 0.15 0.1465 0.15 0.146 0.465;
0.1 0.3 0.56 0.25 0.31 0.4 0.25 0.25 0.24 0.12 0.13 0.14 0.36 0.56 0.456 0.364 0.58 0.546 0.456 0.43 0.156 0.36 0.23 0.25 0.24 0.26 0.28];%权矩阵(3行27列）
M=[-0.5 0 0.5;
-0.2 0 0.2;
-0.3 0 0.3];%隶属度中值
K=[0.9 0.2 0.1;
0.2 0.3 0.5;
0.3 0.32 0.13];%标准差
q=0.5;%学习率
for h=1:1250
B=[0 0 0];
for w=1:3
for a=1:3
for j=1:3
NET2(a,j)=-(U(w,j)-M(a,j))^2/K(a,j)^2;
Y2(a,j)=exp(NET2(a,j));
end
end%隶属度层结束
n=1;
for b=1:3
for j=1:3
for l=1:3
C(n)=Y2(1,b)*Y2(2,j)*Y2(3,l);
n=n+1;
end
end
end%规则层,解模糊层结束
s1=[0 0 0];
for i=1:3
for j=1:27
s1(i)=s1(i)+E(i,j)*C(j);
end
end
s2=0;
for i=1:27
s2=s2+C(i);
end
for i=1:3
G(i)=s1(i)/s2;
end%解模糊层
for i=1:3
dD(i)=D(w,i)-G(i);
end
for i=1:3
for j=1:27
E(i,j)=E(i,j)+q*dD(i)*C(j)/s2;
end
end
B=B+dD;%误差累计
end
Ba(h)=B(1);
Bb(h)=B(2);
Bc(h)=B(3);
end