</pre><pre code_snippet_id="370100" snippet_file_name="blog_20140530_2_4718563" name="code" class="plain">
clc;
clear;
load('sample_alldata.mat');
sample=10000*sample_alldata(:,2:end-1);
wave_s=34;
wave_e=464;
data_train=sample(1:65,wave_s:wave_e);
data_predict=sample(66:end,wave_s:wave_e);
target_out=sample_alldata(1:65,end);%训练集
data_predict=data_predict/max(data_predict);%测试
predict_out=sample_alldata(66:end,end);
num=size(data_train,1);%训练集样本数
for i=1:num
%数据标准化
data_train(i,:)=zscore(data_train(i,:));
end
n=size(data_train,2);%自变量维度
m=size(target_out,2);%因变量维度
%for i=1:num-1
%forj=1:num
% [XL,YL,XS,YS,beta,PCTVAR,MSE]=plsregress(data_train,target_out,i);
%end
%end
handleofwaitbar=waitbar(0);%设置进度条便于观察
%交叉验证
SS=[];
PRESS=[];