chaos_pls3.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

x=Data(1:500, :);
y=Data(2:501,:);


chkData=Data(501:1000, :);ychkData=Data(502:1001,:);

% set parameters
global method;
method=1;

stepnum=10;


for step=5:stepnum
    
    k=1;pre_value=[];
    while k <=size(chkData,1)
        xt=chkData(k,:);
        for j=1:step
        yy=calculatey(x,y,xt); 
        pre_value(k+j-1)=yy(end);
        x=[x;xt];y=[y;yy];
        xt=[xt(2:end) yy(end)];
        end
        k=k+step;
    end;

ape=abs((ychkData(:,4)-pre_value(1:500)')./ychkData(:,4));
mape=mean(ape);
mse1=ychkData(:,4)'-pre_value(1:500);
rmse=sqrt(mse1*mse1'/500);
result2(:,step)=[step;mape;rmse];
step
% result2=[result2 rmse];
end