ブロック行列

数学において、ブロック行列または分割行列とは、ブロックまたは部分行列と呼ばれるセクションに分割されていると解釈される行列である。[ 1 ] [ 2 ]

直感的には、ブロック行列として解釈される行列は、元の行列を水平線と垂直線の集合で分割(または区分)し、小さな行列の集合にしたものとして視覚化できます。[ 3 ] [ 2 ]たとえば、以下に示す 3×4 行列は、水平線と垂直線によって 4 つのブロックに分割されています。左上の 2×3 ブロック、右上の 2×1 ブロック、左下の 1×3 ブロック、右下の 1×1 ブロックです。

[1つの111つの121つの13b11つの211つの221つの23b2c1c2c3d]{\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}メートル{\displaystyle m}M{\displaystyle M}n{\displaystyle n}行グループ{\displaystyle {\text{行グループ}}}メートル{\displaystyle m}コルグループ{\displaystyle {\text{colgroups}}}j{\displaystyle (i,j)}st{\displaystyle (s,t)}×y{\displaystyle (x,y)}×行グループ{\displaystyle x\in {\text{行グループ}}}yコルグループ{\displaystyle y\in {\text{colgroups}}}

ブロック行列代数は一般に行列のカテゴリ二重積から生じる。[ 5 ]

168×168要素のブロック行列。12×12、12×24、24×12、24×24のサブ行列が含まれています。非ゼロ要素は青色、ゼロ要素は灰色で表示されます。

マトリックス

P[1227156233453367]{\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=[1227156233453367].{\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}

P11=[1215],P12=[2762],P21=[3333],P22=[4567].{\displaystyle {\begin{aligned}\mathbf {P} _{11}&={\begin{bmatrix}1&2\\1&5\end{bmatrix}},&\mathbf {P} _{12}&={\begin{bmatrix}2&7\\6&2\end{bmatrix}},\\[1ex]\mathbf {P} _{21}&={\begin{bmatrix}3&3\\3&3\end{bmatrix}},&\mathbf {P} _{22}&={\begin{bmatrix}4&5\\6&7\end{bmatrix}}.\end{aligned}}}

分割された行列は次のように表される[ 8 ]

P=[P11P12P21P22].{\displaystyle \mathbf {P} ={\begin{bmatrix}\mathbf {P} _{11}&\mathbf {P} _{12}\\\mathbf {P} _{21}&\mathbf {P} _{22}\end{bmatrix}}.}

正式な定義

とする。の分割はの 表現であり、ACm×n{\displaystyle A\in \mathbb {C} ^{m\times n}}A{\displaystyle A}A{\displaystyle A}

A=[A11A12A1qA21A22A2qAp1Ap2Apq],{\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 ]AijCmi×nj{\displaystyle A_{ij}\in \mathbb {C} ^{m_{i}\times n_{j}}}i=1pmi=m{\textstyle \sum _{i=1}^{p}m_{i}=m}j=1qnj=n{\textstyle \sum _{j=1}^{q}n_{j}=n}Aij{\displaystyle A_{ij}}

この定義によれば、どの列のブロックも列数が同じでなければなりません。[ 9 ]同様に、どの行のブロックも行数が同じでなければなりません。[ 9 ]

パーティション分割方法

行列は様々な方法で分割することができる。[ 9 ]例えば、行列が次のように表される場合 、その行列は列によって分割されていると言われる。A{\displaystyle A}

A=(a1 a2  an),{\displaystyle A=(a_{1}\ a_{2}\ \cdots \ a_{n}),}

ここではの 番目の列である。[ 9 ]行列は行ごとに分割することもできる。 aj{\displaystyle a_{j}}j{\displaystyle j}A{\displaystyle A}

A=[a1Ta2TamT],{\displaystyle A={\begin{bmatrix}a_{1}^{T}\\a_{2}^{T}\\\vdots \\a_{m}^{T}\end{bmatrix}},}

ここでは の-行目である。[ 9 ]aiT{\displaystyle a_{i}^{T}}i{\displaystyle i}A{\displaystyle A}

共通パーティション

多くの場合、[ 9 ] 2×2分割に遭遇します

A=[A11A12A21A22],{\displaystyle A={\begin{bmatrix}A_{11}&A_{12}\\A_{21}&A_{22}\end{bmatrix}},}

特に の形で表され、 はスカラーである: [ 9 ]A11{\displaystyle A_{11}}

A=[a11a12Ta21A22].{\displaystyle A={\begin{bmatrix}a_{11}&a_{12}^{T}\\a_{21}&A_{22}\end{bmatrix}}.}

ブロック行列演算

転置

させて

A=[A11A12A1qA21A22A2qAp1Ap2Apq]{\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 ] [ 10 ]です。AijCki×j{\displaystyle A_{ij}\in \mathbb {C} ^{k_{i}\times \ell _{j}}}A{\displaystyle A}

AT=[A11TA21TAp1TA12TA22TAp2TA1qTA2qTApqT],{\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 ]

ブロック転置

ブロック行列に対しても、転置の特殊な形式を定義することができます。この形式では、個々のブロックは順序が変更されますが、転置は行われません。 をブロックを持つブロック行列とすると、 のブロック転置はを持つブロック行列です。[ 11 ]従来のトレース演算子と同様に、ブロック転置はとなる線形写像です。[ 10 ]しかし、一般に、と のブロックが可換でない限り、この性質は成り立ちません。 A=(Bij){\displaystyle A=(B_{ij})}k×l{\displaystyle k\times l}m×n{\displaystyle m\times n}Bij{\displaystyle B_{ij}}A{\displaystyle A}l×k{\displaystyle l\times k}AB{\displaystyle A^{\mathcal {B}}}m×n{\displaystyle m\times n}(AB)ij=Bji{\displaystyle \left(A^{\mathcal {B}}\right)_{ij}=B_{ji}}(A+C)B=AB+CB{\displaystyle (A+C)^{\mathcal {B}}=A^{\mathcal {B}}+C^{\mathcal {B}}}(AC)B=CBAB{\displaystyle (AC)^{\mathcal {B}}=C^{\mathcal {B}}A^{\mathcal {B}}}A{\displaystyle A}C{\displaystyle C}

追加

させて

B=[B11B12B1sB21B22B2sBr1Br2Brs],{\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}},}

