CALL EXPMATRIX ルーチン

行列e^tAを返します(入力行列A、乗数tの場合)。

カテゴリ:	行列操作
要件	入力行列と出力行列は正方で、同じ次元が含まれている必要があります。tは、任意のスカラー値です。

構文

CALL EXPMATRIX(X, t, Y);

必須引数

X

次元m x m (つまり、X[m, m])を含む入力行列を指定します。

t

倍精度のスカラー値を指定します。

Y

次のように、次元m x m (つまり、Y[m, m])を含む出力行列を指定します。

table with 1 row and 1 column , row1 column 1 , y equals , epsilon super t x end super , end table. 別の形式を利用するにはイメージをクリックします。

詳細

EXPMATRIX CALLでは、Padé概算アルゴリズム(Golub and van Loan (1989), p. 558に記載)を使用します。このモジュールは、行列の各エントリを累乗しないことに注意してください。詳細については、SAS/IML User's GuideのEXPMATRIXドキュメントを参照してください。

例

次の例では、EXPMATRIX CALLルーチンを使用します。

options pageno=1 nodate;

proc fcmp;
   array mat1[3,3] (0.3, -0.78, -0.82, 0.54, 1.74,
                    1.2, -1.3, 0.25, 1.49);
   array result[3,3];
   call expmatrix(mat1, 3, result);
   put result=;
quit;

options pageno=1 nodate;

proc fcmp;
   array mat1[3,3] (0.3, -0.78, -0.82, 0.54, 1.74,
                    1.2, -1.3, 0.25, 1.49);
   array result[3,3];
   call expmatrix(mat1, 3, result);
   put result=;
quit;

EXPMATRIX CALLルーチンの結果

                                The SAS System                               1

                              The FCMP Procedure

result[1, 1]=365.58043585 result[1, 2]=-589.6358476 result[1, 3]=-897.1034008
result[2, 1]=-507.0874798 result[2, 2]=838.64570481 result[2, 3]=1267.3598426
result[3, 1]=-551.588816 result[3, 2]=858.97629382 result[3, 3]=1324.8187125

最終更新: 2026年9月18日