กลับไปหน้าบทความ

อ่าน 9 นาที

สูตรของจาโคบี

ปัจจัยกำหนด/ทฤษฎีเมทริกซ์

ในแคลคูลัสเมทริกซ์สูตรของ Jacobiแสดงอนุพันธ์ของ ดีเทอร์มิ แนนต์ของเมทริกซ์AในรูปของแอดจูเกตของAและอนุพันธ์ของA

สูตรของจาโคบี

ในแคลคูลัสเมทริกซ์สูตรของ Jacobiแสดงอนุพันธ์ของ ดีเทอร์มิ แนนต์ของเมทริกซ์AในรูปของแอดจูเกตของAและอนุพันธ์ของA [ 1 ]

ถ้าAเป็นฟังก์ชันที่หาอนุพันธ์ได้จากจำนวนจริงไปยัง เมทริกซ์ n  × n แล้ว

ทีเดทเอ(ที)=tr(adj(เอ(ที))เอ(ที)ที)=(เดทเอ(ที))tr(เอ(ที)1เอ(ที)ที){\displaystyle {\frac {d}{dt}}\det A(t)=\operatorname {tr} \left(\operatorname {adj} (A(t))\,{\frac {dA(t)}{dt}}\right)=\left(\det A(t)\right)\cdot \operatorname {tr} \left(A(t)^{-1}\cdot \,{\frac {dA(t)}{dt}}\right)}

โดยที่tr( X )คือร่องรอยของเมทริกซ์Xและadj(X){\displaystyle \operatorname {adj} (X)}คือเมทริกซ์ผกผัน ของมัน (ความเท่าเทียมกันข้อหลังนี้เป็นจริงก็ต่อเมื่อA ( t ) สามารถผกผันได้ )

ในกรณีพิเศษ

เดท(เอ)เอฉันเจ=adj(เอ)เจฉัน=adj(เอ)ฉันเจทีเดท(เอ)เอ=เดท(เอ)เอที ถ้า A เป็นอินเทอร์เฟส{\displaystyle {\partial \det(A) \over \partial A_{ij}}=\operatorname {adj} (A)_{ji}=\operatorname {adj} (A)_{ij}^{T}\Longrightarrow {\partial \det(A) \over \partial A}=\operatorname {det} (A)A^{-T}{\text{ ถ้า A เป็นอินเวอร์สได้}}}

ในทำนองเดียวกัน ถ้าdAแทนอนุพันธ์ของAสูตรทั่วไปคือ

เดท(เอ)=tr(adj(เอ)เอ)=เดท(เอ)tr(เอ1เอ){\displaystyle d\det(A)=\operatorname {tr} (\operatorname {adj} (A)\,dA)=\det(A)\operatorname {tr} \left(A^{-1}dA\right)}

สูตรนี้ตั้งชื่อตามนักคณิตศาสตร์คาร์ล จาโคบี

อนุพันธ์

โดยการคำนวณเมทริกซ์

ทฤษฎีบท (สูตรของจาโคบี) สำหรับฟังก์ชันอนุพันธ์ใดๆAจากจำนวนจริงไปยังเมทริกซ์n  × n 

เดท(เอ)=tr(adj(เอ)เอ).{\displaystyle d\det(A)=\operatorname {tr} (\operatorname {adj} (A)\,dA).}

บทพิสูจน์สูตรของลาปลาซสำหรับดีเทอร์มิแนนต์ของเมทริกซ์Aสามารถเขียนได้ดังนี้

เดท(เอ)=เจเอฉันเจadjที(เอ)ฉันเจ.{\displaystyle \det(A)=\sum _{j}A_{ij}\operatorname {adj} ^{\rm {T}}(A)_{ij}.}

โปรดสังเกตว่าการหาผลรวมนั้นกระทำกับแถวi ใดๆ ของเมทริกซ์

ดีเทอร์มิแนนต์ของAสามารถพิจารณาได้ว่าเป็นฟังก์ชันขององค์ประกอบของA :

เดท(เอ)=เอฟ(เอ11,เอ12,,เอ21,เอ22,,เอnn){\displaystyle \det(A)=F\,(A_{11},A_{12},\ldots ,A_{21},A_{22},\ldots ,A_{nn})}

ดังนั้น โดยใช้กฎลูกโซ่อนุพันธ์ของมันคือ

