plsc.rar_pls

clear all data=xlsread('test'); e=ones(length(data),1); P=data(:,1:3); rank(P); T=data(:,4:end); datamean=mean(P); delta=std(P); e=ones(length(P),1); X2=(P-e*datamean)./(e*delta);%数据标准化 X0=X2; tdatamean=mean(T); tdelta=std(T); e=ones(length(T),1); Y2=(T-e*tdatamean)./(e*tdelta); Y0=Y2; m=rank(X0);% for n=1:m M=X0'*Y0*Y0'*X0; [v,d]=eig(M); [m1,i]=max(diag(d)); w(:,n)=v(:,i); t(:,n)=X0*v(:,i); a(:,n)=X0'*t(:,n)/(t(:,n)'*t(:,n)); X0=X0-t(:,n)*a(:,n)'; ZX(:,:,n)=X0; b(:,n)=Y0'*t(:,n)/(t(:,n)'*t(:,n)); Y0=Y0-t(:,n)*b(:,n)'; ZY(:,:,n)=Y0; end %J(:,:,1)=w(:,1)*b(:,1)' E=eye(size(w(:,1)*a(:,1)')); %J(:,:,2)=(eye(size(w(:,1)*a(:,1)'))-w(:,1)*a(:,1)')*w(:,2)*b(:,2)' %J(:,:,1)+J(:,:,2) t1=w(:,1); t2=E-t1*a(:,1)'; TT(:,:,1)=[1 2 3; 2 3 4; 3 4 5;]; TT(:,:,2)=t2; for j=3:m TT(:,:,j)=TT(:,:,j-1)-TT(:,:,j-1)*w(:,j-1)*a(:,j-1)'; %递推 end Tk(:,:,1)=t1*b(:,1)';; for j=2:m Tk(:,:,j)=TT(:,:,j)*w(:,j)*b(:,j)';%计算tk与X2的关系 end %t3=t2-t2*w(:,2)*a(:,2)'; %t4=t3-t3*w(:,3)*a(:,3)'; R=0; for j=1:m R=R+(Tk(:,:,j)); W(:,:,j)=R; %m个成分的系数 end %(eye(size(w(:,1)*a(:,1)'))-w(:,1)*a(:,1)')*w(:,2)*b(:,2)' %J(:,:,2) %(eye(3)-w(:,1)*a(:,1)')*w(:,2)*b(:,2)'+w(:,1)*b(:,1)'; %还原成原始变量 for j=1:3 for i=1:3 s(:,i,j)=((W(:,i,j)'./delta)*tdelta(i))'; %*tdelta %o(i,:,j)=-s(:,i,j)'; p(j,i)=sum(-s(:,i,j)'.*datamean)+tdatamean(i); %+tdatamean end end %验证效果 for j=1:3 result(:,:,j)=(s(:,:,j)'*P')'+ones(20,1)*p(j,:); end %(s(:,:,2)'*P')'+ones(20,1)*p(2,:); %(s(:,:,3)'*P')'+ones(20,1)*p(3,:); for j=1:3 error(:,:,j)=abs(result(:,:,j)-T); end sum(error)