ここで、 、 を§転置で定義された行列とする。(この行列は§乗算で再利用される。)、、、ならば、[ 9 ]BijCmi×nj{\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}ki=mi{\displaystyle k_{i}=m_{i}}j=nj{\displaystyle \ell _{j}=n_{j}}

A+B=[A11+B11A12+B12A1q+B1qA21+B21A22+B22A2q+B2qAp1+Bp1Ap2+Bp2Apq+Bpq].{\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}}.}

乗算

因子の部分行列の代数のみを扱うブロック分割行列積を用いることも可能です。しかしながら、因子の分割は任意ではなく、2つの行列間の「適合分割」 [ 12 ]と、使用されるすべての部分行列積が定義されるような「適合分割」[ 13 ]が必要です。A{\displaystyle A}B{\displaystyle B}

2 つの行列および が、積について等角的に分割されていると言われるのは、乗算が、部分行列をスカラーであるかのように扱いながら順序を保って実行され、関係する部分行列の積と和がすべて定義されている場合です。A{\displaystyle A}B{\displaystyle B}AB{\displaystyle AB}A{\displaystyle A}B{\displaystyle B}AB{\displaystyle AB}

— アラク・M・マタイとハンス・J・ハウボルド『線形代数:物理学者とエンジニアのためのコース』[ 14 ]

§転置で定義された行列をとし、§加法で定義された行列を とする。このとき、行列積 A{\displaystyle A}B{\displaystyle B}

C=AB{\displaystyle C=AB}

ブロックごとに実行することができ、行列が得られる。結果の行列の各行列は、以下の乗算によって計算される。[ 6 ]C{\displaystyle C}(p×s){\displaystyle (p\times s)}C{\displaystyle C}

Cij=k=1qAikBkj.{\displaystyle C_{ij}=\sum _{k=1}^{q}A_{ik}B_{kj}.}

または、繰り返されるインデックスを暗黙的に合計する アインシュタイン表記法を使用します。

Cij=AikBkj.{\displaystyle C_{ij}=A_{ik}B_{kj}.}

これを行列で表すと[ 9 ]となる。C{\displaystyle C}

C=AB=[i=1qA1iBi1i=1qA1iBi2i=1qA1iBisi=1qA2iBi1i=1qA2iBi2i=1qA2iBisi=1qApiBi1i=1qApiBi2i=1qApiBis].{\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}}.}

反転

行列が 4 つのブロックに分割されている場合、次のようにブロックごとに反転できます。

