Matrix defined using smaller matrices called blocks
数学 において 、 ブロック行列 または 分割行列とは、 ブロック または 部分行列 と呼ばれるセクションに分割されていると解釈される 行列 である 。 [1] [2]
直感的には、ブロック行列として解釈される行列は、元の行列を水平線と垂直線の集合で表し、それが 小さな行列の集合に分割(または 区分)されているものとして視覚化することができます。 [3] [2] たとえば、以下に示す 3x4 行列は、水平線と垂直線によって 4 つのブロックに分割されています。左上の 2x3 ブロック、右上の 2x1 ブロック、左下の 1x3 ブロック、右下の 1x1 ブロックです。
[ a 11 a 12 a 13 b 1 a 21 a 22 a 23 b 2 c 1 c 2 c 3 d ] {\displaystyle \left[{\begin{array}{ccc|c}a_{11}&a_{12}&a_{13}&b_{1}\\a_{21}&a_{22}&a_{23}&b_{2}\\\hline c_{1}&c_{2}&c_{3}&d\end{array}}\right]} 任意の行列は 1 つ以上の方法でブロック行列として解釈することができ、それぞれの解釈は行と列がどのように分割されるかによって定義されます。
この概念は、 を の 集合 に 分割し 、さらに をの集合に分割することで、 の 集合を の集合に分割することで、 の集合を の集合の「合計」としてより正確に表すことができます。元の行列は、元の行列の要素が、 の オフセット 要素と 1対1 で対応する という意味で、これらのグループの「合計」とみなされます。 ここで、 および です 。 [4] n {\displaystyle n} m {\displaystyle m} M {\displaystyle M} n {\displaystyle n} rowgroups {\displaystyle {\text{rowgroups}}} m {\displaystyle m} colgroups {\displaystyle {\text{colgroups}}} ( i , j ) {\displaystyle (i,j)} ( s , t ) {\displaystyle (s,t)} ( x , y ) {\displaystyle (x,y)} x ∈ rowgroups {\displaystyle x\in {\text{rowgroups}}} y ∈ colgroups {\displaystyle y\in {\text{colgroups}}}
ブロック行列代数は一般に行列の カテゴリ の 二重積 から生じる。 [5]
168×168要素のブロック行列。12×12、12×24、24×12、24×24のサブ行列が含まれています。非ゼロ要素は青色、ゼロ要素は灰色で表示されます。
例 マトリックス
P = [ 1 2 2 7 1 5 6 2 3 3 4 5 3 3 6 7 ] {\displaystyle \mathbf {P} ={\begin{bmatrix}1&2&2&7\\1&5&6&2\\3&3&4&5\\3&3&6&7\end{bmatrix}}} 4つのブロックに分割して視覚化することができます。
P = [ 1 2 2 7 1 5 6 2 3 3 4 5 3 3 6 7 ] {\displaystyle \mathbf {P} =\left[{\begin{array}{cc|cc}1&2&2&7\\1&5&6&2\\\hline 3&3&4&5\\3&3&6&7\end{array}}\right]} 。 水平線と垂直線は特別な数学的意味を持たないが [6] [7] 、分割を視覚化するための一般的な方法である。 [6] [7] この分割によって、 は4つの2×2ブロックに分割され、 P {\displaystyle P}
P 11 = [ 1 2 1 5 ] , P 12 = [ 2 7 6 2 ] , P 21 = [ 3 3 3 3 ] , P 22 = [ 4 5 6 7 ] . {\displaystyle \mathbf {P} _{11}={\begin{bmatrix}1&2\\1&5\end{bmatrix}},\quad \mathbf {P} _{12}={\begin{bmatrix}2&7\\6&2\end{bmatrix}},\quad \mathbf {P} _{21}={\begin{bmatrix}3&3\\3&3\end{bmatrix}},\quad \mathbf {P} _{22}={\begin{bmatrix}4&5\\6&7\end{bmatrix}}.} 分割された行列は次のように書ける。
P = [ P 11 P 12 P 21 P 22 ] . {\displaystyle \mathbf {P} ={\begin{bmatrix}\mathbf {P} _{11}&\mathbf {P} _{12}\\\mathbf {P} _{21}&\mathbf {P} _{22}\end{bmatrix}}.} [8]
とする 。 の 分割は の 表現であり、 A ∈ C m × n {\displaystyle A\in \mathbb {C} ^{m\times n}} A {\displaystyle A} A {\displaystyle A}
A = [ A 11 A 12 ⋯ A 1 q A 21 A 22 ⋯ A 2 q ⋮ ⋮ ⋱ ⋮ A p 1 A p 2 ⋯ A p q ] {\displaystyle A={\begin{bmatrix}A_{11}&A_{12}&\cdots &A_{1q}\\A_{21}&A_{22}&\cdots &A_{2q}\\\vdots &\vdots &\ddots &\vdots \\A_{p1}&A_{p2}&\cdots &A_{pq}\end{bmatrix}}} 、 ここで、 は連続する部分行列、 、 である 。 [9] 分割の 要素は ブロック と呼ばれる。 [9] A i j ∈ C m i × n j {\displaystyle A_{ij}\in \mathbb {C} ^{m_{i}\times n_{j}}} ∑ i = 1 p m i = m {\displaystyle \sum _{i=1}^{p}m_{i}=m} ∑ j = 1 q n j = n {\displaystyle \sum _{j=1}^{q}n_{j}=n} A i j {\displaystyle A_{ij}}
この定義によれば、どの列のブロックも列数が同じでなければなりません。 [9] 同様に、どの行のブロックも行数が同じでなければなりません。 [9]
パーティション分割方法 行列は様々な方法で分割できる。 [9] 例えば、行列が 列で分割されて いると言われるのは、 次のように書かれる場合である。 A {\displaystyle A}
A = ( a 1 a 2 ⋯ a n ) {\displaystyle A=(a_{1}\ a_{2}\ \cdots \ a_{n})} 、 ここでは の 番目の列 です 。 [9] 行列は 行ごとに分割する こともできます。 a j {\displaystyle a_{j}} j {\displaystyle j} A {\displaystyle A}
A = [ a 1 T a 2 T ⋮ a m T ] {\displaystyle A={\begin{bmatrix}a_{1}^{T}\\a_{2}^{T}\\\vdots \\a_{m}^{T}\end{bmatrix}}} 、 の 番目の行は どこ ですか 。 [9] a i T {\displaystyle a_{i}^{T}} i {\displaystyle i} A {\displaystyle A}
共通パーティション 多くの場合、 [9] 2x2のパーティションに遭遇します
A = [ A 11 A 12 A 21 A 22 ] {\displaystyle A={\begin{bmatrix}A_{11}&A_{12}\\A_{21}&A_{22}\end{bmatrix}}} , [9] 特に、がスカラーである 形式では、次のようになります。 A 11 {\displaystyle A_{11}}
A = [ a 11 a 12 T a 21 A 22 ] {\displaystyle A={\begin{bmatrix}a_{11}&a_{12}^{T}\\a_{21}&A_{22}\end{bmatrix}}} . [9]
ブロック行列演算
転置 させて
A = [ A 11 A 12 ⋯ A 1 q A 21 A 22 ⋯ A 2 q ⋮ ⋮ ⋱ ⋮ A p 1 A p 2 ⋯ A p q ] {\displaystyle A={\begin{bmatrix}A_{11}&A_{12}&\cdots &A_{1q}\\A_{21}&A_{22}&\cdots &A_{2q}\\\vdots &\vdots &\ddots &\vdots \\A_{p1}&A_{p2}&\cdots &A_{pq}\end{bmatrix}}} ここで である 。(この行列は § 加算と§ 乗算で再利用される。)そしてその転置は A i j ∈ C k i × ℓ j {\displaystyle A_{ij}\in \mathbb {C} ^{k_{i}\times \ell _{j}}} A {\displaystyle A}
A T = [ A 11 T A 21 T ⋯ A p 1 T A 12 T A 22 T ⋯ A p 2 T ⋮ ⋮ ⋱ ⋮ A 1 q T A 2 q T ⋯ A p q T ] {\displaystyle A^{T}={\begin{bmatrix}A_{11}^{T}&A_{21}^{T}&\cdots &A_{p1}^{T}\\A_{12}^{T}&A_{22}^{T}&\cdots &A_{p2}^{T}\\\vdots &\vdots &\ddots &\vdots \\A_{1q}^{T}&A_{2q}^{T}&\cdots &A_{pq}^{T}\end{bmatrix}}} , [9] [10] 転置を共役転置に置き換えても同じ式が成り立つ。 [9]
ブロック転置 ブロック行列に対しても、転置 の特殊な形式 を定義することができます。この形式 では、 個々のブロックは順序が変更されますが、転置は行われません。 が ブロック を持つ ブロック行列であるとすると、 のブロック転置はを 持つ ブロック行列 です 。 [11] 従来のトレース演算子と同様に、ブロック転置は となる 線形写像 です。 [10] しかし、一般に、 と のブロックが可換でない限り、この性質は成り立ちません 。 A = ( B i j ) {\displaystyle A=(B_{ij})} k × l {\displaystyle k\times l} m × n {\displaystyle m\times n} B i j {\displaystyle B_{ij}} A {\displaystyle A} l × k {\displaystyle l\times k} A B {\displaystyle A^{\mathcal {B}}} m × n {\displaystyle m\times n} ( A B ) i j = B j i {\displaystyle \left(A^{\mathcal {B}}\right)_{ij}=B_{ji}} ( A + C ) B = A B + C B {\displaystyle (A+C)^{\mathcal {B}}=A^{\mathcal {B}}+C^{\mathcal {B}}} ( A C ) B = C B A B {\displaystyle (AC)^{\mathcal {B}}=C^{\mathcal {B}}A^{\mathcal {B}}} A {\displaystyle A} C {\displaystyle C}
追加 させて
B = [ B 11 B 12 ⋯ B 1 s B 21 B 22 ⋯ B 2 s ⋮ ⋮ ⋱ ⋮ B r 1 B r 2 ⋯ B r s ] {\displaystyle B={\begin{bmatrix}B_{11}&B_{12}&\cdots &B_{1s}\\B_{21}&B_{22}&\cdots &B_{2s}\\\vdots &\vdots &\ddots &\vdots \\B_{r1}&B_{r2}&\cdots &B_{rs}\end{bmatrix}}} 、 ここで 、 、 を §転置で定義された行列とする。(この行列は§乗算で再利用される。) 、 、 、 の 場合、 B i j ∈ C m i × n j {\displaystyle B_{ij}\in \mathbb {C} ^{m_{i}\times n_{j}}} A {\displaystyle A} B {\displaystyle B} p = r {\displaystyle p=r} q = s {\displaystyle q=s} k i = m i {\displaystyle k_{i}=m_{i}} ℓ j = n j {\displaystyle \ell _{j}=n_{j}}
A + B = [ A 11 + B 11 A 12 + B 12 ⋯ A 1 q + B 1 q A 21 + B 21 A 22 + B 22 ⋯ A 2 q + B 2 q ⋮ ⋮ ⋱ ⋮ A p 1 + B p 1 A p 2 + B p 2 ⋯ A p q + B p q ] {\displaystyle A+B={\begin{bmatrix}A_{11}+B_{11}&A_{12}+B_{12}&\cdots &A_{1q}+B_{1q}\\A_{21}+B_{21}&A_{22}+B_{22}&\cdots &A_{2q}+B_{2q}\\\vdots &\vdots &\ddots &\vdots \\A_{p1}+B_{p1}&A_{p2}+B_{p2}&\cdots &A_{pq}+B_{pq}\end{bmatrix}}} . [9]
乗算 因子の部分行列の代数のみを扱うブロック分割行列積を用いることも可能です。しかしながら、因子の分割は任意ではなく、 2つの行列間の「 適合 分割」 [12] と、 使用されるすべての部分行列積が定義されるような分割が必要です。 [13] A {\displaystyle A} B {\displaystyle B}
2 つの行列 および が、積 について等角 的に分割されていると言われるのは 、 乗算が、 部分行列をスカラーであるかのように扱いながら順序を保って実行され、関係する部分行列の積と和がすべて定義されている場合です。 A {\displaystyle A} B {\displaystyle B} A B {\displaystyle AB} A {\displaystyle A} B {\displaystyle B} A B {\displaystyle AB}
— アラク・M・マタイとハンス・J・ハウボルド『 線形代数:物理学者とエンジニアのためのコース』 [14]
§転置で定義された行列をと し、 §加法で定義された行列を とする。すると、行列積は A {\displaystyle A} B {\displaystyle B}
C = A B {\displaystyle C=AB} ブロック単位で 実行することで行列が得られます 。結果の行列の各行列は、 以下の乗算によって計算されます。 C {\displaystyle C} ( p × s ) {\displaystyle (p\times s)} C {\displaystyle C}
C i j = ∑ k = 1 q A i k B k j . {\displaystyle C_{ij}=\sum _{k=1}^{q}A_{ik}B_{kj}.} [6] または、 繰り返されるインデックスを暗黙的に合計する アインシュタイン表記法を使用します。
C i j = A i k B k j . {\displaystyle C_{ij}=A_{ik}B_{kj}.} 行列として 描くと、 C {\displaystyle C}
C = A B = [ ∑ i = 1 q A 1 i B i 1 ∑ i = 1 q A 1 i B i 2 ⋯ ∑ i = 1 q A 1 i B i s ∑ i = 1 q A 2 i B i 1 ∑ i = 1 q A 2 i B i 2 ⋯ ∑ i = 1 q A 2 i B i s ⋮ ⋮ ⋱ ⋮ ∑ i = 1 q A p i B i 1 ∑ i = 1 q A p i B i 2 ⋯ ∑ i = 1 q A p i B i s ] {\displaystyle C=AB={\begin{bmatrix}\sum _{i=1}^{q}A_{1i}B_{i1}&\sum _{i=1}^{q}A_{1i}B_{i2}&\cdots &\sum _{i=1}^{q}A_{1i}B_{is}\\\sum _{i=1}^{q}A_{2i}B_{i1}&\sum _{i=1}^{q}A_{2i}B_{i2}&\cdots &\sum _{i=1}^{q}A_{2i}B_{is}\\\vdots &\vdots &\ddots &\vdots \\\sum _{i=1}^{q}A_{pi}B_{i1}&\sum _{i=1}^{q}A_{pi}B_{i2}&\cdots &\sum _{i=1}^{q}A_{pi}B_{is}\end{bmatrix}}} . [9]
反転 行列が 4 つのブロックに分割されている場合、次のように ブロックごとに反転 できます。
P = [ A B C D ] − 1 = [ A − 1 + A − 1 B ( D − C A − 1 B ) − 1 C A − 1 − A − 1 B ( D − C A − 1 B ) − 1 − ( D − C A − 1 B ) − 1 C A − 1 ( D − C A − 1 B ) − 1 ] , {\displaystyle {P}={\begin{bmatrix}{A}&{B}\\{C}&{D}\end{bmatrix}}^{-1}={\begin{bmatrix}{A}^{-1}+{A}^{-1}{B}\left({D}-{CA}^{-1}{B}\right)^{-1}{CA}^{-1}&-{A}^{-1}{B}\left({D}-{CA}^{-1}{B}\right)^{-1}\\-\left({D}-{CA}^{-1}{B}\right)^{-1}{CA}^{-1}&\left({D}-{CA}^{-1}{B}\right)^{-1}\end{bmatrix}},} ここで 、A と Dは 任意サイズの正方ブロックであり、 B と C はそれらと分割 可能 である。さらに、 Aと、 P における A のシュアー補集合 : P / A = D − CA −1 B は逆行列でなければならない。 [15]
同様に、ブロックを並べ替えると次のようになります。
P = [ A B C D ] − 1 = [ ( A − B D − 1 C ) − 1 − ( A − B D − 1 C ) − 1 B D − 1 − D − 1 C ( A − B D − 1 C ) − 1 D − 1 + D − 1 C ( A − B D − 1 C ) − 1 B D − 1 ] . {\displaystyle {P}={\begin{bmatrix}{A}&{B}\\{C}&{D}\end{bmatrix}}^{-1}={\begin{bmatrix}\left({A}-{BD}^{-1}{C}\right)^{-1}&-\left({A}-{BD}^{-1}{C}\right)^{-1}{BD}^{-1}\\-{D}^{-1}{C}\left({A}-{BD}^{-1}{C}\right)^{-1}&\quad {D}^{-1}+{D}^{-1}{C}\left({A}-{BD}^{-1}{C}\right)^{-1}{BD}^{-1}\end{bmatrix}}.} [16] ここで、 Dと P における D のシュアー補数 : P / D = A − BD −1 C は逆行列でなければならない。
A と D が 両方とも逆である場合 、次のようになります。
[ A B C D ] − 1 = [ ( A − B D − 1 C ) − 1 0 0 ( D − C A − 1 B ) − 1 ] [ I − B D − 1 − C A − 1 I ] . {\displaystyle {\begin{bmatrix}{A}&{B}\\{C}&{D}\end{bmatrix}}^{-1}={\begin{bmatrix}\left({A}-{B}{D}^{-1}{C}\right)^{-1}&{0}\\{0}&\left({D}-{C}{A}^{-1}{B}\right)^{-1}\end{bmatrix}}{\begin{bmatrix}{I}&-{B}{D}^{-1}\\-{C}{A}^{-1}&{I}\end{bmatrix}}.} ワインスタイン・アロンザイン恒等式 により 、ブロック対角行列の 2 つの行列のうち 1 つが逆行列である場合、もう 1 つも逆行列になります。
完全逆行列から部分逆行列を計算する ブロック反転公式における行列とその逆行列間の対称性により、行列 P とその逆行列 P −1 が共形に分割される場合、
P = [ A B C D ] , P − 1 = [ E F G H ] {\displaystyle P={\begin{bmatrix}{A}&{B}\\{C}&{D}\end{bmatrix}},\quad P^{-1}={\begin{bmatrix}{E}&{F}\\{G}&{H}\end{bmatrix}}} すると、任意の主部分行列の逆行列はP −1 の対応するブロックから計算できる 。
A − 1 = E − F H − 1 G {\displaystyle {A}^{-1}={E}-{FH}^{-1}{G}} D − 1 = H − G E − 1 F {\displaystyle {D}^{-1}={H}-{GE}^{-1}{F}} この関係は、 E −1 = A − BD −1 C (シュアー補集合)であることを認識し、 P と P −1 の役割を逆にして同じブロック反転公式を適用することから導かれる 。 [17] [18]
行列式 上の行列の行列式の公式は、適切な更なる仮定の下で、4つの部分行列 と平方行列 からなる行列に対しても成立する。 ライプニッツの公式 か、 シュール補集合 を含む因数分解 を用いて証明できる最も簡単な公式は、以下の 通りである
。 2 × 2 {\displaystyle 2\times 2} A , B , C , D {\displaystyle A,B,C,D} A {\displaystyle A} D {\displaystyle D}
det [ A 0 C D ] = det ( A ) det ( D ) = det [ A B 0 D ] . {\displaystyle \det {\begin{bmatrix}A&0\\C&D\end{bmatrix}}=\det(A)\det(D)=\det {\begin{bmatrix}A&B\\0&D\end{bmatrix}}.} [16] この式を用いると、と の 特性多項式は と の特性多項式の積と同じである ことが分かります 。さらに、 または が 対角化 可能であれば、 と も 対角化可能です。逆は偽です。 を確認してください 。 [ A 0 C D ] {\displaystyle {\begin{bmatrix}A&0\\C&D\end{bmatrix}}} [ A B 0 D ] {\displaystyle {\begin{bmatrix}A&B\\0&D\end{bmatrix}}} A {\displaystyle A} D {\displaystyle D} [ A 0 C D ] {\displaystyle {\begin{bmatrix}A&0\\C&D\end{bmatrix}}} [ A B 0 D ] {\displaystyle {\begin{bmatrix}A&B\\0&D\end{bmatrix}}} A {\displaystyle A} D {\displaystyle D} [ 1 1 0 1 ] {\displaystyle {\begin{bmatrix}1&1\\0&1\end{bmatrix}}}
が 逆行列を 持つ場合 、 A {\displaystyle A}
det [ A B C D ] = det ( A ) det ( D − C A − 1 B ) , {\displaystyle \det {\begin{bmatrix}A&B\\C&D\end{bmatrix}}=\det(A)\det \left(D-CA^{-1}B\right),} [16] そして 逆行列が成り立つ場合、 D {\displaystyle D}
det [ A B C D ] = det ( D ) det ( A − B D − 1 C ) . {\displaystyle \det {\begin{bmatrix}A&B\\C&D\end{bmatrix}}=\det(D)\det \left(A-BD^{-1}C\right).} [19] [16] ブロックが 同じ 大きさの正方行列である場合、さらに次の式が成り立ちます。例えば、 と が 可換 (つまり )である場合、 C {\displaystyle C} D {\displaystyle D} C D = D C {\displaystyle CD=DC}
det [ A B C D ] = det ( A D − B C ) . {\displaystyle \det {\begin{bmatrix}A&B\\C&D\end{bmatrix}}=\det(AD-BC).} [20] 、、 または の場合も同様のことが言えます 。つまり、 の
場合、 A B = B A {\displaystyle AB=BA} A C = C A {\displaystyle AC=CA} B D = D B {\displaystyle BD=DB} A C = C A {\displaystyle AC=CA}
det [ A B C D ] = det ( A D − C B ) . {\displaystyle \det {\begin{bmatrix}A&B\\C&D\end{bmatrix}}=\det(AD-CB).} と の順序が変わっていることに注目してください ( の代わりに )。同様に、 の場合、 は に置き換えてください (つまり )。また 、 の場合、 となります 。最後の2つの結果では、基底環の可換性を使用する必要がありますが、最初の2つの結果では使用する必要はありません。 C {\displaystyle C} B {\displaystyle B} C B {\displaystyle CB} B C {\displaystyle BC} B D = D B {\displaystyle BD=DB} A D {\displaystyle AD} D A {\displaystyle DA} det ( D A − B C ) {\displaystyle \det(DA-BC)} A B = B A {\displaystyle AB=BA} det ( D A − C B ) {\displaystyle \det(DA-CB)}
この式は、個々のブロック間の適切な交換条件の下で、 個 以上のブロックから構成される行列に一般化されている。 [21] 2 × 2 {\displaystyle 2\times 2}
と については 、次の式が成り立ちます( と が 可換でなくても)。 A = D {\displaystyle A=D} B = C {\displaystyle B=C} A {\displaystyle A} B {\displaystyle B}
det [ A B B A ] = det ( A − B ) det ( A + B ) . {\displaystyle \det {\begin{bmatrix}A&B\\B&A\end{bmatrix}}=\det(A-B)\det(A+B).} [16]
特殊なタイプのブロック行列
直和とブロック対角行列
直和 任意の行列 A (サイズ m × n )と B (サイズ p × q )に対して、 A と B の 直和は A B と表され 、次のように定義されます。 ⊕ {\displaystyle \oplus }
A ⊕ B = [ a 11 ⋯ a 1 n 0 ⋯ 0 ⋮ ⋱ ⋮ ⋮ ⋱ ⋮ a m 1 ⋯ a m n 0 ⋯ 0 0 ⋯ 0 b 11 ⋯ b 1 q ⋮ ⋱ ⋮ ⋮ ⋱ ⋮ 0 ⋯ 0 b p 1 ⋯ b p q ] . {\displaystyle {A}\oplus {B}={\begin{bmatrix}a_{11}&\cdots &a_{1n}&0&\cdots &0\\\vdots &\ddots &\vdots &\vdots &\ddots &\vdots \\a_{m1}&\cdots &a_{mn}&0&\cdots &0\\0&\cdots &0&b_{11}&\cdots &b_{1q}\\\vdots &\ddots &\vdots &\vdots &\ddots &\vdots \\0&\cdots &0&b_{p1}&\cdots &b_{pq}\end{bmatrix}}.} [10] 例えば、
[ 1 3 2 2 3 1 ] ⊕ [ 1 6 0 1 ] = [ 1 3 2 0 0 2 3 1 0 0 0 0 0 1 6 0 0 0 0 1 ] . {\displaystyle {\begin{bmatrix}1&3&2\\2&3&1\end{bmatrix}}\oplus {\begin{bmatrix}1&6\\0&1\end{bmatrix}}={\begin{bmatrix}1&3&2&0&0\\2&3&1&0&0\\0&0&0&1&6\\0&0&0&0&1\end{bmatrix}}.} この操作は、任意の次元の配列に自然に一般化されます ( A と B の 次元数が同じである 場合)。
2 つの行列の ベクトル空間の 直和 内の任意の要素は、 2 つの行列の直和として表すことができることに注意してください。
ブロック対角行列 ブロック 対角行列 とは、主対角 ブロックが 正方行列で、すべての非対角ブロックが零行列であるような正方行列であるブロック行列である。 [16] つまり、ブロック対角行列 A は次の形式を持つ。
A = [ A 1 0 ⋯ 0 0 A 2 ⋯ 0 ⋮ ⋮ ⋱ ⋮ 0 0 ⋯ A n ] {\displaystyle {A}={\begin{bmatrix}{A}_{1}&{0}&\cdots &{0}\\{0}&{A}_{2}&\cdots &{0}\\\vdots &\vdots &\ddots &\vdots \\{0}&{0}&\cdots &{A}_{n}\end{bmatrix}}} ここで、 A k はk = 1, ..., n に対して正方行列である 。言い換えれば、行列 Aは A 1 , ..., A n の 直和 である 。 [16]これは A 1 ⊕ A 2 ⊕ ... ⊕ A n [10] または diag( A 1 , A 2 , ..., A n ) [10] と表記することもできる (後者は 対角行列 に使用されるのと同じ形式である)。任意の正方行列は、1つのブロックのみを持つブロック対角行列と見なすことができる。
行列式 と トレース については 、次の特性が成り立ちます。
det A = det A 1 × ⋯ × det A n , {\displaystyle {\begin{aligned}\det {A}&=\det {A}_{1}\times \cdots \times \det {A}_{n},\end{aligned}}} [22] [23] と tr A = tr A 1 + ⋯ + tr A n . {\displaystyle {\begin{aligned}\operatorname {tr} {A}&=\operatorname {tr} {A}_{1}+\cdots +\operatorname {tr} {A}_{n}.\end{aligned}}} [16] [23] ブロック対角行列が逆行列となるのは、その主対角ブロックのそれぞれが逆行列となる 場合のみであり 、この場合、その逆行列は次式で表される別のブロック対角行列となる。
[ A 1 0 ⋯ 0 0 A 2 ⋯ 0 ⋮ ⋮ ⋱ ⋮ 0 0 ⋯ A n ] − 1 = [ A 1 − 1 0 ⋯ 0 0 A 2 − 1 ⋯ 0 ⋮ ⋮ ⋱ ⋮ 0 0 ⋯ A n − 1 ] . {\displaystyle {\begin{bmatrix}{A}_{1}&{0}&\cdots &{0}\\{0}&{A}_{2}&\cdots &{0}\\\vdots &\vdots &\ddots &\vdots \\{0}&{0}&\cdots &{A}_{n}\end{bmatrix}}^{-1}={\begin{bmatrix}{A}_{1}^{-1}&{0}&\cdots &{0}\\{0}&{A}_{2}^{-1}&\cdots &{0}\\\vdots &\vdots &\ddots &\vdots \\{0}&{0}&\cdots &{A}_{n}^{-1}\end{bmatrix}}.} [24] の固有値 [25] と固有ベクトル は 、単にsの固有値と固有ベクトルを組み合わせたものである 。 [ 23] A {\displaystyle {A}} A k {\displaystyle {A}_{k}}
ブロック三角行列 ブロック 三重対角行列 は、もう一つの特殊なブロック行列です。ブロック対角行列と同様に、 正方 行列であり、下対角、 主対角 、上対角に正方行列(ブロック)を持ち、その他のブロックは零行列です。ブロック三重対角行列は本質的に 三重対角行列 ですが、スカラーの代わりに部分行列を持ちます。ブロック三重対角行列は、 以下の形式を持ちます
。 A {\displaystyle A}
A = [ B 1 C 1 ⋯ 0 A 2 B 2 C 2 ⋱ ⋱ ⋱ ⋮ A k B k C k ⋮ ⋱ ⋱ ⋱ A n − 1 B n − 1 C n − 1 0 ⋯ A n B n ] {\displaystyle {A}={\begin{bmatrix}{B}_{1}&{C}_{1}&&&\cdots &&{0}\\{A}_{2}&{B}_{2}&{C}_{2}&&&&\\&\ddots &\ddots &\ddots &&&\vdots \\&&{A}_{k}&{B}_{k}&{C}_{k}&&\\\vdots &&&\ddots &\ddots &\ddots &\\&&&&{A}_{n-1}&{B}_{n-1}&{C}_{n-1}\\{0}&&\cdots &&&{A}_{n}&{B}_{n}\end{bmatrix}}} ここで 、、 は それぞれ下対角、主対角、上対角の正方部分行列である。 [26] [27] A k {\displaystyle {A}_{k}} B k {\displaystyle {B}_{k}} C k {\displaystyle {C}_{k}}
ブロック三角行列は、工学問題(例えば、数値流体力学 ) の数値解法においてしばしば用いられます。LU 分解のための最適化された数値解法 [28] が利用可能であり 、ブロック三角行列を係数行列とする方程式系に対する効率的な解法アルゴリズムが存在します。 三角行列 を含む方程式系を効率的に解くために使用される トーマスアルゴリズムは、行列演算を用いてブロック三角行列に適用することもできます( ブロックLU分解 も参照 )。
ブロック三角行列 行列 が 上ブロック三角行列 (または ブロック上三角行列 [29] ) であるとは、 および と なる正の整数が存在 する場合である 。ここで、 行列はすべての に対して となる。 [25] [29] 同様に、が 下ブロック三角行列 であるとは、 となる正の整数が 存在する 場合で
ある。 ここで、行列 はすべての に対して となる 。 [25] n × n {\displaystyle n\times n} A {\displaystyle A} n 1 , … , n k {\displaystyle n_{1},\ldots ,n_{k}} n = n 1 + n 2 + … + n k {\displaystyle n=n_{1}+n_{2}+\ldots +n_{k}} A = [ A 11 A 12 ⋯ A 1 k 0 A 22 ⋯ A 2 k ⋮ ⋮ ⋱ ⋮ 0 0 ⋯ A k k ] , {\displaystyle A={\begin{bmatrix}A_{11}&A_{12}&\cdots &A_{1k}\\0&A_{22}&\cdots &A_{2k}\\\vdots &\vdots &\ddots &\vdots \\0&0&\cdots &A_{kk}\end{bmatrix}},} A i j {\displaystyle A_{ij}} n i × n j {\displaystyle n_{i}\times n_{j}} i , j = 1 , … , k {\displaystyle i,j=1,\ldots ,k} A {\displaystyle A} A = [ A 11 0 ⋯ 0 A 21 A 22 ⋯ 0 ⋮ ⋮ ⋱ ⋮ A k 1 A k 2 ⋯ A k k ] , {\displaystyle A={\begin{bmatrix}A_{11}&0&\cdots &0\\A_{21}&A_{22}&\cdots &0\\\vdots &\vdots &\ddots &\vdots \\A_{k1}&A_{k2}&\cdots &A_{kk}\end{bmatrix}},} A i j {\displaystyle A_{ij}} n i × n j {\displaystyle n_{i}\times n_{j}} i , j = 1 , … , k {\displaystyle i,j=1,\ldots ,k}
ブロックテプリッツ行列 ブロック テプリッツ行列は、別の特殊なブロック行列であり、 テプリッツ行列 には対角線に沿って繰り返される 要素があるため、行列の対角線に沿って繰り返されるブロックが含まれます。
行列が ブロックテプリッツ行列 である のは、 すべての に対して成り立つとき である。つまり、 A {\displaystyle A} A ( i , j ) = A ( k , l ) {\displaystyle A_{(i,j)}=A_{(k,l)}} k − i = l − j {\displaystyle k-i=l-j}
A = [ A 1 A 2 A 3 ⋯ A 4 A 1 A 2 ⋯ A 5 A 4 A 1 ⋯ ⋮ ⋮ ⋮ ⋱ ] {\displaystyle A={\begin{bmatrix}A_{1}&A_{2}&A_{3}&\cdots \\A_{4}&A_{1}&A_{2}&\cdots \\A_{5}&A_{4}&A_{1}&\cdots \\\vdots &\vdots &\vdots &\ddots \end{bmatrix}}} 、 ここで [ 25] A i ∈ F n i × m i {\displaystyle A_{i}\in \mathbb {F} ^{n_{i}\times m_{i}}}
ブロックハンケル行列 行列が ブロックハンケル行列 である のは、 すべての に対して成り立つとき である。つまり、 A {\displaystyle A} A ( i , j ) = A ( k , l ) {\displaystyle A_{(i,j)}=A_{(k,l)}} i + j = k + l {\displaystyle i+j=k+l}
A = [ A 1 A 2 A 3 ⋯ A 2 A 3 A 4 ⋯ A 3 A 4 A 5 ⋯ ⋮ ⋮ ⋮ ⋱ ] {\displaystyle A={\begin{bmatrix}A_{1}&A_{2}&A_{3}&\cdots \\A_{2}&A_{3}&A_{4}&\cdots \\A_{3}&A_{4}&A_{5}&\cdots \\\vdots &\vdots &\vdots &\ddots \end{bmatrix}}} 、 ここで [ 25] A i ∈ F n i × m i {\displaystyle A_{i}\in \mathbb {F} ^{n_{i}\times m_{i}}}
参照
注記 ^ イヴス、ハワード (1980). 『初等行列理論』 (復刻版). ニューヨーク: ドーバー. p. 37. ISBN 0-486-63946-0 2013年 4月24日 閲覧 。 行列を要素の長方形ブロックに分割すると便利な場合があることがわかります。そこで、いわゆる 分割 行列 、あるいは ブロック 行列 について考えてみましょう。 ^ ab Dobrushkin, Vladimir. 「Partition Matrices」. Mathematicaによる線形代数. 2024年3月24日 閲覧 。 ^ アントン・ハワード (1994). 初等線形代数 (第7版). ニューヨーク: ジョン・ワイリー. p. 30. ISBN 0-471-58742-7 選択した行と列の間に水平線と垂直線を挿入することで、 行列をより小さな行列に 分割する ことができます。 ^ Indhumathi, D.; Sarala, S. (2014-05-16). 「F値を用いたアダプティブランダムテストおよびパーティションブロックベースのアダプティブランダムテストにおけるフラグメント分析とテストケース生成」 (PDF) . International Journal of Computer Applications . 93 (6): 13. Bibcode :2014IJCA...93f..11I. doi :10.5120/16218-5662. ^ Macedo, HD; Oliveira, JN (2013). 「線形代数の型付け:二積指向アプローチ」. コンピュータプログラミング科学 . 78 (11): 2160– 2191. arXiv : 1312.4818 . doi :10.1016/j.scico.2012.07.012. ^ abc Johnston, Nathaniel (2021). 線形代数と行列代数入門 . シャム、スイス: Springer Nature. pp. 30, 425. ISBN 978-3-030-52811-9 。 ^ ab Johnston, Nathaniel (2021). 上級線形代数と行列代数 . シャム、スイス: Springer Nature. p. 298. ISBN 978-3-030-52814-0 。 ^ ジェフリー・アラン(2010年)『エンジニアと科学者のための行列演算:線形代数の必須ガイド』ドルドレヒト(オランダ);ニューヨーク:シュプリンガー、p.54、 ISBN 978-90-481-9273-1 . OCLC 639165077。 ^ abcdefghijklmn スチュワート、ギルバート・W. (1998). 行列アルゴリズム. 1: 基本分解 . フィラデルフィア、ペンシルバニア州: 産業応用数学協会. pp. 18– 20. ISBN 978-0-89871-414-2 。 ^ abcde Gentle, James E. (2007). 行列代数:統計における理論、計算、応用 . Springer Texts in Statistics. ニューヨーク: Springer New York Springer e-books. pp. 47, 487. ISBN 978-0-387-70873-7 。 ^ Mackey, D. Steven (2006). 行列多項式の構造化線形化 (PDF) (論文). マンチェスター大学. ISSN 1749-9097. OCLC 930686781. ^ イヴス、ハワード (1980). 『初等行列理論』 (復刻版). ニューヨーク: ドーバー. p. 37. ISBN 0-486-63946-0 . 2013年 4月24日 閲覧 。 定理1.9.4のような分割は、 A と B の 適合分割 と呼ばれる。 ^ アントン・ハワード (1994). 初等線形代数 (第7版). ニューヨーク: ジョン・ワイリー. p. 36. ISBN 0-471-58742-7 ... ただし、A と B のサブマトリックスのサイズが、指定された操作を実行できるサイズである必要があります。 ^ Mathai, Arakaparampil M.; Haubold, Hans J. (2017). 線形代数:物理学者とエンジニアのためのコース . De Gruyter 教科書. ベルリン・ボストン: De Gruyter. p. 162. ISBN 978-3-11-056259-0 。 ^ バーンスタイン、デニス (2005). 行列数学 . プリンストン大学出版局. p. 44. ISBN 0-691-11802-7 。 ^ abcdefgh アバディール, カリム・M.; マグナス, ヤン・R. (2005). 行列代数 . ケンブリッジ大学出版局. pp. 97, 100, 106, 111, 114, 118. ISBN 9781139443647 。 ^ 「このブロック逆行列の式は、逆行列全体に関して既知ですか?」 MathOverflow 。 ^ Escalante-B., Alberto N.; Wiskott, Laurenz (2016). 「改良グラフベースSFA:情報保存は遅さの原則を補完する」. 機械学習 . arXiv : 1412.4679 . doi :10.1007/s10994-016-5563-y. ^ Taboga, Marco (2021). 「ブロック行列の行列式」, 行列代数の講義. ^ Silvester, JR (2000). 「ブロック行列の行列式」 (PDF) . Math. Gaz . 84 (501): 460– 467. doi :10.2307/3620776. JSTOR 3620776. 2015年3月18日時点 のオリジナル (PDF)からアーカイブ。 2021年6月25日 閲覧 。 ^ Sothanaphan, Nat (2017年1月). 「非可換ブロックを持つブロック行列の行列式」. 線形代数とその応用 . 512 : 202–218 . arXiv : 1805.06027 . doi :10.1016/j.laa.2016.10.004. S2CID 119272194. ^ Quarteroni, Alfio; Sacco, Riccardo; Saleri, Fausto (2000). 数値数学 . 応用数学テキスト. ニューヨーク: Springer. pp. 10, 13. ISBN 978-0-387-98959-4 。 ^ abc ジョージ、ラジュ K.;アジャヤクマール、アビジット(2024)。 「線形代数のコース」 。 数理科学の大学テキスト : 35, 407. doi :10.1007/978-981-99-8680-4. ISBN 978-981-99-8679-8 . ISSN 2731-9318. ^ プリンス、サイモン JD (2012). 『コンピュータビジョン:モデル、学習、推論 』 ニューヨーク:ケンブリッジ大学出版局. p. 531. ISBN 978-1-107-01179-3 。 ^ abcde Bernstein, Dennis S. (2009). 『行列数学:理論、事実、公式』 (第2版). プリンストン大学出版局, ニュージャージー州, 168, 298頁. ISBN 978-0-691-14039-1 。 ^ Dietl, Guido KE (2007). クリロフ部分空間における線形推定と検出. 信号処理、通信、ネットワークの基礎. ベルリン; ニューヨーク: Springer. pp. 85, 87. ISBN 978-3-540-68478-7 . OCLC 85898525。 ^ Horn, Roger A.; Johnson, Charles R. (2017). マトリックス分析 (第2版、訂正再版). ニューヨーク:ケンブリッジ大学出版局. p. 36. ISBN 978-0-521-83940-2 。 ^ Datta, Biswa Nath (2010). 数値線形代数とその応用 (第2版). フィラデルフィア, ペンシルバニア州: SIAM. p. 168. ISBN 978-0-89871-685-6 。 ^ ab Stewart, Gilbert W. (2001). 行列アルゴリズム. 2: 固有値システム . フィラデルフィア, ペンシルバニア州: 産業応用数学協会. p. 5. ISBN 978-0-89871-503-3 。
参考文献 ストラング、ギルバート (1999). 「講義3:乗算と逆行列」. MITオープンコースウェア. 18:30–21:10.