R*k
2 楼
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行吗？
谢谢！
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行吗？
谢谢！
w*e
3 楼
没人知道么
l*n
4 楼
需要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)]
【在 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)]
相关阅读
请教各位高手有限元的问题?Re: 请问normal disctribution的cdf函数在FORTRAN[转载] 急问: SAS proc iml 里面...[转载] 面试问题！A question about DQS3.3.2请教一个Patran问题求不等式组是否有解的问题:[转载] 我的论文标题的问题谁有hypermesh？We are building new cluster now!measure theory, probably theory, real anFortran code optimizationHelp! reverse Cuthill-McKee reordering..Abaqus material question?请教ansys高手一个泊松比的问题Re: How to save variables in file in MATMS interview quiz不定积分请教一个作图问题菜鸟问个有限元的弱问题