明確なマトリックス

数学において、実数要素 を持つ対称行列が正定値行列であるとは、実数がすべての非零の実列ベクトルに対して正であるときである。ここで、は[ 1 ]転置行ベクトルである。 より一般的には、エルミート行列(つまり、その共役転置に等しい複素行列)が正定値行列であるとは、実数がすべての非零の複素列ベクトルに対して正であるときである。ここで、はの共役転置を表す。M{\displaystyle M}×TM×{\displaystyle \mathbf {x} ^{\mathsf {T}}M\mathbf {x} }×{\displaystyle \mathbf {x} ,}×T{\displaystyle \mathbf {x} ^{\mathsf {T}}}×{\displaystyle \mathbf {x} .}zMz{\displaystyle \mathbf {z} ^{*}M\mathbf {z} }z{\displaystyle \mathbf {z} ,}z{\displaystyle \mathbf {z} ^{*}}z{\displaystyle \mathbf {z} .}

半正定値行列も同様に定義されますが、スカラーとが正またはゼロ(つまり非負)である必要があります。負定値行列半負定値行列も同様に定義されます。半正定値でも半負定値でもない行列は、不定値行列と呼ばれることがあります。 ×TM×{\displaystyle \mathbf {x} ^{\mathsf {T}}M\mathbf {x} }zMz{\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}×TM×>0{\displaystyle \mathbf {x} ^{\mathsf {T}}M\mathbf {x} >0}×{\displaystyle \mathbf {x} }Rn{\displaystyle \mathbb {R} ^{n}.}

M 正定値×TM×>0 すべての人のために ×Rn{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}×TM×0{\displaystyle \mathbf {x} ^{\mathsf {T}}M\mathbf {x} \geq 0}×{\displaystyle \mathbf {x} }Rn{\displaystyle \mathbb {R} ^{n}.}

M 半正定値×TM×0 すべての人のために ×Rn{\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}×TM×<0{\displaystyle \mathbf {x} ^{\mathsf {T}}M\mathbf {x} <0}×{\displaystyle \mathbf {x} }Rn{\displaystyle \mathbb {R} ^{n}.}

M 負定値×TM×<0 すべての人のために ×Rn{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}×TM×0{\displaystyle \mathbf {x} ^{\mathsf {T}}M\mathbf {x} \leq 0}×{\displaystyle \mathbf {x} }Rn{\displaystyle \mathbb {R} ^{n}.}

M 負の半定値×TM×0 すべての人のために ×Rn{\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}

複素行列の定義

以下の定義はすべてこの用語を含んでいます。これは、任意のエルミート正方行列に対して常に実数であることに注意してください。zMz{\displaystyle \mathbf {z} ^{*}M\mathbf {z} .}M{\displaystyle M.}

エルミート複素行列は、すべての非ゼロの値が正定値であるとは、正式には、 n×n{\displaystyle n\times n}M{\displaystyle M}zMz>0{\displaystyle \mathbf {z} ^{*}M\mathbf {z} >0}z{\displaystyle \mathbf {z} }Cn{\displaystyle \mathbb {C} ^{n}.}

M 正定値zMz>0 すべての人のために zCn{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}zMz0{\displaystyle \mathbf {z} ^{*}M\mathbf {z} \geq 0}z{\displaystyle \mathbf {z} }Cn{\displaystyle \mathbb {C} ^{n}.}

M 半正定値zMz0 すべての人のために zCn{\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}zMz<0{\displaystyle \mathbf {z} ^{*}M\mathbf {z} <0}z{\displaystyle \mathbf {z} }Cn{\displaystyle \mathbb {C} ^{n}.}

M 負定値zMz<0 すべての人のために zCn{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}zMz0{\displaystyle \mathbf {z} ^{*}M\mathbf {z} \leq 0}z{\displaystyle \mathbf {z} }Cn{\displaystyle \mathbb {C} ^{n}.}

M 負の半定値zMz0 すべての人のために zCn{\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}zMz{\displaystyle \mathbf {z} ^{*}M\mathbf {z} }z{\displaystyle \mathbf {z} .}M{\displaystyle M}zMz{\displaystyle \mathbf {z} ^{*}M\mathbf {z} }zMz{\displaystyle \mathbf {z} ^{*}M^{*}\mathbf {z} }z{\displaystyle \mathbf {z} ,}MM{\displaystyle M=M^{*}.}

この定義によれば、正定値実数行列はエルミート行列であり、したがって対称行列である。また、すべての非零列ベクトルに対しては正となる 。しかし、最後の条件だけでは、正定値行列であるためには十分ではない。例えば、 M{\displaystyle M}zTMz{\displaystyle \mathbf {z} ^{\mathsf {T}}M\mathbf {z} }z.{\displaystyle \mathbf {z} .}M{\displaystyle M}M=[1111],{\displaystyle M={\begin{bmatrix}1&1\\-1&1\end{bmatrix}},}

とを持つ任意の実ベクトルに対して、が成り立ち、が0でなければ常に正である。しかし、が1を持つ複素ベクトルである場合 、z{\displaystyle \mathbf {z} }a{\displaystyle a}b{\displaystyle b}zTMz=(a+b)a+(a+b)b=a2+b2,{\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} }i{\displaystyle i}

zMz=[1i]M[1i]=[1+i1i][1i]=2+2i.{\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,}zTMz>0{\displaystyle \mathbf {z} ^{\mathsf {T}}M\mathbf {z} >0}z{\displaystyle \mathbf {z} }M{\displaystyle M}

表記

エルミート行列が半正定値行列である場合、 と書くことがあります。また、が正定値行列である場合は と書きます。が半負定値行列であることを示すには と書き、が負定値行列であることを示すには と書きます。M{\displaystyle M}M0{\displaystyle M\succeq 0}M{\displaystyle M}M0.{\displaystyle M\succ 0.}M{\displaystyle M}M0{\displaystyle M\preceq 0}M{\displaystyle M}M0.{\displaystyle M\prec 0.}

