WZgbsv WZgbsv2 WZgesv WZgesv2 WZgtsv WZgtsv2 WZsysv WZsysv2 WZtrtrs WZtrtrs2

WZgesv() Function WZgesv ( N As  Long, A As  Variant, B As  Variant, Optional Nrhs As  Long = 1  ) 連立一次方程式 AX = B の解 (複素行列) (Excel複素数形式) 目的 WZgesvは次の複素連立一次方程式を解く. A * X = B ここで, AはN×N行列, また, XおよびBはN×Nrhs行列である. まず, 行交換によるピボットの部分選択を行うLU分解を用いて, 次のようにAを分解する. A = P * L * U ここで, Pは置換行列, Lは対角要素が1の下三角行列, そして, Uは上三角行列である. 次に, 分解されたAを用いて連立方程式 A * X = B の解を求める. セル中で複素数を表現するためにExcelの複素数形式(例, 2.5+1i)を使用する. 複素数値はComplexワークシート関数を使って入力することができる. 戻り値 N+2 × Nrhs 列1 列2 ・・・ 列Nrhs 行1〜N 解行列 X 行N+1 1/条件数 0 ・・・ 0 行N+2 リターンコード 0 ・・・ 0 リターンコード. = 0: 正常終了. = i > 0: 行列のi番目のピボットがゼロになった. (行列 A は特異) 引数 [in] N 連立方程式の数, すなわち, 行列Aの行および列数. (N >= 1) [in] A (N × N) N×N係数行列 A. [in] B (N × Nrhs) N×Nrhs右辺行列 B. [in] Nrhs (省略可) 右辺行列Bの列数. (Nrhs >= 1) (省略時 = 1) 出典 LAPACK 使用例 連立一次方程式 Ax = B を解く. また, Aの条件数の逆数の推定値(RCond)を求める. ただし, ( 0.2-0.11i -0.93-0.32i 0.81+0.37i ) A = ( -0.8-0.92i -0.29+0.86i 0.64+0.51i ) ( 0.71+0.59i -0.15+0.19i 0.2+0.94i ) ( -0.5853-0.9457i ) B = ( -2.1697-1.0006i ) ( 0.0116-0.5094i ) とする.