コーシー・ビネの公式

提供: testwiki
ナビゲーションに移動 検索に移動

テンプレート:出典の明記 テンプレート:混同 代数学におけるコーシー・ビネの公式(こーしー・びねのこうしき、テンプレート:Lang-en-short)、あるいは、コーシー・ビネの定理コーシー・ビネの展開とは、ジャック・フィリップ・マリー・ビネおよびオーギュスタン=ルイ・コーシーに由来する恒等式で、2つの行列の積から作られる正方行列行列式は、元の行列から取り出せる最大の小行列式の積の和に等しいというものであり[1]、行列の成分が実数複素数などの可換環において成立する。

定理

テンプレート:Mvar自然数とし、集合 テンプレート:Math2テンプレート:Math と表記する。テンプレート:Mvar を非負整数として、テンプレート:Mvarテンプレート:Math2行列テンプレート:Mvarテンプレート:Math2行列とする。テンプレート:Mvar を要素数 テンプレート:Math2テンプレート:Math の部分集合とする。テンプレート:Mvar を、テンプレート:Mvarテンプレート:Mvar個のから テンプレート:Mvar に含まれる添字のを取り出して得られた テンプレート:Mvar次正方行列、テンプレート:Mvar を、テンプレート:Mvarテンプレート:Mvar個のから テンプレート:Mvar に含まれる添字のを取り出して得られた テンプレート:Mvar次正方行列とする。

このとき積 テンプレート:Mvarテンプレート:Math2行列となり、その行列式

det(AB)=S[n]|S|=mdet(AS)det(BS)

となる。ただし、総和について、テンプレート:Mvar は、テンプレート:Math の部分集合で要素数が テンプレート:Mvar のものすべてを取るとする。なお、テンプレート:Math2 の場合は右辺は 0 である。

成分による表示

A:=(a1,1a1,ma1,nam,1am,mam,n),B:=(b1,1b1,mbm,1bm,mbn,1bn,m)

と成分表示すると、公式は

|k=1na1,kbk,1k=1na1,kbk,mk=1nam,kbk,1k=1nam,kbk,m|=1k1<<kmn|a1,k1a1,kmam,k1am,km||bk1,1bk1,mbkm,1bkm,m|

と表現できる。ただし、右辺の総和は、テンプレート:Math2 を満たす整数の組 テンプレート:Math2 の全てに対して取るとする。なお、テンプレート:Math2 の場合は右辺は テンプレート:Math である。

小行列式記法による表示

記法

A(i1ilj1jl)|ai1,j1ai1,jlail,j1ail,jl|

を使えば

(AB)(1m1m)=1k1<<kmnA(1mk1km)B(k1km1m)

となる。ただし、右辺の総和は、テンプレート:Math2 を満たす整数の組 テンプレート:Math2 の全てに対して和を取るとする。なお、テンプレート:Math2 の場合は右辺は テンプレート:Math である。

定理の証明

行列式の多重線型性により

det(AB)=|p1=1na1,p1bp1,1pm=1na1,pmbpm,mp1=1nam,p1bp1,1pm=1nam,pmbpm,m|=p1=1npm=1n|a1,p1a1,pmam,p1am,pm|bp1,1bpm,m=p:[m][n]A(1mp(1)p(m))bp(1),1bp(m),m

が導かれる。最後の式の テンプレート:Mvarテンプレート:Math2 から テンプレート:Math2 への写像である。

行列式の反対称性により テンプレート:Mvar単射の場合のみ行列式は非零なので、テンプレート:Math2 と置き換えられる。ここで、置換 テンプレート:Math2テンプレート:Mvar対称群 𝔖m の元であり、テンプレート:Math2テンプレート:Math2 を満たす関数である。これより

det(AB)=π𝔖mk:[m][n]i<jk(i)<k(j)A(1mk(π(1))k(π(m)))bk(π(1)),1bk(π(m)),m=π𝔖mk:[m][n]i<jk(i)<k(j)sgn(π)A(1mk(1)k(m))bk(π(1)),1bk(π(m)),m=k:[m][n]i<jk(i)<k(j)A(1mk(1)k(m))π𝔖msgn(π)bk(π(1)),1bk(π(m)),m=k:[m][n]i<jk(i)<k(j)A(1mk(1)k(m))B(k(1)k(m)1m)=1k1<<kmnA(1mk1km)B(k1km1m)

が成り立つ。なお、テンプレート:Math は置換 テンプレート:Π符号であり、行列式の反対称性

A(1mk(π(1))k(π(m)))=sgn(π)A(1mk(1)k(m))

および、

B(k(1)k(m)1m)=π𝔖msgn(π)bk(π(1)),1bk(π(m)),m

を用いた。

具体例

コーシー・ビネの公式を具体例で確認してみる。

(例1)テンプレート:Math2 の場合として、行列 A=(112),B=(130) を考える。AB=(4),det(AB)=4である。

S={1},{2},{3} であるから、

S[n]|S|=mdet(AS)det(BS)=det(1)det(1)+det(1)det(3)+det(2)det(0)=4.

となる。

(例2)テンプレート:Math2 の場合として、行列 A=(112311),B=(113102) を考える。AB=(4662),det(AB)=|4662|=28 である。

S={1,2},{2,3},{1,3} であるから、

S[n]|S|=mdet(AS)det(BS)=|1131||1131|+|1211||3102|+|1231||1102|=(2)×(2)+(3)×6+(7)×2=28.

となる。

(例3)テンプレート:Math2 の場合として、行列 A=(112311110),B=(110310021) を考える。AB=(462621200),det(AB)=|462621200|=20 である。