P=[ABCD]1=[A1+A1B(DCA1B)1CA1A1B(DCA1B)1(DCA1B)1CA1(DCA1B)1],{\displaystyle {\begin{aligned}P&={\begin{bmatrix}A&B\\C&D\end{bmatrix}}^{-1}\\[1ex]&={\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}},\end{aligned}}}

ここで、ADは任意サイズの正方ブロックであり、BCはそれらと分割可能である。さらに、 Aと、 PにおけるAのシュアー補集合:P / A = DCA −1 Bは逆行列でなければならない。[ 15 ]

同様に、ブロックを並べ替えると次のようになる。[ 16 ]

P=[ABCD]1=[(ABD1C)1(ABD1C)1BD1D1C(ABD1C)1D1+D1C(ABD1C)1BD1].{\displaystyle {\begin{aligned}P&={\begin{bmatrix}A&B\\C&D\end{bmatrix}}^{-1}\\[1ex]&={\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}&D^{-1}+D^{-1}C\left(A-BD^{-1}C\right)^{-1}BD^{-1}\end{bmatrix}}.\end{aligned}}}

ここで、DとPにおけるDのシュアー補集合:P / D = ABD −1 Cは逆行列でなければならない。

AD が両方とも逆である場合、次のようになります。

[ABCD]1=[(ABD1C)100(DCA1B)1][IBD1CA1I].{\displaystyle {\begin{bmatrix}A&B\\C&D\end{bmatrix}}^{-1}={\begin{bmatrix}\left(A-BD^{-1}C\right)^{-1}&0\\0&\left(D-CA^{-1}B\right)^{-1}\end{bmatrix}}{\begin{bmatrix}I&-BD^{-1}\\-CA^{-1}&I\end{bmatrix}}.}

ワインスタイン・アロンザイン恒等式により、ブロック対角行列の 2 つの行列のうち 1 つが逆行列である場合、もう 1 つも逆行列になります。

完全逆行列から部分逆行列を計算する

ブロック反転公式における行列とその逆行列間の対称性により、行列Pとその逆行列P −1が共形に分割される場合、

P=[ABCD],P1=[EFGH]{\displaystyle P={\begin{bmatrix}{A}&{B}\\{C}&{D}\end{bmatrix}},\quad P^{-1}={\begin{bmatrix}{E}&{F}\\{G}&{H}\end{bmatrix}}}

すると、任意の主部分行列の逆行列はP −1の対応するブロックから計算できる。

A1=EFH1G{\displaystyle {A}^{-1}={E}-{FH}^{-1}{G}}D1=HGE1F{\displaystyle {D}^{-1}={H}-{GE}^{-1}{F}}

この関係は、E −1 = ABD −1 C (シュアー補集合)であることを認識し、 PP −1の役割を逆にして同じブロック反転式を適用することから導かれる。[ 17 ] [ 18 ]

行列式

上の行列の行列式の公式は、適切な更なる仮定の下で、4つの部分行列と正方行列からなる行列に対しても成立する。ライプニッツの公式かシュール補集合を含む因数分解を用いて証明できる最も簡単な公式は[ 16 ]である。2×2{\displaystyle 2\times 2}A,B,C,D{\displaystyle A,B,C,D}A{\displaystyle A}D{\displaystyle D}det[A0CD]=det(A)det(D)=det[AB0D].{\displaystyle \det {\begin{bmatrix}A&0\\C&D\end{bmatrix}}=\det(A)\det(D)=\det {\begin{bmatrix}A&B\\0&D\end{bmatrix}}.}

この式を用いると、との特性多項式はとの特性多項式の積と同じであることが分かります。さらに、またはが対角化可能であれば、とも 対角化可能です。逆は偽です。 を確認してください。 [A0CD]{\displaystyle {\begin{bmatrix}A&0\\C&D\end{bmatrix}}}[AB0D]{\displaystyle {\begin{bmatrix}A&B\\0&D\end{bmatrix}}}A{\displaystyle A}D{\displaystyle D}[A0CD]{\displaystyle {\begin{bmatrix}A&0\\C&D\end{bmatrix}}}[AB0D]{\displaystyle {\begin{bmatrix}A&B\\0&D\end{bmatrix}}}A{\displaystyle A}D{\displaystyle D}[1101]{\displaystyle {\begin{bmatrix}1&1\\0&1\end{bmatrix}}}

逆行列である場合、[ 16 ]A{\displaystyle A}

det[ABCD]=det(A)det(DCA1B),{\displaystyle \det {\begin{bmatrix}A&B\\C&D\end{bmatrix}}=\det(A)\det \left(D-CA^{-1}B\right),}