この概念は関数解析から来ており、半正定値行列は正の作用素を定義します。2つの行列と がを満たす場合、反射的反対称的、推移的な非厳密な半順序を定義できます。ただし、これは全順序ではなく、一般に は不定となる場合があります。 A{\displaystyle A}B{\displaystyle B}BA0,{\displaystyle B-A\succeq 0,}BA{\displaystyle B\succeq A}BA,{\displaystyle B-A,}

一般的な代替表記法として、半正定値行列と正定値行列、半負定値行列と負定値行列をそれぞれ と と表記する。これは、非負行列(それぞれ非正定値行列)もこの表記法で表記されることがあるため、混乱を招く可能性がある。 M0,{\displaystyle M\geq 0,}M>0,{\displaystyle M>0,}M0,{\displaystyle M\leq 0,}M<0{\displaystyle M<0}

影響

上記の定義から、エルミート行列が正定値行列となるのは、それが正定値二次形式またはエルミート形式の行列である場合に限ります。言い換えれば、エルミート行列が正定値行列となるのは、それが内積を定義する場合に限ります。

正定値行列と半正定値行列は様々な方法で特徴付けることができ、数学の様々な分野におけるこの概念の重要性を説明できるかもしれません。エルミート行列Mが正定値行列であるためには、以下の同値な条件のいずれかを満たす必要があります。

行列が半正定値行列であるとは、"正" を "非負" に置き換え、"可逆行列" を "行列" に置き換え、"先頭" という語を削除した同様の条件を満たす場合です。

正定値行列と半正定値実数行列は凸最適化の基礎となる。なぜなら、2回微分可能な複数の実変数の関数が与えられたとき、そのヘッセ行列(その2番目の偏導関数の行列)が点pで正定値であれば関数はpの近くでであり、逆に関数がpの近くで凸であればヘッセ行列はpで半正定値となるからである。p,{\displaystyle p,}p,{\displaystyle p,}p.{\displaystyle p.}

正定値行列の集合は開いた凸錐であり、半正定値行列の集合は閉じた凸錐である。[ 2 ]

  • 単位行列 は正定値行列(したがって半正定値行列でもある)である 。これは実対称行列であり、実数要素abを持つ任意の非零列ベクトルzに対して、I=[1001]{\displaystyle I={\begin{bmatrix}1&0\\0&1\end{bmatrix}}}

    zTIz=[ab][1001][ab]=a2+b2.{\displaystyle \mathbf {z} ^{\mathsf {T}}I\mathbf {z} ={\begin{bmatrix}a&b\end{bmatrix}}{\begin{bmatrix}1&0\\0&1\end{bmatrix}}{\begin{bmatrix}a\\b\end{bmatrix}}=a^{2}+b^{2}.}複素行列として見ると、複素要素abを持つ任意 の非ゼロ列ベクトルzに対して、zIz=[a¯b¯][1001][ab]=a¯a+b¯b=|a|2+|b|2.{\displaystyle \mathbf {z} ^{*}I\mathbf {z} ={\begin{bmatrix}{\overline {a}}&{\overline {b}}\end{bmatrix}}{\begin{bmatrix}1&0\\0&1\end{bmatrix}}{\begin{bmatrix}a\\b\end{bmatrix}}={\overline {a}}a+{\overline {b}}b=|a|^{2}+|b|^{2}.}

    いずれにしても、 はゼロベクトルではないため(つまり、との少なくとも 1 つはゼロではないため) 、結果は正になります。z{\displaystyle \mathbf {z} }a{\displaystyle a}b{\displaystyle b}
  • 実対称行列は、 要素ab 、およびcを持つ 任意の非ゼロ列ベクトルzに対して次が成り立つため、 正定値行列です。 この結果は平方和であるため、非負です。また、 がゼロベクトルのときのみ、 はゼロになります。M=[210121012]{\displaystyle M={\begin{bmatrix}2&-1&0\\-1&2&-1\\0&-1&2\end{bmatrix}}}zTMz=(zTM)z=[(2ab)(a+2bc)(b+2c)][abc]=(2ab)a+(a+2bc)b+(b+2c)c=2a2baab+2b2cbbc+2c2=2a22ab+2b22bc+2c2=a2+a22ab+b2+b22bc+c2+c2=a2+(ab)2+(bc)2+c2{\displaystyle {\begin{aligned}\mathbf {z} ^{\mathsf {T}}M\mathbf {z} =\left(\mathbf {z} ^{\mathsf {T}}M\right)\mathbf {z} &={\begin{bmatrix}(2a-b)&(-a+2b-c)&(-b+2c)\end{bmatrix}}{\begin{bmatrix}a\\b\\c\end{bmatrix}}\\&=(2a-b)a+(-a+2b-c)b+(-b+2c)c\\&=2a^{2}-ba-ab+2b^{2}-cb-bc+2c^{2}\\&=2a^{2}-2ab+2b^{2}-2bc+2c^{2}\\&=a^{2}+a^{2}-2ab+b^{2}+b^{2}-2bc+c^{2}+c^{2}\\&=a^{2}+(a-b)^{2}+(b-c)^{2}+c^{2}\end{aligned}}}a=b=c=0,{\displaystyle a=b=c=0,}z{\displaystyle \mathbf {z} }
  • 任意の実逆行列 に対して、積は正定値行列となる(Aの列の平均が0の場合、これは共分散行列とも呼ばれる)。簡単な証明は、任意の非零ベクトルに対して、行列の逆行列性は次を意味するので、条件が成立することである。A,{\displaystyle A,}ATA{\displaystyle A^{\mathsf {T}}A}z,{\displaystyle \mathbf {z} ,}zTATAz=(Az)T(Az)=Az2>0,{\displaystyle \mathbf {z} ^{\mathsf {T}}A^{\mathsf {T}}A\mathbf {z} =(A\mathbf {z} )^{\mathsf {T}}(A\mathbf {z} )=\|A\mathbf {z} \|^{2}>0,}A{\displaystyle A}Az0.{\displaystyle A\mathbf {z} \neq 0.}
  • 上の例は、一部の要素が負である行列でも正定値行列となる可能性があることを示しています。逆に、すべての要素が正である行列は必ずしも正定値行列とは限りません。例えば 、M{\displaystyle M}N=[1221],{\displaystyle N={\begin{bmatrix}1&2\\2&1\end{bmatrix}},}[11]N[11]T=2<0.{\displaystyle {\begin{bmatrix}-1&1\end{bmatrix}}N{\begin{bmatrix}-1&1\end{bmatrix}}^{\mathsf {T}}=-2<0.}

