在matlab里如何构建幺模矩阵,幺模矩阵精确逆Python
您可以使用^{}。在
例如,下面是一个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相关推荐
- 用matlab自己搭建bp神经网络,怎样在matlab里建立一个BP神经网络模型?
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼 用以下的数据怎样在matlab里建立一个BP神经网络模型?求高手帮忙!!最好是有详细步骤以及代码 年份 WTI(美元/桶) 2007-1 54.26 20 ...
- semilogx 多条曲线_怎么让两个指数在一个坐标,matlab里怎样一个坐标上显示多个曲线,而且横轴要用指数形式的?谢谢...
Q1:matlab里怎样一个坐标上显示多个曲线,而且横轴要用指数形式的?谢谢 多个纵轴数组分别是y1,y2,y3,横轴数组为x 命令为: semilogx(x,y1,x,y2,x,y3) 完了 Q2: ...
- 用matlab跑神经网络模型,怎样在matlab里建立一个BP神经网络模型?
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼 用以下的数据怎样在matlab里建立一个BP神经网络模型?求高手帮忙!!最好是有详细步骤以及代码 年份 WTI(美元/桶) 2007-1 54.26 20 ...
- 怎样用matlab模型建立,怎样在matlab里建立一个BP神经网络模型?
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼 用以下的数据怎样在matlab里建立一个BP神经网络模型?求高手帮忙!!最好是有详细步骤以及代码 年份 WTI(美元/桶) 2007-1 54.26 20 ...
- Matlab里evalin和assignin的用法
原文地址:Matlab里evalin和assignin的用法作者:了凡春秋 assignin与evalin MATLAB通常的基本工作空间是'base'空间.MATLAB在程序运行过程中,将为 ...
- matlab里插入行和列
matlab里插入行和列 >> A=magic(4) A = 16 2 3 13 5 11 10 8 9 7 ...
- MATLAB里面的filter和filtfilt的C语言源代码
MATLAB里面的filter和filtfilt的C语言源代码 嗯,算法非常简单,就是网上搜不到C代码实现.filter是个很万能的数字滤波器函数,只要有滤波器的差分方程系数,IIR呀FIR呀都能通过 ...
- matlab中多元线性回归regress函数精确剖析(附实例代码)
matlab中多元线性回归regress函数精确剖析(附实例代码) 目录 前言 一.何为regress? 二.regress函数中的参数 三.实例分析 总结 前言 regress函数功能十分强大,它可 ...
- Matlab画圆方向图,如何在matlab里画天线的立体方向图?
如何在matlab里画天线的立体方向图?有没有具体的函数? 有的,你看下一般的matlab教程里都有着方面的介绍,查一下就可以了 有数据,Matlab可以画 归根到底,画极坐标方向图是坐标变换的问题. ...
最新文章
- ATS 5.3.0命令行工具traffic_via
- java反射获得注释_Java反射、注解
- Java:多线程之线程池
- 4x4矩阵键盘工作原理及扫描程序_单片机人机交互矩阵按键
- 图解 SpringMVC 五大组件
- 信息系统项目管理01——信息化和信息系统
- thrift实战教程
- 插桩 java_Javassist进行方法插桩
- matlab psf点扩散函数,什么是点扩散函数PSF?-ZEMAX,OpticStudio-光行天下-国内光电行业门户-光学,光电,光机技术及其软件运用专业网站!...
- tomcat 官网下载
- linux samba 漏洞 exp,smaba漏洞总结
- 腹有诗书气自华,新华三交换机大有不同
- npm安装出现npm err FetchError Invalid response body while trying to fetch httpsregistry.npmjs.orgvu
- qq邮箱里面html的图片无法显示,QQ邮箱如何不显示图片? QQ邮箱如何不显示图片?...
- Linux运维部落 DNS服务,Linux运维部落
- C语言实现猜数字小游戏(详细讲解)
- this.className的使用
- 使用 ALT 键输入特殊符号
- 华为服务器bmc snmp协议,华为S5700开启SNMP协议
- Android手势密码view笔记(一)
热门文章
- java 日期calendar_java成神之——date和calendar日期的用法
- mysql嵌套查询效率低,连接查询代替嵌套查询提高select效率
- 如何做好技术布道——用影响影响影响
- 【具体实现过程】百度地图POI区域的边界经纬度坐标数据采集
- python自动投递_利用python如何在前程无忧高效投递简历
- JS 获取星期几的四种写法(转)
- Java开发面试题目,java培训老师第一次讲课
- 【2017 BSUIR Semifinal D】Friends rescue 题解
- GitChat·管理 | 一篇文章读懂项目管理中的精髓
- android重力传感器横竖反,Android重力感应实现方式