SageMath矩阵操作及解线性方程组
矩阵所处环的表示
- ZZ:整数环
- QQ:有理数环
- Zmod(p):p为素数,定义在ZPZ_PZP上
矩阵操作
声明
超完整定义
mt = matrix(QQ, 3, 2, [1, 2, 3, 4, 5, 6])
完整定义
mt = matrix(ZZ, 3, 4) #定义一个在整数环上的3行4列的矩阵
mt = matrix(ZZ, 0, 4) #定义一个在整数环上0行4列的矩阵,后期可动态添加行
部分定义
mt = matrix(QQ, 3, [1, 2, 3, 4, 5, 6]) #定义在有理数环上的3列矩阵,并赋值
添加行
在末尾添加
mt = mt.stack(vector([1, 2, 3]) #矩阵后添加一行
在指定行插入一行,第一个参数为行号,第二个为元素list
mt = mt.insert_row(1, vector([1, 2, 3])) #在第一行插入
修改元素值
注意,矩阵行列均从0开始!
修改单个值
mt[2, 3] = 8 #将第3行第4列的值修改为8
修改整行值
mt[2] = [1, 2, 3, 4] #将第3行的值依次修改为list中的值
求秩
mt.rank()
求行、列值
mt.nrows() #将返回矩阵行数
mt.ncols() #将返回矩阵列数
求行列式的值
mt.determinant()
求是否为可逆矩阵
是可逆矩阵则输出TRUE,否则输出FALSE
mt.is_invertible()
求逆
两种方法
mt.inverse()
mt = mt^(-1)
改变矩阵所处环
mt.change_ring(ZZ) ##将矩阵的环改为ZZ
求解线性方程组和线性相关组
求矩阵X满足XA=B
X = A.solve_left(B)
求矩阵X满足AX=B
X = A.solve_right(B)
求向量x满足Ax=b
两种方法
x = A \ b
x = A.solve_right(b)
找线性相关的行向量
满足XA=0
X = A.left_kernel()
满足AX=0
X = A.right_kernel()
找X中指定行向量
X.gen(0) #解X中的第1行
求最短正交基
LLL格基约简算法
mt.LLL()
SageMath矩阵操作及解线性方程组相关推荐
- [CF917D]Stranger Trees[矩阵树定理+解线性方程组]
题意 给你 \(n\) 个点的无向完全图,指定一棵树 \(S\),问有多少棵生成树和这棵树的公共边数量为 \(k\in[0,n-1]\) \(n\leq 100\) 分析 考虑矩阵树定理,把对应的树边 ...
- 解线性方程组的python实现(2)——矩阵三角分解法
解线性方程组的python实现2--矩阵三角分解法 1. 矩阵三角分解法 实现代码 2. LU分解 2.1 基本步骤 2.2 LU分解的计算公式 2.3 LU分解的结果表示 实现代码 3. 选主元的L ...
- 科学计算机解多元方程组,南京大学计算机科学与技术系 数值计算方法(第3章)3 向量和矩阵的范数 病态方程组 解线性方程组的迭代法.ppt...
南京大学计算机科学与技术系 数值计算方法(第3章)3 向量和矩阵的范数 病态方程组 解线性方程组的迭代法.ppt 上传人:ct****5 文档编号:43116102 上传时间:2020-01-21 格 ...
- 第三章 矩阵的初等变换与线性方程组概念详解
矩阵的初等变换 一.初等变换和初等矩阵及其联系 设A=(aij)m*n,则以下三种变换称为矩阵A的初等行(列)变换 1.交换A的两行(列) 2.用一个非零常数k乘以A的某一行(列) 3.用一个数乘以A ...
- 范德蒙德矩阵在MATLAB中怎么表示,Python 之 Python与MATLAB 矩阵操作总结
Python 之 Python与MATLAB 矩阵操作小结 一.线形代数理论基础 线形代数(linear algebra)是数学的一个分支,研究矩阵理论.向量空间.线性变换和有限维线形方程组等内容. ...
- 高斯消去法解线性方程组C++实现
一.问题分析:高斯消去法解线性方程组主要面对的问题是矩阵的运算,所以可以定义一个矩阵类Matrix类作为基类,然后由矩阵类Matrix类派生出一个线性方程组类LinearEqu类. 二.矩阵功能分析: ...
- 大数据矩阵计算基础(二):R、Excel、Matlab中矩阵操作
一.R语言 构造矩阵 a=matrix(c(1,2,3,4),ncol=2,byrow=T) b=matrix(c(5,6,7,8),ncol=2,byrow=T) c=matrix(c(1,2,3, ...
- 【数值计算】花式解线性方程组
Assignment1 LU分解 codes 参考文献 Assignment2 迭代法 Problem 数据生成 norm slow_Jacobi 清真_Jacobi Gauss-Seidel SOR ...
- OpenCV初学(5)--------矩阵操作
矩阵操作 基本的矩阵和图像算子 计算矩阵中所有元素的绝对值 cv::abs() 计算两个矩阵差值的绝对值 cv::absdiff() 实现两个矩阵逐元素相加 cv::add() 实现两个矩阵逐元素加权 ...
- MATLAB2016笔记(五):进阶矩阵操作
文章目录 一.特殊矩阵创建 (一)希尔伯特矩阵 hilb() (二)托普利兹矩阵 toeplitz() (三)0~1间均匀分布的随机矩阵 rand() (四)标准正态分布随机矩阵 randn() (五 ...
最新文章
- C语言的内存分配模型
- 故障码123401_电力系统规划设计对电力工程设计的应用
- Spring的两种属性注入方式
- itext设置pdf的尺寸_如何获取pdf文档iText 7的页面大小-问答-阿里云开发者社区-阿里云...
- ubuntu下vscode使用cmake编译运行c++配置文件
- es6 __proto__属性,Object.setPrototypeOf(),Object.getPrototypeOf()
- 正确使用招呼语、称呼语_5
- 腾讯广告招人啦,校招优质岗位你不容错过!
- 变量声明提升和函数声明提升
- WIN32_FIND_DATA详解
- window10官方工具在线升级失败,终极解决方案
- html实现点击下载文件
- 写作的必备技能,markdown的基本语法
- EDA程序设计--数字日历电路
- 【蓝桥杯每日一练:木头加工】
- 用虚拟打印机MacroMedia FlashPaper 制作swf文件
- knuth《mathematical analysis of algorithms》学习笔记
- 既然决定了卷论文,那就开卷吧!
- 创客筑梦,第二届机智云智能硬件36小时开发大赛圆满落幕
- 【物理应用】内联全息图外推算法matlab代码
热门文章
- layui树形组件渲染
- 利润分配计算机公式,利润与利润分配表计算公式大全
- Linux初级入门百篇-lsof工具
- Vscode信任问题
- NOI / 1.3编程基础之算术表达式与顺序执行——12:计算球的体积
- FFmpeg 加水印 加马赛克
- ntp版本导致时间同步不成功问题
- set_global_opts在PyEcharts中实现全局配置项
- 如何更改xp计算机用户名和密码错误,XP更改开机密码提示“Windows不能更改密码”如何解决...
- 2021-07-21