偏最小二乘法Matlab源代码

function [k b y T W]=NEWPLS（XYm）
%PLS回归主程序
%输入自变量表X，因变量Y，主成分值m；
%输出拟合参数：系数向量k，截距b，拟合值y，并绘出拟合效果图

E0=zscore（X）;
F0=zscore（Y）;

E=E0;
F=F0;
W=[];
T=[];
P=[];
R=[];
WX=[];
Qh2=0;
R2=0;
for i=1:m
    w=E‘*F/norm（E‘*F）;
    t=E*w;
    p=E‘*t/（t‘*t）;
    r=F‘*t/（t‘*t）;
    E=E-t*p‘;
    F=F-t*r‘;
%    [w t p r E F]=NPLS_1（EF）;
    W=[Ww];
    T=[Tt];
    P=[Pp];
    R=[Rr];
%    Qh2=Qh2+（norm（t））.^2*（norm（p））.^2/（norm（E0‘fro‘））.^2;
%    R2=R2+（norm（t））.^2*（norm（r））.^2/（norm（F0））.^2;
    
%[rW cW]=size（W）;
%wx=eye（rW）;
%for i=1:cW-1
%    wx=wx*（eye（rW）-W（:i）*P（:i）‘）;
%end
%    wx=wx*W（: