那位朋友有用matlab编的偏最小二乘法的程序

function [a,b]=lsline(x,y)
%拟合曲线y=ax+b
xmean=mean(x);ymean=mean(y);
sumx2=(x-xmean)*(x-xmean)';
sumxy=(y-ymean)*(x-xmean)';
a=sumxy/sumx2;
b=ymean-a*xmean;
文件保存为lsline.m，即可供调用