pline.asv

该代码为原创

function yy=pline(x,y)
n=length(x);
x=x(1:n);
y=y(1:n);
x=reshape(x,n,1);
y=reshape(y,n,1);
A=[x ones(n,1)];
b=y;
B=A'*A;
b=A'*b;
yy=B\b;
yy=yy';