predict.m

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

 % 预测序列
 % 输入模糊集的中心
 InCen  = linspace(gMin,gMax,InNum);
 OutCen = linspace(gMin,gMax,OutNum);
 MyOut  = [];
 for i=SLen(1)+1:SLen(2)-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)]; 
     else
         MyOut = [MyOut 0];
     end
end
figure(2);
plot(k(SLen(1)+n+1:SLen(2)),g(SLen(1)+n+1:SLen(2)),'r',k(SLen(1)+n+1:SLen(2)),MyOut,'b--');
legend('原来序列','拟和序列');