C++之 Eigen-3.4.0 全方位教程：Chapter04-操作篇

最新推荐文章于 2025-06-03 08:58:53 发布

Jacob-xyb

最新推荐文章于 2025-06-03 08:58:53 发布

阅读量1.4k

点赞数 2

CC 4.0 BY-SA版权

分类专栏： C++ 文章标签： c++ 线性代数 Eigen

本文链接：https://blog.csdn.net/weixin_44560698/article/details/120287830

这篇博客详细介绍了Eigen库中矩阵块操作的使用，包括动态和静态大小的块表达式，以及左值和右值操作。内容涵盖块的读写、单独行和列的访问、角块操作、向量的块操作、切片和索引等。通过示例代码展示了如何高效地操作矩阵的特定部分，并解释了不同操作对性能的影响。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

块操作

块操作既可以作为右值，还可以作为左值！

基本操作

Block operation	dynamic-size block expression	fixed-size block expression
Block of size `(p,q)`, starting at `(i,j)`	`matrix.block(i,j,p,q);`	`matrix.block<p,q>(i,j);`

两个版本都可用于固定大小和动态大小的矩阵和数组。这两个表达式在语义上是等价的。

唯一的区别是，如果块大小很小，固定大小的版本通常会给你更快的代码，但需要在编译时知道这个大小。

void Eigen_Introduction_Block()
{
   
   
	Eigen::MatrixXf m(4, 4);
	m << 1, 2, 3, 4,
		 5, 6, 7, 8,
		 9, 10, 11, 12,
		 13, 14, 15, 16;
	cout << "Block in the middle" << endl;
	cout << m.block<2, 2>(1, 1) << endl << endl;
	for (int i = 1; i <= 3; ++i)
	{
   
   
		cout << "Block of size " << i << "x" << i << endl;
		cout << m.block(0, 0, i, i) << endl << endl;
	}
}

输出：

Block in the middle
 6  7
10 11

Block of size 1x1
1

Block of size 2x2
1 2
5 6

Block of size 3x3
 1  2  3
 5  6  7
 9 10 11

块的左值操作

void Eigen_Introduction_Block001()
{
   
   
	Array22f m;
	m << 1, 2,
		 3, 4;
	Array44f a = Array44f::Constant(0.6);
	cout << "Here is the array a:" << endl << a << endl << endl;
	a.block<2, 2>(1, 1) = m;
	cout << "Here is now a with m copied into its central 2x2 block:" << endl << a << endl << endl;
	a.block(0, 0, 2, 3) = a.block(2, 1, 2, 3);
	cout << "Here is now a with bottom-right 2x3 block copied into top-left 2x3 block:" << endl << a << endl << endl;
}

输出

Here is the array a:
0.6 0.6 0.6 0.6
0.6 0.6 0.6 0.6
0.6 0.6 0.6 0.6
0.6 0.6 0.6 0.6

Here is now a with m copied into its central 2x2 block:
0.6 0.6 0.6 0.6
0.6   1   2 0.6
0.6   3   4 0.6
0.6 0.6 0.6 0.6

Here is now a with bottom-right 2x3 block copied into top-left 2x3 block:
  3   4 0.6 0.6
0.6 0.6 0.6 0.6
0.6   3   4 0.6
0.6 0.6 0.6 0.6