Eigen 求解线性方程组

#include <iostream>
#include <Eigen/Dense>
#include <Eigen/Cholesky>
#include <Eigen/LU>
#include <Eigen/QR>
#include <Eigen/SVD>  using namespace std;
using namespace Eigen;  int main()
{      //线性方程求解 Ax =B;  //求解四元一次方程组，Matrix4d 为double类型4*4矩阵，Vectpr4d 为double类型4*1向量,Answer=[1;2;3;4]Matrix4d A;      A <<  1    , 2   ,  3   ,  4,1    , 4   ,  3    , 2,1   ,  3  ,   2    , 4,4   ,  1   ,  1   ,  3;      Vector4d B(30,26,29,21);      Vector4d x1 = A.colPivHouseholderQr().solve(B); //                                      right Answer     Vector4d x2 = A.llt().solve(B);      Vector4d x3 = A.ldlt().solve(B);Vector4d x4 = A.ldlt().solve(B);// A sym. p.s.d.    #include <Eigen/Cholesky>           wrong AnswerVector4d x5 = A.llt() .solve(B);  // A sym. p.d.      #include <Eigen/Cholesky>         wrong AnswerVector4d x6 = A.lu()  .solve(B);  // Stable and fast. #include <Eigen/LU>               right Answer//Vector4d x7 = A.qr()  .solve(B);  // No pivoting.     #include <Eigen/QR>             wrong Answer//Vector4d x8 = A.svd() .solve(B);  // Stable, slowest. #include <Eigen/SVD>            wrong Answerstd::cout <<A<<std::endl;std::cout << "The solution is:\n" << x1 <<"\n\n"<<x2<<"\n\n"<<x3 <<"\n\n"<<x4 <<"\n\n"<<x5 <<"\n\n"<<x6 <<std::endl;  //求解三元一次方程组,Matrix3d 为double类型3*3矩阵，Vectpr3d 为double类型3*1向量, Answer=[1;2;3]Matrix3d A2;      A2 <<  1,2,3,3,1,2,1,2,0;      Vector3d B2(14,11,5);      Vector3d x12 = A2.colPivHouseholderQr().solve(B2); //right Answer    Vector3d x62 = A2.lu().solve(B2);                  // right Answerstd::cout <<A2<<std::endl;std::cout << "The solution is:\n" << x12 <<"\n\n"<<x62 <<std::endl;  }

Eigen 求解线性方程组相关推荐

Eigen求解数学问题（一）
Eigen求解数学问题(一) 有这样的一个问题: 求解过程: 利用Eigen 求解: 输出总结: 有这样的一个问题: 三个朋友A,B,C各饲养家禽,A养鸡,B养鸭,C养兔. 他们同意按照下面的比例分 ...
分别用雅可比（Jacobi）迭代法和高斯—塞德尔（Gauss—Seidel）迭代法求解线性方程组(转载)
版权声明:本文为博主原创文章,未经博主允许不得转载. https://blog.csdn.net/yyywww666/article/details/42805071 算法介绍(迭代法介绍): 代码C ...
python QR分解求解线性方程组和矩阵本征值和本征向量
下面的代码提供了两个函数 solve_linear_equ, 利用QR分解求解线性方程组,输入是一个二维的非奇异的系数方阵和一个常数array,输出是该线性方程组的解 eigen, 输入是一个实方阵( ...
雅克比（Jacobi）迭代法求解线性方程组
长博文不利于翻阅,于是又将Jacobi迭代法单独出来了. 这篇博文把高斯-赛德尔迭代法和雅克比迭代法都放到一起了,个人觉得看着有点累.(迭代法求解线性方程组),不过还是要看的,因为它引出了迭代法. 进 ...
迭代法求解线性方程组的收敛问题总结
本讲之前,先将高斯-赛德尔迭代法和雅克比迭代法以及迭代法求解线性方程组贴出来,毕竟收敛问题研究的是迭代方法的收敛问题. 进入主题: 判断迭代法收敛的办法: 1.首先根据方程组的系数矩阵A的特点判断: ...
Gauss-Seidel迭代求解线性方程组
高斯-赛德尔迭代法考试比较多,所以考虑再三,还是单独提取出来独立一篇,方便查阅,突出重点. 首先举例引入: 通过手动求解下面的线性方程组得到精确解: 再用高斯-赛德尔迭代法求解比较: 本人拙见,将每一 ...
Guass列主元消去法求解线性方程组
上篇博文讲到:Guass消去法求解线性方程组,也提到了此方法求解线性方程组存在的问题,因此有如下: 基本思想: 看不懂没关系,直接看下面的例题,会有更直观的理解,然后再回到这里继续看. 选主元的步骤( ...
Guass消去法求解线性方程组
基本思想: 具体消元过程: 回代解方程组: 高斯顺序消去法求解线性方程组的计算公式: 存在的问题: 这种方法不是很实用,但是通常讲求解线性方程组都会先提到这种方法,这是基础,还是需要看看.
【原创】开源Math.NET基础数学类库使用(06)直接求解线性方程组
阅读目录前言 1.数值分析与线性方程 2.Math.NET解线性方程源码分析 3.Math.NET求解线性方程的实例 4.资源本博客所有文章分类的总目录:[总目录]本 ...

Eigen 求解线性方程组

Eigen 求解线性方程组相关推荐

最新文章

热门文章