关于Eigen库ConjugateGradient是否支持Hermitian矩阵的技术咨询

你提的这个问题我之前也碰到过，一开始也摸不着头脑——明明文档标注支持Hermitian矩阵，怎么测试就出问题了？先给你明确结论：Eigen的ConjugateGradient确实支持Hermitian（厄米特）矩阵，大概率是你的测试代码漏掉了关键细节，下面给你拆解几个容易踩的坑：

1. 先确认矩阵是严格Hermitian的 #

厄米特矩阵必须满足 A == A.adjoint()（实部对称，虚部反对称），如果构造矩阵时虚部符号出错，或者存在数值精度偏差，求解器的行为就会异常。建议测试前用 A.isApprox(A.adjoint()) 验证矩阵的厄米特性质。

2. 必须显式指定模板参数的UpLo选项 #

Eigen的ConjugateGradient默认模板参数是针对实对称正定矩阵的，要让它正确处理复厄米特矩阵，你需要显式声明矩阵的存储三角部分：

• 用Eigen::Lower表示矩阵的下三角（含对角线）存储了厄米特矩阵的信息
• 用Eigen::Upper表示上三角部分

给你一个能正常运行的最小示例：

#include <Eigen/Core>
#include <Eigen/IterativeLinearSolvers>
#include <iostream>

int main() {
    // 构造3x3厄米特矩阵
    Eigen::MatrixXcd A(3,3);
    A << 2, 1+0.5i, 0.5-1i,
         1-0.5i, 3, 1i,
         0.5+1i, -1i, 4;

    // 验证厄米特性质
    assert(A.isApprox(A.adjoint()));

    Eigen::VectorXcd b(3);
    b << 1, 2, 3;

    // 显式指定Lower，匹配矩阵存储方式
    Eigen::ConjugateGradient<Eigen::MatrixXcd, Eigen::Lower> cg;
    cg.compute(A);
    Eigen::VectorXcd x = cg.solve(b);

    std::cout << "Solution:\n" << x << std::endl;
    std::cout << "Residual norm: " << (A*x - b).norm() << std::endl;
    return 0;
}

3. 检查收敛设置是否合理 #

如果矩阵条件数很差，求解器可能提前终止迭代，导致结果看起来“不正常”。可以通过以下方式调整：

• cg.setTolerance(1e-8)：设置更严格的收敛阈值（默认是1e-5）
• cg.setMaxIterations(1000)：增加最大迭代次数

文档里标注的“selfadjoint (hermitian)”是准确的，厄米特矩阵就是复数域上的自伴矩阵，Eigen的实现完全支持，避开上面这些坑就能正常运行啦。

内容的提问来源于stack exchange，提问作者mo_cygorek