tsfuzzyp28.m

Fuzzy Control Systems Design and Analysis [TANAKA& HUA O. WANG] CH1

clear;
clc;



a=1;


x(1:2,1)=[0.9 -0.7]';

for t=1:1:30


A1=[1 -0.5 
    1   0];

A2=[-1 -0.5 
     1  0];


% MBF
%  M1                           M2       
%---------- \            / ----------                  
%             \        /    
%               \    /        
%                 \       
%                /  \ 
%------------------------------------
%           -a           a


%M1
if x(2,t)<=-a; M1=1;end
if x(2,t) > -a;M1=(a-x(2,t))/(a-(-a));end
if x(2,t) > 1 ; M1=0;end

%M2
if x(2,t)<=-a;M2=0;end
if x(2,t) > -a;M2=(x(2,t)+a)/(a-(-a));end
if x(2,t) >= 1 ;M2=1;end



%Dfuzzy


x(1:2,t+1)=(M1*A1*x(1:2,t)+M2*A2*x(1:2,t))/(M1+M2);



end

c=1:1:31;



fi1=figure;set(fi1,'color','white');set(fi1,'name',datestr(now,30));
plot(c,x(1,:),'r','LineWidth',2.5)
axis([1 20 -6 6])
xlabel('t'); ylabel('x1(t)');