Mathematica的“精确解”和“数值解”得到的本征矢为何不一样？(3)

Mathematica的“精确解”和“数值解”得到的本征矢为何不一样

通过Mathematica求解矩阵的Pfaffian，得到

Mathematica求解反对称矩阵的Pfaffian在Wolfram函数资源库中，具体见

https://resources.wolframcloud.com/FunctionRepository/resources/Pfaffian

相关的参考文献见

Wimmer, M., “Algorithm 923: Efficient Numerical Computation of the Pfaffian for Dense and Banded Skew-Symmetric Matrices.” ACM Trans. Math. Softw., vol. 38, no. 4, 1–17, 2012. DOI: 10.1145/2331130.2331138.

数值解

通过上边“精确解”的计算，我们看到一路走下来没有问题！接下来，看一下“数值解”，我这里把正格子基矢用N这个函数数值化。

可以看到，后续相关结果都会数值化，如下