数学的行列の性質
数学 において、 実数 要素 を持つ 対称行列が 正定値 行列であるとは、実数が すべての非ゼロの実 列ベクトル に対して正である場合である。ここで、 はの 行ベクトル 転置を 表す。 [1]
より一般的には、 エルミート行列 (つまり、その 共役転置 に等しい 複素行列 )が 正定値行列 であるとは、実数が すべての非ゼロの複素列ベクトルに対して正である場合である。 ここで、 はの共役転置を表す。
ま
{\displaystyle \M\}
x
⊤
ま
x
{\displaystyle \\mathbf {x} ^{\top }M\mathbf {x} \}
x
、
{\displaystyle \\mathbf {x}\,}
x
⊤
{\displaystyle \ \mathbf {x} ^{\top }\ }
x
.
{\displaystyle \ \mathbf {x} ~.}
z
∗
M
z
{\displaystyle \ \mathbf {z} ^{*}M\mathbf {z} \ }
z
,
{\displaystyle \ \mathbf {z} \ ,}
z
∗
{\displaystyle \ \mathbf {z} ^{*}\ }
z
.
{\displaystyle \ \mathbf {z} ~.}
正の半定値 行列も同様に定義されますが、スカラー とが正 またはゼロ (つまり、非負) である必要があります。 負の定値行列 と 負の半定値 行列も同様に定義されます。正の半定値でも負の半定値でもない行列は、 不定値 と呼ばれることがあります。
x
⊤
M
x
{\displaystyle \ \mathbf {x} ^{\top }M\mathbf {x} \ }
z
∗
M
z
{\displaystyle \ \mathbf {z} ^{*}M\mathbf {z} \ }
影響
上記の定義から、行列が正定値となる のは、それが 正定値二次形式 または エルミート形式 の行列である場合のみであることがわかります 。言い換えると、行列が正定値となるのは、内積を定義する場合のみであるということ です 。
正定値行列と半正定値行列はさまざまな方法で特徴付けることができ、数学のさまざまな部分でこの概念が重要である理由を説明できます。行列 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]
一部の著者は、非対称実数行列や非エルミート複素行列を含む、より一般的な定性の定義を使用しています。
定義
以下の定義では、 は の転置、 は の 共役転置 、 は n 次元の ゼロベクトル
を表します。
x
⊤
{\displaystyle \ \mathbf {x} ^{\top }\ }
x
,
{\displaystyle \ \mathbf {x} \ ,}
z
∗
{\displaystyle \ \mathbf {z} ^{*}\ }
z
,
{\displaystyle \ \mathbf {z} \ ,}
0
{\displaystyle \ \mathbf {0} \ }
実数行列の定義
対称実数行列は 、 すべての非ゼロの 行列に対して 正定値行列 であるという 。正式
には、
n
×
n
{\displaystyle n\times n}
M
{\displaystyle \ M\ }
x
⊤
M
x
>
0
{\displaystyle \ \mathbf {x} ^{\top }M\ \mathbf {x} >0\ }
x
{\displaystyle \ \mathbf {x} \ }
R
n
.
{\displaystyle \ \mathbb {R} ^{n}~.}
M
positive-definite
⟺
x
⊤
M
x
>
0
for all
x
∈
R
n
∖
{
0
}
{\displaystyle \ M{\text{ positive-definite}}\quad \iff \quad \mathbf {x} ^{\top }M\ \mathbf {x} >0{\text{ for all }}\mathbf {x} \in \mathbb {R} ^{n}\setminus \{\mathbf {0} \}\ }
対称 実行列は、 すべての に対して となるとき、 半正定値行列 または 非負定値行列 であると言われる 。正式
には、
n
×
n
{\displaystyle \ n\times n\ }
M
{\displaystyle \ M\ }
x
⊤
M
x
≥
0
{\displaystyle \ \mathbf {x} ^{\top }M\ \mathbf {x} \geq 0\ }
x
{\displaystyle \ \mathbf {x} \ }
R
n
.
{\displaystyle \ \mathbb {R} ^{n}~.}
M
positive semi-definite
⟺
x
⊤
M
x
≥
0
for all
x
∈
R
n
{\displaystyle \ M{\text{ positive semi-definite}}\quad \iff \quad \mathbf {x} ^{\top }M\ \mathbf {x} \geq 0{\text{ for all }}\mathbf {x} \in \mathbb {R} ^{n}\ }
対称実数行列は 、 すべての非ゼロに対して 負 定値 であるという 。正式
には、
n
×
n
{\displaystyle \ n\times n\ }
M
{\displaystyle \ M\ }
x
⊤
M
x
<
0
{\displaystyle \ \mathbf {x} ^{\top }M\ \mathbf {x} <0\ }
x
{\displaystyle \ \mathbf {x} \ }
R
n
.
{\displaystyle \ \mathbb {R} ^{n}~.}
M
negative-definite
⟺
x
⊤
M
x
<
0
for all
x
∈
R
n
∖
{
0
}
{\displaystyle \ M{\text{ negative-definite}}\quad \iff \quad \mathbf {x} ^{\top }M\ \mathbf {x} <0{\text{ for all }}\mathbf {x} \in \mathbb {R} ^{n}\setminus \{\mathbf {0} \}\ }
対称実行列は、 すべて のに対して が成り立つとき、 半負定値行列 または 非正定値行列 であると言われる 。正式
には、
n
×
n
{\displaystyle \ n\times n\ }
M
{\displaystyle \ M\ }
x
⊤
M
x
≤
0
{\displaystyle \ \mathbf {x} ^{\top }M\ \mathbf {x} \leq 0\ }
x
{\displaystyle \ \mathbf {x} \ }
R
n
.
{\displaystyle \ \mathbb {R} ^{n}~.}
M
negative semi-definite
⟺
x
⊤
M
x
≤
0
for all
x
∈
R
n
{\displaystyle M{\text{ negative semi-definite}}\quad \iff \quad \mathbf {x} ^{\top }M\ \mathbf {x} \leq 0{\text{ for all }}\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
positive-definite
⟺
z
∗
M
z
>
0
for all
z
∈
C
n
∖
{
0
}
{\displaystyle \ M{\text{ positive-definite}}\quad \iff \quad \mathbf {z} ^{*}M\ \mathbf {z} >0{\text{ for all }}\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
positive semi-definite
⟺
z
∗
M
z
≥
0
for all
z
∈
C
n
{\displaystyle \ M{\text{ positive semi-definite}}\quad \iff \quad \mathbf {z} ^{*}M\ \mathbf {z} \geq 0{\text{ for all }}\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
negative-definite
⟺
z
∗
M
z
<
0
for all
z
∈
C
n
∖
{
0
}
{\displaystyle \ M{\text{ negative-definite}}\quad \iff \quad \mathbf {z} ^{*}M\ \mathbf {z} <0{\text{ for all }}\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
negative semi-definite
⟺
z
∗
M
z
≤
0
for all
z
∈
C
n
{\displaystyle \ M{\text{ negative semi-definite}}\quad \iff \quad \mathbf {z} ^{*}M\ \mathbf {z} \leq 0{\text{ for all }}\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
⊤
M
z
{\displaystyle \ \mathbf {z} ^{\top }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} \ }
a
{\displaystyle \ a\ }
b
{\displaystyle \ b\ }
z
⊤
M
z
=
(
a
+
b
)
a
+
(
−
a
+
b
)
b
=
a
2
+
b
2
,
{\displaystyle \ \mathbf {z} ^{\top }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\ ,}
z
∗
M
z
=
[
1
−
i
]
M
[
1
i
]
=
[
1
+
i
1
−
i
]
[
1
i
]
=
2
+
2
i
.
{\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
⊤
M
z
>
0
{\displaystyle \ \mathbf {z} ^{\top }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 つの行列 が と を満たす場合、 反射的 、 反対称的 、 推移的な 非厳密な半順序 を定義できます。 ただし、これは 全順序 で はなく、 一般に不定である可能性があります。
A
{\displaystyle \ A\ }
B
{\displaystyle \ B\ }
B
−
A
⪰
0
,
{\displaystyle \ B-A\succeq 0\ ,}
B
⪰
A
{\displaystyle \ B\succeq A\ }
B
−
A
,
{\displaystyle \ B-A\ ,}
一般的な代替表記法は、 それぞれ半正定値行列と正定値行列、半負定値行列と負定値行列です。非負行列 (それぞれ非正行列) もこの方法で表記される場合があるため、混乱を招く可能性が あり ます 。
M
≥
0
,
{\displaystyle \ M\geq 0\ ,}
M
>
0
,
{\displaystyle \ M>0\ ,}
M
≤
0
,
{\displaystyle \ M\leq 0\ ,}
M
<
0
{\displaystyle \ M<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}
正の固有値と負の固有値の両方を持つ場合にのみ不定となります。
を の 固有分解 とします 。 ここで は ユニタリ複素行列 で、 その列は の 固有ベクトル の 正規直交基底 から成り 、は 実 対角行列 で、 その 主対角線には対応する 固有値 が含まれます 。この行列は、 (固有ベクトル)基底の座標で再表現された 対角行列と見なすことができます。 言い換えると、 あるベクトルに を適用する ことは、 を用いて 基底を 固有ベクトル座標系に変更すること と同じです。 結果に ストレッチ変換を 適用してとし、次に を 用いて基底を元に戻します。
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} ~.}
これを念頭に置くと、変数の 1 対 1 変換は 、 任意の複素ベクトルに対して が実数かつ正である場合、かつ が任意 の に対して実数かつ正である 場合、 つまり が 正定値である場合に限り、であることを示します。対角行列の場合、これは主対角線の各要素、つまり のすべての固有値 が正である場合にのみ当てはまります。スペクトル 定理 により、エルミート行列のすべての固有値が実数であることが保証されるため、 実対称行列の 特性多項式 が使用可能な場合、 デカルトの符号交替規則 を使用して、固有値の正値性を確認できます。
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
⊤
B
.
{\displaystyle \ M=B^{\top }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
=
B
∗
B
,
{\displaystyle \ M=B^{*}B\ ,}
rank
(
M
)
=
rank
(
B
)
.
{\displaystyle \ \operatorname {rank} (M)=\operatorname {rank} (B)~.}
証拠
ならば は 半 正定値である。さらに が逆ならば に対して 不等式 は厳密であるので は 正定値である。 が 階数 なら ば
M
=
B
∗
B
,
{\displaystyle \ M=B^{*}B\ ,}
x
∗
M
x
=
(
x
∗
B
∗
)
(
B
x
)
=
‖
B
x
‖
2
≥
0
,
{\displaystyle \ x^{*}Mx=(x^{*}B^{*})(Bx)=\|Bx\|^{2}\geq 0\ ,}
M
{\displaystyle \ M\ }
B
{\displaystyle B}
x
≠
0
,
{\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
=
Q
−
1
D
Q
{\displaystyle \ M=Q^{-1}DQ\ }
Q
{\displaystyle \ Q\ }
D
{\displaystyle \ D\ }
M
{\displaystyle \ M\ }
M
{\displaystyle \ M\ }
D
1
2
{\displaystyle \ D^{\frac {1}{2}}\ }
M
=
Q
−
1
D
Q
=
Q
∗
D
Q
=
Q
∗
D
1
2
D
1
2
Q
=
Q
∗
D
1
2
∗
D
1
2
Q
=
B
∗
B
{\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
=
D
1
2
Q
.
{\displaystyle \ B=D^{\frac {1}{2}}Q~.}
M
{\displaystyle M}
D
1
2
{\displaystyle \ D^{\frac {1}{2}}\ }
B
=
D
1
2
Q
{\displaystyle \ B=D^{\frac {1}{2}}Q\ }
M
{\displaystyle \ M\ }
k
,
{\displaystyle \ k\ ,}
k
{\displaystyle k}
B
=
D
1
2
Q
{\displaystyle \ B=D^{\frac {1}{2}}Q\ }
k
{\displaystyle k}
k
×
n
{\displaystyle \ k\times n\ }
B
′
{\displaystyle \ B'\ }
B
′
∗
B
′
=
B
∗
B
=
M
.
{\displaystyle \ B'^{*}B'=B^{*}B=M~.}
の 列は、それぞれ 複素 ベクトル空間または 実ベクトル空間 のベクトルとして見ることができます 。すると、 の要素はこれらのベクトルの 内積 (つまり 、実数の場合は
ドット積 )です。
言い換えると、エルミート行列が 半正定値行列となるのは、それが いくつかのベクトルの グラム行列である場合に限ります。エルミート行列が半正定値行列となるのは、それがいくつかの 線形独立 ベクトル
のグラム行列である場合に限ります 。一般に、ベクトルのグラム行列の階数は、これらのベクトルが 張る 空間の次元に等しくなります 。 [4]
b
1
,
…
,
b
n
{\displaystyle \ b_{1},\dots ,b_{n}\ }
B
{\displaystyle \ B\ }
R
k
,
{\displaystyle \ \mathbb {R} ^{k}\ ,}
M
{\displaystyle M}
M
i
j
=
⟨
b
i
,
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\ }
k
×
n
{\displaystyle \ k\times n\ }
B
{\displaystyle \ B\ }
Q
{\displaystyle \ Q\ }
k
×
k
{\displaystyle k\times k}
Q
∗
Q
=
Q
Q
∗
=
I
{\displaystyle \ Q^{*}Q=QQ^{*}=I\ }
M
=
B
∗
B
=
B
∗
Q
∗
Q
B
=
A
∗
A
{\displaystyle \ M=B^{*}B=B^{*}Q^{*}QB=A^{*}A\ }
A
=
Q
B
.
{\displaystyle \ A=QB~.}
しかし、これは2つの分解が異なる唯一の方法です。分解は ユニタリ変換 を除いて一意です。より正式には、が行列 であり 、が となる行列 である場合、 正規直交列(つまり )を持つ 行列
が存在し、 [5] となる
とき、は ユニタリ であることを 意味します 。
A
{\displaystyle \ A\ }
k
×
n
{\displaystyle \ k\times n\ }
B
{\displaystyle \ B\ }
ℓ
×
n
{\displaystyle \ \ell \times n\ }
A
∗
A
=
B
∗
B
,
{\displaystyle \ A^{*}A=B^{*}B\ ,}
ℓ
×
k
{\displaystyle \ \ell \times k\ }
Q
{\displaystyle Q}
Q
∗
Q
=
I
k
×
k
{\displaystyle \ Q^{*}Q=I_{k\times k}\ }
B
=
Q
A
.
{\displaystyle \ B=QA~.}
ℓ
=
k
{\displaystyle \ell =k}
Q
{\displaystyle Q}
この記述は、実際のケースでは直感的な幾何学的解釈が可能です。 およびの列を の ベクトル およびとします
。 実ユニタリ行列は 直交行列 で、 0 点を保存する 剛体変換 (ユークリッド空間 の等長変換) を記述します(つまり、平行移動なしで 回転 と 反射 )。 したがって、の何らかの剛体変換によって ベクトル が に変換される (および 0 が 0 に変換される) 場合にのみ 、および の内積が等しくなります 。
A
{\displaystyle A}
B
{\displaystyle B}
a
1
,
…
,
a
n
{\displaystyle a_{1},\dots ,a_{n}}
b
1
,
…
,
b
n
{\displaystyle \ b_{1},\dots ,b_{n}\ }
R
k
.
{\displaystyle \ \mathbb {R} ^{k}~.}
R
k
{\displaystyle \mathbb {R} ^{k}}
a
i
⋅
a
j
{\displaystyle a_{i}\cdot a_{j}}
b
i
⋅
b
j
{\displaystyle b_{i}\cdot b_{j}}
R
k
{\displaystyle \mathbb {R} ^{k}}
a
1
,
…
,
a
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
≻
N
≻
0
{\displaystyle \ M\succ N\succ 0\ }
M
1
2
≻
N
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
)
≡
{
x
∈
R
n
:
x
⊤
M
x
≤
1
}
{\displaystyle \ B_{1}(M)\equiv \{\mathbf {x} \in \mathbb {R} ^{n}:\mathbf {x} ^{\top }M\ \mathbf {x} \leq 1\}\ }
M
.
{\displaystyle \ M~.}
B
1
(
v
v
⊤
)
{\displaystyle \ B_{1}(\mathbf {v} \ \mathbf {v} ^{\top })}
間に挟まれた固体の板である
±
{
w
:
⟨
w
,
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)\ }
もし 、 もし、もし、もし 、もし、もし、もし
N
≻
0
,
{\displaystyle \ N\succ 0\ ,}
M
⪰
N
{\displaystyle \ M\succeq N\ }
B
1
(
M
)
⊆
B
1
(
N
)
;
{\displaystyle \ B_{1}(M)\subseteq B_{1}(N)\ ;}
M
≻
N
{\displaystyle \ M\succ N\ }
B
1
(
M
)
⊆
int
(
B
1
(
N
)
)
.
{\displaystyle \ B_{1}(M)\subseteq \operatorname {int} \!{\bigl (}\ B_{1}(N)\ {\bigr )}~.}
ならば、 すべてに対して ならば 、かつその ときのみ と なる。したがって、楕円体の極双対も、同じ主軸を持ち長さが逆である楕円体なので、次式が成り立つ。 すなわち、 が 正定値であれば、 すべてに対して ならば、かつそのときのみ
N
≻
0
,
{\displaystyle \ N\succ 0\ ,}
M
⪰
v
v
⊤
v
⊤
N
v
{\displaystyle \ M\succeq {\frac {\mathbf {v} \ \mathbf {v} ^{\top }}{\ \mathbf {v} ^{\top }N\ \mathbf {v} \ }}\ }
v
≠
0
{\displaystyle v\neq 0}
B
1
(
M
)
⊂
⋂
v
⊤
N
v
=
1
B
1
(
v
v
⊤
)
.
{\textstyle \ B_{1}(M)\subset \bigcap _{\mathbf {v} ^{\top }N\ \mathbf {v} =1}B_{1}(\mathbf {v} \mathbf {v} ^{\top })~.}
B
1
(
N
−
1
)
=
⋂
v
⊤
N
v
=
1
B
1
(
v
v
⊤
)
=
⋂
v
⊤
N
v
=
1
{
w
:
|
⟨
w
,
v
⟩
|
≤
1
}
.
{\displaystyle \ B_{1}(N^{-1})=\bigcap _{\mathbf {v} ^{\top }N\ \mathbf {v} =1}B_{1}(\mathbf {v} \ \mathbf {v} ^{\top })=\bigcap _{\mathbf {v} ^{\top }N\ \mathbf {v} =1}\{\mathbf {w} :|\langle \mathbf {w} ,\mathbf {v} \rangle |\leq 1\}~.}
N
{\displaystyle \ N\ }
M
⪰
v
v
⊤
v
⊤
N
v
{\displaystyle \ M\succeq {\frac {\mathbf {v} \mathbf {v} ^{\top }}{\ \mathbf {v} ^{\top }N\ \mathbf {v} \ }}\ }
v
≠
0
{\displaystyle \ \mathbf {v} \neq \mathbf {0} \ }
M
⪰
N
−
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}\ }
⟨
x
,
y
⟩
≡
y
∗
M
x
{\displaystyle \ \langle \mathbf {x} ,\mathbf {y} \rangle \equiv \mathbf {y} ^{*}M\ \mathbf {x} \ }
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\ ,}
x
{\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\ }
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
:
R
n
→
R
{\displaystyle \ Q:\mathbb {R} ^{n}\to \mathbb {R} \ }
Q
(
x
)
=
x
⊤
M
x
{\displaystyle \ Q(\mathbf {x} )=\mathbf {x} ^{\top }M\mathbf {x} \ }
x
.
{\displaystyle \ \mathbf {x} ~.}
M
{\displaystyle \ M\ }
1
2
(
M
+
M
⊤
)
,
{\displaystyle \ {\tfrac {1}{2}}\left(M+M^{\top }\right)\ ,}
対称行列は、その二次形式が 厳密に凸関数 である場合に限り、正定値行列となります 。
M
{\displaystyle \ M\ }
より一般的には、から へ の任意の 2次関数は 、対称 行列、 実数 n ベクトル、 実 定数として 次の ように表すことができます。この 場合、これは放物線であり、の 場合と同様に、
R
n
{\displaystyle \ \mathbb {R} ^{n}\ }
R
{\displaystyle \mathbb {R} }
x
⊤
M
x
+
b
⊤
x
+
c
{\displaystyle \ \mathbf {x} ^{\top }M\mathbf {x} +\mathbf {b} ^{\top }\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
⊤
M
v
≤
0
,
{\displaystyle \ \mathbf {v} ^{\top }M\mathbf {v} \leq 0\ ,}
f
(
t
)
≡
(
t
v
)
⊤
M
(
t
v
)
+
b
⊤
(
t
v
)
+
c
{\displaystyle \ f(t)\equiv (t\mathbf {v} )^{\top }M(t\mathbf {v} )+b^{\top }(t\mathbf {v} )+c\ }
このため、正定値行列は 最適化 問題において重要な役割を果たします。
同時対角化
1 つの対称行列と、対称かつ正定値である別の行列は、 同時に対角化できます。同時対角化は 相似変換 で必ずしも実行されるわけではありませんが、これは可能です 。この結果は、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} \ }
x
⊤
N
x
=
1
.
{\displaystyle \ \mathbf {x} ^{\top }N\mathbf {x} =1~.}
N
{\displaystyle \ N\ }
Q
⊤
Q
.
{\displaystyle \ Q^{\top }Q~.}
Q
{\displaystyle \ Q\ }
x
=
Q
⊤
y
,
{\displaystyle \ \mathbf {x} =Q^{\top }\mathbf {y} \ ,}
Q
(
M
−
λ
N
)
Q
⊤
y
=
0
,
{\displaystyle \ Q\left(M-\lambda N\right)Q^{\top }\mathbf {y} =0\ ,}
(
Q
M
Q
⊤
)
y
=
λ
y
{\displaystyle \ \left(QMQ^{\top }\right)\mathbf {y} =\lambda \mathbf {y} \ }
y
⊤
y
=
1
.
{\displaystyle \mathbf {y} ^{\top }\mathbf {y} =1~.}
M
X
=
N
X
Λ
{\displaystyle MX=NX\Lambda }
X
{\displaystyle X}
Λ
{\displaystyle \Lambda }
X
⊤
{\displaystyle X^{\top }}
X
⊤
M
X
=
Λ
{\displaystyle \ X^{\top }MX=\Lambda \ }
X
⊤
N
X
=
I
,
{\displaystyle \ X^{\top }NX=I\ ,}
y
⊤
y
=
1
.
{\displaystyle \ \mathbf {y} ^{\top }\mathbf {y} =1~.}
M
{\displaystyle \ M\ }
N
.
{\displaystyle \ N~.}
この結果は、相似変換による同時対角化について言及している記事「対角化可能な行列」 の同時対角化に関する記述と矛盾しないことに注意してください 。ここでの結果は、2 つの二次形式の同時対角化に似ており、一方の形式をもう一方の形式で最適化するのに役立ちます。
プロパティ
誘導された部分順序
任意の正方行列について、 つまり が 半正定値である 場合、 と書きます。これは、 すべての正方行列の集合上の 半順序 を定義します。同様に、厳密な半順序を定義することもできます。この順序は、 Loewner 順序 と呼ばれます。
M
,
{\displaystyle \ M\ ,}
N
{\displaystyle \ N\ }
M
≥
N
{\displaystyle \ M\geq N\ }
M
−
N
≥
0
{\displaystyle \ M-N\geq 0\ }
M
−
N
{\displaystyle \ M-N\ }
M
>
N
.
{\displaystyle \ M>N~.}
正定値行列の逆行列
全ての正定値行列は 逆行列 であり、その逆行列も正定値である。 [9] ならば [10] さらに 、 最小最大定理 により、の k 番目に大きい固有値は、の k 番目に大きい固有値 以上である。
M
≥
N
>
0
{\displaystyle \ M\geq N>0\ }
N
−
1
≥
M
−
1
>
0
.
{\displaystyle \ N^{-1}\geq M^{-1}>0~.}
M
{\displaystyle \ M\ }
N
.
{\displaystyle \ N~.}
スケーリング
が正定値で 実数であれ ば、は 正定値である。 [11]
M
{\displaystyle \ M\ }
r
>
0
{\displaystyle \ r>0\ }
r
M
{\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\ }
M
N
M
{\displaystyle \ MNM\ }
N
M
N
{\displaystyle NMN}
M
N
=
N
M
,
{\displaystyle \ MN=NM\ ,}
M
N
{\displaystyle \ MN\ }
が半正定値行列である場合 、 任意の(おそらく長方形の)行列に対して半正定値行列である。が 正定値行列であり、 列ランクが完全である 場合、 半正定値行列である。 [12]
M
{\displaystyle \ M\ }
A
∗
M
A
{\displaystyle \ A^{*}MA\ }
A
.
{\displaystyle \ A~.}
M
{\displaystyle \ M\ }
A
{\displaystyle A}
A
∗
M
A
{\displaystyle \ A^{*}MA\ }
トレース
半正定値行列の対角要素 は実数かつ非負である。その結果、 トレースは 、 さらに、 [13] すべての主要な部分行列(特に2行2列)が半正定値であるため、
m
i
i
{\displaystyle \ m_{ii}\ }
tr
(
M
)
≥
0
.
{\displaystyle \ \operatorname {tr} (M)\geq 0~.}
|
m
i
j
|
≤
m
i
i
m
j
j
∀
i
,
j
{\displaystyle \ \left|m_{ij}\right|\leq {\sqrt {m_{ii}m_{jj}}}\quad \forall i,j\ }
そして、
n
≥
1
,
{\displaystyle \ n\geq 1\ ,}
max
i
,
j
|
m
i
j
|
≤
max
i
m
i
i
{\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
a
n
d
(
tr
(
M
)
)
2
tr
(
M
2
)
>
n
−
1
.
{\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
(
M
N
)
≥
0
.
{\displaystyle \ \operatorname {tr} (MN)\geq 0~.}
tr
(
M
N
)
=
tr
(
M
1
2
N
M
1
2
)
.
{\displaystyle \ \operatorname {tr} (MN)=\operatorname {tr} (M^{\frac {1}{2}}NM^{\frac {1}{2}})~.}
M
1
2
N
M
1
2
{\displaystyle M^{\frac {1}{2}}NM^{\frac {1}{2}}}
アダマール積
が 必ずしも正半定値ではない 場合、 アダマール積は 正半定値となる (この結果はしばしば シューア積定理 と呼ばれる)。 [15]
M
,
N
≥
0
,
{\displaystyle \ M,N\geq 0\ ,}
M
N
{\displaystyle \ MN\ }
M
∘
N
≥
0
{\displaystyle \ M\circ N\geq 0\ }
2 つの半正定値行列のアダマール積に関しては、 2 つの注目すべき不等式があります。
M
=
(
m
i
j
)
≥
0
,
{\displaystyle \ M=(m_{ij})\geq 0\ ,}
N
≥
0
,
{\displaystyle \ N\geq 0\ ,}
オッペンハイムの不等式: [16]
det
(
M
∘
N
)
≥
det
(
N
)
∏
i
m
i
i
.
{\displaystyle \ \det(M\circ N)\geq \det(N)\prod \nolimits _{i}m_{ii}~.}
det
(
M
∘
N
)
≥
det
(
M
)
det
(
N
)
.
{\displaystyle \ \det(M\circ N)\geq \det(M)\det(N)~.}
[17]
クロネッカー積
が必ずしも半正定値ではない 場合 、 クロネッカー積
M
,
N
≥
0
,
{\displaystyle \ M,N\geq 0\ ,}
M
N
{\displaystyle \ MN\ }
M
⊗
N
≥
0
.
{\displaystyle M\otimes N\geq 0~.}
フロベニウス積
が 必ずしも半正定値ではない 場合、 フロベニウスの内積 (ランカスター・ティスメネツキー著 『行列理論 』218ページ)。
M
,
N
≥
0
,
{\displaystyle \ M,N\geq 0\ ,}
M
N
{\displaystyle \ MN\ }
M
:
N
≥
0
{\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} \ }
x
⊤
(
α
M
+
(
1
−
α
)
N
)
x
=
α
x
⊤
M
x
+
(
1
−
α
)
x
⊤
N
x
≥
0
.
{\displaystyle \ \mathbf {x} ^{\top }\left(\alpha M+\left(1-\alpha \right)N\right)\mathbf {x} =\alpha \mathbf {x} ^{\top }M\mathbf {x} +(1-\alpha )\mathbf {x} ^{\top }N\mathbf {x} \geq 0~.}
この特性により、 半正定値計画 問題がグローバルに最適な解に収束することが保証されます。
コサインとの関係
行列の正定値は、 任意のベクトル とその像の間の 角度 が常に
A
{\displaystyle A}
θ
{\displaystyle \ \theta \ }
x
{\displaystyle \ \mathbf {x} \ }
A
x
{\displaystyle \ A\mathbf {x} \ }
−
π
/
2
<
θ
<
+
π
/
2
:
{\displaystyle \ -\pi /2<\theta <+\pi /2\ :}
cos
θ
=
x
⊤
A
x
‖
x
‖
‖
A
x
‖
=
⟨
x
,
A
x
⟩
‖
x
‖
‖
A
x
‖
,
θ
=
θ
(
x
,
A
x
)
≡
(
x
,
A
x
)
^
≡
{\displaystyle \ \cos \theta ={\frac {\mathbf {x} ^{\top }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} \ }
A
x
.
{\displaystyle \ A\mathbf {x} ~.}
その他の特性
が対称 テプリッツ行列 である場合 、つまり、要素が 絶対インデックス差の関数として与えられ、 厳密 な 不等式 が成り立つ場合、は 厳密に 正定値 です。
M
{\displaystyle M}
m
i
j
{\displaystyle m_{ij}}
m
i
j
=
h
(
|
i
−
j
|
)
,
{\displaystyle \ m_{ij}=h(|i-j|)\ ,}
∑
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}
N
{\displaystyle N}
M
N
+
N
M
≥
0
{\displaystyle MN+NM\geq 0}
M
N
+
N
M
>
0
{\displaystyle MN+NM>0}
N
≥
0
{\displaystyle N\geq 0}
N
>
0
{\displaystyle N>0}
が実数の場合 、 と なるようなが存在し 、 は 単位行列 です 。
M
>
0
{\displaystyle \ M>0\ }
δ
>
0
{\displaystyle \ \delta >0\ }
M
>
δ
I
,
{\displaystyle \ M>\delta I\ ,}
I
{\displaystyle \ I\ }
がリーディングマイナーを表す 場合 、は LU 分解 中の k 番目のピボットです 。
M
k
{\displaystyle M_{k}}
k
×
k
{\displaystyle \ k\times k\ }
det
(
M
k
)
/
det
(
M
k
−
1
)
{\displaystyle \ \det \left(M_{k}\right)/\det \left(M_{k-1}\right)\ }
が奇数のとき、 k 次の 主小 行列項が負になり 、が偶数のとき、正になる場合、 その行列は負定値行列です 。
k
{\displaystyle \ k\ }
k
{\displaystyle \ k\ }
が実正定値行列である 場合、任意の ベクトルに対して
M
{\displaystyle \ M\ }
m
{\displaystyle \ m\ }
v
,
{\displaystyle \ \mathbf {v} \ ,}
v
⊤
M
v
≥
m
‖
v
‖
2
2
.
{\displaystyle \ \mathbf {v} ^{\top }M\ \mathbf {v} \geq m\ \|\mathbf {v} \|_{2}^{\ \!2}~.}
エルミート行列が半正定値行列となるのは、その主小行列がすべて非負である場合に限ります。ただし、 0 と −1 の要素を持つ対角行列で確認されるように、主要な主小行列のみを考慮するだけでは不十分です。
ブロック行列と部分行列
正の行列は ブロック によって定義されることもあります 。
2
n
×
2
n
{\displaystyle \ 2n\times 2n\ }
M
=
[
A
B
C
D
]
{\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^{*}~.}
私たちは すべての複雑なもの 、特に その後
z
∗
M
z
≥
0
{\displaystyle \ \mathbf {z} ^{*}M\mathbf {z} \geq 0\ }
z
,
{\displaystyle \ \mathbf {z} \ ,}
z
=
[
v
,
0
]
⊤
.
{\displaystyle \ \mathbf {z} =[\mathbf {v} ,0]^{\top }~.}
[
v
∗
0
]
[
A
B
B
∗
D
]
[
v
0
]
=
v
∗
A
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\ ,}
A
{\displaystyle \ A\ }
D
{\displaystyle \ D\ }
M
{\displaystyle \ M\ }
逆の結果は、たとえば、 Schur 補集合 を使用して、ブロックに対するより強い条件で証明できます。
局所的極値
実変数 の 一般 二次形式は 常に次のように記述できます。 ここで はそれらの変数の列ベクトルであり、 は 対称実行列です。したがって、行列が正定値であるということは、がゼロ のときに一意の最小値(ゼロ)を持ち 、その他の場合は厳密に正であることを意味します。
f
(
x
)
{\displaystyle f(\mathbf {x} )}
n
{\displaystyle n}
x
1
,
…
,
x
n
{\displaystyle x_{1},\ldots ,x_{n}}
x
⊤
M
x
{\displaystyle \mathbf {x} ^{\top }M\mathbf {x} }
x
{\displaystyle \mathbf {x} }
M
{\displaystyle M}
f
{\displaystyle f}
x
{\displaystyle \mathbf {x} }
x
.
{\displaystyle \ \mathbf {x} ~.}
より一般的には、実変数 の 2 度微分可能な実関数は、 その 勾配 が 0 で、その ヘッセ行列 (すべての 2 次導関数の行列) がその点で半正定値である場合に、引数で局所最小値を持ちます。負定値行列と半正定値行列についても同様のことが言えます。
f
{\displaystyle f}
n
{\displaystyle n}
x
1
,
…
,
x
n
{\displaystyle x_{1},\ldots ,x_{n}}
共分散
統計学 では 、 多変量確率分布 の 共分散行列 は常に半正定値です。また、1 つの変数が他の変数の正確な線形関数でない限り、半正定値です。逆に、すべての半正定値行列は、何らかの多変量分布の共分散行列です。
非エルミート正方行列の拡張
正定値の定義は、任意の複素行列 (例えば、実非対称)を、 すべての非ゼロ複素ベクトルに対して が成り立つときに正定値であると指定することで一般化できる。ここ で、 は 複素数 の実部を表す [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)\ }
x
{\displaystyle \ \mathbf {x} \ }
M
{\displaystyle \ M\ }
x
⊤
M
x
>
0
{\displaystyle \ \mathbf {x} ^{\top }M\mathbf {x} >0\ }
x
{\displaystyle \ \mathbf {x} \ }
1
2
(
M
+
M
⊤
)
{\textstyle \ {\frac {1}{2}}\left(M+M^{\top }\right)\ }
x
⊤
M
x
=
∑
i
j
x
i
M
i
j
x
j
{\textstyle \ \mathbf {x} ^{\top }M\mathbf {x} =\sum _{ij}x_{i}M_{ij}x_{j}\ }
M
.
{\displaystyle \ M~.}
したがって、正の固有値のみを持つ非対称実行列は、正定値である必要はありません。たとえば、行列 は正の固有値を持ちます が、正定値ではありません。特に、 の負の値は ( の対称部分の負の固有値に関連付けられた固有ベクトル ) を選択することで得られます。
M
=
[
4
9
1
4
]
{\displaystyle M=\left[{\begin{smallmatrix}4&9\\1&4\end{smallmatrix}}\right]}
x
⊤
M
x
{\displaystyle \mathbf {x} ^{\top }M\mathbf {x} }
x
=
[
−
1
1
]
{\displaystyle \mathbf {x} =\left[{\begin{smallmatrix}-1\\1\end{smallmatrix}}\right]}
M
{\displaystyle M}
要約すると、実数と複素数の場合の違いは、複素ヒルベルト空間上の 有界正値演算子は必ずエルミート、つまり自己随伴であるという点です。一般的な主張は、 分極恒等式 を使用して論じることができます 。これは実数の場合ではもはや当てはまりません。
アプリケーション
熱伝導率マトリックス
熱伝導に関するフーリエの法則は、 温度勾配によって熱流束を表わし 、異方性媒体では次のように表されます。 ここで、 は 対称 熱伝導 行列です。フーリエの法則に負の値が挿入されているのは、熱は常に熱いところから冷たいところに流れるという予想を反映するためです。言い換えると、温度勾配は 常に冷たいところから熱いところを指しているので、熱流束は との内積が負になるはずなので、 となります 。 フーリエの法則を代入すると、この予想は となり 、伝導行列は正定値になるはずです。
q
{\displaystyle \ \mathbf {q} \ }
g
=
∇
T
{\displaystyle \ \mathbf {g} =\nabla T\ }
q
=
−
K
g
,
{\displaystyle \ \mathbf {q} =-K\mathbf {g} \ ,}
K
{\displaystyle \ K\ }
g
{\displaystyle \ \mathbf {g} \ }
q
{\displaystyle \ \mathbf {q} \ }
g
{\displaystyle \ \mathbf {g} \ }
q
⊤
g
<
0
.
{\displaystyle \ \mathbf {q} ^{\top }\mathbf {g} <0~.}
g
⊤
K
g
>
0
,
{\displaystyle \ \mathbf {g} ^{\top }K\mathbf {g} >0\ ,}
参照
参考文献
^ 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 日)。 凸最適化 。ケンブリッジ大学出版局。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と
k
=
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). Elsevier: 471– 506.
^ ホーン&ジョンソン(2013)、p.479、定理7.5.3
^ ホーン&ジョンソン(2013)、p.509、定理7.8.16
^ Styan, GP (1973). 「アダマール積と多変量統計解析」. 線形代数とその応用 . 6 : 217–240 . 、補題3.6、p. 227
^ Bhatia, Rajendra (2007). 正定値行列 . プリンストン、ニュージャージー: プリンストン大学出版局. p. 8. ISBN 978-0-691-12918-1 。
^ Weisstein, Eric W. 「正定値行列」 。MathWorld。Wolfram Research 。 2012年 7月26日 閲覧 。
出典
バーンスタイン、B.トゥーパン、RA (1962)。 「厳密な凸関数のヘッセ行列のいくつかの特性」。 数学に関するジャーナル 。 210 : 67– 72. 土井 :10.1515/crll.1962.210.65。
外部リンク
「正定値形式」、 数学百科事典 、 EMS Press 、2001 [1994]
「正定値行列」。 Wolfram MathWorld 。 Wolfram Research。