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)]
相关阅读
怎么拟合这样的曲线(ZT)中国科技大学研制成功世界上第一台量子计算机请问matlab plot如何动态更新关于计算机找工作，有经验师兄师姐请帮下小弟如何实现Matlab作图中的双纵轴Ph.D. positions in Computer Science (with a focus on sensor and wireless networks)a question about random number generator怎么设计boxcar的highpass filter？请问insurance company programmer的工资比一般电脑公司来说是高还是低？有没有能在单台机子上调试mpi程序的simulator？请问matlab如何创建结构体符号变量啊？do you use blas/lapack?想买新机器，做计算用，大家推荐一个配置吧Can we get the data from a figure in the paper?如何快速保存matlab 能不能解符号矩阵方程？ (转载) 如何在wgr614v4后面建立一个ftp how to rip a dvdQuick question: BaoZis are rewared.祝大家节日快乐!