固有値

をエルミート行列とする(実対称行列を含む)。 のすべての固有値は実数であり、その符号はその正定性を特徴付ける: 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}正の固有値と負の固有値の両方を持つ場合にのみ不定となります。

を の固有分解とします。ここではユニタリ複素行列であり、その列は の固有ベクトルの直交基底から成り、行列あり、その主対角線には対応する固有値が含まれます。この行列は、 (固有ベクトル)基底の座標で再表現された対角行列とみなすことができます。言い換えれば、あるベクトルに を適用することは、を用いて基底を固有ベクトル座標系に変更することと同じです。結果に を引き伸ばし変換を適用し、 とし、その後 を用いて基底を元に戻します。PDP1{\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} ,}Mz,{\displaystyle M\mathbf {z} ,}P1,{\displaystyle P^{-1},}P1z,{\displaystyle P^{-1}\mathbf {z} ,}D{\displaystyle D}DP1z,{\displaystyle DP^{-1}\mathbf {z} ,}P,{\displaystyle P,}PDP1z.{\displaystyle PDP^{-1}\mathbf {z} .}

これを念頭に置くと、変数の一対一変換は、任意の複素ベクトルに対してが実数かつ正である場合と、任意の が実数かつ正である場合、言い換えればが正定値である場合に限り、ということが分かります。対角行列の場合、これは主対角線の各要素、つまり のすべての固有値が正である場合にのみ当てはまります。スペクトル定理により、エルミート行列のすべての固有値は実数であることが保証されるため、実対称行列の特性多項式が利用可能な場合、デカルトの符号交替則を用いて固有値の正性を確認できます。 y=Pz{\displaystyle \mathbf {y} =P\mathbf {z} }zMz{\displaystyle \mathbf {z} ^{*}M\mathbf {z} }z{\displaystyle \mathbf {z} }yDy{\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=BB{\displaystyle M=B^{*}B}B{\displaystyle B}

が実数の場合、も実数となり、分解は次のように書ける。M{\displaystyle M}B{\displaystyle B}M=BTB.{\displaystyle M=B^{\mathsf {T}}B.}

M{\displaystyle M}が正定値行列となることと、それが逆行列を持つ分解が存在することは同値である。より一般的には、が階数 を持つ半正定値行列となることと、それがフル行階数(すなわち階数 )の行列を持つ分解が存在することは同値である。さらに、任意の分解[ 3 ]に対してB{\displaystyle B}M{\displaystyle M}k{\displaystyle k}k×n{\displaystyle k\times n}B{\displaystyle B}k{\displaystyle k}M=BB,{\displaystyle M=B^{*}B,}rank(M)=rank(B).{\displaystyle \operatorname {rank} (M)=\operatorname {rank} (B).}

証拠

ならば、半正定値である。さらに、が逆であるなら、不等式は正定値である。が階数であるなら、M=BB,{\displaystyle M=B^{*}B,}xMx=(xB)(Bx)=Bx20,{\displaystyle x^{*}Mx=(x^{*}B^{*})(Bx)=\|Bx\|^{2}\geq 0,}M{\displaystyle M}B{\displaystyle B}x0,{\displaystyle x\neq 0,}M{\displaystyle M}B{\displaystyle B}k×n{\displaystyle k\times n}k,{\displaystyle k,}rank(M)=rank(B)=k.{\displaystyle \operatorname {rank} (M)=\operatorname {rank} (B^{*})=k.}

逆に、が半正定値行列であると仮定する。 はエルミート行列なので、はユニタリ行列であり、固有値を要素とする対角行列である。半正定値行列な ので、固有値は非負の実数であるため、を の非負の平方根を要素とする対角行列として定義することができる。するとについてとなる。 さらに が正定値行列である場合、固有値は(厳密に)正であるため、 は可逆であり、したがって も可逆である。が階数を持つ場合、 はちょうど正の固有値を持ち、その他はゼロであるため、 を除くすべての行がゼロになる。ゼロ行をカットすると、次のような行列が得られる。M{\displaystyle M}M{\displaystyle M}M=Q1DQ{\displaystyle M=Q^{-1}DQ}Q{\displaystyle Q}D{\displaystyle D}M{\displaystyle M}M{\displaystyle M}D12{\displaystyle D^{\frac {1}{2}}}M=Q1DQ=QDQ=QD12D12Q=QD12D12Q=BB{\displaystyle M=Q^{-1}DQ=Q^{*}DQ=Q^{*}D^{\frac {1}{2}}D^{\frac {1}{2}}Q=Q^{*}D^{{\frac {1}{2}}*}D^{\frac {1}{2}}Q=B^{*}B}B=D12Q.{\displaystyle B=D^{\frac {1}{2}}Q.}M{\displaystyle M}D12{\displaystyle D^{\frac {1}{2}}}B=D12Q{\displaystyle B=D^{\frac {1}{2}}Q}M{\displaystyle M}k,{\displaystyle k,}k{\displaystyle k}B=D12Q{\displaystyle B=D^{\frac {1}{2}}Q}k{\displaystyle k}k×n{\displaystyle k\times n}B{\displaystyle B'}BB=BB=M.{\displaystyle B'^{*}B'=B^{*}B=M.}

の列は、それぞれ複素ベクトル空間または実ベクトル空間のベクトルとして見ることができます。すると、 の要素はこれらのベクトルの内積(つまり、実数の場合は ドット積)です。 言い換えれば、エルミート行列が半正定値行列となるのは、それがいくつかのベクトルのグラム行列である場合に 限ります。エルミート行列が半正定値行列となるのは、それがいくつかの線型独立ベクトルのグラム行列である場合に限ります。一般に、ベクトルのグラム行列の階数は、これらのベクトルが張る空間の次元に等しいです。[ 4 ]b1,,bn{\displaystyle b_{1},\dots ,b_{n}}B{\displaystyle B}Rk,{\displaystyle \mathbb {R} ^{k},}M{\displaystyle M}Mij=bi,bj.{\displaystyle M_{ij}=\langle b_{i},b_{j}\rangle .}M{\displaystyle M}b1,,bn.{\displaystyle b_{1},\dots ,b_{n}.}b1,,bn{\displaystyle b_{1},\dots ,b_{n}}

ユニタリ変換までの一意性

分解は一意ではない。ある行列に対して が成り立ち、が任意のユニタリ行列(つまり)であるとき、M=BB{\displaystyle M=B^{*}B}k×n{\displaystyle k\times n}B{\displaystyle B}Q{\displaystyle Q}k×k{\displaystyle k\times k}QQ=QQ=I{\displaystyle Q^{*}Q=QQ^{*}=I}M=BB=BQQB=AA{\displaystyle M=B^{*}B=B^{*}Q^{*}QB=A^{*}A}A=QB.{\displaystyle A=QB.}

しかし、これは2つの分解が異なる可能性がある唯一の方法です。分解はユニタリ変換を除いて一意です。より正式には、が行列であり、がとなる行列である場合、直交列(つまり)を持つ行列 が存在し、[ 5 ]となる 場合、これはユニタリであることを意味します。 A{\displaystyle A}k×n{\displaystyle k\times n}B{\displaystyle B}×n{\displaystyle \ell \times n}AA=BB,{\displaystyle A^{*}A=B^{*}B,}×k{\displaystyle \ell \times k}Q{\displaystyle Q}QQ=Ik×k{\displaystyle Q^{*}Q=I_{k\times k}}B=QA.{\displaystyle B=QA.}=k{\displaystyle \ell =k}Q{\displaystyle Q}

この記述は、実数の場合、直感的な幾何学的解釈が可能です。つまり、と の列を のベクトルととします。 実ユニタリ行列は直交行列であり、 0 点を保存する(つまり、 の回転鏡映 を平行移動なしで)剛体変換(ユークリッド空間 の等長変換)を記述します。したがって、の何らかの剛体変換によってベクトルがに変換される(および 0 が 0 に変換される)場合にのみ、と の内積は等しくなります。 A{\displaystyle A}B{\displaystyle B}a1,,an{\displaystyle a_{1},\dots ,a_{n}}b1,,bn{\displaystyle b_{1},\dots ,b_{n}}Rk.{\displaystyle \mathbb {R} ^{k}.}Rk{\displaystyle \mathbb {R} ^{k}}aiaj{\displaystyle a_{i}\cdot a_{j}}bibj{\displaystyle b_{i}\cdot b_{j}}Rk{\displaystyle \mathbb {R} ^{k}}a1,,an{\displaystyle a_{1},\dots ,a_{n}}b1,,bn{\displaystyle b_{1},\dots ,b_{n}}

平方根

エルミート行列が正半定値行列である必要十分条件は、正半定値行列 (特に はエルミート行列なので)が存在し、 を満たすことである。この行列は一意であり、[ 6 ]は の非負平方根と呼ばれ、 で表される。が正定値行列である とき、 は であるので、の正平方根とも呼ばれる。M{\displaystyle M}B{\displaystyle B}B{\displaystyle B}B=B{\displaystyle B^{*}=B}M=BB.{\displaystyle M=BB.}B{\displaystyle B}M,{\displaystyle M,}B=M12.{\displaystyle B=M^{\frac {1}{2}}.}M{\displaystyle M}M12,{\displaystyle M^{\frac {1}{2}},}M.{\displaystyle M.}

非負の平方根を他の分解と混同しないでください。 著者によっては、このような分解すべてに対して、または特にコレスキー分解 、または形式の分解すべてに対して、平方根という名前を使用する人もいます が、非負の平方根に対してのみこの名前を使用する人もいます。 M=BB.{\displaystyle M=B^{*}B.}M12{\displaystyle M^{\frac {1}{2}}}M=BB;{\displaystyle M=BB;}

もしそうならMN0{\displaystyle M\succ N\succ 0}M12N120.{\displaystyle M^{\frac {1}{2}}\succ N^{\frac {1}{2}}\succ 0.}

コレスキー分解

エルミート正半正定値行列は と表すことができます。ここでは非負の対角成分を持つ下三角行列です(同様には上三角行列です)。これはコレスキー分解です。 が正定値行列である場合、 の対角成分は正であり、コレスキー分解は一意です。逆に が非負の対角成分を持つ下三角行列である場合、は半正定値行列です。コレスキー分解は、特に効率的な数値計算に役立ちます。密接に関連する分解はLDL 分解です。ここでは対角成分を持ち、は下一三角です。 M{\displaystyle M}M=LL,{\displaystyle M=LL^{*},}L{\displaystyle L}M=BB{\displaystyle M=B^{*}B}B=L{\displaystyle B=L^{*}}M{\displaystyle M}L{\displaystyle L}L{\displaystyle L}LL{\displaystyle LL^{*}}M=LDL,{\displaystyle M=LDL^{*},}D{\displaystyle D}L{\displaystyle L}

ウィリアムソン定理

任意の正定値エルミート実行列は、シンプレクティック(実)行列を介して対角化できます。より正確には、ウィリアムソンの定理は、を満たすシンプレクティックおよび対角実正行列の存在を保証します。 2n×2n{\displaystyle 2n\times 2n}M{\displaystyle M}SSp(2n,R){\displaystyle S\in \mathbf {Sp} (2n,\mathbb {R} )}DRn×n{\displaystyle D\in \mathbb {R} ^{n\times n}}SMST=DD{\displaystyle SMS^{T}=D\oplus D}

その他の特徴

を実対称行列とし、を「単位球」と定義すると、次の式が得られる。 M{\displaystyle M}n×n{\displaystyle n\times n}B1(M){xRn:xTMx1}{\displaystyle B_{1}(M)\equiv \{\mathbf {x} \in \mathbb {R} ^{n}:\mathbf {x} ^{\mathsf {T}}M\mathbf {x} \leq 1\}}M.{\displaystyle M.}

  • B1(vvT){\displaystyle B_{1}(\mathbf {v} \mathbf {v} ^{\mathsf {T}})}間に挟まれた固体の板±{w:w,v=1}.{\displaystyle \pm \{\mathbf {w} :\langle \mathbf {w} ,\mathbf {v} \rangle =1\}.}
  • M0{\displaystyle M\succeq 0}が楕円体、または楕円柱である場合に限ります。B1(M){\displaystyle B_{1}(M)}
  • M0{\displaystyle M\succ 0}が有界である場合、つまり楕円体である場合に限ります。B1(M){\displaystyle B_{1}(M)}
  • もし、もし、もし、もし、もし、もし、もしN0,{\displaystyle N\succ 0,}MN{\displaystyle M\succeq N}B1(M)B1(N);{\displaystyle B_{1}(M)\subseteq B_{1}(N);}MN{\displaystyle M\succ N}B1(M)int(B1(N)).{\displaystyle B_{1}(M)\subseteq \operatorname {int} {\bigl (}B_{1}(N){\bigr )}.}
  • ならば、すべてに対して が成り立つこと、そしてその場合のみ成り立つ。楕円体の極双対も、同じ主軸を持ち長さが逆である楕円体であるので、次の関係が成り立つ。すなわち、が正定値ならば、すべてに対して が成り立つこと、そしてその場合のみ成り立つ。N0,{\displaystyle N\succ 0,}MvvTvTNv{\displaystyle M\succeq {\frac {\mathbf {v} \mathbf {v} ^{\mathsf {T}}}{\mathbf {v} ^{\mathsf {T}}N\mathbf {v} }}}v0{\displaystyle v\neq 0}B1(M)vTNv=1B1(vvT).{\textstyle B_{1}(M)\subset \bigcap _{\mathbf {v} ^{\mathsf {T}}N\mathbf {v} =1}B_{1}(\mathbf {v} \mathbf {v} ^{\mathsf {T}}).}B1(N1)=vTNv=1B1(vvT)=vTNv=1{w:|w,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\}.}N{\displaystyle N}MvvTvTNv{\displaystyle M\succeq {\frac {\mathbf {v} \mathbf {v} ^{\mathsf {T}}}{\mathbf {v} ^{\mathsf {T}}N\mathbf {v} }}}v0{\displaystyle \mathbf {v} \neq \mathbf {0} }MN1.{\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 }Cn×Cn{\displaystyle \mathbb {C} ^{n}\times \mathbb {C} ^{n}}Cn{\displaystyle \mathbb {C} ^{n}}x,yyMx{\displaystyle \langle \mathbf {x} ,\mathbf {y} \rangle \equiv \mathbf {y} ^{*}M\mathbf {x} }x{\displaystyle \mathbf {x} }y{\displaystyle \mathbf {y} }Cn,{\displaystyle \mathbb {C} ^{n},}y{\displaystyle \mathbf {y} ^{*}}y.{\displaystyle \mathbf {y} .}M,{\displaystyle M,}x{\displaystyle x}y.{\displaystyle \mathbf {y} .}Cn{\displaystyle \mathbb {C} ^{n}}z,z{\displaystyle \langle \mathbf {z} ,\mathbf {z} \rangle }z;{\displaystyle \mathbf {z} ;}M{\displaystyle M}Cn{\displaystyle \mathbb {C} ^{n}}
主要マイナーはすべてプラス
行列のk番目主要小行列式は、その左上部分行列の行列式です。行列が正定値行列であるためには、これらすべての行列式が正である必要があります。この条件はシルベスターの基準として知られており、対称実行列の正定値を効率的にテストできます。つまり、ガウス消去法の最初の部分と同様に、基本的な行操作を使用して、行列を上三角行列に簡約し、ピボット処理中に行列式の符号が保持されるように注意します。三角行列のk番目の主要小行列式は、行 までの対角要素の積であるため、シルベスターの基準は、対角要素がすべて正であるかどうかを確認することと同じです。この条件は、三角行列の新しい行が取得されるたびに確認できます。M{\displaystyle M}k×k{\displaystyle k\times k}k,{\displaystyle k,}k{\displaystyle k}

半正定値行列は、逆行列である場合に限り、正定値行列である。[ 7 ]半正定値 行列が負定値行列である場合に限り、半正定値行列である。 M{\displaystyle M}M{\displaystyle -M}

二次形式

実行行列に関連付けられた(純粋な)二次形式は、すべての に対してとなる関数です。これは、両面積で非対称部分は すべてゼロになるため、 と置き換えることで対称であると想定できます。n×n{\displaystyle n\times n}M{\displaystyle M}Q:RnR{\displaystyle Q:\mathbb {R} ^{n}\to \mathbb {R} }Q(x)=xTMx{\displaystyle Q(\mathbf {x} )=\mathbf {x} ^{\mathsf {T}}M\mathbf {x} }x.{\displaystyle \mathbf {x} .}M{\displaystyle M}12(M+MT),{\displaystyle {\tfrac {1}{2}}\left(M+M^{\mathsf {T}}\right),}

対称行列は、その二次形式が厳密に凸関数である場合に限り、正定値行列となります。 M{\displaystyle M}

より一般的には、からへの任意の二次関数は、対称行列、実数ベクトル 、実定数で表すことができます。この場合、これは放物線であり、の場合と同様に、 Rn{\displaystyle \mathbb {R} ^{n}}R{\displaystyle \mathbb {R} }xTMx+bTx+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}M1b,{\displaystyle M^{-1}\mathbf {b} ,}M{\displaystyle M}v{\displaystyle \mathbf {v} }vTMv0,{\displaystyle \mathbf {v} ^{\mathsf {T}}M\mathbf {v} \leq 0,}f(t)(tv)TM(tv)+bT(tv)+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}N{\displaystyle N}(MλN)x=0{\displaystyle \left(M-\lambda N\right)\mathbf {x} =0}x{\displaystyle \mathbf {x} }xTNx=1.{\displaystyle \mathbf {x} ^{\mathsf {T}}N\mathbf {x} =1.}N{\displaystyle N}QTQ.{\displaystyle Q^{\mathsf {T}}Q.}Q{\displaystyle Q}x=QTy,{\displaystyle \mathbf {x} =Q^{\mathsf {T}}\mathbf {y} ,}Q(MλN)QTy=0,{\displaystyle Q\left(M-\lambda N\right)Q^{\mathsf {T}}\mathbf {y} =0,}(QMQT)y=λy{\displaystyle \left(QMQ^{\mathsf {T}}\right)\mathbf {y} =\lambda \mathbf {y} }yTy=1.{\displaystyle \mathbf {y} ^{\mathsf {T}}\mathbf {y} =1.}MX=NXΛ{\displaystyle MX=NX\Lambda }X{\displaystyle X}Λ{\displaystyle \Lambda }XT{\displaystyle X^{\mathsf {T}}}XTMX=Λ{\displaystyle X^{\mathsf {T}}MX=\Lambda }XTNX=I,{\displaystyle X^{\mathsf {T}}NX=I,}yTy=1.{\displaystyle \mathbf {y} ^{\mathsf {T}}\mathbf {y} =1.}M{\displaystyle M}N.{\displaystyle N.}

