predict.asv

建立Fuzzy模型－序列预测matlab代码

 % 预测序列
 % 输入模糊集的中心
 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];
     end
end
figure(2);
plot(x(56:101),g(56:101),'r',x(56:101),MyOut(1:46),'b');
legend('原来序列','拟和序列');