|
|
◆ WZtrtrs2()
| Function WZtrtrs2 |
( |
Uplo As |
String, |
|
|
N As |
Long, |
|
|
A As |
Variant, |
|
|
B As |
Variant, |
|
|
Optional Trans As |
String = "N", |
|
|
Optional Nrhs As |
Long = 1 |
|
) |
| |
連立一次方程式 AX = B, ATX = B または AHX = B の解 (複素三角行列)
- 目的
- WZtrtrs2は三角行列の連立一次方程式
A * X = B, A^T * X = B または A^H * X = B
複素数を表現するために実数部と虚数部を隣り合ったセルに格納する(左が実数部, 右が虚数部). 得られた解も実数部と虚数部が隣り合った別々のセルに出力される.
- 戻り値
- N+1 × 2Nrhs
| 列1 | 列2 | ・・・ | 列2Nrhs |
| 行1〜N | 解行列 X (実数部と虚数部を隣り合った列に格納(左が実数部, 右が虚数部)) |
| 行N+1 | 1/条件数 | リターンコード | ・・・ | 0 |
リターンコード.
= 0: 正常終了.
= i > 0: 行列のi番目のピボットがゼロになった. (行列 A は特異)
- 引数
-
| [in] | Uplo | = "U": Aは上三角行列.
= "L": Aは下三角行列. |
| [in] | N | 連立方程式の数, すなわち, 行列Aの行および列数. (N >= 1) |
| [in] | A | (N × 2N) N×N係数行列 A. (Uplo = 0のとき下三角要素のみ使用, Uplo = 1のとき上三角要素のみ使用) |
| [in] | B | (N × 2Nrhs) N×Nrhs右辺行列 B. |
| [in] | Trans | (省略可)
連立方程式の形を指定. (省略時 = "N")
= "N": A * X = B. (転置なし)
= "T": A^T * X = B. (転置あり)
= "C": A^H * X = B. (共役転置あり) |
| [in] | Nrhs | (省略可)
右辺行列Bの列数. (Nrhs >= 1) (省略時 = 1) |
- 出典
- LAPACK
- 使用例
- 連立一次方程式 Ax = B を解く. また, Aの条件数の逆数の推定値(RCond)を求める. ただし, Aは三角行列で
( 0.20-0.11i 0 0 )
A = ( -0.93-0.32i 0.81+0.37i 0 )
( -0.80-0.92i -0.29+0.86i 0.64+0.51i )
( 0.2069+0.0399i )
B = ( -0.6633-0.6775i )
( -0.4965-0.6057i )
|
1.9.7