シューア補行列

数値解析 (特に数値線形代数) や統計学、行列解析の分野での主要な道具

線型代数学関連分野におけるシューア補行列(シューアほぎょうれつ、: Schur complement; シューア補元)は区分行列に対して定義される。名称はイサイ・シューアシューアの補題の証明に用いたことに由来するが、それ以前からの使用が認められる[1]。これを Schur complement と呼び始めたのはエミリー・ヘインズワースである[2]。シューア補行列は数値解析 (特に数値線形代数) や統計学行列解析の分野では主要な道具の一つとなっている。

定義

編集

行列 A, B, C, D のサイズをそれぞれ p × p, p × q, q × p, q × q として区分行列   を考える。全体として M(p + q) × (p + q) 行列になっている。以下本項で M と書けば断りなくこの区分行列を意味するものとする。

D正則であるとき、区分行列 M の区画 D に関するシューア補行列とは   で定義される p × p 行列を言う。同様に A が正則であるとき、MA に関するシューア補行列とは   で定義される q × q 行列を言う。

AD が正則でない場合にも、逆行列の代わりに一般化逆行列を用いることにより、一般化シューア補行列を定義することはできる。

背景

編集

シューア補行列は上記の行列 M にブロック下半三角行列  Ipp × p 単位行列)を右から掛けるという形でガウス消去法を施した結果として生じる(L を掛ければ、上側の p × p 行列としてシューア補行列が現れる)。実際、行列の積は   となっている。これはLDU分解のブロック行列版となるもので、実際 M について解けば   であり、M の逆行列は D が正則かつ D に関するシューア補行列の逆行列が(存在するならば)既知のときには、D, M/D の逆行列だけから   と計算できる。

行列の反転補題英語版の項では、上記の式と A, D の役割を入れ替えて同様の導出をした式との間の関係性が詳しく述べられる。

性質

編集
  • 区分行列 M正定値対称行列ならば、シューア補行列 M/D もそうである。
  • A, B, C, D が全てスカラー(p = q = 1)のとき、2 × 2 行列の逆行列の公式   はよく知られている(ただし AD − BC は零でないものとする)。
  • 一般のサイズの場合、A が正則とすれば   が逆行列の存在する限りにおいて成立する。
  • 区分行列 M の行列式は明らかに   で計算できる。これは 2 × 2 行列の行列式の定義式を区分行列版に一般化したものと見なせる。
  • ガットマンの階数加法定理: 区分行列 M階数  で与えられる。
  • ヘインズワースの慣性加法定理英語版: 区分行列 M慣性指数A の慣性指数と M/A の慣性指数との和に等しい。

線型方程式の解法への応用

編集

x, ap-次元列ベクトルy, bq-次元列ベクトルで、区分行列 M が上記の如く与えられているとき、線型方程式系   の解法にシューア補行列は自然に表れる。D が可逆のとき、下の式に BD−1 を掛けて上の式から引けば   を得るから、シューア補行列 M/D も可逆ならば、x について解ける(さらに   から y も分かる)。これにより、もともとのサイズ (p + q) × (p + q) の係数行列の逆行列を計算する問題が、それぞれのサイズが p × pq × q のふたつの行列の逆行列を計算することに帰着される。実用上は、このアルゴリズムが数値的に良い評価を与えるようにするために、D が十分素性が良いものとなるような条件を課す。

電気工学においては、このことをしばしばノード除去 (node elimination) やクロン縮約英語版などと言う。

確率論・統計学への応用

編集

確率列ベクトル X および Y はそれぞれ Rn および Rm を動くものとし、ベクトル (X, Y) ∈ Rn+m共分散が正定値対称行列   で与えられる多変量正規分布に従うものとする。ただし、 X の共分散行列、 Y の共分散行列、 XY の間の共分散行列である。

このとき、Y が既知であるときの X条件付き共分散英語版 Cov(X | Y)C に関する Σ のシューア補行列によって   と与えられる[3](条件付き期待値は   となる)。

上記の如く Σ を(しかし確率ベクトルの共分散としてではなく)標本共分散として与えたならば、ウィッシャート分布に従う。この場合、シューア補行列 Σ/C もまたウィッシャート分布に従う[要出典]

定値性の判定条件

編集

対称行列 X  で与えられるものとする。このとき、XA および C に関するシューア補行列は   と書ける。

  1. X が正定値となるための必要十分条件は、A および X/A がともに正定値となることである:  
  2. X が正定値となるための必要十分条件は C および X/C がともに正定値となることである:  
  3. A が正定値のとき、X が半正定値となるための必要十分条件は X/A が半正定値となることである。
  4. C が正定値のとき、X が半正定値となるための必要十分条件は X/C が半正定値となることである。

1. および 3. は u を止めて v の函数とみた量   の最小化を考えることで導出できる[4]。さらに、  である(半正定値でも同様のことが言える)から、1. および 3. からそれぞれ 2. および 4. が直ちに得られる。

同じように、一般化シューア補行列を用いても X の半正定値性を判定する必要十分条件を述べることができる[1]。つまり、AgA一般化逆行列とすれば   および   が成り立つ。

関連項目

編集

参考文献

編集
  1. ^ a b Zhang, Fuzhen (2005). The Schur Complement and Its Applications. Springer. doi:10.1007/b105056. ISBN 0-387-24271-6 
  2. ^ Haynsworth, E. V., "On the Schur Complement", Basel Mathematical Notes, #BNB 20, 17 pages, June 1968.
  3. ^ von Mises, Richard (1964). “Chapter VIII.9.3”. Mathematical theory of probability and statistics. Academic Press. ISBN 978-1483255385 
  4. ^ Boyd, S. and Vandenberghe, L. (2004), "Convex Optimization", Cambridge University Press (Appendix A.5.5)

外部リンク

編集