が可逆であれば、[ 19 ] [ 16 ]D{\displaystyle D}

det[ABCD]=det(D)det(ABD1C).{\displaystyle \det {\begin{bmatrix}A&B\\C&D\end{bmatrix}}=\det(D)\det \left(A-BD^{-1}C\right).}

ブロックが同じサイズの正方行列である場合、さらに式が成り立ちます。例えば、と が可換である場合(つまり、)、[ 20 ] 同様のステートメントが、、の場合にも成り立ちます。つまり、 の場合、 と の順序の変更に注意してください(の代わりに)。同様に、の場合、 はに置き換える必要があります(つまり になります)。また、 の場合、 になります。最後の 2 つの結果では、基礎となるリングの可換性を使用する必要がありますが、最初の 2 つの結果では必要ありません。 C{\displaystyle C}D{\displaystyle D}CD=DC{\displaystyle CD=DC}det[ABCD]=det(ADBC).{\displaystyle \det {\begin{bmatrix}A&B\\C&D\end{bmatrix}}=\det(AD-BC).}AB=BA{\displaystyle AB=BA}AC=CA{\displaystyle AC=CA}BD=DB{\displaystyle BD=DB}AC=CA{\displaystyle AC=CA}det[ABCD]=det(ADCB).{\displaystyle \det {\begin{bmatrix}A&B\\C&D\end{bmatrix}}=\det(AD-CB).}C{\displaystyle C}B{\displaystyle B}CB{\displaystyle CB}BC{\displaystyle BC}BD=DB{\displaystyle BD=DB}AD{\displaystyle AD}DA{\displaystyle DA}det(DABC){\displaystyle \det(DA-BC)}AB=BA{\displaystyle AB=BA}det(DACB){\displaystyle \det(DA-CB)}

この式は、個々のブロック間に適切な可換条件が満たされる場合、 個以上のブロックから構成される行列に一般化されている。 [ 21 ]2×2{\displaystyle 2\times 2}

とについては、次の式が成り立つ(と が可換でなくても)[ 16 ]A=D{\displaystyle A=D}B=C{\displaystyle B=C}A{\displaystyle A}B{\displaystyle B}det[ABBA]=det(AB)det(A+B).{\displaystyle \det {\begin{bmatrix}A&B\\B&A\end{bmatrix}}=\det(A-B)\det(A+B).}

特殊なタイプのブロック行列

直和とブロック対角行列

直和

任意の行列A(サイズm  ×  n)とB(サイズp  ×  q)に対して、AB直和はA⊕Bと表され、のように定義される[ 10 ]  。 

AB=[a11a1n00am1amn0000b11b1q00bp1bpq].{\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}}.}

例えば、

[132231][1601]=[13200231000001600001].{\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}}.}

この操作は、任意の次元の配列に自然に一般化されます ( AB の次元数が同じである 場合)。

2 つの行列のベクトル空間の直和内の任意の要素は、2 つの行列の直和として表すことができることに注意してください。

ブロック対角行列

ブロック対角行列とは、主対角ブロックが正方行列で、対角外ブロックがすべて零行列であるような正方行列であるブロック行列である。 [ 16 ]つまり、ブロック対角行列Aは次の形式を持つ。

A=[A1000A2000An]{\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つのブロックのみを持つブロック対角行列と見なすことができる。

行列式トレースについては、次の特性が成り立ちます。

detA=detA1××detAn,{\displaystyle {\begin{aligned}\det {A}&=\det {A}_{1}\times \cdots \times \det {A}_{n},\end{aligned}}}[ 22 ] [ 23 ]および
trA=trA1++trAn.{\displaystyle {\begin{aligned}\operatorname {tr} {A}&=\operatorname {tr} {A}_{1}+\cdots +\operatorname {tr} {A}_{n}.\end{aligned}}}[ 16 ] [ 23 ]

ブロック対角行列が逆行列となるのは、その主対角ブロックのそれぞれが逆行列となる場合のみであり、この場合、その逆行列は[ 24 ]で与えられる別のブロック対角行列となる。[A1000A2000An]1=[A11000A21000An1].{\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}}.}

の固有値[ 25 ]と固有ベクトル単にそれらを結合したものである。[ 23 ]A{\displaystyle {A}}Ak{\displaystyle {A}_{k}}

ブロック三角行列

