1

次のマトリックスを構築したいと思います

              A(3x3)              B(3x3N) 
    F     = [|1 0 0|  |0 0 0 0 0 0 0 0 0 0 0 0 ... 0 0 0|;
(6x3+3N)     |0 1 0|  |0 0 0 0 0 0 0 0 0 0 0 0 ... 0 0 0|;
             |0 0 1|  |0 0 0 0 0 0 0 0 0 0 0 0 ... 0 0 0|;
             -------  -----------------------------------
             |0 0 0|  |0 0 0 1 0 0 0 0 0 0 0 0 ... 0 0 0|;
             |0 0 0|  |0 0 0 0 1 0 0 0 0 0 0 0 ... 0 0 0|;
             |0 0 0|  |0 0 0 0 0 1 0 0 0 0 0 0 ... 0 0 0|];
              C(3x3)              D(3x3N)

B & C常にゼロです。A単位行列です。Dトリッキーです。それらはインデックスに基づいて指定されます。たとえば、インデックスが の場合0、D

                   |1 0 0 0 0 0 0 0 0 0 0 0 ... 0 0 0|;
                   |0 1 0 0 0 0 0 0 0 0 0 0 ... 0 0 0|;
                   |0 0 1 0 0 0 0 0 0 0 0 0 ... 0 0 0|;
                                D(3x3N)

インデックスが1_D

                   |0 0 0 1 0 0 0 0 0 0 0 0 ... 0 0 0|;
                   |0 0 0 0 1 0 0 0 0 0 0 0 ... 0 0 0|;
                   |0 0 0 0 0 1 0 0 0 0 0 0 ... 0 0 0|;
                                D(3x3N)

Eigen Library でこの手順を実行するにはどうすればよいですか? Eigen で行列を作成する方法は知っていますが、それらを 1 つの行列として構築する方法がわかりません。

4

1 に答える 1

2

Eigen での高度な初期化のページをご覧ください。次の（テストされていない）コードは、必要なマトリックスを構築すると思います：

MatrixXd F(6, 3 + 3 * N);  // you need to specify the size before doing F << ...
F << MatrixXd::Identity(3, 3),   // matrix A
     MatrixXd::Zero(3, 3 * N),   // matrix B
     MatrixXd::Zero(3, 3 + 3 * index),   // matrix C plus left zero block in D
     MatrixXd::Identity(3, 3),   // indentity block in D
     MatrixXd::Zero(3, 3 * (N - index - 1));  // right zero block in D
于 2014-08-07T09:17:27.607 に答える