この結果は、記事「対角化可能行列」における同時対角化に関する記述と矛盾しないことに注意してください。記事「対角化可能行列」では、相似変換による同時対角化について言及されています。ここでの私たちの結果は、2つの二次形式の同時対角化に近く、一方の形式をもう一方の形式に与えられた条件下で最適化する際に役立ちます。

プロパティ

誘導半順序

任意の正方行列について、すなわちが半正定値行列であるとき、 と書きます。これは、すべての正方行列の集合における半順序を定義します。同様に、厳密な半順序を定義することもできます。この順序はLoewner順序と呼ばれます。 M,{\displaystyle M,}N{\displaystyle N}MN{\displaystyle M\geq N}MN0{\displaystyle M-N\geq 0}MN{\displaystyle M-N}M>N.{\displaystyle M>N.}

正定値行列の逆行列

すべての正定値行列は逆行列であり、その逆行列も正定値である。[ 9 ]ならば[ 10 ]さらに、最小最大定理により、のk番目に大きい固有値はのk番目に大きい固有値以上である。MN>0{\displaystyle M\geq N>0}N1M1>0.{\displaystyle N^{-1}\geq M^{-1}>0.}M{\displaystyle M}N.{\displaystyle N.}

スケーリング

が正定値で実数であれば、は正定値である。[ 11 ]M{\displaystyle M}r>0{\displaystyle r>0}rM{\displaystyle rM}

