我想用符号而不是数字来计算线性代数运算。例如,有以下矩阵:
A = matrix(QQ,[
[ 2, 1, 2, -6],
[-1, 2, 1, 7],
[ 3, -1, -3, -1],
[ 1, 5, 6, 0],
[ 2, 2, 1, 1]
])
我想对带有符号变量的向量进行乘法运算,如下所示:
t = 'real'
var('x1')
assume(x1, t)
var('x2')
assume(x2, t)
var('x3')
assume(x3, t)
var('x4')
assume(x4, t)
xx = vector(QQ, [x1, x2, x3, x4])
A * xx.transpose()
不幸的是,构建xx
向量不成功,产生以下错误消息:
TypeError: Cannot evaluate symbolic expression to a numeric value.
这不起作用,那么我如何在Sage的线性代数框架中使用符号?
一种解决方案是处理多项式环中的项。
sage: R.<x1,x2,x3,x4> = PolynomialRing(QQ)
sage: R
Multivariate Polynomial Ring in x1, x2, x3, x4 over Rational Field
然后用R
中的坐标定义你的向量。
sage: xx = vector(R,[x1,x2,x3,x4])
sage: A * xx
结果是另一个在R
中有条目的向量。
另一个解决方案是在符号环SR
中工作,正如@kcrisman在ask-sage上指出的。