eigen矩阵默认不初始化，元素为随机值；应显式调用.setzero()、.setones()或使用matrixxd::zero(3,4)等静态方法初始化，避免未定义行为。
如何用 Eigen 初始化矩阵并避免内存未初始化问题
Eigen 默认不初始化矩阵元素，直接声明的矩阵值是随机的，这在调试时极难发现，但会导致数值计算结果不可复现。必须显式调用 .setZero()、.setOnes() 或带初始化列表的构造函数。
-
MatrixXd A(3, 4);—— 元素全为垃圾值，不能直接参与计算 -
MatrixXd A = MatrixXd::Zero(3, 4);—— 安全，推荐用于需要零初值的场景 VectorXd b(5); b —— 使用逗号初始化语法前，必须先分配大小，否则运行时报错 <code>Assertion failed: cols() == x.cols()- 动态矩阵（
MatrixXd,VectorXd）比固定尺寸（Matrix3d,Vector4f）稍慢，但更灵活；若尺寸已知且小，优先用固定尺寸类型提升性能和编译期检查
矩阵乘法写错顺序或混用转置会直接导致维度错误
Eigen 的乘法是严格按数学定义的：左操作数列数必须等于右操作数行数。常见错误不是语法报错，而是运行时断言失败，比如 assertion failed: rows() == other.rows()——这其实是你误把 A * B.transpose() 写成了 A * B，而 B 是列向量。
-
A * B要求A.cols() == B.rows()；若想做“每行点积”，常需A * B.transpose()或A.rowwise().dot(B) -
.transpose()返回的是表达式对象，不是新矩阵；链式调用如A.transpose() * B.transpose()是合法的，但A.transpose().eval() * B才真正触发拷贝——多数情况不需要.eval()，除非你明确要打断表达式模板优化 - 对称矩阵运算（如
A.transpose() * A）建议改用A.adjoint() * A（复数安全），即使当前是实数，也避免未来扩展时出错
解线性方程组时别直接用 .inverse()
A.inverse() * b 看起来直观，但效率低、数值不稳定，且当 A 接近奇异时可能返回毫无意义的结果而不报错。Eigen 提供更优的分解接口，应根据矩阵性质选：
C++ Code Review Master
下载
组合式C++代码评审方案，融合静态分析、AI推理、多轮迭代评审和C++专项检查，适用于PR审查、增量代码审查、全项目评审和代码质量评分，触发词包括review cpp、cpp代码评审、C++review、代码审查。
- 一般稠密矩阵 → 用
A.colPivHouseholderQr().solve(b)（稳定、快、自动处理秩亏） - 对称正定矩阵 → 用
A.llt().solve(b)（Cholesky，速度最快，但不验证正定性） - 超定最小二乘 → 用
A.jacobiSvd(ComputeThinU | ComputeThinV).solve(b)（最鲁棒，适合病态系统） - 所有这些方法返回
VectorXd或VectorXd::NullaryExpr类型，可直接赋值；若解不存在，solve()不抛异常，而是返回一个“尽力而为”的结果——务必用info() == Success检查，例如：auto x = A.colPivHouseholderQr().solve(b); if (x.size() == 0 || A.colPivHouseholderQr().info() != Success) { /* 处理失败 */ }
编译时链接和头文件包含的硬性要求
Eigen 是纯头文件库，不需要链接 .a/.so，但有两处极易出错：
- 必须 #include
<eigen></eigen>（基础矩阵/向量）或<eigen></eigen>（稀疏矩阵），而不是只写<eigen3></eigen3>或<eigen></eigen>—— 后者不包含MatrixXd等常用类型定义 - CMake 中只需
find_package(Eigen3 REQUIRED)+target_include_directories(my_target PRIVATE ${EIGEN3_INCLUDE_DIR})；如果手动指定路径，请确认${EIGEN3_INCLUDE_DIR}指向的是 eigen3 目录的父目录（即包含eigen3/子目录的那个路径），否则#include <eigen></eigen>会找不到 - 启用
-DEIGEN_MPL2_ONLY编译宏可禁用部分 Boost 兼容代码，减小模板膨胀，但非必需；不用它也不会出错
事情说清了就结束。Eigen 的坑不在语法多难，而在默认行为反直觉、错误不报在编译期、运行时断言又容易被忽略——尤其 .inverse() 和未初始化矩阵，这两个点几乎每个新人都踩两次。
C++免费学习笔记（深入）：立即使用
在学习笔记中，你将探索 C++ 的入门与实战技巧！
