plsfactor.m

MATLAB科学计算与工程分析源代码源程序3

function [omega,t,pp,XXX,r,YYY]=plsfactor(X0,Y0)

XX=X0'*Y0*Y0'*X0;
[V,D]=eig(XX);
Lamda=max(D);
[MAXLamda,I]=max(Lamda);
omega=V(:,I);          %zui da te zheng  zhi  dui  ying  de  te  zheng  xiang  liang  

%%%di yi  ge  zhu  yuan 
t=X0*omega;
pp=X0'*t/(t'*t);
XXX=X0-t*pp';
r=Y0'*t/(t'*t);
YYY=Y0-t*r';