% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%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:1118, Data(t-117,:)=[x(t-18) x(t-12) x(t-6) x(t) ]; endglobal method;method=1;start=701;stepnum=20;yest=[];for step=1:stepnumx=Data(1:start-1, :);y=Data(step+1:start-1+step,:);chkData=Data(start:1000, :);ychkData=Data(start+step:end,:);xt=chkData(1,:);% set parameters        yy=calculatey(x,y,xt);         yest=[yest;yy(end)];end;ape=abs((Data(start+1:start+stepnum,4)-yest)./Data(start+1:start+stepnum,4));mape=mean(ape)