首先看两个个结论:
结论一:方程组Ax=b的最小二乘解的通式为x=Gb+(I-GA)y, 其中G\in A\{1, 3\}, y是\mathbb C^n中的任意向量....结论二:只有A是满秩时, 矛盾方程组Ax=b 的最小二乘解才是唯一的, 且为x_0=(A^HA)^{-1}A^Hb. 否则, 便有无穷多个最小二乘解....下面看一个实例:
求矛盾方程组
\begin{cases}x_1+2x_2=1, \\2x_1+x_2=0, \\x_1+x_2=0\end{cases}的最小二乘解。...numpy矩阵运算
import numpy as np
A = np.mat([[1, 2], [2, 1], [1, 1]])
A13 = (A.H * A).I * A.H
print(A13)
利用最小二乘法做线性拟合...\\kx_n+b=y_n\end{cases}
这里的k和b为变量,使用上述公式求解出k和b的值,则可以得到变量的最小二乘线性拟合方程。