eigen的简单用法汇总
Eigen帮助文档的地址:http://eigen.tuxfamily.org/dox/pages.html
Eigen的论坛:http://forum.kde.org/viewforum.php?f=74
1.一些基本运算
#include <iostream>
using namespace std;
#include <ctime>
//核心部分
#include <Eigen/Core>
//稠密矩阵的运算
#include <Eigen/Dense>
using namespace Eigen;
#define MATRIX_SIZE 50
int main() {
//eigen中的所有向量和矩阵都是Eigen::Matrix,三个参数为数据类型,行,列
//声明一个2*3的float矩阵
Matrix<float,2,3> matrix_23;
Matrix3d m1; //旋转矩阵3*3 双精度,也可改为f
AngleAxisd m2; //旋转向量 3*1
Vector3d m3; //欧拉角 3*1
Quaterniond m4; //四元数 4*1
Isometry3d m5; //欧氏变换矩阵 4*4
Affine3d m6; //仿射变换4*4
Projective3d m7; //射影变换4*4
//Vector3d 本质上还是Eigen::Matrix<double,3,1>即三维向量
Vector3d v_3d;
Matrix<float,3,1> vd_3d; //类似
//本质上是Eigen::Matrix<double,3,3>
Matrix3d matrix_33 = Matrix3d::Zero(); //初始化为0
//不确定矩阵大小,使用动态矩阵
Matrix<double, Dynamic,Dynamic> matrix_dynamic;
//更简单的:
MatrixXd matrix_x;
matrix_23 << 1,2,3,4,5,6;
cout << matrix_23 << endl;
v_3d << 3,2,1;
vd_3d << 4,5,6;
//基本操作
cout << matrix_23.row(2); //取第二行的元素
cout << matrix_23.col(1); //取第一列的元素
matrix_23.dot(matrix_23); //点积,对应元素相乘后相加,结果为一个数
matrix_23.cross(matrix_23); //叉乘,结果为一个矩阵
//乘法,不同类型需要显性的转换
Matrix<double,2,1> result = matrix_23.cast<double>() * v_3d;
cout << "[1,2,3;4,5,6]*[3,2,1]=\n" << result << endl;
/*******矩阵运算*********/
matrix_33 = Matrix3d::Random();
cout << "random matrix33:\n" << matrix_33 << endl;
cout << "transpose:\n" << matrix_33.transpose() << endl; //转置
cout << "sum:" << matrix_33.sum() << endl; //各元素求和
cout << "trace:" << matrix_33.trace() << endl; //迹
cout << "times 10:\n" << matrix_33 * 10 << endl; //数乘
cout << "inverse:\n" << matrix_33.inverse() << endl; //逆
cout << "det:" << matrix_33.determinant() << endl; //行列式
/***********************/
//特征值
//实对称矩阵可以保证对角化成功
SelfAdjointEigenSolver<Matrix3d> eigen_solver(matrix_33.transpose() * matrix_33);
cout << "eigen values = \n" << eigen_solver.eigenvalues() << endl;
cout << "Eigen vectors = \n" << eigen_solver.eigenvectors() << endl;
//解方程
//求解matrix_nn * x = v_Nd这个方程
//直接求逆最直接,但是运算较大
Matrix<double,MATRIX_SIZE,MATRIX_SIZE> matrix_NN
= MatrixXd::Random(MATRIX_SIZE,MATRIX_SIZE);
matrix_NN = matrix_NN * matrix_NN.transpose(); //保证半正定
Matrix<double,MATRIX_SIZE,1> v_Nd = MatrixXd::Random(MATRIX_SIZE,1);
clock_t time_str = clock();
//直接求逆
Matrix<double,MATRIX_SIZE,1> x = matrix_NN.inverse() * v_Nd;
cout << "time is:" << 1000*(clock() - time_str) / (double) CLOCKS_PER_SEC << "ms" << endl;
cout << "x=" << x.transpose() << endl;
//QR分解,速度快很多
time_str = clock();
x = matrix_NN.colPivHouseholderQr().solve(v_Nd);
cout << "time is:" << 1000*(clock() - time_str) / (double) CLOCKS_PER_SEC << "ms" << endl;
cout << "x=" << x.transpose() << endl;
//对于正定矩阵,还可以用cholesky分解来解方程
time_str = clock();
x = matrix_NN.ldlt().solve(v_Nd);
cout << "time is:" << 1000*(clock() - time_str) / (double) CLOCKS_PER_SEC << "ms" << endl;
cout << "x=" << x.transpose() << endl;
//旋转
double theta = n * 2 * M_PI / (poseNums * 4); // 1/4 圆
R = Eigen::AngleAxisd(theta, Eigen::Vector3d::UnitZ());
return 0;
}
CMakeLists.txt
cmake_minimum_required(VERSION 2.8.3)
project (main)
set(CMAKE_BUILD_TYPE "Release")
set(CMAKE_CXX_FLAGS "-O3")
add_definitions(-std=c++11)
include_directories(inc)
aux_source_directory(src DIR_SRCS)
SET(SOUR_FILE ${DIR_SRCS})
include_directories("/usr/include/eigen3")
add_executable(main ${SOUR_FILE})
eigen的简单用法汇总的更多相关文章
- CATransition(os开发之画面切换) 的简单用法
CATransition 的简单用法 //引进CATransition 时要添加包“QuartzCore.framework”,然后引进“#import <QuartzCore/QuartzCo ...
- jquery.validate.js 表单验证简单用法
引入jquery.validate.js插件以及Jquery,在最后加上这个插件的方法名来引用.$('form').validate(); <!DOCTYPE html PUBLIC " ...
- NSCharacterSet 简单用法
NSCharacterSet 简单用法 NSCharacterSet其实是许多字符或者数字或者符号的组合,在网络处理的时候会用到 NSMutableCharacterSet *base = [NSMu ...
- [转]Valgrind简单用法
[转]Valgrind简单用法 http://www.cnblogs.com/sunyubo/archive/2010/05/05/2282170.html Valgrind的主要作者Julian S ...
- Oracle的substr函数简单用法
substr(字符串,截取开始位置,截取长度) //返回截取的字 substr('Hello World',0,1) //返回结果为 'H' *从字符串第一个字符开始截取长度为1的字符串 subst ...
- Ext.Net学习笔记19:Ext.Net FormPanel 简单用法
Ext.Net学习笔记19:Ext.Net FormPanel 简单用法 FormPanel是一个常用的控件,Ext.Net中的FormPanel控件同样具有非常丰富的功能,在接下来的笔记中我们将一起 ...
- TransactionScope简单用法
记录TransactionScope简单用法,示例如下: void Test() { using (TransactionScope scope = new TransactionScope()) { ...
- WPF之Treeview控件简单用法
TreeView:表示显示在树结构中分层数据具有项目可展开和折叠的控件 TreeView 的内容是可以包含丰富内容的 TreeViewItem 控件,如 Button 和 Image 控件.TreeV ...
- listActivity和ExpandableListActivity的简单用法
http://www.cnblogs.com/limingblogs/archive/2011/10/09/2204866.html 今天自己简单的总结了listActivity和Expandable ...
- SQL*Plus break与compute的简单用法
SQL*Plus break与compute的简单用法在SQL*Plus提示符下输出求和报表,我们可以借助break与compute两个命令来实现.这个两个命令简单易用,可满足日常需求,其实质也相当于 ...
随机推荐
- 【iOS逆向】某车之家sign签名分析
阅读此文档的过程中遇到任何问题,请关注公众号[移动端Android和iOS开发技术分享]或加QQ群[309580013] 1.目标 分析某车之家sign签名算法的实现 2.操作环境 frida mac ...
- SpringBoot使用@Async的总结!
一些业务场景我们需要使用多线程异步执行任务,加快任务执行速度. 之前有写过一篇文章叫做: 异步编程利器:CompletableFuture 在实际工作中也更加推荐使用CompletableFuture ...
- 命令查询职责分离 - CQRS
概念 CQRS是一种与领域驱动设计和事件溯源相关的架构模式, 它的全称是Command Query Responsibility Segregation, 又叫命令查询职责分离, Greg Young ...
- 在实际应用中联合体union的妙用
关键字union,又称为联合体.共用体,联合体的声明和结构体类似,但是它的行为方式又和结构体不同,这里的行为方式主要指的是其在内存中的体现,结构体中的成员每一个占据不同的内存空间,而联合体中的所有成员 ...
- 【kafka】connect的timestamp模式无法同一秒插入多条记录问题解决
一.现在问题 同时插入多条时间戳相同的记录 INSERT INTO "ABANK" VALUES ('1', 'CH', '00211', 'UBS Switzerland AG' ...
- 【每日一题】【排序sort重载】【工具类】2021年12月23日-31. 下一个排列
实现获取 下一个排列 的函数,算法需要将给定数字序列重新排列成字典序中下一个更大的排列(即,组合出下一个更大的整数). 如果不存在下一个更大的排列,则将数字重新排列成最小的排列(即升序排列). 必须 ...
- PW6276是一颗高效同步升压转换芯片,内部集成低阻抗功率 MOS
概述PW6276是一颗高效同步升压转换芯片,内部集成低阻抗功率 MOS. 具有短路保护功能内部集成软启动电路,无需外部补偿电容,外部反馈网络.PW6276采用 SOP8-EP 封装配合较少的外围原件使 ...
- 【JUC】交换器Exchanger详解
欢迎关注专栏[JAVA并发] 前言 JDK中提供了不少的同步工具,现在分享一个相对比较冷门的同步工具--交换器(Exchanger).你知道Exchanger的作用是什么吗?实现机制是什么?可以用来做 ...
- java (String)强制转换与toString()方法
1. Object.toString()介绍 Object中是自带有toString()方法的,也就是说java中的所有类的对象都是可以转换为字符串的. 首先,先看看Object.toString() ...
- CORS与CSRF在Spring Security中的使用
背景 在项目使用了Spring Security之后,很多接口无法访问了,从浏览器的网络调试窗看到的是CORS的报错和403的报错 分析 我们先来看一下CORS是什么,和它很相似的CSRF是什么,在S ...