zycci.m

fuzzy control rule design

function Y=zycci(n,m)
%%%%E=-6and Ed=-6
%%%%%%%%%%%%%%%%%%%%语言变量X的隶属度%%%%%%%%%%%%%%%%%%%%%
A=[1.0 0.8 0.7 0.4 0.1 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0;%NB
    0.2 0.7 1.0 0.7 0.3 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0;%NM
    0.0 0.1 0.3 0.7 1.0 0.7 0.2 0.0 0.0 0.0 0.0 0.0 0.0;%NS
    0.0 0.0 0.0 0.0 0.1 0.6 1.0 0.0 0.0 0.0 0.0 0.0 0.0;%NZ
    0.0 0.0 0.0 0.0 0.0 0.0 1.0 0.6 0.1 0.0 0.0 0.0 0.0;%PZ
    0.0 0.0 0.0 0.0 0.0 0.0 0.2 0.7 1.0 0.3 0.1 0.0 0.0;%PS
    0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.2 0.7 1.0 0.7 0.3;%PM
    0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.1 0.4 0.7 0.8 1.0];%PB
%%%%%%%%%%%%%%%%%误差，误差的变化及U的隶属函数%%%%%%%%%%%%
B=[1.0 0.7 0.3 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0;%NB
    0.3 0.7 1.0 0.7 0.3 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0;%NM
    0.0 0.0 0.3 0.7 1.0 0.7 0.3 0.0 0.0 0.0 0.0 0.0 0.0;%NS
    0.0 0.0 0.0 0.0 0.3 0.7 1.0 0.7 0.3 0.0 0.0 0.0 0.0;%Z
    0.0 0.0 0.0 0.0 0.0 0.0 0.3 0.7 1.0 0.7 0.3 0.0 0.0;%PS
    0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.3 0.7 1.0 0.7 0.3;%PM
    0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.3 0.7 1.0];%PB
C=[1.0 0.7 0.3 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0;%NB
    0.3 0.7 1.0 0.7 0.3 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0;%NM
    0.0 0.0 0.3 0.7 1.0 0.7 0.3 0.0 0.0 0.0 0.0 0.0 0.0;%NS
    0.0 0.0 0.0 0.0 0.3 0.7 1.0 0.7 0.3 0.0 0.0 0.0 0.0;%Z
    0.0 0.0 0.0 0.0 0.0 0.0 0.3 0.7 1.0 0.7 0.3 0.0 0.0;%PS
    0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.3 0.7 1.0 0.7 0.3;%PM
    0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.3 0.7 1.0];%PB
%模糊控制规则表的计算
for a=1:1:8;
    for b=1:1:7;
        cc(a,b)=fsbiao2(a,b);%函数fsbiao2.m的调用
    end
end
MO=cc;
L1=size(MO,1);%8行数的求取
L2=size(MO,2);%7列数的求取
CC=[0 0 0 0 0 0 0 0 0 0 0 0 0];%累加初始化
for i=1:1:L1;
    for j=1:1:L2;
        %计算每一条规则的值C(k),C(k) U C(k+1).
        k=MO(i,j);%K的返回值表示Cxx
        c=C(k,:);
        a=A(i,:);
        b=B(j,:);       
        %取Ap和Bp，计算C和C1，并求C U C1的值
        Ap=[0 0 0 0 0 0 0 0 0 0 0 0 0];
        Bp=[0 0 0 0 0 0 0 0 0 0 0 0 0];
        Ap(n)=A(i,n);
        Bp(m)=B(j,m);%关于对Ap和Bp矩阵中元素的处理
        %调用合成以及取小的函数
        d=hec(a,c);
        d1=hec1(b,c);
        c1=qumin(Ap,d);%取小函数的调用
        c2=qumin(Bp,d1);%取小函数的调用
        R1=min(c1,c2);
        CC=max(R1,CC);
    end
end
CC;
% d=hec(a,b);
% d1=hec1(b,c);
% Ap=[1 0 0 0 0 0 0 0 0 0 0 0 0];
% Bp=[1 0 0 0 0 0 0 0 0 0 0 0 0];
% c1=qumin(Ap,d);%取小函数的调用
% c2=qumin(Bp,d1);%取小函数的调用
% R=max(c1,c2)
for i=1:1:length(CC);
    X(i)=CC(i)*(i-7);
end
Z=X;
Y=sum(Z)/sum(CC);