Eigen 学习之块操作

　　Eigen 为 Matrix 、Array 和 Vector提供了块操作方法。块区域可以被用作左值和右值。在Eigen中最常用的块操作函数是 .block() 。

block() 方法的定义如下：

　　block of size (p,q) ,starting at (i,j)。matrix.block(i,j,p,q); matrix.block<p,q>(i,j);

上述两种形式都可以被用在固定大小和动态大小的矩阵中。

　　举例如下：

 #include <iostream>

 #include <eigen3/Eigen/Dense>

 using namespace Eigen;

 using namespace std;

 int main(int argc ,char** argv)

 {

     MatrixXf m(,);

     m << ,,,,

          ,,,,

          ,,,,

          ,,,;

     cout<<"Block in the middle"<<endl;

     cout<<m.block<,>(,)<<endl<<endl;

     for(int i = ; i <= ; ++i)

     {

         cout<<"Block of size "<<i<<" x "<<i<<endl;

         cout<<m.block(,,i,i)<<endl<<endl;

     }

     return ;

 }

　　block也可以被用作左值，即block可以进行赋值操作。

　　举例如下：

 #include <iostream>

 #include <eigen3/Eigen/Dense>

 using namespace Eigen;

 using namespace std;

 int main(int argc ,char** argv)

 {

     Array22f m;

     m << ,,

          ,;

     Array44f a = Array44f::Constant(0.6);

     cout<<"Here is the array a"<<endl<<a<<endl<<endl;

     a.block<,>(,) = m;

     cout<<"Here is now a with m copoed into its central 2x2 block"<<endl<<a<<endl<<endl;

     a.block(,,,) = a.block(,,,);

     cout<<"Here is now a with bottom-right 2x3 copied into top-left 2x3 block:"<<endl<<a<<endl<<endl;

     return ;

 }

　　运行结果如下：

　　特殊情况下的块操作，比如取整行或者整列，取上面的若干行或者底部的若干行。

　　取整行和整列的操作如下：

　　　　matrix.row(i);

matrix.col(j);

　　访问矩阵的行和列的操作如下：

 #include <iostream>

 #include <eigen3/Eigen/Dense>

 using namespace Eigen;

 using namespace std;

 int main(int argc ,char** argv)

 {

     MatrixXf m(,);

     m << ,,,

          ,,,

          ,,;

     cout<<"Here is the matrix m:"<<endl<<m<<endl;

     cout<<"2nd Row: "<<m.row()<<endl;

     m.col() += *m.col();

     cout<<"After adding 4 times the first column into the third column,the matrix m is:\n";

     cout<<m<<endl;

     return ;

 }

Eigen 还提供了以下对角点等特殊块操作方法：

　　矩阵的块操作：

　　举例如下;

 #include <iostream>

 #include <eigen3/Eigen/Dense>

 using namespace Eigen;

 using namespace std;

 int main(int argc ,char** argv)

 {

     ArrayXf v();

     v << ,,,,,;

     cout<<"v.head(3) = "<<endl<<v.head()<<endl<<endl;

     cout<<"v.tail<>() = "<<endl<<v.tail<>()<<endl<<endl;

     v.segemnt(,) *= ;

     cout<<"after v.segemt(1,4) *= 2，v="<<endl<<v<<endl;

     return ;

 }

　　运行结果如下：

巴特西

Eigen 学习之块操作

最新文章

热门文章