一元线性回归函数程序.m

来自「对“曲线的校正”是很多行业和领域都要面临的问题。本人编写并总结了变量校正的源代码」· M 代码 · 共 20 行

function [a,b,r]=linear(x,y,n)
x1=0;
x2=0;
y1=0;
y2=0;
xy=0;
for i=1:n						
    x1=x1+x(i);						%计算Σxi
    x2=x2+x(i)*x(i);						%！计算Σxi2
    y1=y1+y(i);						%！计算Σyi
    y2=y2+y(i)*y(i);						%！计算Σyi2
    xy=xy+x(i)*y(i);						%！计算Σxiyi
end
lxy=xy-x1*y1/n;							%！计算Lxy
lxx=x2-x1*x1/n;							%！计算Lxx
lyy=y2-y1*y1/n;							%！计算Lyy
b=lxy/lxx;							%！计算b
a=(y1-b*x1)/n;							%！计算a
r=lxy/sqrt(lxx*lyy);						%！计算R