% 預測序列 % 輸入模糊集的中心 InCen  = linspace(gMin,gMax,InNum); OutCen = linspace(gMin,gMax,OutNum); MyOut  = []; for i=52:101-n     MFVal = [];     for j =1 : InNum         MFVal = [MFVal;evalmf(g(i),InMFP(j,:),'trimf'),evalmf(g(i+1),InMFP(j,:),'trimf'),...            evalmf(g(i+2),InMFP(j,:),'trimf'),evalmf(g(i+3),InMFP(j,:),'trimf')];     end     [MaxMF Index] = max(MFVal,[],1);     for j =1 :size(Rule,1)         if Rule(j,1:n)==Index             OIndex = Rule(j,n+1);             break;           else             OIndex =0;         end     end     if OIndex>0         MyOut = [MyOut,OutCen(OIndex)*prod(MaxMF)/];      else         MyOut = [MyOut 0];     endendfigure(2);plot(x(56:101),g(56:101),'r',x(56:101),MyOut(1:46),'b');legend('原來序列','擬和序列');