计算对称矩阵的实特征值和特征向量

计算对称矩阵的实特征值和特征向量

在线性代数中，特征值和特征向量是矩阵的一种基本性质。求解对称矩阵的特征值和特征向量是一个非常重要且普遍的问题，有很多种方法可以解决这个问题，其中一种是利用Jacobi方法来求解。

Jacobi方法，是一种通过正交相似变换将任意实矩阵对角化的方法，其关键步骤是通过相似变换将矩阵转化为对角矩阵。在矩阵相似的过程中，特征值不变，特征向量也随之发生变换，因此可以通过这个方法求解矩阵的特征值和特征向量。

下面是使用C++编写的求解对称矩阵的特征值和特征向量的程序，该程序采用了Jacobi方法。其中包括了完整的源代码和相应的注释。

#include <iostream>
#include <cmath>

const double eps = 1e-8; // 终止条件，当元素绝对值小于该值时认为为0

template <class T>
class SymmetricMatrix {
private:
    int size; // 矩阵大小
    T *data; // 存储矩阵元素的数组
public:
    SymmetricMatrix() : size(0), data(NULL) {}
    SymmetricMatrix(int sz) : size(sz) {
        data = new T[size*(size+1)/2];
    }
    ~SymmetricMatrix() { delete[] data; }

    T& operator () (int i, int j) { // 重载()运算符，用来取矩阵元素