matlab 编写的模糊控制程序

[neumoon]

我想编写一个模糊控制程序
输出模糊控制表
代码如下
clear all;

PB=[0 0 0 0 0 0 0 0 0 0 0 0.5 1];
PM=[0 0 0 0 0 0 0 0 0 0.5 1 0.5 0];
PS=[0 0 0 0 0 0 0 0.5 1 0.5 0 0 0];
ZO=[0 0 0 0 0 0.5 1 0.5 0 0 0 0 0];
NS=[0 0 0 0.5 1 0.5 0 0 0 0 0 0 0];
NM=[0 0.5 1 0.5 0 0 0 0 0 0 0 0 0];
NB=[1 0.5 0 0 0 0 0 0 0 0 0 0 0];

Z=[-6 -5 -4 -3 -2 -1 0 1 2 3 4 5 6 ];

for i=1:size (NB,2)
Aa= zeros (1,size(NB,2));
Aa(1,i)=1;
for j=1:size(NB,2)
Bb=zeros (1,size(NB,2));
Bb(1,j)=1;

A=NB';B=NB';C=PB;
[R1]=fun(A,B,C,Aa,Bb);
R=R1;
A=NB';B=NM';C=PB;
[R1]=fun(A,B,C,Aa,Bb);
R=max(R,R1);
A=NB';B=NS';C=PB;
[R1]=fun(A,B,C,Aa,Bb);
R=max(R,R1);
A=NB';B=ZO';C=PB;
[R1]=fun(A,B,C,Aa,Bb);
R=max(R,R1);
A=NB';B=PS';C=PM;
[R1]=fun(A,B,C,Aa,Bb);
R=max(R,R1);
A=NB';B=PM';C=ZO;
[R1]=fun(A,B,C,Aa,Bb);
R=max(R,R1);
A=NB';B=PB';C=ZO;
[R1]=fun(A,B,C,Aa,Bb);
R=max(R,R1);
A=NM';B=NB';C=PB;
[R1]=fun(A,B,C,Aa,Bb);
R=max(R,R1);
A=NM';B=NM';C=PB;
[R1]=fun(A,B,C,Aa,Bb);
R=max(R,R1);
A=NM';B=NS';C=PB;
[R1]=fun(A,B,C,Aa,Bb);
R=max(R,R1);
A=NM';B=ZO';C=PB;
[R1]=fun(A,B,C,Aa,Bb);
R=max(R,R1);
A=NM';B=PS';C=PM;
[R1]=fun(A,B,C,Aa,Bb);
R=max(R,R1);
A=NM';B=PM';C=ZO;
[R1]=fun(A,B,C,Aa,Bb);
R=max(R,R1);
A=NM';B=PB';C=ZO;
[R1]=fun(A,B,C,Aa,Bb);
R=max(R,R1);
A=NS';B=NB';C=PM;
[R1]=fun(A,B,C,Aa,Bb);
R=max(R,R1);
A=NS';B=NM';C=PM;
[R1]=fun(A,B,C,Aa,Bb);
R=max(R,R1);
A=NS';B=NS';C=PM;
[R1]=fun(A,B,C,Aa,Bb);
R=max(R,R1);
A=NS';B=ZO';C=PM;
[R1]=fun(A,B,C,Aa,Bb);
R=max(R,R1);
A=NS';B=PS';C=ZO;
[R1]=fun(A,B,C,Aa,Bb);
R=max(R,R1);
A=NS';B=PM';C=NS;
[R1]=fun(A,B,C,Aa,Bb);
R=max(R,R1);
A=NS';B=PB';C=NS;
[R1]=fun(A,B,C,Aa,Bb);
R=max(R,R1);
A=ZO';B=NB';C=PM;
[R1]=fun(A,B,C,Aa,Bb);
R=max(R,R1);
A=ZO';B=NM';C=PM;
[R1]=fun(A,B,C,Aa,Bb);
R=max(R,R1);
A=ZO';B=NS';C=PS;
[R1]=fun(A,B,C,Aa,Bb);
R=max(R,R1);
A=ZO';B=ZO';C=ZO;
[R1]=fun(A,B,C,Aa,Bb);
R=max(R,R1);
A=ZO';B=PS';C=NS;
[R1]=fun(A,B,C,Aa,Bb);
R=max(R,R1);
A=ZO';B=PM';C=NM;
[R1]=fun(A,B,C,Aa,Bb);
R=max(R,R1);
A=ZO';B=PB';C=NM;
[R1]=fun(A,B,C,Aa,Bb);
R=max(R,R1);
A=PS';B=NB';C=PS;
[R1]=fun(A,B,C,Aa,Bb);
R=max(R,R1);
A=PS';B=NM';C=PS;
[R1]=fun(A,B,C,Aa,Bb);
R=max(R,R1);
A=PS';B=NS';C=ZO;
[R1]=fun(A,B,C,Aa,Bb);
R=max(R,R1);
A=PS';B=ZO';C=NM;
[R1]=fun(A,B,C,Aa,Bb);
R=max(R,R1);
A=PS';B=PS';C=NM;
[R1]=fun(A,B,C,Aa,Bb);
R=max(R,R1);
A=PS';B=PM';C=NM;
[R1]=fun(A,B,C,Aa,Bb);
R=max(R,R1);
A=PS';B=PB';C=NM;
[R1]=fun(A,B,C,Aa,Bb);
R=max(R,R1);
A=PM';B=NB';C=ZO;
[R1]=fun(A,B,C,Aa,Bb);
R=max(R,R1);
A=PM';B=NM';C=ZO;
[R1]=fun(A,B,C,Aa,Bb);
R=max(R,R1);
A=PM';B=NS';C=NM;
[R1]=fun(A,B,C,Aa,Bb);
R=max(R,R1);
A=PM';B=ZO';C=NB;
[R1]=fun(A,B,C,Aa,Bb);
R=max(R,R1);
A=PM';B=PS';C=NB;
[R1]=fun(A,B,C,Aa,Bb);
R=max(R,R1);
A=PM';B=PM';C=NB;
[R1]=fun(A,B,C,Aa,Bb);
R=max(R,R1);
A=PM';B=PB';C=NB;
[R1]=fun(A,B,C,Aa,Bb);
R=max(R,R1);
A=PB';B=NB';C=ZO;
[R1]=fun(A,B,C,Aa,Bb);
R=max(R,R1);
A=PB';B=NM';C=ZO;
[R1]=fun(A,B,C,Aa,Bb);
R=max(R,R1);
A=PB';B=NS';C=NM;
[R1]=fun(A,B,C,Aa,Bb);
R=max(R,R1);
A=PB';B=ZO';C=NB;
[R1]=fun(A,B,C,Aa,Bb);
R=max(R,R1);
A=PB';B=PS';C=NB;
[R1]=fun(A,B,C,Aa,Bb);
R=max(R,R1);
A=PB';B=PM';C=NB;
[R1]=fun(A,B,C,Aa,Bb);
R=max(R,R1);
A=PB';B=PB';C=NB;
[R1]=fun(A,B,C,Aa,Bb);
R=max(R,R1);

Zz(i,j)=(R*Z')/sum(R)
end
end


function[R1]=fun(A,B,C,Aa,Bb)

RA = zeros(size(A,1),size(C,2));
for i=1 : size(A,1)
for j =1 : size(C,2)
RA(i,j) = max (min(A(i,:),C(:,j)'));
end
end

CA = zeros (size(Aa,1),size(RA,2));
for i =1:size(Aa,1)
for j=1 : size(RA,2)
CA= max (min(Aa(i,:),RA(:,j)'));
end
end

RB = zeros(size(B,1),size(C,2));
for i=1 : size(B,1)
for j =1 : size(C,2)
RB(i,j) = max (min(B(i,:),C(:,j)'));
end
end

CB = zeros (size(Bb,1),size(RB,2));
for i =1:size(Bb,1)
for j=1 : size(RB,2)
CB= max (min(Bb(i,:),RB(:,j)'));
end
end

R1 = min (CA,CB);


运行后提示
??? Subscripted assignment dimension mismatch.

Error in ==> mohu at 168
Zz(i,j)=(R*Z')/sum(R)

可能是R 的问题 R应该是一个1*13的矩阵 可是运行完好像输出的是一个数

很急 请高手帮忙
万分感谢 ！~~~~~