sage中的矩阵运算

sage中的矩阵运算

sage 是基于python的数学工具，在没有类似mathematics这样的软件环境时可以直接拿来用。

下面简单介绍一下其中的矩阵运算：

矩阵定义

a=matrix([[1,1], [0,4]])
a

结果为：

[1 1]
[0 4]

矩阵加法和乘法

m1=matrix([[1,2], [1,3]])
m2=matrix([[3,4], [5,6]])
m1+m2
m1*m2

加法结果为：

[4 6]
[6 9]

乘法结果为：

[13 16]
[18 22]

矩阵的转置和逆

a.transpose()
a.inverse()

结果为：

[1 0]
[1 4]

和：

[   1 -1/4]
[   0  1/4]

矩阵的特征系统

a.eigenvalues()
a.eigenvectors_right()

结果为：

[4, 1]

和：

[(4, [
  (1, 3)
  ], 1), (1, [
  (1, 0)
  ], 1)]

矩阵的行列式

a.det()

结果为：

4

线性方程求解

b=vector(RDF,[1,2])
a.solve_right(b)

结果为：

(1/2, 1/2)

小结

通过对矩阵运算的初步尝试，掌握其用法，那么也可以类似mathematics那样运算，特别是在jupyter中运行时。

参考文献

1. sage 的帮助文档：matrices.PDF