4.有无限多解的非齐次线性方程组AX=b
一个有无限多解的非齐次线性方程组的解有通解和特解组成
通解只需要解AX=0即可
在这一段中我们要找出求特解的方法.
1)对于行满秩的长方阵
我们需要得到它的右逆
这里
当然
2)对于列满秩的长方阵
我们需要得到它的左逆,方法与上面求右逆的方法类似
3)对于秩亏损的长方阵
这时候我们需要得到他的伪逆,这里将在下期作解释
4)Python 里numpy功能十分强大,用pinv()方法可自动求出上述三个逆阵
但是,首先我们得对目标做行简化,具体代码将在下期给出
行满秩的时候,特解就是他的右逆点乘以b
列满秩的时候,特解就是他的左逆点乘以b
代码如下
np.linalg.pinv(A).dot(b)
5.自动行简化程序
见下期