追加

  • とが正定値であれば、その和も正定値である。[ 11 ]M{\displaystyle M}N{\displaystyle N}M+N{\displaystyle M+N}
  • および が半正定値である場合、和も半正定値です。M{\displaystyle M}N{\displaystyle N}M+N{\displaystyle M+N}
  • が正定値で が半正定値である場合、その合計も正定値になります。M{\displaystyle M}N{\displaystyle N}M+N{\displaystyle M+N}

乗算

  • と が正定値であれば、との積も正定値です。 であれば、 も正定値です。M{\displaystyle M}N{\displaystyle N}MNM{\displaystyle MNM}NMN{\displaystyle NMN}MN=NM,{\displaystyle MN=NM,}MN{\displaystyle MN}
  • が半正定値行列であれば、任意の(おそらく長方形の)行列に対して半正定値行列となる。が正定値行列であり、かつフル列ランクを持つ場合、は正定値行列となる。[ 12 ]M{\displaystyle M}AMA{\displaystyle A^{*}MA}A.{\displaystyle A.}M{\displaystyle M}A{\displaystyle A}AMA{\displaystyle A^{*}MA}

トレース

半正定値行列の対角成分は実数かつ非負である。その結果、トレースは、さらに[ 13 ] 、すべての主部分行列(特に2行2列)が半正定値行列であるため 、mii{\displaystyle m_{ii}}tr(M)0.{\displaystyle \operatorname {tr} (M)\geq 0.}|mij|miimjji,j{\displaystyle \left|m_{ij}\right|\leq {\sqrt {m_{ii}m_{jj}}}\quad \forall i,j}n1,{\displaystyle n\geq 1,}maxi,j|mij|maximii{\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)>0and(tr(M))2tr(M2)>n1.{\displaystyle \operatorname {tr} (M)>0\quad \mathrm {and} \quad {\frac {(\operatorname {tr} (M))^{2}}{\operatorname {tr} (M^{2})}}>n-1.}

