Eigen库解二元一次方程(不懂矩阵操作的范例

https://eigen.tuxfamily.org/index.php?title=Main_Page
下载Eigen库

包含目录操作
在这里插入图片描述

下载完解压,打开到如图所示的路径,复制它,丢编译器里
在这里插入图片描述
首先,列出要解的方程,一定要写成这种格式的:
2x - y =1
-x + 3y = 3

等式左边的常数装进Matrix2f A
等式右边的常数装进Vector2f b
结果是Vector2f x

#include <iostream>
#include<Eigen>
using namespace std;
using namespace Eigen;
int main()
{
   Matrix2f A;
   Vector2f b;
    A << 2, -1, -1, 3;
    b << 1,  3;
    cout  <<"Here is the matrix A:\n" << A << endl;
    cout << "Here is the right hand side b:\n" << b << endl;
    Vector2f x = A.colPivHouseholderQr().solve(b);
    // Vector2f x =A.ldlt().solve(b);也可以
    cout << "The solution is:\n" << x << endl;
    return 0;

}

在这里插入图片描述

解三元一次方程的接口是一样的,但是要把容器类型改一下,Matrix2f 改 Matrix3f,Vector2f 改 Vector3f

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值