在Mathematica中数值地找到给定范围内的大符号矩阵的能量本征值

Question

我有一个136x136哈密顿量（矩阵），需要找到它的特征值。

由于需要求解136阶多项式，因此无法用特征值[H]进行分析求解。

我需要在计算其特征值之前用值替换符号项来以数字方式求解。然而，需要针对符号术语的值的范围绘制它，例如-1 <1。 x＆lt; 1.

是否有一种方法可以数值求解并绘制一系列值？

{
 {10.1358 - 6.72029 x, 0., 0.},...
 {0., 10.1358 - 6.72029 x, 0.},
 {0., 0., 10.1358 - 6.72029 x},
 {0., 0., 0.},
 {0., 0., 0.},
 {0., 0., 0.},
 {0., 0., 0.},
 {0., 0., 0.},
 {0., 0.204252, 0.},
 {0., 0., 0.267429}
...

矩阵的角落作为例子。矩阵是真实的，对称的。