もう 1 つの重要な結果は、任意の半正定値行列と半正定値行列に対して、次の式が成り立つことです。 行列は半正定値であるため、非負の固有値を持ち、その和であるトレースも非負になります。 M{\displaystyle M}N{\displaystyle N}tr(MN)0.{\displaystyle \operatorname {tr} (MN)\geq 0.}tr(MN)=tr(M12NM12).{\displaystyle \operatorname {tr} (MN)=\operatorname {tr} (M^{\frac {1}{2}}NM^{\frac {1}{2}}).}M12NM12{\displaystyle M^{\frac {1}{2}}NM^{\frac {1}{2}}}

アダマール積

が必ずしも正半定値ではない場合、アダマール積は正半定値となる(この結果はしばしばシュアー積定理と呼ばれる)。[ 15 ]M,N0,{\displaystyle M,N\geq 0,}MN{\displaystyle MN}MN0{\displaystyle M\circ N\geq 0}

2つの半正定値行列のアダマール積に関しては、2つの注目すべき不等式があります。 M=(mij)0,{\displaystyle M=(m_{ij})\geq 0,}N0,{\displaystyle N\geq 0,}

  • オッペンハイムの不等式: [ 16 ]det(MN)det(N)imii.{\displaystyle \det(M\circ N)\geq \det(N)\prod \nolimits _{i}m_{ii}.}
  • det(MN)det(M)det(N).{\displaystyle \det(M\circ N)\geq \det(M)\det(N).}[ 17 ]