เดท(เอ)=ฉันเจเอฟเอฉันเจเอฉันเจ.{\displaystyle d\det(A)=\sum _{i}\sum _{j}{\partial F \over \partial A_{ij}}\,dA_{ij}.}

การหาผลรวมนี้กระทำกับองค์ประกอบทั้งหมดขนาด n × nของเมทริกซ์

ในการหา ∂ F /∂ A ให้พิจารณาว่าทางด้านขวามือของสูตรลาปลาส ดัชนีiสามารถเลือกได้ตามต้องการ (เพื่อเพิ่มประสิทธิภาพในการคำนวณ: การเลือกแบบอื่นจะให้ผลลัพธ์เดียวกันในที่สุด แต่การคำนวณอาจยากขึ้นมาก) โดยเฉพาะอย่างยิ่ง สามารถเลือกให้ตรงกับดัชนีแรกของ ∂  /  A ได้:

เดท(เอ)เอฉันเจ=เคเอฉันเคadjที(เอ)ฉันเคเอฉันเจ=เค(เอฉันเคadjที(เอ)ฉันเค)เอฉันเจ{\displaystyle {\partial \det(A) \over \partial A_{ij}}={\partial \sum _{k}A_{ik}\operatorname {adj} ^{\rm {T}}(A)_{ik} \over \partial A_{ij}}=\sum _{k}{\partial (A_{ik}\operatorname {adj} ^{\rm {T}}(A)_{ik}) \over \partial A_{ij}}}

ดังนั้น ตามกฎการคูณ

เดท(เอ)เอฉันเจ=เคเอฉันเคเอฉันเจadjที(เอ)ฉันเค+เคเอฉันเคadjที(เอ)ฉันเคเอฉันเจ.{\displaystyle {\partial \det(A) \over \partial A_{ij}}=\sum _{k}{\partial A_{ik} \over \partial A_{ij}}\operatorname {adj} ^{\rm {T}}(A)_{ik}+\sum _{k}A_{ik}{\partial \operatorname {adj} ^{\rm {T}}(A)_{ik} \over \partial A_{ij}}.}

ทีนี้ ถ้าหากองค์ประกอบของเมทริกซ์A และโคแฟกเตอร์ adj T ( A ) ขององค์ประกอบA อยู่ในแถวเดียวกัน (หรือคอลัมน์เดียวกัน) โคแฟกเตอร์นั้นจะไม่ใช่ฟังก์ชันของA เพราะโคแฟกเตอร์ของA ถูกแสดงในรูปขององค์ประกอบที่ไม่ได้อยู่ในแถวเดียวกัน (หรือคอลัมน์เดียวกัน) ดังนั้น

adjที(เอ)ฉันเคเอฉันเจ=0,{\displaystyle {\partial \operatorname {adj} ^{\rm {T}}(A)_{ik} \over \partial A_{ij}}=0,}

ดังนั้น

เดท(เอ)เอฉันเจ=เคadjที(เอ)ฉันเคเอฉันเคเอฉันเจ.{\displaystyle {\partial \det(A) \over \partial A_{ij}}=\sum _{k}\operatorname {adj} ^{\rm {T}}(A)_{ik}{\partial A_{ik} \over \partial A_{ij}}.}

องค์ประกอบทั้งหมดของAเป็นอิสระต่อกัน กล่าวคือ

เอฉันเคเอฉันเจ=δเจเค,{\displaystyle {\partial A_{ik} \over \partial A_{ij}}=\delta _{jk},}

โดยที่δคือเดลต้าโครเนกเกอร์ดังนั้น

เดท(เอ)เอฉันเจ=เคadjที(เอ)ฉันเคδเจเค=adjที(เอ)ฉันเจ.{\displaystyle {\partial \det(A) \over \partial A_{ij}}=\sum _{k}\operatorname {adj} ^{\rm {T}}(A)_{ik}\delta _{jk}=\operatorname {adj} ^{\rm {T}}(A)_{ij}.}

ดังนั้น,

