Numpy与线性方程组2(Python与高等代数02)

本文探讨了非齐次线性方程组AX=b的解法,包括行满秩、列满秩和秩亏损情况下的处理。对于行满秩或列满秩矩阵,可以通过右逆或左逆求得特解;而对于秩亏损矩阵,需要用到伪逆。Python的numpy库中的pinv()函数可用于计算逆或伪逆。特解的计算公式为右逆或左逆与b的点乘。下期将介绍如何进行自动行简化程序及具体代码实现。
摘要由CSDN通过智能技术生成

4.有无限多解的非齐次线性方程组AX=b

一个有无限多解的非齐次线性方程组的解有通解和特解组成

通解只需要解AX=0即可

在这一段中我们要找出求特解的方法.

1)对于行满秩的长方阵

我们需要得到它的右逆\left ( A^{-1} \right )_{R}^{}

这里\left ( A^{-1} \right )_{R}=A_{}^{T}\left ( AA^{T} \right )^{-1}

当然A\left ( A^{-1} \right )_{R}=E

2)对于列满秩的长方阵

我们需要得到它的左逆,方法与上面求右逆的方法类似

3)对于秩亏损的长方阵

这时候我们需要得到他的伪逆,这里将在下期作解释

4)Python 里numpy功能十分强大,用pinv()方法可自动求出上述三个逆阵

但是,首先我们得对目标做行简化,具体代码将在下期给出

行满秩的时候,特解就是他的右逆点乘以b

列满秩的时候,特解就是他的左逆点乘以b

代码如下

np.linalg.pinv(A).dot(b)

5.自动行简化程序

见下期

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

River Chandler

谢谢,我会更努力学习工作的!!

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值