文章来源:http://pocore.com/blog/article_529.html
原始数据 1 x y 最终拟合为y=a*x+1*b
1.000000 0.067732 3.176513
...
xArr: [[1.0, 0.067732], [1.0, 0.42781], ...]
yArr: [3.176513, 3.816464, 4.550095,...]
矩阵转换
xMat:n*2
[[1. 0.067732]
[1. 0.42781 ]
....]
yMat:n*1
[[3.176513]
[3.816464]
...]
平方误差
普通最小二乘法讲解OLS
链接普通最小二乘法
系数求法公式:
矩阵的逆可能并不存在,因此必须要在代码中对此作出判断
xMat.T是xMat的转置 2*n
矩阵行列式 xTx=xMat.T*xMat 2*n * n*2 =2*2
xTx:
[[200. 97.6638 ]