クロネッカー積

必ずしも半正定値ではないが、クロネッカー積M,N0,{\displaystyle M,N\geq 0,}MN{\displaystyle MN}MN0.{\displaystyle M\otimes N\geq 0.}

フロベニウス積

が必ずしも半正定値ではない場合、フロベニウスの内積(ランカスター・ティスメネツキー著『行列理論』218 ページ)。 M,N0,{\displaystyle M,N\geq 0,}MN{\displaystyle MN}M:N0{\displaystyle M:N\geq 0}

凸状性

半正定値対称行列の集合は凸行列である。つまり、と が半正定値行列であるならば、0から1までの任意の に対しても半正定値行列となる。任意のベクトル に対して、 M{\displaystyle M}N{\displaystyle N}α{\displaystyle \alpha }αM+(1α)N{\displaystyle \alpha M+\left(1-\alpha \right)N}x{\displaystyle \mathbf {x} }xT(αM+(1α)N)x=αxTMx+(1α)xTNx0.{\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.}

この特性により、半正定値計画問題がグローバルに最適な解に収束することが保証されます。

コサインとの関係

行列の正定値は、任意のベクトルとその像の間の角度が常にA{\displaystyle A}θ{\displaystyle \theta }x{\displaystyle \mathbf {x} }Ax{\displaystyle A\mathbf {x} }π/2<θ<+π/2:{\displaystyle -\pi /2<\theta <+\pi /2:}

cosθ=xTAxxAx=x,AxxAx,θ=θ(x,Ax)(x,Ax)^{\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 }と間の角度x{\displaystyle \mathbf {x} }Ax.{\displaystyle A\mathbf {x} .}

その他の特性

  1. が対称テプリッツ行列である場合、つまり、エントリが絶対インデックス差の関数として与えられ、厳密不等式が成り立つ場合、は厳密に正定値です。M{\displaystyle M}mij{\displaystyle m_{ij}}mij=h(|ij|),{\displaystyle m_{ij}=h(|i-j|),}j0|h(j)|<h(0){\textstyle \sum _{j\neq 0}\left|h(j)\right|<h(0)}M{\displaystyle M}
  2. とエルミートとする。(それぞれ)ならば(それぞれ)。[ 18 ]M>0{\displaystyle M>0}N{\displaystyle N}MN+NM0{\displaystyle MN+NM\geq 0}MN+NM>0{\displaystyle MN+NM>0}N0{\displaystyle N\geq 0}N>0{\displaystyle N>0}
  3. が実数の場合、 が存在し、ここでは単位行列です。M>0{\displaystyle M>0}δ>0{\displaystyle \delta >0}M>δI,{\displaystyle M>\delta I,}I{\displaystyle I}
  4. がリーディングマイナーを表す場合、はLU 分解中のk番目のピボットです。Mk{\displaystyle M_{k}}k×k{\displaystyle k\times k}det(Mk)/det(Mk1){\displaystyle \det \left(M_{k}\right)/\det \left(M_{k-1}\right)}
  5. が奇数のとき、そのk次の主小行列項が負であり、が偶数のとき、正である場合、その行列は負定値行列です。k{\displaystyle k}k{\displaystyle k}
  6. が実正定値行列である場合、任意のベクトルに対して、M{\displaystyle M}m{\displaystyle m}v,{\displaystyle \mathbf {v} ,}vTMvmv22.{\displaystyle \mathbf {v} ^{\mathsf {T}}M\mathbf {v} \geq m\|\mathbf {v} \|_{2}^{2}.}
  7. エルミート行列が半正定値行列であるための必要十分条件は、そのすべての主小行列が非負であることです。しかし、 0と-1の要素を持つ対角行列で確認されているように、主小行列のみを考慮するだけでは不十分です

ブロック行列と部分行列

正の行列はブロックによって定義されることもあります。 2n×2n{\displaystyle 2n\times 2n}M=[ABCD]{\displaystyle M={\begin{bmatrix}A&B\\C&D\end{bmatrix}}}

ここで各ブロックは正値条件を適用することにより、およびはエルミートであることが直ちに示され、n×n,{\displaystyle n\times n,}A{\displaystyle A}D{\displaystyle D}C=B.{\displaystyle C=B^{*}.}

私たちはすべての複雑なもの、特にThen についてそれを持っていますzMz0{\displaystyle \mathbf {z} ^{*}M\mathbf {z} \geq 0}z,{\displaystyle \mathbf {z} ,}z=[v,0]T.{\displaystyle \mathbf {z} =[\mathbf {v} ,0]^{\mathsf {T}}.}[v0][ABBD][v0]=vAv0.{\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,}A{\displaystyle A}D{\displaystyle D}M{\displaystyle M}

逆の結果は、たとえば、シュアー補集合を使用して、ブロックに対するより強い条件で証明できます。

局所的極値

実変数の一般二次形式は 常に次のように表すことができます。ここで、はそれらの変数を含む列ベクトルであり、は対称実行列です。したがって、行列が正定値であるということは、がゼロのときに唯一の最小値(ゼロ)を持ち、それ以外の場合には正であることを意味します。f(x){\displaystyle f(\mathbf {x} )}n{\displaystyle n}x1,,xn{\displaystyle x_{1},\ldots ,x_{n}}xTMx{\displaystyle \mathbf {x} ^{\mathsf {T}}M\mathbf {x} }x{\displaystyle \mathbf {x} }M{\displaystyle M}f{\displaystyle f}x{\displaystyle \mathbf {x} }x.{\displaystyle \mathbf {x} .}

より一般的には、実変数の二回微分可能な実関数は、その勾配がゼロで、そのヘッセ行列(すべての二階微分行列)がその点で半正定値であるとき、引数において局所最小値を持つ。負定値行列と半正定値行列についても同様のことが言える。 f{\displaystyle f}n{\displaystyle n}x1,,xn{\displaystyle x_{1},\ldots ,x_{n}}

共分散

統計学において、多変量確率分布共分散行列は常に半正定値行列である。また、ある変数が他の変数の正確な線形関数でない限り、半正定値行列である。逆に、すべての半正定値行列は、何らかの多変量分布の共分散行列である。

