数学的行列の性質
数学 において、 実数 要素 を持つ 対称行列が 正定値行列 であるとは、実数が すべての非零の実 列ベクトル に対して正であるときである。ここで、 は行ベクトル の 転置 である [1]。 より一般的には、 エルミート行列 (つまり、その 共役転置 に等しい 複素行列 )が 正定値行列 であるとは、実数が すべての非零の複素列ベクトルに対して正であるときである。 ここで、 はの共役転置を表す。 M {\displaystyle M} × T M × {\displaystyle \mathbf {x} ^{\mathsf {T}}M\mathbf {x} } × 、 {\displaystyle \mathbf {x} ,} × T {\displaystyle \mathbf {x} ^{\mathsf {T}}} × 。 {\displaystyle \mathbf {x} .} z ∗ M z {\displaystyle \mathbf {z} ^{*}M\mathbf {z} } z 、 {\displaystyle \mathbf {z} ,} z ∗ {\displaystyle \mathbf {z} ^{*}} z 。 {\displaystyle \mathbf {z} .}
半正定値 行列も同様に定義されますが、スカラー とが正 またはゼロ (つまり非負) である必要があります。 負定値行列 と 半負定値 行列も同様に定義されます。半正定値でも半負定値でもない行列は、 不定値行列 と呼ばれることがあります。 × T M × {\displaystyle \mathbf {x} ^{\mathsf {T}}M\mathbf {x} } z ∗ M z {\displaystyle \mathbf {z} ^{*}M\mathbf {z} }
一部の著者は、より一般的な正定性の定義を用いて、行列が非対称行列や非エルミート行列であっても許容する。これらの一般化された正定値行列の性質については、後述の「非エルミート正方行列への拡張」で考察するが、本稿の主題ではない。
定義 以下の定義では、 は の転置、 は の 共役転置 、 は n 次元の 零ベクトル を表します。 × T {\displaystyle \mathbf {x} ^{\mathsf {T}}} × 、 {\displaystyle \mathbf {x} ,} z ∗ {\displaystyle \mathbf {z} ^{*}} z 、 {\displaystyle \mathbf {z} ,} 0 {\displaystyle \mathbf {0} }
実数行列の定義 対称 実行列は、 すべての非ゼロの 値が、 正式 に は 、 n × n {\displaystyle n\times n} M {\displaystyle M} × T M × > 0 {\displaystyle \mathbf {x} ^{\mathsf {T}}M\mathbf {x} >0} × {\displaystyle \mathbf {x} } R n 。 {\displaystyle \mathbb {R} ^{n}.}
M 正定値 ⟺ × T M × > 0 すべての人のために × ∈ R n ∖ { 0 } {\displaystyle M{\text{ 正定値}}\quad \iff \quad \mathbf {x} ^{\mathsf {T}}M\mathbf {x} >0{\text{ すべての }}\mathbf {x} \in \mathbb {R} ^{n}\setminus \{\mathbf {0} \}} に対して
対称 実行列は、 すべてのに対して が成り立つ とき 、半正定値 行列または 非負定値行列 であると言われる 。 n × n {\displaystyle n\times n} M {\displaystyle M} × T M × ≥ 0 {\displaystyle \mathbf {x} ^{\mathsf {T}}M\mathbf {x} \geq 0} × {\displaystyle \mathbf {x} } R n 。 {\displaystyle \mathbb {R} ^{n}.}
M 半正定値 ⟺ × T M × ≥ 0 すべての人のために × ∈ R n {\displaystyle M{\text{ 半正定値}}\quad \iff \quad \mathbf {x} ^{\mathsf {T}}M\mathbf {x} \geq 0{\text{ すべての }}\mathbf {x} \in \mathbb {R} ^{n}} に対して
対称 実行列は、 すべての非ゼロの値が 負定値 である とは、正式 には 、 n × n {\displaystyle n\times n} M {\displaystyle M} × T M × < 0 {\displaystyle \mathbf {x} ^{\mathsf {T}}M\mathbf {x} <0} × {\displaystyle \mathbf {x} } R n 。 {\displaystyle \mathbb {R} ^{n}.}
M 負定値 ⟺ × T M × < 0 すべての人のために × ∈ R n ∖ { 0 } {\displaystyle M{\text{ 負定値}}\quad \iff \quad \mathbf {x} ^{\mathsf {T}}M\mathbf {x} <0{\text{ 全ての }}\mathbf {x} \in \mathbb {R} ^{n}\setminus \{\mathbf {0} \}} に対して
対称 実行列は、 すべての に対して で あるとき 、半負定値行列 または 非正定値行列 であるといわれる 。 n × n {\displaystyle n\times n} M {\displaystyle M} × T M × ≤ 0 {\displaystyle \mathbf {x} ^{\mathsf {T}}M\mathbf {x} \leq 0} × {\displaystyle \mathbf {x} } R n 。 {\displaystyle \mathbb {R} ^{n}.}
M 負の半定値 ⟺ × T M × ≤ 0 すべての人のために × ∈ R n {\displaystyle M{\text{ 負の半定値}}\quad \iff \quad \mathbf {x} ^{\mathsf {T}}M\mathbf {x} \leq 0{\text{ すべての }}\mathbf {x} \in \mathbb {R} ^{n}} に対して
半正定値でも半負定値でもない対称 実数行列は 不定行列 と呼ばれます。 n × n {\displaystyle n\times n}
複素行列の定義 以下の定義はすべてこの用語を含んでいます。 これは、任意のエルミート正方行列に対して常に実数であることに注意してください。 z ∗ M z 。 {\displaystyle \mathbf {z} ^{*}M\mathbf {z} .} M 。 {\displaystyle M.}
エルミート 複素行列は、 すべての非ゼロの値が 正定値 である とは、正式 には 、 n × n {\displaystyle n\times n} M {\displaystyle M} z ∗ M z > 0 {\displaystyle \mathbf {z} ^{*}M\mathbf {z} >0} z {\displaystyle \mathbf {z} } C n 。 {\displaystyle \mathbb {C} ^{n}.}
M 正定値 ⟺ z ∗ M z > 0 すべての人のために z ∈ C n ∖ { 0 } {\displaystyle M{\text{ 正定値}}\quad \iff \quad \mathbf {z} ^{*}M\mathbf {z} >0{\text{ 全ての }}\mathbf {z} \in \mathbb {C} ^{n}\setminus \{\mathbf {0} \}} に対して
エルミート 複素行列は、 すべての に対して で あるとき、 半正定値行列 または 非負定値行列 であると言われる 。 n × n {\displaystyle n\times n} M {\displaystyle M} z ∗ M z ≥ 0 {\displaystyle \mathbf {z} ^{*}M\mathbf {z} \geq 0} z {\displaystyle \mathbf {z} } C n 。 {\displaystyle \mathbb {C} ^{n}.}
M 半正定値 ⟺ z ∗ M z ≥ 0 すべての人のために z ∈ C n {\displaystyle M{\text{ 半正定値}}\quad \iff \quad \mathbf {z} ^{*}M\mathbf {z} \geq 0{\text{ すべての }}\mathbf {z} \in \mathbb {C} ^{n}} に対して
エルミート 複素行列は、 すべての非ゼロの 値が 負定値 であるとされる 。正式 には、 n × n {\displaystyle n\times n} M {\displaystyle M} z ∗ M z < 0 {\displaystyle \mathbf {z} ^{*}M\mathbf {z} <0} z {\displaystyle \mathbf {z} } C n 。 {\displaystyle \mathbb {C} ^{n}.}
M 負定値 ⟺ z ∗ M z < 0 すべての人のために z ∈ C n ∖ { 0 } {\displaystyle M{\text{ 負定値}}\quad \iff \quad \mathbf {z} ^{*}M\mathbf {z} <0{\text{ 全ての }}\mathbf {z} \in \mathbb {C} ^{n}\setminus \{\mathbf {0} \}} に対して
エルミート 複素行列は、 すべての に対して で あるとき、 負の半定値行列 または 非正定値行列 であると言われる。正式に は 、 n × n {\displaystyle n\times n} M {\displaystyle M} z ∗ M z ≤ 0 {\displaystyle \mathbf {z} ^{*}M\mathbf {z} \leq 0} z {\displaystyle \mathbf {z} } C n 。 {\displaystyle \mathbb {C} ^{n}.}
M 負の半定値 ⟺ z ∗ M z ≤ 0 すべての人のために z ∈ C n {\displaystyle M{\text{ 負の半定値}}\quad \iff \quad \mathbf {z} ^{*}M\mathbf {z} \leq 0{\text{ すべての }}\mathbf {z} \in \mathbb {C} ^{n}} に対して
半正定値でも半負定値でもないエルミート 複素行列は 不定値 と呼ばれます。 n × n {\displaystyle n\times n}
実数定義と複素数定義の一貫性 すべての実行列は複素行列でもあるため、2 つのクラスの「定性」の定義は一致する必要があります。
複素行列の場合、最も一般的な定義は、 任意の非零複素列ベクトルに対して が実数かつ正である 場合に限り、 が正定値行列である、というものである。 この条件は、 がエルミート行列(すなわち、その転置が共役行列に等しい)であることを意味する。なぜなら、 が実数であるので、 任意の に対して が共役転置に等しいからである。 M {\displaystyle M} z ∗ M z {\displaystyle \mathbf {z} ^{*}M\mathbf {z} } z 。 {\displaystyle \mathbf {z} .} M {\displaystyle M} z ∗ M z {\displaystyle \mathbf {z} ^{*}M\mathbf {z} } z ∗ M ∗ z {\displaystyle \mathbf {z} ^{*}M^{*}\mathbf {z} } z 、 {\displaystyle \mathbf {z} ,} M = M ∗ 。 {\displaystyle M=M^{*}.}
この定義によれば、正定値 実数 行列 はエルミート行列であり、したがって対称行列である。また、 すべての非零 実 列ベクトルに対しては正となる 。しかし、最後の条件だけでは、 正定値行列であるためには十分ではない。例えば、 M {\displaystyle M} z T M z {\displaystyle \mathbf {z} ^{\mathsf {T}}M\mathbf {z} } z 。 {\displaystyle \mathbf {z} .} M {\displaystyle M} M = [ 1 1 − 1 1 ] 、 {\displaystyle M={\begin{bmatrix}1&1\\-1&1\end{bmatrix}},}
とを持つ任意の 実ベクトルに対して、 が成り立ち 、が0でなければ 常に正である 。しかし、が 1 と を持つ複素ベクトルの場合 、 z {\displaystyle \mathbf {z} } 1つの {\displaystyle a} b {\displaystyle b} z T M z = ( 1つの + b ) 1つの + ( − 1つの + b ) b = 1つの 2 + b 2 、 {\displaystyle \mathbf {z} ^{\mathsf {T}}M\mathbf {z} =\left(a+b\right)a+\left(-a+b\right)b=a^{2}+b^{2},} z {\displaystyle \mathbf {z} } z {\displaystyle \mathbf {z} } 私 {\displaystyle i}
z ∗ M z = [ 1 − 私 ] M [ 1 私 ] = [ 1 + 私 1 − 私 ] [ 1 私 ] = 2 + 2 私 。 {\displaystyle \mathbf {z} ^{*}M\mathbf {z} ={\begin{bmatrix}1&-i\end{bmatrix}}M{\begin{bmatrix}1\\i\end{bmatrix}}={\begin{bmatrix}1+i&1-i\end{bmatrix}}{\begin{bmatrix}1\\i\end{bmatrix}}=2+2i.}
これは実数ではない。したがって、 正定値ではない。 M {\displaystyle M}
一方、 対称 実数行列の場合、「 すべての非ゼロ実数ベクトルに対して 」 という条件は、 複素数の意味で正定値である ことを意味 します。 M 、 {\displaystyle M,} z T M z > 0 {\displaystyle \mathbf {z} ^{\mathsf {T}}M\mathbf {z} >0} z {\displaystyle \mathbf {z} } M {\displaystyle M}
表記 エルミート行列が 半正定値行列である場合、 と書くことがあります。 また、 が正定値行列である場合は と書きます。 が半負定値行列である ことを示すには と書き、 が負定値行列である ことを示すには と書きます。 M {\displaystyle M} M ⪰ 0 {\displaystyle M\succeq 0} M {\displaystyle M} M ≻ 0。 {\displaystyle M\succ 0.} M {\displaystyle M} M ⪯ 0 {\displaystyle M\preceq 0} M {\displaystyle M} M ≺ 0。 {\displaystyle M\prec 0.}
この概念は 関数解析 から来ており、半正定値行列は 正の作用素 を定義します。2つの行列 と が を満たす場合、 反射的 、 反対称的 、 推移 的な 非厳密な半順序 を定義できます 。ただし、これは 全順序 ではなく、 一般に は不定となる場合があります。 あ {\displaystyle A} B {\displaystyle B} B − あ ⪰ 0 、 {\displaystyle BA\succeq 0,} B ⪰ あ {\displaystyle B\succeq A} B − あ 、 {\displaystyle BA,}
一般的な代替表記法 として、 半正定値行列と正定値行列、半負定値行列と負定値行列をそれぞれ と と表記する。これは、 非負行列 (それぞれ非正定値行列)もこの表記法で表記されることがあるため、混乱を招く可能性がある。 M ≥ 0 、 {\displaystyle M\geq 0,} M > 0 、 {\displaystyle M>0,} M ≤ 0 、 {\displaystyle M\leq 0,} M < 0 {\displaystyle M<0}
影響 上記の定義から、エルミート行列が正定値行列となる のは、それが 正定値二次形式 または エルミート形式 の行列である場合に限ります 。言い換えれば、エルミート行列が正定値行列となるのは、それが 内積 を定義する場合に限ります。
正定値行列と半正定値行列は様々な方法で特徴付けることができ、数学の様々な分野におけるこの概念の重要性を説明できるかもしれません。エルミート行列 M が正定値行列であるためには、以下の同値な条件のいずれかを満たす必要があります。
M {\displaystyle M} は、正の実数要素を持つ対角行列 と 合同 です 。 M {\displaystyle M} はエルミート行列であり、その 固有値は すべて実数かつ正です。 M {\displaystyle M} はエルミート行列であり、その主要な 主マイナー行列 はすべて正です。 共役転置を持つ 可逆行列 が存在し 、 B {\displaystyle B} B ∗ {\displaystyle B^{*}} M = B ∗ B 。 {\displaystyle M=B^{*}B.} 行列が半正定値行列であるとは、"正" を "非負" に置き換え、"可逆行列" を "行列" に置き換え、"先頭" という語を削除した同様の条件を満たす場合です。
正定値行列と半正定値実数行列は凸最適化 の基礎となる 。なぜなら、 2回 微分可能な 複数の実変数の関数 が与えられたとき、その ヘッセ行列 (その2番目の偏導関数の行列)が点pで正定値で あれば関数は pの 近くで 凸 であり、逆に関数がpの近くで凸であれば ヘッセ行列はpで半正定値となるからである。 p 、 {\displaystyle p,} p 、 {\displaystyle p,} p 。 {\displaystyle p.}
正定値行列の集合は 開いた 凸錐 であり、半正定値行列の集合は 閉じた 凸錐である。 [2]
例
固有値 をエルミート行列 とする ( 実 対称行列 を含む)。 のすべての固有値は 実数であり、その符号はその正定性を特徴付ける: M {\displaystyle M} n × n {\displaystyle n\times n} M {\displaystyle M}
M {\displaystyle M} は、その固有値がすべて正である場合にのみ正定値です。 M {\displaystyle M} 全ての固有値が非負である場合にのみ、半正定値となります。 M {\displaystyle M} は、その固有値がすべて負である場合にのみ負定値です。 M {\displaystyle M} すべての固有値が非正である場合に限り、半正定値となります。 M {\displaystyle M} 正の固有値と負の固有値の両方を持つ場合にのみ不定となります。 を の 固有分解 とします。 ここ では ユニタリ複素行列 であり、 その列は の固有ベクトルの直交基底から成り 、 は 実 対 角 行列 で あり、 その 主対角線 には対応する 固有値 が含まれます。この行列は、 (固有ベクトル)基底の座標で再表現された 対角行列とみなすことができます。 言い換えれば、 あるベクトルに を 適用することは、 を用いて 基底を 固有ベクトル座標系に変更する ことと同じです。 結果に を 引き伸ばし変換 を適用し、 とし 、その後 を用いて基底を元に戻します 。 P D P − 1 {\displaystyle PDP^{-1}} M 、 {\displaystyle M,} P {\displaystyle P} M 、 {\displaystyle M,} D {\displaystyle D} M {\displaystyle M} D {\displaystyle D} P 。 {\displaystyle P.} M {\displaystyle M} z 、 {\displaystyle \mathbf {z} ,} M z 、 {\displaystyle M\mathbf {z} ,} P − 1 、 {\displaystyle P^{-1},} P − 1 z 、 {\displaystyle P^{-1}\mathbf {z} ,} D {\displaystyle D} D P − 1 z 、 {\displaystyle DP^{-1}\mathbf {z} ,} P 、 {\displaystyle P,} P D P − 1 z 。 {\displaystyle PDP^{-1}\mathbf {z} .}
これを念頭に置くと、変数の一対一変換は、任意の複素ベクトルに対して が 実数かつ正である 場合と、任意 の が実数かつ正である場合 、言い換えれば が正定値である場合に限り、ということが分かります。対角行列の場合、これは主対角線の各要素、つまり のすべての固有値 が正である場合にのみ当てはまります。 スペクトル定理 により、エルミート行列のすべての固有値は実数であることが保証されるため、実対称行列の 特性多項式 が利用可能な場合、 デカルトの符号交替則 を用いて固有値の正性を確認できます 。 y = P z {\displaystyle \mathbf {y} =P\mathbf {z} } z ∗ M z {\displaystyle \mathbf {z} ^{*}M\mathbf {z} } z {\displaystyle \mathbf {z} } y ∗ D y {\displaystyle \mathbf {y} ^{*}D\mathbf {y} } y ; {\displaystyle y;} D {\displaystyle D} M {\displaystyle M} M {\displaystyle M}
分解 をエルミート行列 とし ます 。
が半正定値行列となるのは、 共役転置行列 と の積 として分解できる場合のみです 。 M {\displaystyle M} n × n {\displaystyle n\times n} M {\displaystyle M} M = B ∗ B {\displaystyle M=B^{*}B} B {\displaystyle B}
が実数の場合 、 も実数となり、分解は次のように書ける。 M {\displaystyle M} B {\displaystyle B} M = B T B 。 {\displaystyle M=B^{\mathsf {T}}B.}
M {\displaystyle M} が正定値であることと、逆行列 を 持つ分解が存在することは同値である 。より一般的には、が階数 を持つ 半正 定値 であることと、階数 を持つ分解が存在することは同値である。さらに、 任意 の分解 [3]に対して B {\displaystyle B} M {\displaystyle M} け {\displaystyle k} け × n {\displaystyle k\times n} B {\displaystyle B} け {\displaystyle k} M = B ∗ B 、 {\displaystyle M=B^{*}B,} ランク ( M ) = ランク ( B ) 。 {\displaystyle \operatorname {rank} (M)=\operatorname {rank} (B).}
の 列は、それぞれ 複素 ベクトル空間または 実ベクトル空間 のベクトルとして見ることができます 。すると、 の要素はこれらのベクトルの 内積 (つまり 、実数の場合は ドット積 )です。 言い換えれば、エルミート行列 が半正定値行列となるのは、それが いくつかのベクトルの グラム行列 である場合に
限ります。エルミート行列が半正定値行列となるのは、それがいくつかの 線形独立 ベクトルのグラム行列である場合に限ります。一般に、ベクトルのグラム行列の階数は、これらのベクトルが 張る 空間の次元に等しいです 。 [4] b 1 、 … 、 b n {\displaystyle b_{1},\dots ,b_{n}} B {\displaystyle B} R け 、 {\displaystyle \mathbb {R} ^{k},} M {\displaystyle M} M 私 j = ⟨ b 私 、 b j ⟩ 。 {\displaystyle M_{ij}=\langle b_{i},b_{j}\rangle .} M {\displaystyle M} b 1 、 … 、 b n 。 {\displaystyle b_{1},\dots ,b_{n}.} b 1 、 … 、 b n {\displaystyle b_{1},\dots ,b_{n}}
分解は一意ではない。 ある 行列に対して が成り立ち 、が任意の ユニタリ 行列(つまり ) である とき、 M = B ∗ B {\displaystyle M=B^{*}B} け × n {\displaystyle k\times n} B {\displaystyle B} 質問 {\displaystyle Q} け × け {\displaystyle k\times k} 質問 ∗ 質問 = 質問 質問 ∗ = 私 {\displaystyle Q^{*}Q=QQ^{*}=I} M = B ∗ B = B ∗ 質問 ∗ 質問 B = あ ∗ あ {\displaystyle M=B^{*}B=B^{*}Q^{*}QB=A^{*}A} あ = 質問 B 。 {\displaystyle A=QB.}
しかし、これは2つの分解が異なる唯一の方法である。分解は ユニタリ変換 を除いて一意である。より正式には、が 行列 であり、が行列 であって 、となる場合 、直交列(つまり )を持つ 行列 が存在し、 [5] となる。この とき、 は ユニタリで あることを意味する。 あ {\displaystyle A} け × n {\displaystyle k\times n} B {\displaystyle B} ℓ × n {\displaystyle \ell \times n} あ ∗ あ = B ∗ B 、 {\displaystyle A^{*}A=B^{*}B,} ℓ × け {\displaystyle \ell \times k} 質問 {\displaystyle Q} 質問 ∗ 質問 = 私 け × け {\displaystyle Q^{*}Q=I_{k\times k}} B = 質問 あ 。 {\displaystyle B=QA.} ℓ = け {\displaystyle \ell =k} 質問 {\displaystyle Q}
この記述は、実数の場合、直感的な幾何学的解釈が可能です。つまり、 と の列を の ベクトル ととし ます。 実ユニタリ行列は 直交行列 であり、 0 点を保存する(つまり、 の 回転 と 鏡映 を平行移動なしで) 剛体変換 (ユークリッド空間 の等長変換) を記述します 。したがって、の何らかの剛体変換によって ベクトルが に変換される(および 0 が 0 に変換される) 場合にのみ、 と の内積は等しくなります 。 あ {\displaystyle A} B {\displaystyle B} 1つの 1 、 … 、 1つの n {\displaystyle a_{1},\dots ,a_{n}} b 1 、 … 、 b n {\displaystyle b_{1},\dots ,b_{n}} R け 。 {\displaystyle \mathbb {R} ^{k}.} R け {\displaystyle \mathbb {R} ^{k}} 1つの 私 ⋅ 1つの j {\displaystyle a_{i}\cdot a_{j}} b 私 ⋅ b j {\displaystyle b_{i}\cdot b_{j}} R け {\displaystyle \mathbb {R} ^{k}} 1つの 1 、 … 、 1つの n {\displaystyle a_{1},\dots ,a_{n}} b 1 、 … 、 b n {\displaystyle b_{1},\dots ,b_{n}}
平方根 エルミート行列 が正半定値行列である必要十分条件は、次を満たす正半定値行列 (特に は エルミート行列なので )が存在することである。 この行列は 一意であり、 [6] は の 非負 平方根 と呼ばれ 、 で表される。 が正定値行列である とき、 は であるので、 の 正平方根 とも呼ばれる。 M {\displaystyle M} B {\displaystyle B} B {\displaystyle B} B ∗ = B {\displaystyle B^{*}=B} M = B B 。 {\displaystyle M=BB.} B {\displaystyle B} M 、 {\displaystyle M,} B = M 1 2 。 {\displaystyle B=M^{\frac {1}{2}}.} M {\displaystyle M} M 1 2 、 {\displaystyle M^{\frac {1}{2}},} M 。 {\displaystyle M.}
非負の平方根を他の分解と混同しないでください。 著者によっては 、このような分解すべてに対して、または特にコレスキー分解 、または形式の分解すべてに対して 、 平方根 という 名前 を使用する人もいます
が、非負の平方根に対してのみこの名前を使用する人もいます。 M = B ∗ B 。 {\displaystyle M=B^{*}B.} M 1 2 {\displaystyle M^{\frac {1}{2}}} M = B B ; {\displaystyle M=BB;}
もし そうなら M ≻ 北 ≻ 0 {\displaystyle M\succ N\succ 0} M 1 2 ≻ 北 1 2 ≻ 0。 {\displaystyle M^{\frac {1}{2}}\succ N^{\frac {1}{2}}\succ 0.}
コレスキー分解 エルミート正半正定値行列 は と表すことができます。 ここで は非負の対角成分を持つ下三角行列です(同様に は 上三角行列です)。これは コレスキー分解 です。 が正定値行列である場合 、 の対角成分は 正であり、コレスキー分解は一意です。逆に が 非負の対角成分を持つ下三角行列である場合、 は半正定値行列です。コレスキー分解は、特に効率的な数値計算に役立ちます。密接に関連する分解は LDL 分解 です。 ここで は対角成分を持ち、は 下一三角 です 。 M {\displaystyle M} M = L L ∗ 、 {\displaystyle M=LL^{*},} L {\displaystyle L} M = B ∗ B {\displaystyle M=B^{*}B} B = L ∗ {\displaystyle B=L^{*}} M {\displaystyle M} L {\displaystyle L} L {\displaystyle L} L L ∗ {\displaystyle LL^{*}} M = L D L ∗ 、 {\displaystyle M=LDL^{*},} D {\displaystyle D} L {\displaystyle L}
ウィリアムソン定理 任意の 正定値エルミート実行列は、 シンプレクティック(実)行列を介して対角化できます。より正確には、 ウィリアムソンの定理は、 を満たすシンプレクティック および対角実正行列 の存在を保証します 。 2 n × 2 n {\displaystyle 2n\times 2n} M {\displaystyle M} S ∈ S p ( 2 n 、 R ) {\displaystyle S\in \mathbf {Sp} (2n,\mathbb {R} )} D ∈ R n × n {\displaystyle D\in \mathbb {R} ^{n\times n}} S M S T = D ⊕ D {\displaystyle SMS^{T}=D\oplus D}
その他の特徴 を実対称行列 とし 、 を 「単位球」と定義すると、 次の式が得られる。 M {\displaystyle M} n × n {\displaystyle n\times n} B 1 ( M ) ≡ { × ∈ R n : × T M × ≤ 1 } {\displaystyle B_{1}(M)\equiv \{\mathbf {x} \in \mathbb {R} ^{n}:\mathbf {x} ^{\mathsf {T}}M\mathbf {x} \leq 1\}} M 。 {\displaystyle M.}
B 1 ( v v T ) {\displaystyle B_{1}(\mathbf {v} \mathbf {v} ^{\mathsf {T}})} 間に挟まれた固体の板 ± { わ : ⟨ わ 、 v ⟩ = 1 } 。 {\displaystyle \pm \{\mathbf {w} :\langle \mathbf {w} ,\mathbf {v} \rangle =1\}。} M ⪰ 0 {\displaystyle M\succeq 0} が楕円体、または楕円柱である 場合に限ります。 B 1 ( M ) {\displaystyle B_{1}(M)} M ≻ 0 {\displaystyle M\succ 0} が有界である場合 、つまり楕円体である場合に限ります。 B 1 ( M ) {\displaystyle B_{1}(M)} もし、 もし 、もし、もし 、もし、もし、もし 北 ≻ 0 、 {\displaystyle N\succ 0,} M ⪰ 北 {\displaystyle M\succeq N} B 1 ( M ) ⊆ B 1 ( 北 ) ; {\displaystyle B_{1}(M)\subseteq B_{1}(N);} M ≻ 北 {\displaystyle M\succ N} B 1 ( M ) ⊆ 整数 ( B 1 ( 北 ) ) 。 {\displaystyle B_{1}(M)\subseteq \operatorname {int} {\bigl (}B_{1}(N){\bigr )}.} ならば 、 すべてに対して が 成り立つこと、そしてその場合のみ成り立つ。 楕円体の極双対も、同じ主軸を持ち長さが逆である楕円体であるので、 次の関係が成り立つ。すなわち、 が正定値ならば、 すべてに対して が成り立つこと、 そしてその場合のみ成り立つ。 北 ≻ 0 、 {\displaystyle N\succ 0,} M ⪰ v v T v T 北 v {\displaystyle M\succeq {\frac {\mathbf {v} \mathbf {v} ^{\mathsf {T}}}{\mathbf {v} ^{\mathsf {T}}N\mathbf {v} }}} v ≠ 0 {\displaystyle v\neq 0} B 1 ( M ) ⊂ ⋂ v T 北 v = 1 B 1 ( v v T ) 。 {\textstyle B_{1}(M)\subset \bigcap _{\mathbf {v} ^{\mathsf {T}}N\mathbf {v} =1}B_{1}(\mathbf {v} \mathbf {v} ^{\mathsf {T}})。 B 1 ( 北 − 1 ) = ⋂ v T 北 v = 1 B 1 ( v v T ) = ⋂ v T 北 v = 1 { わ : | ⟨ わ 、 v ⟩ | ≤ 1 } 。 {\displaystyle B_{1}(N^{-1})=\bigcap _{\mathbf {v} ^{\mathsf {T}}N\mathbf {v} =1}B_{1}(\mathbf {v} \mathbf {v} ^{\mathsf {T}})=\bigcap _{\mathbf {v} ^{\mathsf {T}}N\mathbf {v} =1}\{\mathbf {w} :|\langle \mathbf {w} ,\mathbf {v} \rangle |\leq 1\}.} 北 {\displaystyle N} M ⪰ v v T v T 北 v {\displaystyle M\succeq {\frac {\mathbf {v} \mathbf {v} ^{\mathsf {T}}}{\mathbf {v} ^{\mathsf {T}}N\mathbf {v} }}} v ≠ 0 {\displaystyle \mathbf {v} \neq \mathbf {0} } M ⪰ 北 − 1 。 {\displaystyle M\succeq N^{-1}.} をエルミート行列 とする 。 以下の性質は 正定値行列であることと同等である。 M {\displaystyle M} n × n {\displaystyle n\times n} M {\displaystyle M}
関連するセクスティリニア形式は内積である によって定義される セクスティ 線形形式は 、 から への 関数で、 すべての および において と なるものです。 ここで は の共役転置です。 任意の複素行列に対して、 この形式は で線形 、 で半線形です。 したがって、 が の 内積 となるのは 、 が実数かつ正である 場合に限ります。すべての非ゼロに対して となるのは、 が正定値である場合に限ります 。(実際、 のすべての内積は、 エルミート正定値行列からこのように生じます。) M {\displaystyle M} ⟨ ⋅ 、 ⋅ ⟩ {\displaystyle \langle \cdot ,\cdot \rangle } C n × C n {\displaystyle \mathbb {C} ^{n}\times \mathbb {C} ^{n}} C n {\displaystyle \mathbb {C} ^{n}} ⟨ × 、 y ⟩ ≡ y ∗ M × {\displaystyle \langle \mathbf {x} ,\mathbf {y} \rangle \equiv \mathbf {y} ^{*}M\mathbf {x} } × {\displaystyle \mathbf {x} } y {\displaystyle \mathbf {y} } C n 、 {\displaystyle \mathbb {C} ^{n},} y ∗ {\displaystyle \mathbf {y} ^{*}} y 。 {\displaystyle \mathbf {y} .} M 、 {\displaystyle M,} × {\displaystyle x} y 。 {\displaystyle \mathbf {y} .} C n {\displaystyle \mathbb {C} ^{n}} ⟨ z 、 z ⟩ {\displaystyle \langle \mathbf {z} ,\mathbf {z} \rangle } z ; {\displaystyle \mathbf {z} ;} M {\displaystyle M} C n {\displaystyle \mathbb {C} ^{n}} 主要マイナーはすべてプラス 行列の k 番目 の 主要小行列 式は、その左上部分行列の 行列式 です 。行列が正定値行列であるためには、これらすべての行列式が正である必要があります。この条件は シルベスターの基準 として知られており、対称実行列の正定値を効率的にテストできます。つまり、 ガウス消去 法の最初の部分と同様に、 基本的な行操作 を使用して、行列 を上三角行列に簡約し、 ピボット 処理中に行列式の符号が保持されるように注意します 。三角行列の k 番目の主要小行列式は、行 までの対角要素の積であるため、シルベスターの基準は、対角要素がすべて正であるかどうかを確認することと同じです。この条件は 、三角行列の新しい行が取得される たびに確認できます。 M {\displaystyle M} け × け {\displaystyle k\times k} け 、 {\displaystyle k,} け {\displaystyle k} 半正定値行列が正定値行列となるのは、それが 逆行列で ある場合に限ります。 [7] 半正定値 行列が負定値行列となるのは 、それが正定値行列である場合に限ります。 M {\displaystyle M} − M {\displaystyle -M}
実行 行列に関連付けられた (純粋な) 二次形式は 、すべての に対してとなる 関数です。これ は、両面積で非対称部分は すべてゼロになるため 、 と置き換えることで対称であると想定できます。 n × n {\displaystyle n\times n} M {\displaystyle M} 質問 : R n → R {\displaystyle Q:\mathbb {R} ^{n}\to \mathbb {R} } 質問 ( × ) = × T M × {\displaystyle Q(\mathbf {x} )=\mathbf {x} ^{\mathsf {T}}M\mathbf {x} } × 。 {\displaystyle \mathbf {x} .} M {\displaystyle M} 1 2 ( M + M T ) 、 {\displaystyle {\tfrac {1}{2}}\left(M+M^{\mathsf {T}}\right),}
対称行列は、その二次形式が 厳密に凸関数 である場合に限り、正定値行列となります 。 M {\displaystyle M}
より一般的には、 から への任意の 二次関数は 、対称 行列 、 実数 ベクトル 、 実定数 で表すことができます。この 場合、これは放物線であり、の 場合と同様に、 R n {\displaystyle \mathbb {R} ^{n}} R {\displaystyle \mathbb {R} } × T M × + b T × + c {\displaystyle \mathbf {x} ^{\mathsf {T}}M\mathbf {x} +\mathbf {b} ^{\mathsf {T}}\mathbf {x} +c} M {\displaystyle M} n × n {\displaystyle n\times n} b {\displaystyle \mathbf {b} } c {\displaystyle c} n = 1 {\displaystyle n=1} n = 1 {\displaystyle n=1}
定理: この二次関数は厳密に凸関数であるため、正定値関数である場合に限り、唯一の有限大域最小値を持ちます 。 M {\displaystyle M}
証明: が正定値関数である場合 、関数は厳密に凸である。 関数が厳密に凸であるため、その勾配は唯一の点においてゼロとなり、その点が大域的最小値となる。 が正定値関数でない場合、 関数が直線または下向き放物線となるような ベクトルが存在する ため 、厳密に凸ではなく、大域的最小値を持たない。 M {\displaystyle M} M − 1 b 、 {\displaystyle M^{-1}\mathbf {b} ,} M {\displaystyle M} v {\displaystyle \mathbf {v} } v T M v ≤ 0 、 {\displaystyle \mathbf {v} ^{\mathsf {T}}M\mathbf {v} \leq 0,} f ( t ) ≡ ( t v ) T M ( t v ) + b T ( t v ) + c {\displaystyle f(t)\equiv (t\mathbf {v} )^{\mathsf {T}}M(t\mathbf {v} )+b^{\mathsf {T}}(t\mathbf {v} )+c}
このため、正定値行列は 最適化 問題において重要な役割を果たします。
同時対角化 対称行列と、対称かつ正定値行列である別の行列は、 同時に対角化すること ができます。同時対角化は必ずしも 相似変換 によって行われるわけではありませんが、これは可能です。この結果は、3つ以上の行列の場合には適用されません。この節では実数の場合について記述します。複素数の場合への拡張は直ちに可能です。
を対称行列と 対称正定値行列とする。一般化固有値方程式を次のように書く。 ここで 、正規化、すなわち を課す。ここで、 コレスキー分解を 用いて の逆行列を と書く。を掛け て とする と 、 となり、 これは次のように書き直すことができる。 ここ で、 の操作により が 得られる。 ここで 、 は一般化固有ベクトルを列として持つ行列であり、 は一般化固有値の対角行列である。ここで を前置すると、 最終結果は となり、 ただし 、これはもはや内積 に関して直交対角化ではないことに注意する。 実際には、 [8] によって誘導される内積に関して対角化している。 M {\displaystyle M} 北 {\displaystyle N} ( M − λ 北 ) × = 0 {\displaystyle \left(M-\lambda N\right)\mathbf {x} =0} × {\displaystyle \mathbf {x} } × T 北 × = 1. {\displaystyle \mathbf {x} ^{\mathsf {T}}N\mathbf {x} =1.} 北 {\displaystyle N} 質問 T 質問 。 {\displaystyle Q^{\mathsf {T}}Q.} 質問 {\displaystyle Q} × = 質問 T y 、 {\displaystyle \mathbf {x} =Q^{\mathsf {T}}\mathbf {y} ,} 質問 ( M − λ 北 ) 質問 T y = 0 、 {\displaystyle Q\left(M-\lambda N\right)Q^{\mathsf {T}}\mathbf {y} =0,} ( 質問 M 質問 T ) y = λ y {\displaystyle \left(QMQ^{\mathsf {T}}\right)\mathbf {y} =\lambda \mathbf {y} } y T y = 1. {\displaystyle \mathbf {y} ^{\mathsf {T}}\mathbf {y} =1.} M X = 北 X Λ {\displaystyle MX=NX\Lambda } X {\displaystyle X} Λ {\displaystyle \Lambda } X T {\displaystyle X^{\mathsf {T}}} X T M X = Λ {\displaystyle X^{\mathsf {T}}MX=\Lambda } X T 北 X = 私 、 {\displaystyle X^{\mathsf {T}}NX=I,} y T y = 1. {\displaystyle \mathbf {y} ^{\mathsf {T}}\mathbf {y} =1.} M {\displaystyle M} 北 。 {\displaystyle N.}
この結果は、記事「対角化可能行列」 における同時対角化に関する記述と矛盾しないことに注意してください。記事「対角化可能行列 」では、相似変換による同時対角化について言及されています。ここでの私たちの結果は、2つの二次形式の同時対角化に近く、一方の形式をもう一方の形式に与えられた条件下で最適化する際に役立ちます。
プロパティ
誘導半順序 任意の正方行列について、 すなわち が半正定値行列である とき 、 と書きます。これは 、すべての正方行列の集合における 半順序を定義します。同様に、厳密な半順序を定義することもできます。この順序は Loewner順序 と呼ばれます 。 M 、 {\displaystyle M,} 北 {\displaystyle N} M ≥ 北 {\displaystyle M\geq N} M − 北 ≥ 0 {\displaystyle MN\geq 0} M − 北 {\displaystyle MN} M > 北 。 {\displaystyle M>N.}
正定値行列の逆行列 すべての正定値行列は 逆行列 であり、その逆行列も正定値である。 [9] ならば [10] さらに、 最小最大定理 により 、の k 番目に大きい固有値は、の k 番目に大きい固有値 以上である。 M ≥ 北 > 0 {\displaystyle M\geq N>0} 北 − 1 ≥ M − 1 > 0。 {\displaystyle N^{-1}\geq M^{-1}>0.} M {\displaystyle M} 北 。 {\displaystyle N.}
スケーリング が正定値で 実数であれ ば、 は正定値である。 [11] M {\displaystyle M} r > 0 {\displaystyle r>0} r M {\displaystyle rM}
追加 とが正定値で あれば 、和 も正定値である。 [11] M {\displaystyle M} 北 {\displaystyle N} M + 北 {\displaystyle M+N} および が半正定値である 場合 、和 も半正定値です。 M {\displaystyle M} 北 {\displaystyle N} M + 北 {\displaystyle M+N} が正定値で が 半正定値である 場合、その合計 も正定値になります。 M {\displaystyle M} 北 {\displaystyle N} M + 北 {\displaystyle M+N}
乗算 と が 正定値であれ ば、 積 も正定値です。 で あれば 、 も正定値です。 M {\displaystyle M} 北 {\displaystyle N} M 北 M {\displaystyle MNM} 北 M 北 {\displaystyle NMN} M 北 = 北 M 、 {\displaystyle MN=NM,} M 北 {\displaystyle MN} が半正定値行列であれ ば、 任意の(おそらく長方形の)行列に対して半正定値行列となる。 が正定値行列であり 、かつフル列ランクを持つ 場合、 半正定値行列となる。 [12] M {\displaystyle M} あ ∗ M あ {\displaystyle A^{*}MA} あ 。 {\displaystyle A.} M {\displaystyle M} あ {\displaystyle A} あ ∗ M あ {\displaystyle A^{*}MA}
トレース 半正定値行列の対角成分は 実数かつ非負である。したがって、 トレースは 、 さらに [13] 、すべての主部分行列(特に2行2列)が半正定値行列であるため 、 メートル 私 私 {\displaystyle m_{ii}} tr ( M ) ≥ 0。 {\displaystyle \operatorname {tr} (M)\geq 0.} | メートル 私 j | ≤ メートル 私 私 メートル j j た 私 、 j {\displaystyle \left|m_{ij}\right|\leq {\sqrt {m_{ii}m_{jj}}}\quad \forall i,j} n ≥ 1 、 {\displaystyle n\geq 1,} 最大 私 、 j | メートル 私 j | ≤ 最大 私 メートル 私 私 {\displaystyle \max_{i,j}\left|m_{ij}\right|\leq \max_{i}m_{ii}}
エルミート 行列 が正定値行列であるとは、次のトレース不等式を満たす場合である。 [14] n × n {\displaystyle n\times n} M {\displaystyle M} tr ( M ) > 0 1つの n d ( tr ( M ) ) 2 tr ( M 2 ) > n − 1. {\displaystyle \operatorname {tr} (M)>0\quad \mathrm {および} \quad {\frac {(\operatorname {tr} (M))^{2}}{\operatorname {tr} (M^{2})}}>n-1.}
もう 1 つの重要な結果は、任意 の 半正定値行列と半正定値行列に対して、 次の式が成り立つことです。 行列は半正定値であるため、非負の固有値を持ち、その和である トレース も非負になります。 M {\displaystyle M} 北 {\displaystyle N} tr ( M 北 ) ≥ 0。 {\displaystyle \operatorname {tr} (MN)\geq 0.} tr ( M 北 ) = tr ( M 1 2 北 M 1 2 ) 。 {\displaystyle \operatorname {tr} (MN)=\operatorname {tr} (M^{\frac {1}{2}}NM^{\frac {1}{2}}).} M 1 2 北 M 1 2 {\displaystyle M^{\frac {1}{2}}NM^{\frac {1}{2}}}
アダマール積 が必ずしも半正定値ではない としても 、 アダマール積 は、 (この結果はしばしば シュアー積定理 と呼ばれる)。 [15] M 、 北 ≥ 0 、 {\displaystyle M,N\geq 0,} M 北 {\displaystyle MN} M ∘ 北 ≥ 0 {\displaystyle M\circ N\geq 0}
2つの半正定値行列のアダマール積に関しては、 2つの注目すべき不等式があります。 M = ( メートル 私 j ) ≥ 0 、 {\displaystyle M=(m_{ij})\geq 0,} 北 ≥ 0 、 {\displaystyle N\geq 0,}
オッペンハイムの不等式: [16] 詳細 ( M ∘ 北 ) ≥ 詳細 ( 北 ) ∏ 私 メートル 私 私 。 {\displaystyle \det(M\circ N)\geq \det(N)\prod \nolimits _{i}m_{ii}.} 詳細 ( M ∘ 北 ) ≥ 詳細 ( M ) 詳細 ( 北 ) 。 {\displaystyle \det(M\circ N)\geq \det(M)\det(N).} [17]
クロネッカー積 必ずしも半正定値ではない が 、 クロネッカー積 M 、 北 ≥ 0 、 {\displaystyle M,N\geq 0,} M 北 {\displaystyle MN} M ⊗ 北 ≥ 0。 {\displaystyle M\otimes N\geq 0.}
フロベニウス積 が 必ずしも半正定値ではない 場合、 フロベニウスの内積 (ランカスター・ティスメネツキー著 『行列理論 』218 ページ)。 M 、 北 ≥ 0 、 {\displaystyle M,N\geq 0,} M 北 {\displaystyle MN} M : 北 ≥ 0 {\displaystyle M:N\geq 0}
凸状性 半正定値対称行列の集合は 凸行列で ある。つまり、 と が半正定値行列であるならば 、 0 から 1 までの 任意の に対して も半正定値行列となる。任意のベクトル に対して 、 M {\displaystyle M} 北 {\displaystyle N} α {\displaystyle \alpha} α M + ( 1 − α ) 北 {\displaystyle \alpha M+\left(1-\alpha \right)N} × {\displaystyle \mathbf {x} } × T ( α M + ( 1 − α ) 北 ) × = α × T M × + ( 1 − α ) × T 北 × ≥ 0。 {\displaystyle \mathbf {x} ^{\mathsf {T}}\left(\alpha M+\left(1-\alpha \right)N\right)\mathbf {x} =\alpha \mathbf {x} ^{\mathsf {T}}M\mathbf {x} +(1-\alpha )\mathbf {x} ^{\mathsf {T}}N\mathbf {x} \geq 0.}
この特性により、 半正定値計画 問題がグローバルに最適な解に収束することが保証されます。
コサインとの関係 行列の正定値は、 任意のベクトル とその像の間の 角度 が常に あ {\displaystyle A} θ {\displaystyle \theta} × {\displaystyle \mathbf {x} } あ × {\displaystyle A\mathbf {x} } − π / 2 < θ < + π / 2 : {\displaystyle -\pi /2<\theta <+\pi /2:}
コス θ = × T あ × ‖ × ‖ ‖ あ × ‖ = ⟨ × 、 あ × ⟩ ‖ × ‖ ‖ あ × ‖ 、 θ = θ ( × 、 あ × ) ≡ ( × 、 あ × ) ^ ≡ {\displaystyle \cos \theta ={\frac {\mathbf {x} ^{\mathsf {T}}A\mathbf {x} }{\lVert \mathbf {x} \rVert \lVert A\mathbf {x} \rVert }}={\frac {\langle \mathbf {x} ,A\mathbf {x} \rangle }{\lVert \mathbf {x} \rVert \lVert A\mathbf {x} \rVert }},\theta =\theta (\mathbf {x} ,A\mathbf {x} )\equiv {\widehat {\left(\mathbf {x} ,A\mathbf {x} \right)}}\equiv } と 間の角度 × {\displaystyle \mathbf {x} } あ × 。 {\displaystyle A\mathbf {x} .}
その他の特性 が対称 テプリッツ行列 である場合 、つまり、エントリが 絶対インデックス差の関数として与えられ、 厳密 な 不等式 が成り立つ場合、は 厳密に 正定値 です。 M {\displaystyle M} メートル 私 j {\displaystyle m_{ij}} メートル 私 j = h ( | 私 − j | ) 、 {\displaystyle m_{ij}=h(|ij|),} ∑ j ≠ 0 | h ( j ) | < h ( 0 ) {\textstyle \sum _{j\neq 0}\left|h(j)\right|<h(0)} M {\displaystyle M} とエルミート とする 。 (それぞれ )ならば (それぞれ )。 [18] M > 0 {\displaystyle M>0} 北 {\displaystyle N} M 北 + 北 M ≥ 0 {\displaystyle MN+NM\geq 0} M 北 + 北 M > 0 {\displaystyle MN+NM>0} 北 ≥ 0 {\displaystyle N\geq 0} 北 > 0 {\displaystyle N>0} が実数の 場合、 が存在し、 ここで は 単位行列 です 。 M > 0 {\displaystyle M>0} δ > 0 {\displaystyle \delta >0} M > δ 私 、 {\displaystyle M>\delta I,} 私 {\displaystyle I} がリーディングマイナーを表す 場合 、は LU 分解中の k 番目のピボット です 。 M け {\displaystyle M_{k}} け × け {\displaystyle k\times k} 詳細 ( M け ) / 詳細 ( M け − 1 ) {\displaystyle \det \left(M_{k}\right)/\det \left(M_{k-1}\right)} が奇数のとき、その k 次の 主小 行列項が負であり、が偶数のとき、正である場合、 その行列は負定値行列です 。 け {\displaystyle k} け {\displaystyle k} が実正定値行列である 場合、任意の ベクトルに対して、 M {\displaystyle M} メートル {\displaystyle m} v 、 {\displaystyle \mathbf {v} ,} v T M v ≥ メートル ‖ v ‖ 2 2 。 {\displaystyle \mathbf {v}^{\mathsf {T}}M\mathbf {v} \geq m\|\mathbf {v} \|_{2}^{2}.} エルミート行列が半正定値行列であるための必要十分条件は、そのすべての主小行列が非負であることです。しかし、 0 と-1 の要素を持つ対角行列で確認されているように、主小行列のみを考慮するだけでは不十分です 。
ブロック行列と部分行列 正の行列は ブロック によって定義されることもあります 。 2 n × 2 n {\displaystyle 2n\times 2n} M = [ あ B C D ] {\displaystyle M={\begin{bmatrix}A&B\\C&D\end{bmatrix}}}
ここで各ブロックは 正値条件を適用することにより、 および はエルミートであることが直ちに示され、 n × n 、 {\displaystyle n\times n,} あ {\displaystyle A} D {\displaystyle D} C = B ∗ 。 {\displaystyle C=B^{*}.}
私たちはすべての複雑なもの 、特に Then
について それを持っています z ∗ M z ≥ 0 {\displaystyle \mathbf {z} ^{*}M\mathbf {z} \geq 0} z 、 {\displaystyle \mathbf {z} ,} z = [ v 、 0 ] T 。 {\displaystyle \mathbf {z} =[\mathbf {v},0]^{\mathsf {T}}.} [ v ∗ 0 ] [ あ B B ∗ D ] [ v 0 ] = v ∗ あ v ≥ 0。 {\displaystyle {\begin{bmatrix}\mathbf {v} ^{*}&0\end{bmatrix}}{\begin{bmatrix}A&B\\B^{*}&D\end{bmatrix}}{\begin{bmatrix}\mathbf {v} \\0\end{bmatrix}}=\mathbf {v} ^{*}A\mathbf {v} \geq 0.}
同様の議論は にも適用でき、したがって と は両方 とも正定値行列である と結論付けられます 。この議論を拡張すると、 の任意の 主部分行列は それ自体が正定値行列であることが示されます。 D 、 {\displaystyle D,} あ {\displaystyle A} D {\displaystyle D} M {\displaystyle M}
逆の結果は、たとえば、シュアー補集合 を使用して、ブロックに対するより強い条件で証明できます 。
局所的極値 実変数 の 一般 二次形式は 常に次のように表すことができます。 ここで 、はそれらの変数を含む列ベクトルであり、 は対称実行列です。したがって、行列が正定値であるということは、がゼロの ときに唯一の最小値(ゼロ)を持ち 、それ以外の場合には正であることを意味します。 f ( × ) {\displaystyle f(\mathbf {x} )} n {\displaystyle n} × 1 、 … 、 × n {\displaystyle x_{1},\ldots,x_{n}} × T M × {\displaystyle \mathbf {x} ^{\mathsf {T}}M\mathbf {x} } × {\displaystyle \mathbf {x} } M {\displaystyle M} f {\displaystyle f} × {\displaystyle \mathbf {x} } × 。 {\displaystyle \mathbf {x} .}
より一般的には、実変数 の 二回微分可能な実関数は、 その 勾配 がゼロで、その ヘッセ行列 (すべての二階微分行列)がその点で半正定値であるとき、引数において局所最小値を持つ。負定値行列と半正定値行列についても同様のことが言える。 f {\displaystyle f} n {\displaystyle n} × 1 、 … 、 × n {\displaystyle x_{1},\ldots,x_{n}}
共分散 統計学 において 、 多変量確率分布 の 共分散行列 は常に半正定値行列である。また、ある変数が他の変数の正確な線形関数でない限り、半正定値行列である。逆に、すべての半正定値行列は、何らかの多変量分布の共分散行列である。
非エルミート正方行列への拡張 正定値の定義は、任意の複素行列 (例えば実非対称)を、 すべての非零複素ベクトルに対して が成り立つとき、正定値であると定義することで一般化できる。 ここで は 複素数 の実部を表す [19] 行列が正定値であるかどうかはエルミート部によってのみ 決定され、これは上記の狭義の意味で評価される。同様に、 と が実数であるとき、 すべての非零実ベクトルに対して が成り立つ ことと、対称部が狭義の意味で正定値であるとき、かつその場合に限る。 が の転置に対して無反応である ことはすぐに明らかである。 M {\displaystyle M} R e { z ∗ M z } > 0 {\displaystyle {\mathcal {R_{e}}}\left\{\mathbf {z} ^{*}M\mathbf {z} \right\}>0} z 、 {\displaystyle \mathbf {z} ,} R e { c } {\displaystyle {\mathcal {R_{e}}}\{c\}} c 。 {\displaystyle c.} 1 2 ( M + M ∗ ) {\textstyle {\frac {1}{2}}\left(M+M^{*}\right)} × {\displaystyle \mathbf {x} } M {\displaystyle M} × T M × > 0 {\displaystyle \mathbf {x} ^{\mathsf {T}}M\mathbf {x} >0} × {\displaystyle \mathbf {x} } 1 2 ( M + M T ) {\textstyle {\frac {1}{2}}\left(M+M^{\mathsf {T}}\right)} × T M × = ∑ 私 j × 私 M 私 j × j {\textstyle \mathbf {x} ^{\mathsf {T}}M\mathbf {x} =\sum _{ij}x_{i}M_{ij}x_{j}} M 。 {\displaystyle M.}
正の固有値のみを持つ非対称実行列は、負の固有値を持つ対称部分を持つ場合があります。その場合、行列は正(半)定値行列にはなりません。例えば、行列は 正の固有値1と7を持ちますが、 と選択します 。 M = [ 4 9 1 4 ] {\textstyle M=\left[{\begin{smallmatrix}4&9\\1&4\end{smallmatrix}}\right]} × T M × = − 2 {\displaystyle \mathbf {x} ^{\mathsf {T}}M\mathbf {x} =-2} × = [ − 1 1 ] {\displaystyle \mathbf {x} =\left[{\begin{smallmatrix}-1\\1\end{smallmatrix}}\right]}
まとめると、実数と複素数の場合の違いは、複素ヒルベルト空間上の 有界正作用素は必然的にエルミート、つまり自己随伴であるという点です。この一般的な主張は 分極恒等式 を用いて論じることができます 。しかし、実数の場合、これはもはや成り立ちません。
アプリケーション
熱伝導率マトリックス 熱伝導に関するフーリエの法則は、熱流束を 温度勾配で 表し、異方性媒体では と表されます。 ここで は 熱伝導率 行列です 。フーリエの法則に負の値が挿入されているのは、熱は常に高温から低温へ流れるという期待を反映しているためです。言い換えれば、温度勾配は 常に低温から高温に向かうため、熱流束は との内積が負になると考えられます。 したがって、 となり ます。フーリエの法則を に代入すると、この期待は となり 、伝導率行列は正定値行列となるはずです。通常は は対称ですが、 熱ホール効果 のように磁場が存在する場合は非対称になります 。 q {\displaystyle \mathbf {q} } グラム = ∇ T {\displaystyle \mathbf {g} =\nabla T} q = − K グラム 、 {\displaystyle \mathbf {q} =-K\mathbf {g} ,} K {\displaystyle K} グラム {\displaystyle \mathbf {g} } q {\displaystyle \mathbf {q} } グラム {\displaystyle \mathbf {g} } q T グラム < 0。 {\displaystyle \mathbf {q} ^{\mathsf {T}}\mathbf {g} <0.} グラム T K グラム > 0 、 {\displaystyle \mathbf {g} ^{\mathsf {T}}K\mathbf {g} >0,} K {\displaystyle K}
より一般的には、熱力学では、熱と粒子の流れは、 オンサガーの逆関係 で説明される完全に結合したシステムであり、エントロピー生成が非負になるためには、結合行列は半正定値(非対称の場合もある)である必要があります。
参照
参考文献 ^ van den Bos, Adriaan (2007年3月). 「付録C: 半正定値行列と正定値行列」. 科学者とエンジニアのためのパラメータ推定 (.pdf) (オンライン版). John Wiley & Sons. pp. 259– 263. doi :10.1002/9780470173862. ISBN 978-047-017386-2 。 印刷版 ISBN 9780470147818 ^ Boyd, Stephen; Vandenberghe, Lieven (2004年3月8日). 凸最適化 . Cambridge University Press. doi :10.1017/cbo9780511804441. ISBN 978-0-521-83378-3 。 ^ ホーン&ジョンソン(2013)、p.440、定理7.2.7 ^ ホーン&ジョンソン(2013)、p.441、定理7.2.10 ^ ホーン&ジョンソン(2013)、p.452、定理7.3.11 ^ Horn & Johnson (2013)、p. 439、定理7.2.6と け = 2 {\displaystyle k=2} ^ ホーン&ジョンソン(2013)、p.431、補論7.1.7 ^ ホーン&ジョンソン(2013)、p.485、定理7.6.1 ^ ホーン&ジョンソン(2013)、p.438、定理7.2.1 ^ ホーン&ジョンソン(2013)、p.495、補論7.7.4(a) ^ ab Horn & Johnson (2013)、p. 430、観察7.1.3 ^ ホーン&ジョンソン(2013)、p.431、観察7.1.8 ^ ホーン&ジョンソン(2013)、430ページ ^ Wolkowicz, Henry; Styan, George PH (1980). 「トレースを用いた固有値の境界」. 線形代数とその応用 . 29 (29). Elsevier: 471– 506. doi :10.1016/0024-3795(80)90258-X. ^ ホーン&ジョンソン(2013)、p.479、定理7.5.3 ^ ホーン&ジョンソン(2013)、p.509、定理7.8.16 ^ Styan, GP (1973). 「アダマール積と多変量統計解析」. 線形代数とその応用 . 6 : 217–240 . doi :10.1016/0024-3795(73)90023-2. 、系3.6、p. 227 ^ Bhatia, Rajendra (2007). 正定値行列 . プリンストン大学出版局, ニュージャージー州, 8頁. ISBN 978-0-691-12918-1 。 ^ Weisstein, Eric W. 「正定値行列」. MathWorld . Wolfram Research . 2012年 7月26日 閲覧 。
出典
外部リンク 「正定値形式」、 数学百科事典 、 EMS Press 、2001 [1994] 「正定値行列」。Wolfram MathWorld。Wolfram Research。