ブロック三重対角行列は、もう一つの特殊なブロック行列です。ブロック対角行列と同様に、正方行列であり、下対角、主対角、上対角に正方行列(ブロック)を持ち、その他のブロックは零行列です。ブロック三重対角行列は本質的に三重対角行列ですが、スカラーの代わりに部分行列を持ちます。ブロック三重対角行列は、以下の形式を持ちます 。A{\displaystyle A}

A=[B1C10A2B2C2AkBkCkAn1Bn1Cn10AnBn]{\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 ]Ak{\displaystyle {A}_{k}}Bk{\displaystyle {B}_{k}}Ck{\displaystyle {C}_{k}}

ブロック三角行列は、工学問題(例えば、数値流体力学)の数値解法においてしばしば用いられる。LU分解のための最適化された数値解法[ 28 ]が利用可能であり、ブロック三角行列を係数行列として持つ方程式系に対する効率的な解法アルゴリズムが存在する。三角行列を含む方程式系を効率的に解くために使用されるトーマスアルゴリズムは、行列演算を用いてブロック三角行列に適用することもできる(ブロックLU分解も参照)。

ブロック三角行列

行列が上ブロック三角行列(またはブロック上三角行列[ 29 ])であるとは、およびと なる正の整数が存在する場合である。 ここで、 行列はすべての に対して となる。 [ 25 ] [ 29 ] 同様に、が下ブロック三角行列である場合は、 および となる正の整数が存在する場合である。[ 25 ]n×n{\displaystyle n\times n}A{\displaystyle A}n1,,nk{\displaystyle n_{1},\ldots ,n_{k}}n=n1+n2++nk{\displaystyle n=n_{1}+n_{2}+\ldots +n_{k}}A=[A11A12A1k0A22A2k00Akk],{\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}},}Aij{\displaystyle A_{ij}}ni×nj{\displaystyle n_{i}\times n_{j}}i,j=1,,k{\displaystyle i,j=1,\ldots ,k}A{\displaystyle A}A=[A1100A21A220Ak1Ak2Akk],{\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}},}Aij{\displaystyle A_{ij}}ni×nj{\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)}}ki=lj{\displaystyle k-i=l-j}