非エルミート正方行列への拡張

正定値の定義は、任意の複素行列(例えば実非対称)を、すべての非零複素ベクトルに対して が成り立つとき、正定値であると定義することで一般化できる。ここで は複素数の実部を表す[ 19 ] 。行列が正定値であるかどうかはエルミート部によってのみ決定され、これは上記の狭義の意味で評価される。同様に、と が実数であるとき、すべての非零実ベクトルに対してが成り立つことと、対称部が狭義の意味で正定値であるときに限る。が の転置に対して無反応であることはすぐに明らかである。M{\displaystyle M}Re{zMz}>0{\displaystyle {\mathcal {R_{e}}}\left\{\mathbf {z} ^{*}M\mathbf {z} \right\}>0}z,{\displaystyle \mathbf {z} ,}Re{c}{\displaystyle {\mathcal {R_{e}}}\{c\}}c.{\displaystyle c.}12(M+M){\textstyle {\frac {1}{2}}\left(M+M^{*}\right)}x{\displaystyle \mathbf {x} }M{\displaystyle M}xTMx>0{\displaystyle \mathbf {x} ^{\mathsf {T}}M\mathbf {x} >0}x{\displaystyle \mathbf {x} }12(M+MT){\textstyle {\frac {1}{2}}\left(M+M^{\mathsf {T}}\right)}xTMx=ijxiMijxj{\textstyle \mathbf {x} ^{\mathsf {T}}M\mathbf {x} =\sum _{ij}x_{i}M_{ij}x_{j}}M.{\displaystyle M.}

正の固有値のみを持つ非対称実行列は、負の固有値を持つ対称部分を持つ場合があります。その場合、行列は正(半)定値行列にはなりません。例えば、行列は正の固有値1と7を持ちますが、と選択します。 M=[4914]{\textstyle M=\left[{\begin{smallmatrix}4&9\\1&4\end{smallmatrix}}\right]}xTMx=2{\displaystyle \mathbf {x} ^{\mathsf {T}}M\mathbf {x} =-2}x=[11]{\displaystyle \mathbf {x} =\left[{\begin{smallmatrix}-1\\1\end{smallmatrix}}\right]}

まとめると、実数と複素数の場合の違いは、複素ヒルベルト空間上の有界正作用素は必然的にエルミート、つまり自己随伴であるという点です。この一般的な主張は分極恒等式を用いて論じることができます。しかし、実数の場合、これはもはや成り立ちません。

アプリケーション

熱伝導率マトリックス

熱伝導に関するフーリエの法則は、熱流束を温度勾配で表し、異方性媒体では と表されます。ここで は熱伝導率行列です。フーリエの法則に負の値が挿入されているのは、熱は常に高温から低温へ流れるという期待を反映しているためです。言い換えれば、温度勾配は常に低温から高温に向かうため、熱流束はとの内積が負になると考えられます。したがって、 となります。フーリエの法則を に代入すると、この期待は となり、伝導率行列は正定値行列となるはずです。通常は は対称ですが、熱ホール効果のように磁場が存在する場合は非対称になります。 q{\displaystyle \mathbf {q} }g=T{\displaystyle \mathbf {g} =\nabla T}q=Kg,{\displaystyle \mathbf {q} =-K\mathbf {g} ,}K{\displaystyle K}g{\displaystyle \mathbf {g} }q{\displaystyle \mathbf {q} }g{\displaystyle \mathbf {g} }qTg<0.{\displaystyle \mathbf {q} ^{\mathsf {T}}\mathbf {g} <0.}gTKg>0,{\displaystyle \mathbf {g} ^{\mathsf {T}}K\mathbf {g} >0,}K{\displaystyle K}

より一般的には、熱力学では、熱と粒子の流れは、オンサガーの逆関係で説明される完全に結合したシステムであり、エントロピー生成が非負になるためには、結合行列は半正定値(非対称の場合もある)である必要があります。

参照

参考文献

  1. ^ 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
  2. ^ Boyd, Stephen; Vandenberghe, Lieven (2004年3月8日).凸最適化. Cambridge University Press. doi : 10.1017/cbo9780511804441 . ISBN 978-0-521-83378-3
  3. ^ Horn & Johnson (2013)、p. 440、定理7.2.7
  4. ^ Horn & Johnson (2013)、p. 441、定理7.2.10
  5. ^ Horn & Johnson (2013)、p. 452、定理7.3.11
  6. ^ Horn & Johnson (2013)、p. 439、定理7.2.6とk=2{\displaystyle k=2}
  7. ^ホーン&ジョンソン(2013)、p. 431、補論7.1.7
  8. ^ Horn & Johnson (2013)、p. 485、定理7.6.1
  9. ^ Horn & Johnson (2013)、p. 438、定理7.2.1
  10. ^ホーン&ジョンソン(2013)、p. 495、補論7.7.4(a)
  11. ^ a bホーン&ジョンソン(2013)、p. 430、観察7.1.3
  12. ^ホーン&ジョンソン(2013)、p.431、観察7.1.8
  13. ^ホーン&ジョンソン(2013)、430ページ
  14. ^ Wolkowicz, Henry; Styan, George PH (1980). 「トレースを用いた固有値の境界」.線形代数とその応用. 29 (29). Elsevier: 471– 506. doi : 10.1016/0024-3795(80)90258-X .
  15. ^ Horn & Johnson (2013)、p. 479、定理7.5.3
  16. ^ Horn & Johnson (2013)、p. 509、定理7.8.16
  17. ^ Styan, GP (1973). 「アダマール積と多変量統計解析」.線形代数とその応用. 6 : 217–240 . doi : 10.1016/0024-3795(73)90023-2 .、系3.6、p. 227
  18. ^ Bhatia, Rajendra (2007).正定値行列. プリンストン大学出版局, ニュージャージー州, 8頁. ISBN 978-0-691-12918-1
  19. ^ Weisstein, Eric W. 「正定値行列」 . MathWorld . Wolfram Research . 2012年7月26日閲覧

出典