fuzzyc.m

matlab环境下开发的模糊控制原程序。

NL=[1,0.34,0,0,0,0,0,0,0];
NS=[0,0.5,1,0.5,0,0,0,0,0];
O=[0,0,0,0.24,1,0.24,0,0,0];
PS=[0,0,0,0,0,0.5,1,0.5,0];
PL=[0,0,0,0,0,0,0,0.34,1];
CONTAB=[2,2,2,1,-1;2,1,1,0,-1;2,1,0,-1,-2;2,0,-1,-1,-2;2,-1,-2,-2,-2];
SUBTAB=[NL;NS;O;PS;PL];
Uout=[];
C_Sub=[];
E_Sub=[];
for c=-4:4
	for e=-4:4
        U=[];
        UW=[];
		C_Sub=SUBTAB(:,c+5);
		E_Sub=SUBTAB(:,e+5);
		CONTR=C_Sub*E_Sub';
		for i=1:5
		    for j=1:5
			if CONTR(i,j)~=0
			    U=[U,CONTAB(i,j)];                              %查询控制规律表格
			    CONTR(i,j)=min(C_Sub(i),E_Sub(j));              %隶属度取小
			    UW=[UW,CONTR(i,j)];                             %记录控制权重
			end
		    end
		end
		lu=length(U);
		Udeltab=zeros(size(U));                                                 %标记重复的控制量
		Utem=[];
		UWtem=[];
        if lu>1
    	    for i=1:lu                                                      %合并相同的控制量权重取大    
        		for j=i+1:lu
               	    if Udeltab(j)~=1
    				    if U(i)==U(j)
				            UW(i)=max(UW(i),UW(j));
    				        Udeltab(j)=1;           
    				    end
        				Utem(i)=U(i);
				        UWtem(i)=UW(i);
    			    end
    			end
	        end
        else
            Utem=U;
            UWtem=UW;
        end
		Uend=[];
		UWend=[];
		for i=1:length(Utem)
		    Uend=[Uend;SUBTAB(Utem(i)+3,:)];
		    for j=1:9
			if Uend(i,j)~=0
			    Uend(i,j)=Uend(i,j)*UWtem(i);
			end
		    end
		end
		a=size(Uend);
		for j=1:9
		    for i=1:a(1)
			Uend(1,j)=max(Uend(1,j),Uend(i,j));
		    end
		end
		Uout(c+5,e+5)=0;
		sum=0;
		for i=1:9
		    Uout(c+5,e+5)=Uout(c+5,e+5)+Uend(1,i)*(i-5);
		    sum=sum+Uend(1,i);
		end
		Uout(c+5,e+5)=Uout(c+5,e+5)/sum;
        Uout(c+5,e+5)=round(Uout(c+5,e+5));
    end
end