我们利用QR分解来计算矩阵的特征值。这一方法是迭代的,并生成一个上三角矩阵,特征值即为该矩阵的对角元素。我们的发现显示,这些特征值与Matlab内置函数eig计算结果一致。您可以在以下链接中找到类似用Mathematica实现的程序:http://library.wolfram.com/infocenter/MathSource/6612/