3
大学数学基礎解説
文献あり

det(I+AB)の対称的な表示

347
0
$$\newcommand{adari}[0]{\mathrm{adari}} \newcommand{adari}[0]{\mathrm{adari}} \newcommand{adgari}[0]{\mathrm{adgari}} \newcommand{al}[0]{\mathrm{al}} \newcommand{amit}[0]{\mathrm{amit}} \newcommand{anit}[0]{\boldsymbol{anit}} \newcommand{anit}[0]{\mathrm{anit}} \newcommand{answamu}[0]{\mathrm{answamu}} \newcommand{anti}[0]{\mathrm{anti}} \newcommand{ari}[0]{\mathrm{ari}} \newcommand{ARI}[0]{\mathrm{ARI}} \newcommand{ARI}[0]{\mathrm{ARI}} \newcommand{arit}[0]{\mathrm{arit}} \newcommand{as}[0]{\mathrm{as}} \newcommand{axi}[0]{\mathrm{axi}} \newcommand{axit}[0]{\mathrm{axit}} \newcommand{ba}[0]{\boldsymbol{a}} \newcommand{bb}[0]{\boldsymbol{b}} \newcommand{bc}[0]{\boldsymbol{c}} \newcommand{bd}[0]{\boldsymbol{d}} \newcommand{be}[0]{\boldsymbol{e}} \newcommand{bk}[0]{\boldsymbol{k}} \newcommand{bl}[0]{\boldsymbol{l}} \newcommand{BQ}[5]{{}_{#1}\psi_{#2}\left[\begin{matrix}#3\\#4\end{matrix};#5\right]} \newcommand{bw}[0]{\boldsymbol{w}} \newcommand{bx}[0]{\boldsymbol{x}} \newcommand{by}[0]{\boldsymbol{y}} \newcommand{calA}[0]{\mathcal{A}} \newcommand{calS}[0]{\mathcal{S}} \newcommand{CC}[0]{\mathbb{C}} \newcommand{crash}[0]{\mathrm{crash}} \newcommand{der}[0]{\mathrm{der}} \newcommand{DIFF}[0]{\mathrm{DIFF}} \newcommand{EE}[0]{\mathfrak{E}} \newcommand{Eneg}[0]{\mathfrak{E}\text{-}\mathrm{neg}} \newcommand{Enegpush}[0]{\mathfrak{E}\text{-}\mathrm{negpush}} \newcommand{Epush}[0]{\mathfrak{E}\text{-}\mathrm{push}} \newcommand{es}[0]{\mathfrak{es}} \newcommand{Esena}[0]{\mathfrak{E}\text{-}\mathrm{sena}} \newcommand{ess}[0]{\mathfrak{ess}} \newcommand{Eswap}[0]{\mathfrak{E}\text{-}\swap} \newcommand{Eter}[0]{\mathfrak{E}\text{-}\mathrm{ter}} \newcommand{expari}[0]{\mathrm{expari}} \newcommand{ez}[0]{\mathfrak{ez}} \newcommand{F}[5]{{}_{#1}F_{#2}\left[\begin{matrix}#3\\#4\end{matrix};#5\right]} \newcommand{fragari}[0]{\mathrm{fragari}} \newcommand{fragira}[0]{\mathrm{fragira}} \newcommand{gami}[0]{\mathrm{gami}} \newcommand{gamit}[0]{\mathrm{gamit}} \newcommand{gani}[0]{\mathrm{gani}} \newcommand{ganit}[0]{\mathrm{ganit}} \newcommand{gantar}[0]{\mathrm{gantar}} \newcommand{gari}[0]{\mathrm{gari}} \newcommand{GARI}[0]{\mathrm{GARI}} \newcommand{GARI}[0]{\mathrm{GARI}} \newcommand{garit}[0]{\mathrm{garit}} \newcommand{gaxi}[0]{\mathrm{gaxi}} \newcommand{gaxit}[0]{\mathrm{gaxit}} \newcommand{gepar}[0]{\mathrm{gepar}} \newcommand{GIFF}[0]{\mathrm{GIFF}} \newcommand{gira}[0]{\mathrm{gira}} \newcommand{girat}[0]{\mathrm{girat}} \newcommand{gush}[0]{\mathrm{gush}} \newcommand{H}[5]{{}_{#1}H_{#2}\left[\begin{matrix}#3\\#4\end{matrix};#5\right]} \newcommand{He}[0]{\mathfrak{He}} \newcommand{inv}[0]{\mathrm{inv}} \newcommand{invgami}[0]{\mathrm{invgami}} \newcommand{invgani}[0]{\mathrm{invgani}} \newcommand{invgari}[0]{\mathrm{invgari}} \newcommand{invgaxi}[0]{\mathrm{invgaxi}} \newcommand{invgira}[0]{\mathrm{invgira}} \newcommand{invmu}[0]{\mathrm{invmu}} \newcommand{ira}[0]{\mathrm{ira}} \newcommand{irat}[0]{\mathrm{irat}} \newcommand{iwat}[0]{\mathrm{iwat}} \newcommand{lu}[0]{\mathrm{lu}} \newcommand{LU}[0]{\mathrm{LU}} \newcommand{maj}[0]{\mathrm{maj}} \newcommand{mantar}[0]{\mathrm{mantar}} \newcommand{MU}[0]{\mathrm{MU}} \newcommand{neg}[0]{\mathrm{neg}} \newcommand{ol}[0]{\overline} \newcommand{Omantar}[0]{\mathfrak{O}\text{-}\mathrm{mantar}} \newcommand{OO}[0]{\mathfrak{O}} \newcommand{os}[0]{\mathfrak{os}} \newcommand{oss}[0]{\mathfrak{oss}} \newcommand{oz}[0]{\mathfrak{oz}} \newcommand{pari}[0]{\mathrm{pari}} \newcommand{preari}[0]{\mathrm{preari}} \newcommand{preira}[0]{\mathrm{preira}} \newcommand{pus}[0]{\mathrm{pus}} \newcommand{push}[0]{\mathrm{push}} \newcommand{pusnu}[0]{\mathrm{pusnu}} \newcommand{Q}[5]{{}_{#1}\phi_{#2}\left[\begin{matrix}#3\\#4\end{matrix};#5\right]} \newcommand{QQ}[0]{\mathbb{Q}} \newcommand{ras}[0]{\mathrm{ras}} \newcommand{rash}[0]{\mathrm{rash}} \newcommand{re}[0]{\mathfrak{re}} \newcommand{ro}[0]{\mathfrak{r\ddot{o}}} \newcommand{Se}[0]{\mathfrak{Se}} \newcommand{sgn}[0]{\mathrm{sgn}} \newcommand{sh}[0]{\,\text{ш}\,} \newcommand{So}[0]{\mathfrak{S\ddot{o}}} \newcommand{swamu}[0]{\mathrm{swamu}} \newcommand{swap}[0]{\mathrm{swap}} \newcommand{To}[0]{\mathfrak{T\ddot{o}}} \newcommand{ZZ}[0]{\mathbb{Z}} $$

$m\times n$行列のことを$(m,n)$型行列ということにする. $I_n$$n$次の単位行列とするとき, 以下が成り立つことが知られている.

$(m,n)$型行列$A$$(n,m)$型行列$B$に対して,
\begin{align} \det(I_m+AB)=\det(I_n+BA) \end{align}
が成り立つ.

これは, Sylvesterの行列式恒等式, Weinstein–Aronszajn identityと呼ばれているようであり, この等式自体はブロック行列$\left(\begin{matrix}I_m&A\\-B&I_n\end{matrix}\right)$を用いることによって簡潔に示すことができるようである( Wikipediaの記事 を参照). 今回は別のアプローチとして, $A,B$に関して対称的な表示を与えることによって定理1の証明を与えたいと思う.

Cauchy-Binetの公式

$[N]:=\{1,\dots,N\}$と表す. $A$$(m,n)$型行列とする. $[n]$の部分集合$S=\{s_1,\dots,s_l\}$に対し, $A$$s_1,\dots,s_l$列目以外を除いた行列を$A_S$として, $[m]$の部分集合$T=\{t_1,\dots,t_l\}$に対し, $t_1,\dots,t_l$行目以外を除いた行列を$A^T$と書くことにする(これは転置行列の記号と紛らわしいが, この記事では転置行列を用いない).

このとき, 以下が成り立つことが知られている.

Cauchy-Binetの公式

$(m,n)$型行列$A$$(n,m)$型行列$B$に対して,
\begin{align} \det(AB)&=\sum_{\substack{S\subset[n]\\|S|=m}}\det(A_S)\det(B^S) \end{align}
が成り立つ.

証明は探せばすぐに出てくるので, ここでは省略する.

$A_S^T:=(A_S)^T$と書くことにする. 1つ補題を用意する.

$m$次正方行列$C$に対し
\begin{align} \det(I_m+C)=\sum_{S\subset[m]}C_S^S \end{align}
が成り立つ.

$C=(c_{ij})_{1\leq i,j\leq m}$として, $\delta_{ij}$をKroneckerのデルタ, $\mathfrak{S}_m$$m$次対称群とすると,
\begin{align} \det(I_m+C)&=\sum_{\sigma\in\mathfrak{S}_m}\sgn(\sigma)(c_{1\sigma(1)}+\delta_{1\sigma(1)})\cdots(c_{m\sigma(m)}+\delta_{m\sigma(m)})\\ &=\sum_{\sigma\in\mathfrak{S}_m}\sgn(\sigma)\sum_{S\subset \{i\in[m];i=\sigma(i)\}}\prod_{j\in [m]\setminus S}c_{j\sigma(j)}\\ &=\sum_{S\subset [m]}\sum_{\substack{\sigma\in\mathfrak{S}_m\\\sigma(i)=i\,\mathrm{if}\,i\in S}}\sgn(\sigma)\prod_{j\in [m]\setminus S}c_{j\sigma(j)} \end{align}
となる. ここで, 内側の和は$S$の元を固定するような$\mathfrak{S}_m$の元を渡る. それは$[m]\setminus S$の置換群$\mathfrak{S}([m]\setminus S)$の元と見なすことができ,
\begin{align} &\sum_{S\subset [m]}\sum_{\substack{\sigma\in\mathfrak{S}_m\\\sigma(i)=i\,\mathrm{if}\,i\in S}}\sgn(\sigma)\prod_{j\in [m]\setminus S}c_{j\sigma(j)}\\ &=\sum_{S\subset [m]}\sum_{\substack{\sigma\in\mathfrak{S}([m]\setminus S)}}\sgn(\sigma)\prod_{j\in [m]\setminus S}c_{j\sigma(j)}\\ &=\sum_{S\subset [m]}\sum_{\substack{\sigma\in\mathfrak{S}(S)}}\sgn(\sigma)\prod_{j\in S}c_{j\sigma(j)} &S\mapsto [m]\setminus S\\ &=\sum_{S\subset [m]}\det(C_S^S) \end{align}
となるので, 示すべき等式を得る.

対称的な表示

これらを用いることによって, 以下を示すことができる.

$(m,n)$型行列$A$$(n,m)$型行列$B$に対して,
\begin{align} \det(I_m+AB)&=\sum_{\substack{S\subset [m]\\T\subset [n]\\|S|=|T|}}\det(A_T^S)\det(B^T_S) \end{align}
が成り立つ.

補題3より,
\begin{align} \det(I_m+AB)=\sum_{S\subset [m]}\det((AB)_S^S) \end{align}
である. ここで, 行列の積の定義から$(AB)_S^S=A^SB_S$となることから, Cauchy-Binetの公式より
\begin{align} \sum_{S\subset [m]}\det((AB)_S^S)&=\sum_{S\subset [m]}\det(A^SB_S)\\ &=\sum_{\substack{S\subset [m]\\T\subset [n]\\|S|=|T|}}\det(A^S_T)\det(B_S^T) \end{align}
となって示すべき等式を得る.

定理4は$A,B$に関して対称的な表示であることから, 特に
\begin{align} \det(I_m+AB)=\det(I_n+BA) \end{align}
を得ることができる.

参考文献

投稿日:6日前
更新日:6日前
数学の力で現場を変える アルゴリズムエンジニア募集 - Mathlog served by OptHub

この記事を高評価した人

高評価したユーザはいません

この記事に送られたバッジ

バッジはありません。

投稿者

Wataru
Wataru
1154
84396
超幾何関数, 直交関数, 多重ゼータ値などに興味があります

コメント

他の人のコメント

コメントはありません。
読み込み中...
読み込み中