線形代数 の特定の用途では、 ランダム行列 の有限 和 の 最大 固有値の 確率分布 の特性を知っておくと便利です 。 はランダム行列の有限シーケンスであるとします。よく知られているスカラーの和の チェルノフ境界 と同様に、与えられたパラメータ t に対して次の境界が求められます 。
{
バツ
け
}
{\displaystyle \{\mathbf {X} _{k}\}}
広報
{
λ
最大
(
∑
け
バツ
け
)
≥
t
}
{\displaystyle \Pr \left\{\lambda _{\max }\left(\sum _{k}\mathbf {X} _{k}\right)\geq t\right\}}
以下の定理は、さまざまな仮定の下でこの一般的な質問に答えます。これらの仮定は、古典的なスカラーの対応物との類推によって以下で名前が付けられています。これらの定理はすべて、以下で導出される一般的な結果の特定の応用として (Tropp 2010) に記載されています。関連する研究の概要を示します。
行列ガウスとラーデマッハ級数
自己随伴行列の場合
次元 の固定された自己随伴行列の 有限列を考え 、 を独立した 標準 正規変数 または独立した ラデマッハ 乱数変数の有限列とします 。
{
あ
け
}
{\displaystyle \{\mathbf {A} _{k}\}}
d
{\displaystyle d}
{
ξ
け
}
{\displaystyle \{\xi _{k}\}}
そして、すべての に対して 、
t
≥
0
{\displaystyle t\geq 0}
広報
{
λ
最大
(
∑
け
ξ
け
あ
け
)
≥
t
}
≤
d
⋅
e
−
t
2
/
2
σ
2
{\displaystyle \Pr \left\{\lambda _{\text{max}}\left(\sum _{k}\xi _{k}\mathbf {A} _{k}\right)\geq t\right\}\leq d\cdot e^{-t^{2}/2\sigma ^{2}}}
どこ
σ
2
=
‖
∑
け
あ
け
2
‖
。
{\displaystyle \sigma^{2}={\bigg \Vert }\sum _{k}\mathbf {A} _{k}^{2}{\bigg \Vert }.}
長方形ケース
次元の固定行列の 有限列を考え 、 を独立した標準正規分布または独立したラデマッハ乱数変数の有限列とする。分散パラメータを定義する。
{
B
け
}
{\displaystyle \{\mathbf {B} _{k}\}}
d
1
×
d
2
{\displaystyle d_{1}\times d_{2}}
{
ξ
け
}
{\displaystyle \{\xi _{k}\}}
σ
2
=
最大
{
‖
∑
け
B
け
B
け
∗
‖
、
‖
∑
け
B
け
∗
B
け
‖
}
。
{\displaystyle \sigma ^{2}=\max \left\{{\bigg \Vert }\sum _{k}\mathbf {B} _{k}\mathbf {B} _{k}^{*}{\bigg \Vert },{\bigg \Vert }\sum _{k}\mathbf {B} _{k}^{*}\mathbf {B} _{k}{\bigg \Vert }\right\}.}
そして、すべての に対して 、
t
≥
0
{\displaystyle t\geq 0}
広報
{
‖
∑
け
ξ
け
B
け
‖
≥
t
}
≤
(
d
1
+
d
2
)
⋅
e
−
t
2
/
2
σ
2
。
{\displaystyle \Pr \left\{{\bigg \Vert }\sum _{k}\xi _{k}\mathbf {B} _{k}{\bigg \Vert }\geq t\right\}\leq (d_{1}+d_{2})\cdot e^{-t^{2}/2\sigma ^{2}}.}
行列チェルノフ不等式
古典的な チェルノフ境界は 、独立で非負で一様に境界付けられたランダム変数の合計に関係します。行列設定では、類似の定理は、一様な固有値境界に従う 半正定値 ランダム行列の合計に関係します。
マトリックス チェルノフ I
次元を持つ、独立したランダムな自己随伴行列の 有限列を考える 。各ランダム行列が次式を満たすと仮定する。
{
バツ
け
}
{\displaystyle \{\mathbf {X} _{k}\}}
d
{\displaystyle d}
バツ
け
⪰
0
そして
λ
最大
(
バツ
け
)
≤
R
{\displaystyle \mathbf {X} _{k}\succeq \mathbf {0} \quad {\text{and}}\quad \lambda _{\text{max}}(\mathbf {X} _{k})\leq R}
ほぼ確実です。
定義する
μ
分
=
λ
分
(
∑
け
え
バツ
け
)
そして
μ
最大
=
λ
最大
(
∑
け
え
バツ
け
)
。
{\displaystyle \mu _{\text{min}}=\lambda _{\text{min}}\left(\sum _{k}\mathbb {E} \,\mathbf {X} _{k}\right)\quad {\text{および}}\quad \mu _{\text{max}}=\lambda _{\text{max}}\left(\sum _{k}\mathbb {E} \,\mathbf {X} _{k}\right).}
それから
広報
{
λ
分
(
∑
け
バツ
け
)
≤
(
1
−
δ
)
μ
分
}
≤
d
⋅
[
e
−
δ
(
1
−
δ
)
1
−
δ
]
μ
分
/
R
のために
δ
∈
[
0
、
1
)
、 そして
{\displaystyle \Pr \left\{\lambda _{\text{min}}\left(\sum _{k}\mathbf {X} _{k}\right)\leq (1-\delta )\mu _{\text{min}}\right\}\leq d\cdot \left[{\frac {e^{-\delta }}{(1-\delta )^{1-\delta }}}\right]^{\mu _{\text{min}}/R}\quad {\text{for }}\delta \in [0,1){\text{, かつ}}}
広報
{
λ
最大
(
∑
け
バツ
け
)
≥
(
1
+
δ
)
μ
最大
}
≤
d
⋅
[
e
δ
(
1
+
δ
)
1
+
δ
]
μ
最大
/
R
のために
δ
≥
0.
{\displaystyle \Pr \left\{\lambda _{\text{max}}\left(\sum _{k}\mathbf {X} _{k}\right)\geq (1+\delta )\mu _{\text{max}}\right\}\leq d\cdot \left[{\frac {e^{\delta }}{(1+\delta )^{1+\delta }}}\right]^{\mu _{\text{max}}/R}\quad {\text{for }}\delta \geq 0.}
マトリックス チェルノフ II
次式を満たす、独立でランダムな自己随伴行列の
列を考える。
{
バツ
け
:
け
=
1
、
2
、
…
、
ん
}
{\displaystyle \{\mathbf {X} _{k}:k=1,2,\ldots ,n\}}
バツ
け
⪰
0
そして
λ
最大
(
バツ
け
)
≤
1
{\displaystyle \mathbf {X} _{k}\succeq \mathbf {0} \quad {\text{and}}\quad \lambda _{\text{max}}(\mathbf {X} _{k})\leq 1}
ほぼ確実です。
平均期待値の最小および最大固有値を計算する。
μ
¯
分
=
λ
分
(
1
ん
∑
け
=
1
ん
え
バツ
け
)
そして
μ
¯
最大
=
λ
最大
(
1
ん
∑
け
=
1
ん
え
バツ
け
)
。
{\displaystyle {\bar {\mu}}_{\text{min}}=\lambda _{\text{min}}\left({\frac {1}{n}}\sum _{k=1}^{n}\mathbb {E} \,\mathbf {X} _{k}\right)\quad {\text{および}}\quad {\bar {\mu}}_{\text{max}}=\lambda _{\text{max}}\left({\frac {1}{n}}\sum _{k=1}^{n}\mathbb {E} \,\mathbf {X} _{k}\right).}
それから
広報
{
λ
分
(
1
ん
∑
け
=
1
ん
バツ
け
)
≤
α
}
≤
d
⋅
e
−
ん
だ
(
α
‖
μ
¯
分
)
のために
0
≤
α
≤
μ
¯
分
、 そして
{\displaystyle \Pr \left\{\lambda _{\text{min}}\left({\frac {1}{n}}\sum _{k=1}^{n}\mathbf {X} _{k}\right)\leq \alpha \right\}\leq d\cdot e^{-nD(\alpha \Vert {\bar {\mu }}_{\text{min}})}\quad {\text{for }}0\leq \alpha \leq {\bar {\mu }}_{\text{min}}{\text{, and}}}
Pr
{
λ
max
(
1
n
∑
k
=
1
n
X
k
)
≥
α
}
≤
d
⋅
e
−
n
D
(
α
‖
μ
¯
max
)
for
μ
¯
max
≤
α
≤
1.
{\displaystyle \Pr \left\{\lambda _{\text{max}}\left({\frac {1}{n}}\sum _{k=1}^{n}\mathbf {X} _{k}\right)\geq \alpha \right\}\leq d\cdot e^{-nD(\alpha \Vert {\bar {\mu }}_{\text{max}})}\quad {\text{for }}{\bar {\mu }}_{\text{max}}\leq \alpha \leq 1.}
バイナリ情報分散は次のように定義される。
D
(
a
‖
u
)
=
a
(
log
a
−
log
u
)
+
(
1
−
a
)
(
log
(
1
−
a
)
−
log
(
1
−
u
)
)
{\displaystyle D(a\Vert u)=a\left(\log a-\log u\right)+(1-a)\left(\log(1-a)-\log(1-u)\right)}
のために 。
a
,
u
∈
[
0
,
1
]
{\displaystyle a,u\in [0,1]}
行列ベネット不等式とバーンスタイン不等式
スカラー設定では、 ベネット不等式とバーンスタイン不等式は、有界または 指数以下 の独立したゼロ平均ランダム変数の合計の上端を表します 。行列の場合、同様の結果はゼロ平均ランダム行列の合計に関係します。
境界ケース
次元を持つ、独立したランダムな自己随伴行列の 有限列を考える 。各ランダム行列が次式を満たすと仮定する。
{
X
k
}
{\displaystyle \{\mathbf {X} _{k}\}}
d
{\displaystyle d}
E
X
k
=
0
and
λ
max
(
X
k
)
≤
R
{\displaystyle \mathbb {E} \mathbf {X} _{k}=\mathbf {0} \quad {\text{and}}\quad \lambda _{\text{max}}(\mathbf {X} _{k})\leq R}
ほぼ確実です。
総分散のノルムを計算する。
σ
2
=
‖
∑
k
E
(
X
k
2
)
‖
.
{\displaystyle \sigma ^{2}={\bigg \Vert }\sum _{k}\mathbb {E} \,(\mathbf {X} _{k}^{2}){\bigg \Vert }.}
すると、次の不等式連鎖がすべての に対して成立します 。
t
≥
0
{\displaystyle t\geq 0}
Pr
{
λ
max
(
∑
k
X
k
)
≥
t
}
≤
d
⋅
exp
(
−
σ
2
R
2
⋅
h
(
R
t
σ
2
)
)
≤
d
⋅
exp
(
−
t
2
σ
2
+
R
t
/
3
)
≤
{
d
⋅
exp
(
−
3
t
2
/
8
σ
2
)
for
t
≤
σ
2
/
R
;
d
⋅
exp
(
−
3
t
/
8
R
)
for
t
≥
σ
2
/
R
.
{\displaystyle {\begin{aligned}\Pr \left\{\lambda _{\text{max}}\left(\sum _{k}\mathbf {X} _{k}\right)\geq t\right\}&\leq d\cdot \exp \left(-{\frac {\sigma ^{2}}{R^{2}}}\cdot h\left({\frac {Rt}{\sigma ^{2}}}\right)\right)\\&\leq d\cdot \exp \left({\frac {-t^{2}}{\sigma ^{2}+Rt/3}}\right)\\&\leq {\begin{cases}d\cdot \exp(-3t^{2}/8\sigma ^{2})\quad &{\text{for }}t\leq \sigma ^{2}/R;\\d\cdot \exp(-3t/8R)\quad &{\text{for }}t\geq \sigma ^{2}/R.\\\end{cases}}\end{aligned}}}
関数は に対して として定義されます 。
h
(
u
)
{\displaystyle h(u)}
h
(
u
)
=
(
1
+
u
)
log
(
1
+
u
)
−
u
{\displaystyle h(u)=(1+u)\log(1+u)-u}
u
≥
0
{\displaystyle u\geq 0}
指数関数的ではないケース
次元を持つ、独立したランダムな自己随伴行列の 有限列を考える 。
{
X
k
}
{\displaystyle \{\mathbf {X} _{k}\}}
d
{\displaystyle d}
E
X
k
=
0
and
E
(
X
k
p
)
⪯
p
!
2
⋅
R
p
−
2
A
k
2
{\displaystyle \mathbb {E} \,\mathbf {X} _{k}=\mathbf {0} \quad {\text{and}}\quad \mathbb {E} \,(\mathbf {X} _{k}^{p})\preceq {\frac {p!}{2}}\cdot R^{p-2}\mathbf {A} _{k}^{2}}
のために 。
p
=
2
,
3
,
4
,
…
{\displaystyle p=2,3,4,\ldots }
分散パラメータを計算する。
σ
2
=
‖
∑
k
A
k
2
‖
.
{\displaystyle \sigma ^{2}={\bigg \Vert }\sum _{k}\mathbf {A} _{k}^{2}{\bigg \Vert }.}
すると、次の不等式連鎖がすべての に対して成立します 。
t
≥
0
{\displaystyle t\geq 0}
Pr
{
λ
max
(
∑
k
X
k
)
≥
t
}
≤
d
⋅
exp
(
−
t
2
/
2
σ
2
+
R
t
)
≤
{
d
⋅
exp
(
−
t
2
/
4
σ
2
)
for
t
≤
σ
2
/
R
;
d
⋅
exp
(
−
t
/
4
R
)
for
t
≥
σ
2
/
R
.
{\displaystyle {\begin{aligned}\Pr \left\{\lambda _{\text{max}}\left(\sum _{k}\mathbf {X} _{k}\right)\geq t\right\}&\leq d\cdot \exp \left({\frac {-t^{2}/2}{\sigma ^{2}+Rt}}\right)\\&\leq {\begin{cases}d\cdot \exp(-t^{2}/4\sigma ^{2})\quad &{\text{for }}t\leq \sigma ^{2}/R;\\d\cdot \exp(-t/4R)\quad &{\text{for }}t\geq \sigma ^{2}/R.\\\end{cases}}\end{aligned}}}
長方形ケース
次元 の独立したランダム行列の 有限列を考える 。各ランダム行列が次式を満たすと仮定する。
{
Z
k
}
{\displaystyle \{\mathbf {Z} _{k}\}}
d
1
×
d
2
{\displaystyle d_{1}\times d_{2}}
E
Z
k
=
0
and
‖
Z
k
‖
≤
R
{\displaystyle \mathbb {E} \,\mathbf {Z} _{k}=\mathbf {0} \quad {\text{and}}\quad \Vert \mathbf {Z} _{k}\Vert \leq R}
ほぼ確実です。分散パラメータを定義する
σ
2
=
max
{
‖
∑
k
E
(
Z
k
Z
k
∗
)
‖
,
‖
∑
k
E
(
Z
k
∗
Z
k
)
‖
}
.
{\displaystyle \sigma ^{2}=\max \left\{{\bigg \Vert }\sum _{k}\mathbb {E} \,(\mathbf {Z} _{k}\mathbf {Z} _{k}^{*}){\bigg \Vert },{\bigg \Vert }\sum _{k}\mathbb {E} \,(\mathbf {Z} _{k}^{*}\mathbf {Z} _{k}){\bigg \Vert }\right\}.}
そして、すべての
t
≥
0
{\displaystyle t\geq 0}
Pr
{
‖
∑
k
Z
k
‖
≥
t
}
≤
(
d
1
+
d
2
)
⋅
exp
(
−
t
2
/
2
σ
2
+
R
t
/
3
)
{\displaystyle \Pr \left\{{\bigg \Vert }\sum _{k}\mathbf {Z} _{k}{\bigg \Vert }\geq t\right\}\leq (d_{1}+d_{2})\cdot \exp \left({\frac {-t^{2}/2}{\sigma ^{2}+Rt/3}}\right)}
保持する。 [1]
東、ホエフディング、マクダーミドの行列不等式
マトリックスあずま
アズマの不等式 のスカラー版は 、スカラー マルチンゲールが 平均値を中心に正規集中を示し、偏差のスケールは差分シーケンスの合計最大二乗範囲によって制御されると述べています。以下は、行列設定での拡張です。
次元 の自己随伴行列の 有限適応列と 、次式 を満たす自己随伴行列の
固定列を考える。
{
X
k
}
{\displaystyle \{\mathbf {X} _{k}\}}
d
{\displaystyle d}
{
A
k
}
{\displaystyle \{\mathbf {A} _{k}\}}
E
k
−
1
X
k
=
0
and
X
k
2
⪯
A
k
2
{\displaystyle \mathbb {E} _{k-1}\,\mathbf {X} _{k}=\mathbf {0} \quad {\text{and}}\quad \mathbf {X} _{k}^{2}\preceq \mathbf {A} _{k}^{2}}
ほぼ確実です。
分散パラメータを計算する
σ
2
=
‖
∑
k
A
k
2
‖
.
{\displaystyle \sigma ^{2}={\bigg \Vert }\sum _{k}\mathbf {A} _{k}^{2}{\bigg \Vert }.}
そして、すべての
t
≥
0
{\displaystyle t\geq 0}
Pr
{
λ
max
(
∑
k
X
k
)
≥
t
}
≤
d
⋅
e
−
t
2
/
8
σ
2
{\displaystyle \Pr \left\{\lambda _{\text{max}}\left(\sum _{k}\mathbf {X} _{k}\right)\geq t\right\}\leq d\cdot e^{-t^{2}/8\sigma ^{2}}}
追加情報が利用できる場合、定数 1/8 は 1/2 に改善できます。1 つのケースは、各加数が条件付きで対称である場合に発生します。別の例では、 が とほぼ確実に交換する という仮定が必要です 。
X
k
{\displaystyle \mathbf {X} _{k}}
X
k
{\displaystyle \mathbf {X} _{k}}
A
k
{\displaystyle \mathbf {A} _{k}}
マトリックス・ホーフディング
東行列の加数が独立であるという加法仮定を置くと、 Hoeffding の不等式 の行列拡張が得られます 。
次元 の独立したランダムな自己随伴行列の 有限列を考え 、 を 固定された自己随伴行列の列とする。各ランダム行列が次式を満たすと仮定する。
{
X
k
}
{\displaystyle \{\mathbf {X} _{k}\}}
d
{\displaystyle d}
{
A
k
}
{\displaystyle \{\mathbf {A} _{k}\}}
E
X
k
=
0
and
X
k
2
⪯
A
k
2
{\displaystyle \mathbb {E} \,\mathbf {X} _{k}=\mathbf {0} \quad {\text{and}}\quad \mathbf {X} _{k}^{2}\preceq \mathbf {A} _{k}^{2}}
ほぼ確実です。
そして、すべての
t
≥
0
{\displaystyle t\geq 0}
Pr
{
λ
max
(
∑
k
X
k
)
≥
t
}
≤
d
⋅
e
−
t
2
/
8
σ
2
{\displaystyle \Pr \left\{\lambda _{\text{max}}\left(\sum _{k}\mathbf {X} _{k}\right)\geq t\right\}\leq d\cdot e^{-t^{2}/8\sigma ^{2}}}
どこ
σ
2
=
‖
∑
k
A
k
2
‖
.
{\displaystyle \sigma ^{2}={\bigg \Vert }\sum _{k}\mathbf {A} _{k}^{2}{\bigg \Vert }.}
この結果の改良は(Mackey et al. 2012)で確立されました:すべての
t
≥
0
{\displaystyle t\geq 0}
Pr
{
λ
max
(
∑
k
X
k
)
≥
t
}
≤
d
⋅
e
−
t
2
/
2
σ
2
{\displaystyle \Pr \left\{\lambda _{\text{max}}\left(\sum _{k}\mathbf {X} _{k}\right)\geq t\right\}\leq d\cdot e^{-t^{2}/2\sigma ^{2}}}
どこ
σ
2
=
1
2
‖
∑
k
A
k
2
+
E
X
k
2
‖
≤
‖
∑
k
A
k
2
‖
.
{\displaystyle \sigma ^{2}={\frac {1}{2}}{\bigg \Vert }\sum _{k}\mathbf {A} _{k}^{2}+\mathbb {E} \,\mathbf {X} _{k}^{2}{\bigg \Vert }\leq {\bigg \Vert }\sum _{k}\mathbf {A} _{k}^{2}{\bigg \Vert }.}
行列の境界差 (マクダーミッド)
スカラー設定では、 マクダーミッド不等式は、 アズマ不等式を Doob マルチンゲール に 適用することにより、差を境界付ける一般的な方法の 1 つを提供します 。境界付き差不等式のバージョンは、行列設定でも保持されます。
を独立したランダム変数の族とし 、を 変数を次元 の自己随伴行列に 写像する関数とする。 次式を満たす固定自己随伴行列の
列を考える。
{
Z
k
:
k
=
1
,
2
,
…
,
n
}
{\displaystyle \{Z_{k}:k=1,2,\ldots ,n\}}
H
{\displaystyle \mathbf {H} }
n
{\displaystyle n}
d
{\displaystyle d}
{
A
k
}
{\displaystyle \{\mathbf {A} _{k}\}}
(
H
(
z
1
,
…
,
z
k
,
…
,
z
n
)
−
H
(
z
1
,
…
,
z
k
′
,
…
,
z
n
)
)
2
⪯
A
k
2
,
{\displaystyle \left(\mathbf {H} (z_{1},\ldots ,z_{k},\ldots ,z_{n})-\mathbf {H} (z_{1},\ldots ,z'_{k},\ldots ,z_{n})\right)^{2}\preceq \mathbf {A} _{k}^{2},}
ここで 、 各インデックスに対する の すべての可能な値の範囲 。分散パラメータを計算する。
z
i
{\displaystyle z_{i}}
z
i
′
{\displaystyle z'_{i}}
Z
i
{\displaystyle Z_{i}}
i
{\displaystyle i}
σ
2
=
‖
∑
k
A
k
2
‖
.
{\displaystyle \sigma ^{2}={\bigg \Vert }\sum _{k}\mathbf {A} _{k}^{2}{\bigg \Vert }.}
そして、すべての
t
≥
0
{\displaystyle t\geq 0}
Pr
{
λ
max
(
H
(
z
)
−
E
H
(
z
)
)
≥
t
}
≤
d
⋅
e
−
t
2
/
8
σ
2
,
{\displaystyle \Pr \left\{\lambda _{\text{max}}\left(\mathbf {H} (\mathbf {z} )-\mathbb {E} \,\mathbf {H} (\mathbf {z} )\right)\geq t\right\}\leq d\cdot e^{-t^{2}/8\sigma ^{2}},}
どこ 。
z
=
(
Z
1
,
…
,
Z
n
)
{\displaystyle \mathbf {z} =(Z_{1},\ldots ,Z_{n})}
この結果の改良は(Paulin、Mackey & Tropp 2013)で確立されました((Paulin、Mackey & Tropp 2016)も参照)。
t
≥
0
{\displaystyle t\geq 0}
Pr
{
λ
max
(
H
(
z
)
−
E
H
(
z
)
)
≥
t
}
≤
d
⋅
e
−
t
2
/
σ
2
,
{\displaystyle \Pr \left\{\lambda _{\text{max}}\left(\mathbf {H} (\mathbf {z} )-\mathbb {E} \,\mathbf {H} (\mathbf {z} )\right)\geq t\right\}\leq d\cdot e^{-t^{2}/\sigma ^{2}},}
どこで そして
z
=
(
Z
1
,
…
,
Z
n
)
{\displaystyle \mathbf {z} =(Z_{1},\ldots ,Z_{n})}
σ
2
=
‖
∑
k
A
k
2
‖
.
{\displaystyle \sigma ^{2}={\bigg \Vert }\sum _{k}\mathbf {A} _{k}^{2}{\bigg \Vert }.}
このタイプの境界は、(Ahlswede & Winter 2003)によって初めて導かれました。自己随伴行列のガウス境界とラデマッハ境界に関する上記の定理を思い出してください。次元を 持つ固定された自己随伴行列の有限列と、 独立した 標準 正規 または独立した ラデマッハ ランダム変数の有限列 に対して 、
{
A
k
}
{\displaystyle \{\mathbf {A} _{k}\}}
d
{\displaystyle d}
{
ξ
k
}
{\displaystyle \{\xi _{k}\}}
Pr
{
λ
max
(
∑
k
ξ
k
A
k
)
≥
t
}
≤
d
⋅
e
−
t
2
/
2
σ
2
{\displaystyle \Pr \left\{\lambda _{\text{max}}\left(\sum _{k}\xi _{k}\mathbf {A} _{k}\right)\geq t\right\}\leq d\cdot e^{-t^{2}/2\sigma ^{2}}}
どこ
σ
2
=
‖
∑
k
A
k
2
‖
.
{\displaystyle \sigma ^{2}={\bigg \Vert }\sum _{k}\mathbf {A} _{k}^{2}{\bigg \Vert }.}
アールスヴェーデとウィンターは、同じ結果を出すだろうが、
σ
A
W
2
=
∑
k
λ
max
(
A
k
2
)
{\displaystyle \sigma _{AW}^{2}=\sum _{k}\lambda _{\max }\left(\mathbf {A} _{k}^{2}\right)}
。
比較すると、 上の定理の は と と可換です 。つまり、これは最大固有値の和ではなく、和の最大固有値です。これは、アールスウェーデ・ウィンター値 ( ノルム 三角不等式 より) よりも大きくなることはありませんが、はるかに小さくなる可能性があります。したがって、上の定理は、アールスウェーデ・ウィンターの結果よりも厳しい境界を与えます。
σ
2
{\displaystyle \sigma ^{2}}
Σ
{\displaystyle \Sigma }
λ
max
{\displaystyle \lambda _{\max }}
(Ahlswede & Winter 2003) の主な貢献は、スカラー チェルノフ境界 ( チェルノフ境界#加法形式 (絶対誤差) を 参照) を証明するために使用されるラプラス変換法を自己随伴行列の場合に拡張したことです。以下の導出で手順を示します。このトピックに関する最近の研究はすべてこの同じ手順に従っており、主な違いは後続のステップから生じます。Ahlswede & Winter は ゴールデン–トムソン不等式を 使用して処理を進めますが、Tropp (Tropp 2010) は リーブの定理 を使用します。
右辺をほぼ一定に保ちながら、 級数の長さ ( n ) と行列の次元 ( d ) を変化させたいとします。この場合、 n は d の対数にほぼ比例して変化する必要があります。いくつかの論文では、次元に依存せずに境界を確立しようと試みられています。Rudelson と Vershynin (Rudelson & Vershynin 2007) は、2 つのベクトルの外積である行列の結果を示しています。 (Magen & Zouzias 2010) は、 低ランク行列について次元に依存しない結果を示しています 。元の結果は、Ahlswede–Winter アプローチとは独立して導き出されましたが、(Oliveira 2010b) は、Ahlswede–Winter アプローチを使用して同様の結果を証明しています。
最後に、Oliveira (Oliveira 2010a) は、Ahlswede–Winter フレームワークとは独立して、行列マルチンゲールの結果を証明しています。Tropp (Tropp 2011) は、Ahlswede–Winter フレームワークを使用して結果をわずかに改善しています。どちらの結果もこの記事では示されていません。
導出と証明
アールスヴェーデと冬
(Ahlswede & Winter 2003) で発見されたラプラス変換の議論は、それ自体が重要な結果である。 ランダムな自己随伴行列をとろう。すると、
Y
{\displaystyle \mathbf {Y} }
Pr
{
λ
max
(
Y
)
≥
t
}
≤
inf
θ
>
0
{
e
−
θ
t
⋅
E
[
tr
e
θ
Y
]
}
.
{\displaystyle \Pr \left\{\lambda _{\max }(Y)\geq t\right\}\leq \inf _{\theta >0}\left\{e^{-\theta t}\cdot \operatorname {E} \left[\operatorname {tr} e^{\theta \mathbf {Y} }\right]\right\}.}
これを証明するには、 を固定します 。すると
θ
>
0
{\displaystyle \theta >0}
Pr
{
λ
max
(
Y
)
≥
t
}
=
Pr
{
λ
max
(
θ
Y
)
≥
θ
t
}
=
Pr
{
e
λ
max
(
θ
Y
)
≥
e
θ
t
}
≤
e
−
θ
t
E
e
λ
max
(
θ
Y
)
≤
e
−
θ
t
E
tr
e
(
θ
Y
)
{\displaystyle {\begin{aligned}\Pr \left\{\lambda _{\max }(\mathbf {Y} )\geq t\right\}&=\Pr \left\{\lambda _{\max }(\mathbf {\theta Y} )\geq \theta t\right\}\\&=\Pr \left\{e^{\lambda _{\max }(\theta \mathbf {Y} )}\geq e^{\theta t}\right\}\\&\leq e^{-\theta t}\operatorname {E} e^{\lambda _{\max }(\theta \mathbf {Y} )}\\&\leq e^{-\theta t}\operatorname {E} \operatorname {tr} e^{(\theta \mathbf {Y} )}\end{aligned}}}
最後から 2 番目の不等式は マルコフの不等式 です。最後の不等式は であるため成り立ちます 。左端の量は とは独立しているため 、 上の最小値は の上限のままです。
e
λ
max
(
θ
Y
)
=
λ
max
(
e
θ
Y
)
≤
tr
(
e
θ
Y
)
{\displaystyle e^{\lambda _{\max }(\theta \mathbf {Y} )}=\lambda _{\max }(e^{\theta \mathbf {Y} })\leq \operatorname {tr} (e^{\theta \mathbf {Y} })}
θ
{\displaystyle \theta }
θ
>
0
{\displaystyle \theta >0}
したがって、私たちの課題は、 を理解することです 。ただし、トレースと期待値はどちらも線形なので、それらを交換することができ、 を考慮するだけで十分です 。これを行列生成関数と呼びます。ここが、(Ahlswede & Winter 2003) と (Tropp 2010) の方法が分岐するところです。次のプレゼンテーションは、(Ahlswede & Winter 2003) に続きます。
E
[
tr
(
e
θ
Y
)
]
{\displaystyle \operatorname {E} [\operatorname {tr} (e^{\theta \mathbf {Y} })]}
E
e
θ
Y
:=
M
Y
(
θ
)
{\displaystyle \operatorname {E} e^{\theta \mathbf {Y} }:=\mathbf {M} _{\mathbf {Y} }(\theta )}
ゴールデン ・トムソン不等式 は、
tr
M
X
1
+
X
2
(
θ
)
≤
tr
[
(
E
e
θ
X
1
)
(
E
e
θ
X
2
)
]
=
tr
M
X
1
(
θ
)
M
X
2
(
θ
)
{\displaystyle \operatorname {tr} \mathbf {M} _{\mathbf {X} _{1}+\mathbf {X} _{2}}(\theta )\leq \operatorname {tr} \left[\left(\operatorname {E} e^{\theta \mathbf {X} _{1}}\right)\left(\operatorname {E} e^{\theta \mathbf {X} _{2}}\right)\right]=\operatorname {tr} \mathbf {M} _{\mathbf {X} _{1}}(\theta )\mathbf {M} _{\mathbf {X} _{2}}(\theta )}
ここでは、期待値の線形性を何度か使用しました。
と仮定します 。この結果を繰り返すことで の上限を見つけることができます 。 であることに注目すると、
Y
=
∑
k
X
k
{\displaystyle \mathbf {Y} =\sum _{k}\mathbf {X} _{k}}
tr
M
Y
(
θ
)
{\displaystyle \operatorname {tr} \mathbf {M} _{\mathbf {Y} }(\theta )}
tr
(
A
B
)
≤
tr
(
A
)
λ
max
(
B
)
{\displaystyle \operatorname {tr} (\mathbf {AB} )\leq \operatorname {tr} (\mathbf {A} )\lambda _{\max }(\mathbf {B} )}
tr
M
Y
(
θ
)
≤
tr
[
(
E
e
∑
k
=
1
n
−
1
θ
X
k
)
(
E
e
θ
X
n
)
]
≤
tr
(
E
e
∑
k
=
1
n
−
1
θ
X
k
)
λ
max
(
E
e
θ
X
n
)
.
{\displaystyle \operatorname {tr} \mathbf {M} _{\mathbf {Y} }(\theta )\leq \operatorname {tr} \left[\left(\operatorname {E} e^{\sum _{k=1}^{n-1}\theta \mathbf {X} _{k}}\right)\left(\operatorname {E} e^{\theta \mathbf {X} _{n}}\right)\right]\leq \operatorname {tr} \left(\operatorname {E} e^{\sum _{k=1}^{n-1}\theta \mathbf {X} _{k}}\right)\lambda _{\max }(\operatorname {E} e^{\theta \mathbf {X} _{n}}).}
これを繰り返すと、
tr
M
Y
(
θ
)
≤
(
tr
I
)
[
Π
k
λ
max
(
E
e
θ
X
k
)
]
=
d
e
∑
k
λ
max
(
log
E
e
θ
X
k
)
{\displaystyle \operatorname {tr} \mathbf {M} _{\mathbf {Y} }(\theta )\leq (\operatorname {tr} \mathbf {I} )\left[\Pi _{k}\lambda _{\max }(\operatorname {E} e^{\theta \mathbf {X} _{k}})\right]=de^{\sum _{k}\lambda _{\max }\left(\log \operatorname {E} e^{\theta \mathbf {X} _{k}}\right)}}
これまで、 上の下限値を持つ境界を見つけました 。今度は、これを境界で囲むことができます。いずれにせよ、アールスウェーデ-ウィンター境界が最大固有値の合計としてどのように生じるかがわかります。
θ
{\displaystyle \theta }
トロップ
(Tropp 2010) の主な貢献は、 (Ahlswede & Winter 2003) が ゴールデン・トムソン不等式 を適用したのに対し、 リープの定理 を適用したことです。Tropp の系は次のようになります。 が固定された自己随伴行列で が ランダムな自己随伴行列である場合、
H
{\displaystyle H}
X
{\displaystyle X}
E
tr
e
H
+
X
≤
tr
e
H
+
log
(
E
e
X
)
{\displaystyle \operatorname {E} \operatorname {tr} e^{\mathbf {H} +\mathbf {X} }\leq \operatorname {tr} e^{\mathbf {H} +\log(\operatorname {E} e^{\mathbf {X} })}}
証明: とする 。するとリープの定理は次を言う。
Y
=
e
X
{\displaystyle \mathbf {Y} =e^{\mathbf {X} }}
f
(
Y
)
=
tr
e
H
+
log
(
Y
)
{\displaystyle f(\mathbf {Y} )=\operatorname {tr} e^{\mathbf {H} +\log(\mathbf {Y} )}}
は凹です。最後のステップは、 ジェンセンの不等式 を使用して関数内の期待値を移動することです。
E
tr
e
H
+
log
(
Y
)
≤
tr
e
H
+
log
(
E
Y
)
.
{\displaystyle \operatorname {E} \operatorname {tr} e^{\mathbf {H} +\log(\mathbf {Y} )}\leq \operatorname {tr} e^{\mathbf {H} +\log(\operatorname {E} \mathbf {Y} )}.}
これにより、本論文の主要な結果、すなわち行列生成関数の対数の劣加法性が得られます。
log mgf の劣加法性
を独立したランダム自己随伴行列の有限列と します。すると、すべての に対して 、
X
k
{\displaystyle \mathbf {X} _{k}}
θ
∈
R
{\displaystyle \theta \in \mathbb {R} }
tr
M
∑
k
X
k
(
θ
)
≤
tr
e
∑
k
log
M
X
k
(
θ
)
{\displaystyle \operatorname {tr} \mathbf {M} _{\sum _{k}\mathbf {X} _{k}}(\theta )\leq \operatorname {tr} e^{\sum _{k}\log \mathbf {M} _{\mathbf {X} _{k}}(\theta )}}
証明: とすれば十分である 。定義を拡張して、次のことを示す必要がある。
θ
=
1
{\displaystyle \theta =1}
E
tr
e
∑
k
θ
X
k
≤
tr
e
∑
k
log
E
e
θ
X
k
.
{\displaystyle \operatorname {E} \operatorname {tr} e^{\sum _{k}\theta \mathbf {X} _{k}}\leq \operatorname {tr} e^{\sum _{k}\log \operatorname {E} e^{\theta \mathbf {X} _{k}}}.}
証明を完了するために、 全期待値の法則 を使用します。 を を条件とする期待値とします 。 はすべて 独立していると仮定するため、
E
k
{\displaystyle \operatorname {E} _{k}}
X
1
,
…
,
X
k
{\displaystyle \mathbf {X} _{1},\ldots ,\mathbf {X} _{k}}
X
i
{\displaystyle \mathbf {X} _{i}}
E
k
−
1
e
X
k
=
E
e
X
k
.
{\displaystyle \operatorname {E} _{k-1}e^{\mathbf {X} _{k}}=\operatorname {E} e^{\mathbf {X} _{k}}.}
定義する 。
Ξ
k
=
log
E
k
−
1
e
X
k
=
log
M
X
k
(
θ
)
{\displaystyle \mathbf {\Xi } _{k}=\log \operatorname {E} _{k-1}e^{\mathbf {X} _{k}}=\log \mathbf {M} _{\mathbf {X} _{k}}(\theta )}
最後に、
E
tr
e
∑
k
=
1
n
X
k
=
E
0
⋯
E
n
−
1
tr
e
∑
k
=
1
n
−
1
X
k
+
X
n
≤
E
0
⋯
E
n
−
2
tr
e
∑
k
=
1
n
−
1
X
k
+
log
(
E
n
−
1
e
X
n
)
=
E
0
⋯
E
n
−
2
tr
e
∑
k
=
1
n
−
2
X
k
+
X
n
−
1
+
Ξ
n
⋮
=
tr
e
∑
k
=
1
n
Ξ
k
{\displaystyle {\begin{aligned}\operatorname {E} \operatorname {tr} e^{\sum _{k=1}^{n}\mathbf {X} _{k}}&=\operatorname {E} _{0}\cdots \operatorname {E} _{n-1}\operatorname {tr} e^{\sum _{k=1}^{n-1}\mathbf {X} _{k}+\mathbf {X} _{n}}\\&\leq \operatorname {E} _{0}\cdots \operatorname {E} _{n-2}\operatorname {tr} e^{\sum _{k=1}^{n-1}\mathbf {X} _{k}+\log(\operatorname {E} _{n-1}e^{\mathbf {X} _{n}})}\\&=\operatorname {E} _{0}\cdots \operatorname {E} _{n-2}\operatorname {tr} e^{\sum _{k=1}^{n-2}\mathbf {X} _{k}+\mathbf {X} _{n-1}+\mathbf {\Xi } _{n}}\\&\vdots \\&=\operatorname {tr} e^{\sum _{k=1}^{n}\mathbf {\Xi } _{k}}\end{aligned}}}
ここで、各ステップmにおいて、Troppの系を用いて、
H
m
=
∑
k
=
1
m
−
1
X
k
+
∑
k
=
m
+
1
n
Ξ
k
{\displaystyle \mathbf {H} _{m}=\sum _{k=1}^{m-1}\mathbf {X} _{k}+\sum _{k=m+1}^{n}\mathbf {\Xi } _{k}}
マスターテールバウンド
前回の結果から次のことが分かります。
Pr
{
λ
max
(
∑
k
X
k
)
≥
t
}
≤
inf
θ
>
0
{
e
−
θ
t
tr
e
∑
k
log
M
X
k
(
θ
)
}
{\displaystyle \Pr \left\{\lambda _{\max }\left(\sum _{k}\mathbf {X} _{k}\right)\geq t\right\}\leq \inf _{\theta >0}\left\{e^{-\theta t}\operatorname {tr} e^{\sum _{k}\log \mathbf {M} _{\mathbf {X} _{k}}(\theta )}\right\}}
上で示した定理はすべてこの境界から導き出されたもので、定理はさまざまな方法で下限値を制限します。これらの手順は、示した証明よりもはるかに簡単です。
参考文献
^ ランダム行列の和に対するユーザーフレンドリーな末尾境界
Ahlswede, R.; Winter, A. (2003). 「量子チャネルによる識別のための強い逆」. IEEE Transactions on Information Theory . 48 (3): 569–579. arXiv : quant-ph/0012127 . doi :10.1109/18.985947. S2CID 523176.
Mackey, L.; Jordan, MI; Chen, RY; Farrell, B .; Tropp, JA (2012 ) . 「 交換 可能 ペア法によるマトリックス濃度不等式」。 確率年報 。42 (3): 906–945。arXiv : 1201.6002。doi :10.1214/13-AOP892。S2CID 9635314。
Magen, A. ; Zouzias, A. (2010). 「低ランク行列値チェルノフ境界と近似行列乗算」. arXiv : 1005.2724 [cs.DS].
Oliveira, RI (2010a). 「独立エッジを持つランダムグラフにおける隣接行列とラプラシアンの集中」 arXiv : 0911.0600 [math.CO].
Oliveira, RI (2010b). 「Rudelson によるランダム エルミート行列の和と不等式」. arXiv : 1004.3821 [math.PR].
Paulin, D.; Mackey, L.; Tropp, JA (2013). 「カーネル結合からの行列濃度不等式の導出」. arXiv : 1305.0612 [math.PR].
Paulin, D.; Mackey, L.; Tropp, JA (2016). 「ランダム行列のエフロン–スタイン不等式」. The Annals of Probability . 44 (5): 3431–3473. arXiv : 1408.3470 . doi :10.1214/15-AOP1054. S2CID 16263460.
Rudelson, M.; Vershynin, R. (2007). 「大規模行列からのサンプリング: 幾何関数解析によるアプローチ」 Journal of the Association for Computing Machinery . 54 (4 ed.). arXiv : math/9608208 . Bibcode :1996math......8208R. doi :10.1145/1255443.1255449. S2CID 6054789.
Tropp, J. (2011). 「行列マルチンゲールに対するフリードマンの不等式」. arXiv : 1101.3039 [math.PR].
Tropp, J. (2010). 「ランダム行列の合計に対するユーザーフレンドリーな末尾境界」. 計算数学の基礎 . 12 (4): 389–434. arXiv : 1004.4389 . doi :10.1007/s10208-011-9099-z. S2CID 17735965.