% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%function load_main(iteration)%close all;% clear;% clear;% load mgdata.dat% t = mgdata(:, 1); x = mgdata(:, 2); % plot(t, x);clear;load('G:\200510\Special\data\xdata.mat');% result2=[];% noiselev=[0,0.001,0.01,0.1];% for ii=1:4 % % x(1:617)=x(1:617)+sqrt(noiselev(ii))*randn(1,length(x(1:617)));for t=118:1117, Data(t-117,:)=[x(t-18) x(t-12) x(t-6) x(t) x(t+6)]; end% yData=[];% for ii=1:1000%     if (ii+6)>1000%         break;%     end;%     yData=[yData;Data(ii+6,:)];%   % endx=Data(1:500, 1:4);% y=yData(1:500,:);y=Data(1:500, 5);x=x(300:500,:);y=y(300:500,:);chkData=Data(501:end, 1:4);ychkData=Data(501:end,5);% set parametersglobal method;method=1;% step=20;result2=zeros(3,20);for step=1:5k=0;pre_value=[];while k < size(chkData,1)    for j=1:6        if (k+j>500)             break;         end        yy=calculatekery(x,y,chkData(k+j,:));                pre_value(k+j)=yy(end);           input=chkData(k+j,:);           result=pre_value(k+j);%            x=[x;chkData(k+j,:)];y=[y;input(2:end) result];        for i=1:step-1           input=[input(2:end) result];           yy=calculatekery(x,y,input);                       pre_value(k+j+i*6)=yy(end);           result=pre_value(k+j+i*6);%            x=[x;input];y=[y;input(2:end) result];       end;    end    k=k+6*step;endape=abs((ychkData(1:500)-pre_value(1:500)')./ychkData(1:500));mape=mean(ape);mse1=ychkData(1:500)'-pre_value(1:500);rmse=sqrt(mse1*mse1'/500);result2(:,step)=[step;mape;rmse];step% result2=[result2 rmse];end