(เดท(เอ))=ฉันเจadjที(เอ)ฉันเจเอฉันเจ=เจฉันadj(เอ)เจฉันเอฉันเจ=เจ(adj(เอ)เอ)เจเจ=tr(adj(เอ)เอ). {\displaystyle d(\det(A))=\sum _{i}\sum _{j}\operatorname {adj} ^{\rm {T}}(A)_{ij}\,dA_{ij}=\sum _{j}\sum _{i}\operatorname {adj} (A)_{ji}\,dA_{ij}=\sum _{j}(\operatorname {adj} (A)\,dA)_{jj}=\operatorname {tr} (\operatorname {adj} (A)\,dA).\ \square }

ผ่านกฎลูกโซ่

ทฤษฎีบทเสริมข้อที่ 1.เดท(ฉัน)=ที{\displaystyle \det '(I)=\mathrm {tr} }, ที่ไหนเดท{\displaystyle \det '}คือผลต่างของเดท{\displaystyle \det }.

สมการนี้หมายความว่าอนุพันธ์ของเดท{\displaystyle \det }เมื่อประเมินที่เมทริกซ์เอกลักษณ์จะมีค่าเท่ากับร่องรอย (trace) อนุพันธ์เดท(ฉัน){\displaystyle \det '(I)}เป็นตัวดำเนินการเชิงเส้นที่แปลง เมทริกซ์ขนาด n × nไปเป็นจำนวนจริง

บทพิสูจน์โดยใช้คำนิยามของอนุพันธ์ทิศทางร่วมกับคุณสมบัติพื้นฐานข้อหนึ่งของอนุพันธ์ทิศทางสำหรับฟังก์ชันที่หาอนุพันธ์ได้ เราจะได้ว่า

เดท(ฉัน)(ที)=ทีเดท(ฉัน)=ลิมε0เดท(ฉัน+εที)เดทฉันε{\displaystyle \det '(I)(T)=\nabla _{T}\det(I)=\lim _{\varepsilon \to 0}{\frac {\det(I+\varepsilon T)-\det I}{\varepsilon }}}

เดท(ฉัน+εที){\displaystyle \det(I+\varepsilon T)}เป็นพหุนามในε{\displaystyle \varepsilon }มีอันดับnมีความสัมพันธ์อย่างใกล้ชิดกับพหุนามลักษณะเฉพาะของที{\displaystyle T}พจน์คงที่ในพหุนามนั้น (พจน์ที่มีε=0{\displaystyle \varepsilon =0}) คือ 1 ในขณะที่พจน์เชิงเส้นในε{\displaystyle \varepsilon }เป็นที ที{\displaystyle \mathrm {tr} \ T}ดังนั้นลิมิตจึงเท่ากับที ที{\displaystyle \mathrm {tr} \ T}ซึ่งเป็นข้อกล่าวอ้างนั้น

ทฤษฎีบทเสริมข้อที่ 2สำหรับเมทริกซ์ผกผันAเราจะได้ว่า:เดท(เอ)(ที)=เดทเอที(เอ1ที){\displaystyle \det '(A)(T)=\det A\;\mathrm {tr} (A^{-1}T)}.

บทพิสูจน์พิจารณาฟังก์ชันต่อไปนี้ของX :

เดทX=เดท(เอเอ1X)=เดท(เอ) เดท(เอ1X){\displaystyle \det X=\det(AA^{-1}X)=\det(A)\ \det(A^{-1}X)}

เราคำนวณอนุพันธ์ของเดทX{\displaystyle \det X}และประเมินผลที่X=เอ{\displaystyle X=A}โดยใช้ Lemma 1 สมการข้างต้น และกฎลูกโซ่:

เดท(เอ)(ที)=เดทเอ เดท(ฉัน)(เอ1ที)=เดทเอ ที(เอ1ที){\displaystyle \det '(A)(T)=\det A\ \det '(I)(A^{-1}T)=\det A\ \mathrm {tr} (A^{-1}T)}

ทฤษฎีบท (สูตรของจาโคบี) ทีเดทเอ=ที(เอเจ เอเอที){\displaystyle {\frac {d}{dt}}\det A=\mathrm {tr} \left(\mathrm {adj} \ A{\frac {dA}{dt}}\right)}

หลักฐาน.ถ้าเอ{\displaystyle A}สามารถผกผันได้ โดยใช้ Lemma 2 โดยมีที=เอ/ที{\displaystyle T=dA/dt}

ทีเดทเอ=เดทเอที(เอ1เอที)=ที(เอเจ เอเอที){\displaystyle {\frac {d}{dt}}\det A=\det A\;\mathrm {tr} \left(A^{-1}{\frac {dA}{dt}}\right)=\mathrm {tr} \left(\mathrm {adj} \ A\;{\frac {dA}{dt}}\right)}

โดยใช้สมการที่เชื่อมโยงค่าผกผันของเอ{\displaystyle A}ถึงเอ1{\displaystyle A^{-1}}สูตรนี้ใช้ได้กับเมทริกซ์ทุกตัว เนื่องจากเซตของเมทริกซ์เชิงเส้นที่ผกผันได้นั้นมีความหนาแน่นในปริภูมิของเมทริกซ์

โดยการหาค่าแนวทแยง

ทั้งสองข้างของสูตร Jacobi เป็นพหุนามในสัมประสิทธิ์เมทริกซ์ของAและA'ดังนั้นจึงเพียงพอที่จะตรวจสอบเอกลักษณ์ของพหุนามบนเซตย่อยหนาแน่นซึ่งค่าลักษณะเฉพาะของAแตกต่างกันและไม่เป็นศูนย์

ถ้าAแยกตัวประกอบได้ต่างชนิดกันดังนี้เอ=บีซี{\displaystyle A=BC}, แล้ว

ที(เอ1เอ)=ที((บีซี)1(บีซี))=ที(บี1บี)+ที(ซี1ซี).{\displaystyle \mathrm {tr} (A^{-1}A')=\mathrm {tr} ((BC)^{-1}(BC)')=\mathrm {tr} (B^{-1}B')+\mathrm {tr} (C^{-1}C').}

โดยเฉพาะอย่างยิ่ง ถ้าLเป็นอินเทอร์เฟสแล้วฉัน=แอล1แอล{\displaystyle I=L^{-1}L}และ

0=ที(ฉัน1ฉัน)=ที(แอล(แอล1))+ที(แอล1แอล).{\displaystyle 0=\mathrm {tr} (I^{-1}I')=\mathrm {tr} (L(L^{-1})')+\mathrm {tr} (L^{-1}L').}

เนื่องจาก เมทริก ซ์ Aมีค่าลักษณะเฉพาะที่แตกต่างกัน จึงมีเมทริกซ์เชิงซ้อนผกผันได้L ที่สามารถหาอนุพันธ์ได้ ซึ่งทำให้ เอ=แอล1ดีแอล{\displaystyle A=L^{-1}DL}และDเป็นแนวทแยงมุม จากนั้น

ที(เอ1เอ)=ที(แอล(แอล1))+ที(ดี1ดี)+ที(แอล1แอล)=ที(ดี1ดี).{\displaystyle \mathrm {tr} (A^{-1}A')=\mathrm {tr} (L(L^{-1})')+\mathrm {tr} (D^{-1}D')+\mathrm {tr} (L^{-1}L')=\mathrm {tr} (D^{-1}D').}

อนุญาตλฉัน{\displaystyle \lambda _{i}},ฉัน=1,,n{\displaystyle i=1,\ldots ,n} ให้ค่าไอเกนของเมทริกซ์ A เป็นค่าไอเกน จากนั้น

(lnเดทเอ)=(ฉัน=1nlnλฉัน)=ฉัน=1nλฉัน/λฉัน=ที(ดี1ดี)=ที(เอ1เอ),{\displaystyle \left(\ln \det A\right)'=\left(\sum _{i=1}^{n}\ln \lambda _{i}\right)'=\sum _{i=1}^{n}\lambda _{i}'/\lambda _{i}=\mathrm {tr} (D^{-1}D')=\mathrm {tr} (A^{-1}A'),}

ซึ่งเป็นสูตรของจาโคบีสำหรับเมทริกซ์Aที่มีค่าไอเกนที่ไม่เป็นศูนย์และแตกต่างกัน

บทสรุป

ต่อไปนี้เป็นความสัมพันธ์ที่มีประโยชน์ซึ่งเชื่อมโยงร่องรอย (trace)กับดีเทอร์มิแนนต์ของเมทริกซ์เอกซ์โพเนนเชียล ที่เกี่ยวข้อง :

เดทอีบี=อีtr(บี){\displaystyle \det e^{B}=e^{\operatorname {tr} \left(B\right)}}

ข้อความนี้ชัดเจนสำหรับเมทริกซ์แนวทแยง และการพิสูจน์ข้ออ้างทั่วไปมีดังต่อไปนี้

สำหรับเมทริกซ์ผกผัน ใดๆเอ(ที){\displaystyle A(t)}ในหัวข้อก่อนหน้านี้"โดยใช้กฎลูกโซ่"เราได้แสดงให้เห็นแล้วว่า

ทีเดทเอ(ที)=เดทเอ(ที)tr(เอ(ที)1ทีเอ(ที)){\displaystyle {\frac {d}{dt}}\det A(t)=\det A(t)\;\operatorname {tr} \left(A(t)^{-1}\,{\frac {d}{dt}}A(t)\right)}

เมื่อพิจารณาแล้วเอ(ที)=เอ็กซ์(ทีบี){\displaystyle A(t)=\exp(tB)}จากสมการนี้จะได้ผลลัพธ์ดังนี้:

ทีเดทอีทีบี=tr(บี)เดทอีทีบี{\displaystyle {\frac {d}{dt}}\det e^{tB}=\operatorname {tr} (B)\det e^{tB}}

ผลลัพธ์ที่ต้องการได้มาจากการแก้สมการเชิงอนุพันธ์สามัญ นี้

แอปพลิเคชัน

สูตรหลายรูปแบบเป็นพื้นฐานของอัลกอริทึม Faddeev–LeVerrierสำหรับการคำนวณพหุนามลักษณะเฉพาะและการประยุกต์ใช้ทฤษฎีบท Cayley–Hamilton อย่างชัดเจน ตัวอย่างเช่น เริ่มจากสมการต่อไปนี้ ซึ่งได้รับการพิสูจน์แล้วข้างต้น:

ทีเดทเอ(ที)=เดทเอ(ที) tr(เอ(ที)1ทีเอ(ที)){\displaystyle {\frac {d}{dt}}\det A(t)=\det A(t)\ \operatorname {tr} \left(A(t)^{-1}\,{\frac {d}{dt}}A(t)\right)}

และการใช้เอ(ที)=ทีฉันบี{\displaystyle A(t)=tI-B}เราจึงได้ผลลัพธ์ดังนี้:

ทีเดท(ทีฉันบี)=เดท(ทีฉันบี)tr[(ทีฉันบี)1]=tr[adj(ทีฉันบี)]{\displaystyle {\frac {d}{dt}}\det(tI-B)=\det(tI-B)\operatorname {tr} [(tI-B)^{-1}]=\operatorname {tr} [\operatorname {adj} (tI-B)]}

โดยที่ adj หมายถึงเมทริกซ์ผกผัน

หมายเหตุ

  1. Magnus & Neudecker (1999 , หน้า149–150) , ส่วนที่สาม, หัวข้อ 8.3 
ดึงข้อมูลมาจาก " https://en.wikipedia.org/w/index.php?title=Jacobi%27s_formula&oldid=1358847198 "

สรุปเนื้อหา

ข้อมูลสำคัญจากบทความ

ข้อมูลสำคัญเกี่ยวกับ สูตรของจาโคบี

ในแคลคูลัสเมทริกซ์สูตรของ Jacobiแสดงอนุพันธ์ของ ดีเทอร์มิ แนนต์ของเมทริกซ์AในรูปของแอดจูเกตของAและอนุพันธ์ของA

โดยการคำนวณเมทริกซ์

ทฤษฎีบท (สูตรของจาโคบี) สำหรับฟังก์ชันอนุพันธ์ใดๆ A จากจำนวนจริงไปยังเมทริกซ์ n × n

ผ่านกฎลูกโซ่

ทฤษฎีบทเสริมข้อที่ 1. เดท ′ ( ฉัน ) = ที ร {\displaystyle \det '(I)=\mathrm {tr} } , ที่ไหน เดท ′ {\displaystyle \det '} คือผลต่างของ เดท {\displaystyle \det } .

โดยการหาค่าแนวทแยง

ทั้งสองข้างของสูตร Jacobi เป็นพหุนามในสัมประสิทธิ์เมทริกซ์ของ A และ A' ดังนั้นจึงเพียงพอที่จะตรวจสอบเอกลักษณ์ของพหุนามบนเซตย่อยหนาแน่นซึ่งค่าลักษณะเฉพาะของ A แตกต่างกันและไม่เป็นศูนย์