A=[A1A2A3A4A1A2A5A4A1],{\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 ]AiFni×mi{\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=[A1A2A3A2A3A4A3A4A5],{\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 ]AiFni×mi{\displaystyle A_{i}\in \mathbb {F} ^{n_{i}\times m_{i}}}

参照

注記

  1. ^イヴスハワード(1980). 『初等行列理論』(復刻版). ニューヨーク: ドーバー. p.  37. ISBN 0-486-63946-02013年4月24日閲覧行列を要素の長方形ブロックに分割すると便利な場合があることがわかります。そこで、いわゆる分割行列、あるいはブロック行列について考えてみましょう。
  2. ^ a b Dobrushkin, Vladimir. 「Partition Matrices」 . Mathematicaによる線形代数. 2024年3月24日閲覧
  3. ^アントン・ハワード (1994).初等線形代数(第7版). ニューヨーク: ジョン・ワイリー. p. 30. ISBN 0-471-58742-7選択した行と列の間に水平線と垂直線を挿入することで、行列をより小さな行列に分割することができます。
  4. ^ 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 .
  5. ^ Macedo, HD; Oliveira, JN (2013). 「線形代数の型付け:二積指向アプローチ」.コンピュータプログラミング科学. 78 (11): 2160– 2191. arXiv : 1312.4818 . doi : 10.1016/j.scico.2012.07.012 .
  6. ^ a b cジョンストン、ナサニエル (2021).線形代数と行列代数入門. シャム、スイス: シュプリンガー・ネイチャー. pp. 30, 425. ISBN 978-3-030-52811-9
  7. ^ a bジョンストン、ナサニエル (2021).上級線形代数と行列代数. シャム、スイス: シュプリンガー・ネイチャー. p. 298. ISBN 978-3-030-52814-0
  8. ^ジェフリー・アラン (2010).エンジニアと科学者のための行列演算:線形代数の必須ガイド. ドルドレヒト [オランダ] ; ニューヨーク: シュプリンガー. p. 54. ISBN 978-90-481-9273-1. OCLC  639165077 .
  9. ^ a b c d e f g h i j k l mスチュワート、ギルバート・W. (1998).行列アルゴリズム. 1: 基本分解. フィラデルフィア、ペンシルバニア州: 産業応用数学協会. pp.  18– 20. ISBN 978-0-89871-414-2
  10. ^ a b c d e Gentle, James E. (2007). Matrix Algebra: Theory, Computations, and Applications in Statistics . Springer Texts in Statistics. New York: Springer New York Springer e-books. pp. 47, 487. ISBN 978-0-387-70873-7
  11. ^ Mackey, D. Steven (2006).行列多項式の構造化線形化(PDF) (論文). マンチェスター大学. ISSN 1749-9097 . OCLC 930686781 .  
  12. ^イヴスハワード(1980). 『初等行列理論』(復刻版). ニューヨーク: ドーバー. p.  37. ISBN 0-486-63946-0. 2013年4月24日閲覧定理1.9.4のような分割は、AB適合分割と呼ばれる。
  13. ^アントン・ハワード (1994).初等線形代数(第7版). ニューヨーク: ジョン・ワイリー. p. 36. ISBN 0-471-58742-7...ただし、A と B のサブマトリックスのサイズが、指定された操作を実行できるサイズである必要があります。
  14. ^ Mathai, Arakaparampil M.; Haubold, Hans J. (2017).線形代数:物理学者とエンジニアのためのコース. De Gruyter 教科書. ベルリン・ボストン: De Gruyter. p. 162. ISBN 978-3-11-056259-0
  15. ^バーンスタイン、デニス (2005).行列数学. プリンストン大学出版局. p. 44. ISBN 0-691-11802-7
  16. ^ a b c d e f g hアバディール、カリム M.;マグナス、ジャン R. (2005)。行列代数。ケンブリッジ大学出版局。 97、100、106、111、114、118ページ。ISBN 9781139443647
  17. ^ 「このブロック逆行列の式は、逆行列全体に関して既知ですか?」MathOverflow
  18. ^ Escalante-B., Alberto N.; Wiskott, Laurenz (2016). 「改良されたグラフベースSFA:情報保存は遅さの原則を補完する」 .機械学習. arXiv : 1412.4679 . doi : 10.1007/s10994-016-5563-y .
  19. ^ Taboga, Marco (2021). 「ブロック行列の行列式」, 行列代数の講義.
  20. ^ Silvester, JR (2000). 「ブロック行列の行列式」(PDF) . Math. Gaz . 84 (501): 460– 467. doi : 10.2307/3620776 . JSTOR 3620776. 2015年3月18日時点のオリジナル(PDF)からアーカイブ。 2021年6月25日閲覧 
  21. ^ Sothanaphan, Nat (2017年1月). 「非可換ブロックを持つブロック行列の行列式」.線形代数とその応用. 512 : 202–218 . arXiv : 1805.06027 . doi : 10.1016/j.laa.2016.10.004 . S2CID 119272194 . 
  22. ^ Quarteroni, Alfio; Sacco, Riccardo; Saleri, Fausto (2000).数値数学. 応用数学テキスト. ニューヨーク: Springer. pp. 10, 13. ISBN 978-0-387-98959-4
  23. ^ a b cジョージ、ラジュ K.;アジャヤクマール、アビジット(2024)。「線形代数のコース」数理科学の大学テキスト: 35, 407. doi : 10.1007/978-981-99-8680-4ISBN 978-981-99-8679-8. ISSN  2731-9318 .
  24. ^プリンス、サイモン JD (2012). 『コンピュータビジョン:モデル、学習、推論』 ニューヨーク:ケンブリッジ大学出版局. p. 531. ISBN 978-1-107-01179-3
  25. ^ a b c d eバーンスタイン、デニス・S. (2009).行列数学:理論、事実、公式(第2版). プリンストン、ニュージャージー:プリンストン大学出版局. pp. 168, 298. ISBN 978-0-691-14039-1
  26. ^ Dietl, Guido KE (2007).クリロフ部分空間における線形推定と検出. 信号処理、通信、ネットワークの基礎 . ベルリン; ニューヨーク: Springer. pp. 85, 87. ISBN 978-3-540-68478-7. OCLC  85898525 .
  27. ^ Horn, Roger A.; Johnson, Charles R. (2017).マトリックス分析(第2版、訂正再版). ニューヨーク:ケンブリッジ大学出版局. p. 36. ISBN 978-0-521-83940-2
  28. ^ Datta, Biswa Nath (2010).数値線形代数とその応用(第2版). フィラデルフィア, ペンシルバニア州: SIAM. p. 168. ISBN 978-0-89871-685-6
  29. ^ a bスチュワート、ギルバート・W. (2001).行列アルゴリズム. 2: 固有値システム. フィラデルフィア、ペンシルバニア州: 産業応用数学協会. p. 5. ISBN 978-0-89871-503-3

参考文献