R*k
发帖数: 264 | 1 1. 有两个矩阵 X 和 Y，如果求迹
Tr[X.Y.Y.X-X.X.Y.Y]
由于求迹时矩阵的可置换性，上面的结果应该为0。但是Mathematica却只能原样输出。
2. 如果 X 和 Y 做变换到矩阵 A 和 B,如
X=A+B
Y=A-B
然后将结果展开，用 A 和 B 输出：
Expand[Tr[X.Y.Y.X-X.X.Y.Y]]
但是输出结果却为：
Tr[(A + B).(A - B).(A + B).(A - B) - (A + B).(A + B).(A - B).(A - B)]
完全没有做简化。
如果X，Y，A，B不是矩阵，而是一般的变量，就没有问题。
难道mathematica没法做到这一点吗？如果是这样，maple或matlab行吗？
谢谢！ | l******n
发帖数: 9344 | 2 需要simplify
【在 R*k 的大作中提到】
: 1. 有两个矩阵 X 和 Y，如果求迹
: Tr[X.Y.Y.X-X.X.Y.Y]
: 由于求迹时矩阵的可置换性，上面的结果应该为0。但是Mathematica却只能原样输出。
: 2. 如果 X 和 Y 做变换到矩阵 A 和 B,如
: X=A+B
: Y=A-B
: 然后将结果展开，用 A 和 B 输出：
: Expand[Tr[X.Y.Y.X-X.X.Y.Y]]
: 但是输出结果却为：
: Tr[(A + B).(A - B).(A + B).(A - B) - (A + B).(A + B).(A - B).(A - B)]
|
|
