mg_chaos_pls4.m

一个简单的pls程序吧

% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%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) ]; 
end

global method;
method=1;
start=701;
stepnum=20;
yest=[];
for step=1:stepnum
x=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)