对于图中这种厄米矩阵，所有的矩阵元都是实数，
，
我直接用Eigenvalues这个命令进行对角化，结果是Root形式的，命令是：
Hk = -{{0, 2 t*c1, 2 t*c2}, {2 t*c1, 0, 2 t*c12}, {2 t*c2, 2 t*c12,
0}};
Eigenvalues[Hk]
然后运行结果是这样的
加个Cubics命令就变成
然后用FullSimplify命令都无法继续化简下去，我也尝试了Assuming所有矩阵元都是实数等等技巧也都不行。
事实上这个矩阵的本征值应该是很简洁的形式，我在一份讲义中找到了答案是
所以很奇怪为什么mathematica给不出这种形式的解，或者是我没有进行适当的操作？
此外我用了matlab进行尝试居然发现可以得到正确答案。
，我直接用Eigenvalues这个命令进行对角化，结果是Root形式的，命令是：
Hk = -{{0, 2 t*c1, 2 t*c2}, {2 t*c1, 0, 2 t*c12}, {2 t*c2, 2 t*c12,
0}};
Eigenvalues[Hk]
然后运行结果是这样的
加个Cubics命令就变成
然后用FullSimplify命令都无法继续化简下去，我也尝试了Assuming所有矩阵元都是实数等等技巧也都不行。
事实上这个矩阵的本征值应该是很简洁的形式，我在一份讲义中找到了答案是
所以很奇怪为什么mathematica给不出这种形式的解，或者是我没有进行适当的操作？
此外我用了matlab进行尝试居然发现可以得到正确答案。
