PLS的matlab实现主程序

%PLS的matlab实现代码
x=x';%转置，波长一列转为一行
 X=X';%样本数作为行 X=XX';
plot(x,XX); %
%预处理
[Xsm11]=smooth(Xs,11);%对Xs取11点平滑，只能是奇数可以修改，看哪个最好 Xs=Xtr; Xs=Xte;   平滑变成一列 
             % [Xsm11]= smooth(X,11); 对未排序 的X进行平滑 Xs=Xte(1,:); Xte=Xte';
Xsm11r=reshape(Xsm11,[16 12287]);

xsnv=SNV(Xtr);%数据预处理改到对应路径，G:\Matlab 学习\数据处理全部相关内容\预处理，  xsnv=SNV(Xsm11r); %44个样全用来建模
xtnv=SNV(Xte);%   标准正态变量变换 standardization of Matlab process variables 
% xtnv=SNV(Xsm11r);%对平滑后的一条光谱进行SNV处理，将平滑后的全作为未知样
[xmsc,me]=MSC(Xtr,1,16);%307和Xtr变量对应，307个波长变量，多元散射校正预处理（MSC）Multiplicative scatter correction ，
%以消除光谱平移带来的误差
[xtmsc,me]=MSC(Xte,1,16);

d1Xtr=deriv(Xtr,1);%一阶导数，训练集
d1Xte=deriv(Xte,1);%预测集
%  d1Xte=deriv(Xsm11r,1); %对平滑后的一条光谱进行一阶导,%format rat 使输出为分数
d2Xtr=deriv(Xtr,2);%二阶导数 33x4096(样品数x波长点数）下同
d2Xte=deriv(Xte,2);
%%%中心化，还要对Y也进行 中心化%%%%%
[xcen,me,ctest]=center(Xtr,1,Xtr);%中心化预处理
[xtcen,me,ctest]=center(Xte,1,Xte);%至此处
%对Y也进行中心化
[Ytrcen,me,ctest]=center(Ytr,1,Ytr);%中心化预处理
[Ytecen,me,ctest]=center(Yte,1,Yte);

%定量分析