您可以使用^{}。在

例如,下面是一个numpy整数数组:In [148]: m = np.array([[10, 11, 0, 0], [9, 10, 0, 0], [100, 0, 20, 3], [-10000, 200, 133, 20]])

In [149]: m

Out[149]:

array([[ 10, 11, 0, 0],

[ 9, 10, 0, 0],

[ 100, 0, 20, 3],

[-10000, 200, 133, 20]])

精确的行列式是1,但是np.linalg.det()和np.linalg.inv()都会引入浮点错误:

^{pr2}$

从numpy数组创建一个sypy ^{}对象:In [156]: import sympy

In [157]: M = sympy.Matrix(m)

In [158]: M

Out[158]:

Matrix([

[ 10, 11, 0, 0],

[ 9, 10, 0, 0],

[ 100, 0, 20, 3],

[-10000, 200, 133, 20]])

sympy计算是精确的:In [159]: M.det()

Out[159]: 1

In [160]: M.inv()

Out[160]:

Matrix([

[ 10, -11, 0, 0],

[ -9, 10, 0, 0],

[-325400, 358000, 20, -3],

[2169000, -2386300, -133, 20]])

要将逆矩阵转换回numpy数组,可以执行以下操作:In [185]: Minv = M.inv()

In [186]: minv = np.asarray(Minv).astype(int)

In [187]: minv

Out[187]:

array([[ 10, -11, 0, 0],

[ -9, 10, 0, 0],

[ -325400, 358000, 20, -3],

[ 2169000, -2386300, -133, 20]])

必须使用astype(int)方法,因为numpy数组中没有object的数据类型:In [188]: np.asarray(Minv)

Out[188]:

array([[10, -11, 0, 0],

[-9, 10, 0, 0],

[-325400, 358000, 20, -3],

[2169000, -2386300, -133, 20]], dtype=object)

这个结果是一个整数的整数数组。在

在转换回numpy时要小心:逆整数可能比用64位整数表示的大。在

在matlab里如何构建幺模矩阵,幺模矩阵精确逆Python相关推荐

  1. 用matlab自己搭建bp神经网络,怎样在matlab里建立一个BP神经网络模型?

    该楼层疑似违规已被系统折叠 隐藏此楼查看此楼 用以下的数据怎样在matlab里建立一个BP神经网络模型?求高手帮忙!!最好是有详细步骤以及代码 年份 WTI(美元/桶) 2007-1 54.26 20 ...

  2. semilogx 多条曲线_怎么让两个指数在一个坐标,matlab里怎样一个坐标上显示多个曲线,而且横轴要用指数形式的?谢谢...

    Q1:matlab里怎样一个坐标上显示多个曲线,而且横轴要用指数形式的?谢谢 多个纵轴数组分别是y1,y2,y3,横轴数组为x 命令为: semilogx(x,y1,x,y2,x,y3) 完了 Q2: ...

  3. 用matlab跑神经网络模型,怎样在matlab里建立一个BP神经网络模型?

    该楼层疑似违规已被系统折叠 隐藏此楼查看此楼 用以下的数据怎样在matlab里建立一个BP神经网络模型?求高手帮忙!!最好是有详细步骤以及代码 年份 WTI(美元/桶) 2007-1 54.26 20 ...

  4. 怎样用matlab模型建立,怎样在matlab里建立一个BP神经网络模型?

    该楼层疑似违规已被系统折叠 隐藏此楼查看此楼 用以下的数据怎样在matlab里建立一个BP神经网络模型?求高手帮忙!!最好是有详细步骤以及代码 年份 WTI(美元/桶) 2007-1 54.26 20 ...

  5. Matlab里evalin和assignin的用法

    原文地址:Matlab里evalin和assignin的用法作者:了凡春秋 assignin与evalin     MATLAB通常的基本工作空间是'base'空间.MATLAB在程序运行过程中,将为 ...

  6. matlab里插入行和列

    matlab里插入行和列 >> A=magic(4) A =     16     2     3    13      5    11    10     8      9     7 ...

  7. MATLAB里面的filter和filtfilt的C语言源代码

    MATLAB里面的filter和filtfilt的C语言源代码 嗯,算法非常简单,就是网上搜不到C代码实现.filter是个很万能的数字滤波器函数,只要有滤波器的差分方程系数,IIR呀FIR呀都能通过 ...

  8. matlab中多元线性回归regress函数精确剖析(附实例代码)

    matlab中多元线性回归regress函数精确剖析(附实例代码) 目录 前言 一.何为regress? 二.regress函数中的参数 三.实例分析 总结 前言 regress函数功能十分强大,它可 ...

  9. Matlab画圆方向图,如何在matlab里画天线的立体方向图?

    如何在matlab里画天线的立体方向图?有没有具体的函数? 有的,你看下一般的matlab教程里都有着方面的介绍,查一下就可以了 有数据,Matlab可以画 归根到底,画极坐标方向图是坐标变换的问题. ...

最新文章

  1. ATS 5.3.0命令行工具traffic_via
  2. java反射获得注释_Java反射、注解
  3. Java:多线程之线程池
  4. 4x4矩阵键盘工作原理及扫描程序_单片机人机交互矩阵按键
  5. 图解 SpringMVC 五大组件
  6. 信息系统项目管理01——信息化和信息系统
  7. thrift实战教程
  8. 插桩 java_Javassist进行方法插桩
  9. matlab psf点扩散函数,什么是点扩散函数PSF?-ZEMAX,OpticStudio-光行天下-国内光电行业门户-光学,光电,光机技术及其软件运用专业网站!...
  10. tomcat 官网下载
  11. linux samba 漏洞 exp,smaba漏洞总结
  12. 腹有诗书气自华,新华三交换机大有不同
  13. npm安装出现npm err FetchError Invalid response body while trying to fetch httpsregistry.npmjs.orgvu
  14. qq邮箱里面html的图片无法显示,QQ邮箱如何不显示图片? QQ邮箱如何不显示图片?...
  15. Linux运维部落 DNS服务,Linux运维部落
  16. C语言实现猜数字小游戏(详细讲解)
  17. this.className的使用
  18. 使用 ALT 键输入特殊符号
  19. 华为服务器bmc snmp协议,华为S5700开启SNMP协议
  20. Android手势密码view笔记(一)

热门文章

  1. java 日期calendar_java成神之——date和calendar日期的用法
  2. mysql嵌套查询效率低,连接查询代替嵌套查询提高select效率
  3. 如何做好技术布道——用影响影响影响
  4. 【具体实现过程】百度地图POI区域的边界经纬度坐标数据采集
  5. python自动投递_利用python如何在前程无忧高效投递简历
  6. JS 获取星期几的四种写法(转)
  7. Java开发面试题目,java培训老师第一次讲课
  8. 【2017 BSUIR Semifinal D】Friends rescue 题解
  9. GitChat·管理 | 一篇文章读懂项目管理中的精髓
  10. android重力传感器横竖反,Android重力感应实现方式