0

Jacobiの補小行列公式

23
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{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}} $$

前の記事 の記法を用いる. $V$を内積を備えたベクトル空間とすると, $v,w\in V$に対し, 正則行列$A\in\mathrm{GL}(V)$
\begin{align} \langle v,Aw\rangle=\langle A^{\top}v,w\rangle \end{align}
を満たす. よって, $n$ベクトル$v=v_1\wedge\cdots\wedge v_n,w=w_1\wedge\cdots\wedge w_n$に対して
\begin{align} Av&=Av_1\wedge\cdots\wedge Av_n\\ Aw&=Aw_1\wedge\cdots\wedge Aw_n \end{align}
のように表すとすると, $n$ベクトルに対しても
\begin{align} \langle v,Aw\rangle=\langle A^{\top}v,w\rangle \end{align}
が成り立つことが分かる. これと 前の記事 の命題2を合わせると
\begin{align} v\wedge \star(Aw)=A^{\top}v\wedge \star w \end{align}
であることが分かる. ここで, 左辺は
\begin{align} A^T(v\wedge (A^{\top})^{-1}(\star w))&=\det(A)(v\wedge (A^{\top})^{-1}(\star w))\\ &=v\wedge \det(A)((A^{\top})^{-1}(\star w)) \end{align}
となるから, 以下が得られる.

$V$を内積を備えたベクトル空間とする. 正則行列$A\in\mathrm{GL}(V)$$n$ベクトル$w$に対し,
\begin{align} \star(Aw)=\det(A)((A^{\top})^{-1}(\star w)) \end{align}
が成り立つ.

Jacobiの補小行列公式

$m\times n$行列$A$$S\subset[m]$に対し, $A$$S$に含まれる添字の行を以外を除いた行列を$A^S$と表す. また, $A_S^T=(A_S)^T$のように表すとする.

Jacobiの補小行列公式

$r$次の正則行列$A$$S,T\subset[r],|S|=|T|=n$に対し,
\begin{align} \varepsilon(S,S^c)\varepsilon(T,T^c)\det(A^{T}_{S})=\det(A)\det((A^{-1})^{S^c}_{T^c}) \end{align}
が成り立つ.

これは余因子行列と逆行列の関係の一般化を与えている.

$e_1,\dots,e_r$$V$の標準基底として, 命題1において$w=e_T$とすると, 左辺は
\begin{align} \star(Aw)&=\star\left(\sum_{\substack{S\subset [r]\\|S|=n}}\det(A_S^T)e_S\right)\\ &=\sum_{\substack{S\subset [r]\\|S|=n}}\det(A_S^T)\varepsilon(S,S^c)e_{S^c} \end{align}
となる. 一方右辺は
\begin{align} &\det(A)((A^{\top})^{-1}(\star w))\\ &=\varepsilon(T,T^c)\det(A)((A^{\top})^{-1}e_{T^c})\\ &=\varepsilon(T,T^c)\det(A)\sum_{S\subset[m]}\det(((A^{\top})^{-1})^{T^c}_S)e_S\\ &=\varepsilon(T,T^c)\det(A)\sum_{S\subset[m]}\det((A^{-1})^S_{T^c})e_S\\ \end{align}
となる. よって, 両辺の$e_{S^c}$の係数を比較して示すべき等式を得る.

これを用いると以下のような公式を示すことができる.

$r$次正則行列$A,B$
\begin{align} AB^{\top}=C \end{align}
であり, ある$S,T\subset [r],|S|+|T|=r$に対し, $C_S^T=0$となるものとする. このとき, $R\subset [r],|R|=|T|$に対し,
\begin{align} \det(A_R^T)&=\frac{\varepsilon(S,S^c)\varepsilon(R^c,R)\det(C^T_{S^c})\det(B_{R^c}^S)}{\det(B)} \end{align}
が成り立つ.

\begin{align} A_R^T&=(C(B^\top)^{-1})_{R}^T\\ &=C^T((B^{\top})^{-1})_R \end{align}
あり, $C^T$$S$に含まれる添字の列は$0$であるから, Cauchy-Binetの公式と定理2を用いて,
\begin{align} \det(A_R^T)&=\det(C^T_{S^c})\det(((B^\top)^{-1})_R^{S^c})\\ &=\det(C^T_{S^c})\det(B)^{-1}\varepsilon(R^c,R)\varepsilon(S,S^c)\det((B^\top)_{S}^{R^c})\\ &=\frac{\varepsilon(S,S^c)\varepsilon(R^c,R)\det(C^T_{S^c})\det(B_{R^c}^S)}{\det(B)} \end{align}
と示される.

条件から
\begin{align} 0=C_S^T=(AB^{\top})_S^T=A^T(B^S)^{\top} \end{align}
であるから, $A^T$の行ベクトルと$B^S$の行ベクトルは直交している. よって, 前の記事 の定理3から
\begin{align} \det(A^T_R)=c\varepsilon(R^c,R)\det(B^S_{R^c}) \end{align}
を満たす$R$によらない$c$が存在することが分かるが, 上の公式は$C$が分かっている場合に$c$を明示的に表すことができることを意味している.

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

この記事を高評価した人

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

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

バッジはありません。

投稿者

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

コメント

他の人のコメント

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