S={1,2,3} であるから、

S[n]|S|=mdet(AS)det(BS)=|112311110||110310021|=(10)×(2)=20.

となる。

(例4)テンプレート:Math2 の場合として、行列 A=(112311110100),B=(110131000210) を考える。AB=(4621621320011101),det(AB)=|4621621320011101|=0 である。

S は存在しないから、

S[n]|S|=mdet(AS)det(BS)=0.

となる。

一般化されたクロネッカーのデルタとの関係

A(δ1i1δmi1δni1δ1imδmimδnim),B(δj11δjm1δj1mδjmmδj1nδjmn)(ik[n],jk[n]fork[m])

とする。ただし、δはクロネッカーのデルタ

δνμ={1(μ=ν)0(μν)

である。 これをコーシー・ビネの公式に代入し、一般化されたクロネッカーのデルタ

δj1jmi1im|δj1i1δjmi1δj1imδjmim|

を使えば、

δj1jmi1im=1k1<<kmnδk1kmi1imδj1jmk1km(ip[n],jp[n]forp[m])

が得られる。逆にこの式からコーシー・ビネの公式を導くこともできる。

これは単位行列の基本的性質

δji=k=1nδkiδjk(i[n],j[n])

の一般化である。

特別な場合

テンプレート:Math2 の場合、テンプレート:Math2 となる整数の組 テンプレート:Math2 は存在しないから、公式の右辺は 0 となる。このとき テンプレート:Math2階数はこの場合高々 テンプレート:Mvar だから、テンプレート:Math2 行列 テンプレート:Mvar の階数も高々 テンプレート:Math2 であるので、公式の左辺 テンプレート:Math は 0 となり、公式が成り立つ。

テンプレート:Math2 のとき、テンプレート:Math2 は正方行列である。テンプレート:Math2 となる整数の組 テンプレート:Math2テンプレート:Math に等しいから、公式は

(AB)(1m1m)=A(1m1n)B(1n1m)

すなわち テンプレート:Math2 となる。

テンプレート:Math2 のとき、テンプレート:Math2 そして テンプレート:Mvar は空行列(ただし、テンプレート:Math2 なら行列の型は異なる)であり、空行列の行列式は定義により 1 だから、公式は 1 = 1 を述べているに過ぎない。

テンプレート:Math2 のとき、公式は detj=1nA1,jBj,1=k=1ndet(Ak1)det(B1k) となるが、1×1行列 A に対して det(A) = A だから、自明の式を述べているに過ぎない。

テンプレート:Math2 は非自明な公式を与える最小の テンプレート:Mvar であり、そのときの公式

|k=1na1,kbk,1k=1na1,kbk,2k=1na2,kbk,1k=1na2,kbk,2|=1k1<k2n|a1,k1a1,k2a2,k1a2,k2||bk1,1bk1,2bk2,1bk2,2|

ビネ・コーシーの恒等式と呼ばれる。

n = 3 の場合の具体例

𝒂,𝒙 などは3次元ベクトルとする。 1=1(m=0)𝒂𝒙=a1x1+a2x2+a3x3(m=1)|𝒂𝒙𝒂𝒚𝒃𝒙𝒃𝒚|=|a2a3b2b3||x2y2x3y3|+|a3a1b3b1||x3y3x1y1|+|a1a2b1b2||x1y1x2y2|=(𝒂×𝒃)(𝒙×𝒚)(m=2)|𝒂𝒙𝒂𝒚𝒂𝒛𝒃𝒙𝒃𝒚𝒃𝒛𝒄𝒙𝒄𝒚𝒄𝒛|=|a1a2a3b1b2b3c1c2c3||x1y1z1x2y2z2x3y3z3|={𝒂(𝒃×𝒄)}{𝒙(𝒚×𝒛)}(m=3)|𝒂𝒙𝒂𝒚𝒂𝒛𝒂𝒘𝒃𝒙𝒃𝒚𝒃𝒛𝒃𝒘𝒄𝒙𝒄𝒚𝒄𝒛𝒄𝒘𝒅𝒙𝒅𝒚𝒅𝒛𝒅𝒘|=0(m=4)

m > 3 の場合、右辺は常に0である。なお、

m = 2 の式はスカラー四重積に対するビネ・コーシーの恒等式
m = 3 の式はスカラー三重積の積に対する公式

であり、m = 4 の式より四重積 (ベクトル解析) の公式:

[𝒂,𝒃,𝒄]𝒅=[𝒅,𝒃,𝒄]𝒂+[𝒂,𝒅,𝒄]𝒃+[𝒂,𝒃,𝒅]𝒄([𝒂,𝒃,𝒄]=𝒂(𝒃×𝒄))

が導かれる。

脚注

テンプレート:脚注ヘルプ テンプレート:Reflist

参考文献

  • テンプレート:Cite book
  • Joel G. Broida & S. Gill Williamson (1989) A Comprehensive Introduction to Linear Algebra, §4.6 Cauchy- Binet theorem, pp.208-14, Addison-Wesley ISBN 0-201-50065-5.
  • Jin Ho Kwak & Sungpyo Hong (2004) Linear Algebra 2nd edition, Example 2.15 Binet-Cauchy formula, pp.66,7, Birkhäuser ISBN 0-8176-4294-3.
  • Igor R. Shafarevich & Alexey O. Remizov (2012) Linear Algebra and Geometry, §2.9 (p.68) & §10.5 (p. 377), Springer ISBN 978-3-642-30993-9.
  • Aaron Lauve (2004) A short combinatoric proof of Cauchy–Binet formula from Université du Québec à Montréal.

関連項目

外部リンク

テンプレート:線形代数