z444_579的博客

最近需要用 C++ 做一些数值计算，之前一直采用Matlab 混合编程的方式处理矩阵运算，非常麻烦，直到发现了 Eigen 库，简直相见恨晚，好用哭了。 Eigen 是一个基于C++模板的线性代数库，直接将库下载后放在项目目录下，然后包含头文件就能使用，非常方便。此外，Eigen的接口清晰，稳定高效。唯一的问题是之前一直用 Matlab，对 Eigen 的 API 接口不太熟悉，如果能有 Ei

double dbX[9] = {1,2,3,4,5,6,7,8,9};//这里是列优先...MatrixXd T2 = Map(dbX,3,3);//Matrix中的数据存在一维数组中，默认是行优先的格式，即一行行的存  //data()返回Matrix中的指针 dymMat.data();//静态矩阵，编译时确定维数 Matrix Eigen

// 利用Eigen 求解非线性最小二乘;// 示例：y = 10*(x0+3)^2 + (x1-5)^2#include "math.h"#include "iostream"#include "vector"#include "list"using namespace std;#include "Eigen/Dense"#include "Eigen/Core"#inc

// 目标函数：f(x1, x2, x3, x4) = (x1 + 10*x2)^2 + 5*(x3 - x4)^2 + (x2 - 2*x3)^4 + 10.0*(x1 - x4)^4; 求取输入向量x?// 牛顿法： x(k+1) = x(k) - inv(F(x(k))) * g(x(k)), F是黑塞矩阵，g是导数;double CPowell::ComputePowellFunc(f

#include "math.h"#include "iostream"#include "vector"#include "list"using namespace std;#include "Eigen/Dense"#include "Eigen/Core"#include "Eigen/QR"#include "Eigen/SVD"#include "Eigen/Ge