linefit.m

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

function [k,b]=linefit(x,y)
n=length(x);
x=reshape(x,n,1);
y=reshape(y,n,1);
A=[x,ones(n,1)];
bb=y;
B=A'*A;
bb=A'*bb;
yy=B\bb;
size(B)
size(bb)
k=yy(1);
b=yy(2);