Python 如何求矩阵的逆
更新时间：2021年03月11日 10:21:28 作者：Kenn7
这篇文章主要介绍了Python 如何求矩阵的逆案例，具有很好的参考价值，希望对大家有所帮助。一起跟随小编过来看看吧
我就废话不多说了，大家还是直接看代码吧~
import numpy as np kernel = np.array([1, 1, 1, 2]).reshape((2, 2)) print(kernel) print(np.linalg.inv(kernel))
注意，Singular matrix奇异矩阵不可求逆
补充：python+numpy中矩阵的逆和伪逆的区别
定义：
对于矩阵A，如果存在一个矩阵B，使得AB=BA=E，其中E为与A,B同维数的单位阵，就称A为可逆矩阵（或者称A可逆），并称B是A的逆矩阵，简称逆阵。（此时的逆称为凯利逆）
矩阵A可逆的充分必要条件是|A|≠0。
伪逆矩阵是逆矩阵的广义形式。由于奇异矩阵或非方阵的矩阵不存在逆矩阵，但可以用函数pinv(A)求其伪逆矩阵。
基本语法为X=pinv(A),X=pinv(A,tol),其中tol为误差,pinv为pseudo-inverse的缩写：max(size(A))*norm(A)*eps。
函数返回一个与A的转置矩阵A' 同型的矩阵X，并且满足：AXA=A,XAX=X.此时，称矩阵X为矩阵A的伪逆，也称为广义逆矩阵。
pinv(A)具有inv(A)的部分特性，但不与inv(A)完全等同。
如果A为非奇异方阵，pinv(A)=inv(A)，但却会耗费大量的计算时间，相比较而言，inv(A)花费更少的时间。
代码如下：
1.矩阵求逆
import numpy as np a = np.array([[1, 2], [3, 4]]) # 初始化一个非奇异矩阵(数组) print(np.linalg.inv(a)) # 对应于MATLAB中 inv() 函数 # 矩阵对象可以通过 .I 求逆,但必须先使用matirx转化 A = np.matrix(a) print(A.I)
2.矩阵求伪逆
import numpy as np # 定义一个奇异阵 A A = np.zeros((4, 4)) A[0, -1] = 1 A[-1, 0] = -1 A = np.matrix(A) print(A) # print(A.I) 将报错，矩阵 A 为奇异矩阵，不可逆 print(np.linalg.pinv(A)) # 求矩阵 A 的伪逆（广义逆矩阵），对应于MATLAB中 pinv() 函数
这就是矩阵的逆和伪逆的区别
截至2020/10/4，matrix函数还可以使用，但已经过时，应该是mat函数这种。
以上为个人经验，希望能给大家一个参考，也希望大家多多支持脚本之家。如有错误或未考虑完全的地方，望不吝赐教。
相关文章
这篇文章主要介绍了Python 爬虫多线程详解及实例代码的相关资料,需要的朋友可以参考下2016-10-10
这篇文章主要介绍了一文带你了解Python 四种常见基础爬虫方法介绍，文中通过示例代码介绍的非常详细，对大家的学习或者工作具有一定的参考学习价值，需要的朋友们下面随着小编来一起学习学习吧2020-12-12
Python的强大之处在于他有非常丰富和强大的标准库和第三方库，几乎你想实现的任何功能都有相应的Python库支持。下面通过本文给大家介绍python入门基础之用户输入与模块初认识，一起看看吧2016-11-11
这篇文章主要介绍了TensorFlow2基本操作之合并分割与统计，本文给大家介绍的非常详细，对大家的学习或工作具有一定的参考借鉴价值，需要的朋友可以参考下2021-09-09
这篇文章主要介绍了OpenCV基本操作之图像轮廓，本文给大家介绍的非常详细，对大家的学习或工作具有一定的参考借鉴价值，需要的朋友可以参考下2021-09-09
这篇文章主要介绍了Python实现简易Web爬虫详解，具有一定借鉴价值,需要的朋友可以参考下2018-01-01
这篇文章主要为大家详细介绍了python实现简单聊天程序的相关代码，包括客户端和服务端，具有一定的参考价值，感兴趣的小伙伴们可以参考一下2018-04-04
这篇文章主要介绍了Python标准库与第三方库,需要的朋友可以参考下2014-07-07
这篇文章主要介绍了python调用Delphi写的Dll代码示例，具有一定参考价值，需要的朋友可以了解下。2017-12-12
这篇文章主要介绍了Python关于抽奖系统的思考与设计思路,本文通过一些简单的例子来说一说抽奖系统背后的逻辑,看看究竟是你运气不好还是系统逻辑在作怪,需要的朋友可以参考下2023-03-03
最新评论