sage之线性代数

#定义矩阵，默认定义在实数域
A = Matrix([[1,2,3,5],[3,2,1,2],[1,1,1,0],[3,7,2,2]])
A^-1
#定义在其他域上的矩阵，如有限域
A = Matrix(GF(13),[[1,2,3,5],[3,2,1,2],[1,1,1,0],[3,7,2,2]])
A^-1
#可以看到两个逆矩阵不一样

#定义向量，定义在有限域，默认定义在实数域
w = vector(GF(13),[1,1,4,3])
Y=A*w;Y
Z=w*A;Z

#解线性方程组
X = A.solve_right(Y);X
#也可以使用符号\
A\Y

X = A.solve_left(Z);X

#格基约减
A = Matrix([[1,2,3,5],[3,2,1,2],[1,1,1,0],[3,7,2,2]])
A.LLL()