ในพีชคณิตเชิงเส้น เวกเตอร์ลักษณะทั่วไป ของเมทริกซ์ คือเวกเตอร์ ที่ตรงตามเกณฑ์บางประการซึ่งผ่อนปรนกว่าเกณฑ์สำหรับเวกเตอร์ลักษณะ ( ทั่วไป) [ 1 ] n × n {\displaystyle n\times n} เอ {\displaystyle A}
ให้เป็นปริภูมิเวกเตอร์ มิติและให้เป็นเมทริกซ์แทน ของแผนที่เชิงเส้นจากไปโดยสัมพันธ์กับฐาน เรียง ลำดับบางฐานวี {\displaystyle V} n {\displaystyle n} เอ {\displaystyle A} วี {\displaystyle V} วี {\displaystyle V}
อาจไม่มีชุด เวกเตอร์ลักษณะ เฉพาะที่เป็นอิสระเชิงเส้นที่ สมบูรณ์ ซึ่งประกอบเป็นฐานที่สมบูรณ์สำหรับ เสมอไป กล่าวคือ เมทริกซ์อาจไม่สามารถทำให้เป็นแนวทแยงได้ [ 2 ] [ 3 ] สิ่ง นี้เกิดขึ้นเมื่อความซ้ำเชิงพีชคณิต ของ ค่าลักษณะเฉพาะ อย่างน้อยหนึ่งค่ามากกว่าความซ้ำเชิงเรขาคณิต ( มิติ ของปริภูมิว่างของเมทริกซ์หรือมิติ ของปริภูมิว่าง ) ในกรณีนี้เรียกว่าค่าลักษณะเฉพาะที่บกพร่อง และเรียกว่า เมทริก ซ์ที่บกพร่อง [ 4 ] n {\displaystyle n} เอ {\displaystyle A} วี {\displaystyle V} เอ {\displaystyle A} λ ฉัน {\displaystyle \lambda _{i}} ( เอ − λ ฉัน ฉัน ) {\displaystyle (A-\lambda _{i}I)} λ ฉัน {\displaystyle \lambda _{i}} เอ {\displaystyle A}
เวกเตอร์ลักษณะทั่วไปที่ สอดคล้องกับพร้อมกับเมทริกซ์จะสร้างโซ่จอร์แดน ของเวกเตอร์ลักษณะทั่วไปที่เป็นอิสระเชิงเส้น ซึ่งเป็นฐานสำหรับปริภูมิย่อยไม่แปรเปลี่ยน ของ[ 5 ] [ 6 ] [ 7 ] x ฉัน {\displaystyle x_{i}} λ ฉัน {\displaystyle \lambda _{i}} ( เอ − λ ฉัน ฉัน ) {\displaystyle (A-\lambda _{i}I)} วี {\displaystyle V}
การใช้เวกเตอร์ลักษณะเฉพาะทั่วไป ชุดของเวกเตอร์ลักษณะเฉพาะที่เป็นอิสระเชิงเส้นของสามารถขยายได้ หากจำเป็น ไปสู่ฐานที่สมบูรณ์สำหรับ[ 8 ] ฐาน นี้สามารถใช้เพื่อกำหนด "เมทริกซ์เกือบทแยงมุม" ในรูปแบบปกติของจอร์แดน คล้ายกับ ซึ่งมีประโยชน์ในการคำนวณฟังก์ชันเมทริกซ์บางอย่างของ [ 9 ] เมท ริก ซ์นี้ยังมีประโยชน์ในการแก้ระบบสมการเชิงอนุพันธ์เชิงเส้น โดยที่ไม่จำเป็นต้องทำให้เป็นเมทริกซ์ทแยงมุม[ 10 ] [ 11 ] เอ {\displaystyle A} วี {\displaystyle V} เจ {\displaystyle J} เอ {\displaystyle A} เอ {\displaystyle A} เจ {\displaystyle J} x ′ = เอ x , {\displaystyle \mathbf {x} '=A\mathbf {x} ,} เอ {\displaystyle A}
มิติของปริภูมิไอเกนทั่วไปที่สอดคล้องกับค่าไอเกนที่กำหนดคือความหลากหลายทางพีชคณิตของ[ 12 ] λ {\displaystyle \lambda } λ {\displaystyle \lambda }
ภาพรวมและคำจำกัดความ มีหลายวิธีเทียบเท่ากันในการกำหนด เวก เตอร์ลักษณะเฉพาะทั่วไป [ 13 ] [ 14 ] [ 15 ] [ 16 ] [ 17 ] [ 18 ] [ 19 ] [ 20 ] สำหรับ วัตถุประสงค์ ของ เรา เวกเตอร์ลักษณะเฉพาะที่เกี่ยวข้องกับค่าลักษณะเฉพาะของเมทริกซ์× คือเวกเตอร์ที่ไม่เป็น ศูนย์ ซึ่งโดยที่คือเมทริกซ์เอกลักษณ์ × และคือเวกเตอร์ศูนย์ ที่มีความยาว[ 21 ] นั่น คืออยู่ในเคอร์เนล ของการแปลง ถ้ามีเวกเตอร์ลักษณะเฉพาะที่เป็นอิสระเชิงเส้น แล้วจะคล้ายกับเมทริกซ์แนวทแยง นั่นคือ มีเมทริกซ์ผกผันได้ ซึ่ง สามารถทำให้เป็นเมทริกซ์แนวทแยงได้ผ่านการแปลงความคล้ายคลึง [ 22 ] [ 23 ] เมทริกซ์เรียกว่าเมท ริก ซ์ สเปกตรัม สำหรับ เมท ริกซ์เรียกว่าเมทริกซ์ โมดอ ล สำหรับ[ 24 ] เมทริกซ์ที่สามารถทำให้เป็น เมทริกซ์ แนวทแยง ได้ มีความน่าสนใจเป็นพิเศษ เนื่องจากสามารถคำนวณฟังก์ชันเมทริกซ์ของเมทริกซ์เหล่านั้นได้ง่าย[ 25 ] คุณ {\displaystyle \mathbf {u} } λ {\displaystyle \lambda } n {\displaystyle n} n {\displaystyle n} เอ {\displaystyle A} ( เอ − λ ฉัน ) คุณ = 0 {\displaystyle (A-\lambda I)\mathbf {u} =\mathbf {0} } ฉัน {\displaystyle I} n {\displaystyle n} n {\displaystyle n} 0 {\displaystyle \mathbf {0} } n {\displaystyle n} คุณ {\displaystyle \mathbf {u} } ( เอ − λ ฉัน ) {\displaystyle (A-\lambda I)} A {\displaystyle A} n {\displaystyle n} A {\displaystyle A} D {\displaystyle D} M {\displaystyle M} A {\displaystyle A} D = M − 1 A M {\displaystyle D=M^{-1}AM} D {\displaystyle D} A {\displaystyle A} M {\displaystyle M} A {\displaystyle A}
ในทางกลับกัน หากไม่มีเวกเตอร์ลักษณะเฉพาะที่เป็นอิสระเชิงเส้นที่เกี่ยวข้องกับมัน ก็จะไม่สามารถทำให้เป็นแนวทแยงได้[ 26 ] [ 27 ] A {\displaystyle A} n {\displaystyle n} A {\displaystyle A}
นิยาม: เวกเตอร์คือเวกเตอร์ลักษณะเฉพาะทั่วไปที่มีอันดับ m ของเมทริกซ์และสอดคล้องกับค่าลักษณะเฉพาะถ้า x m {\displaystyle \mathbf {x} _{m}} A {\displaystyle A} λ {\displaystyle \lambda }
( A − λ I ) m x m = 0 {\displaystyle (A-\lambda I)^{m}\mathbf {x} _{m}=\mathbf {0} } แต่
( A − λ I ) m − 1 x m ≠ 0 . {\displaystyle (A-\lambda I)^{m-1}\mathbf {x} _{m}\neq \mathbf {0} .} [ 28 ] เห็นได้ชัดว่า เวกเตอร์ลักษณะเฉพาะทั่วไปที่มีอันดับ 1 คือ เวกเตอร์ลักษณะเฉพาะธรรมดา[ 29 ] เมทริกซ์ ทุก เมทริกซ์ × มีเวกเตอร์ลักษณะเฉพาะทั่วไปที่เป็นอิสระเชิงเส้นที่เกี่ยวข้องกับมัน และสามารถแสดงให้เห็นว่าคล้ายกับเมทริกซ์ "เกือบทแยงมุม" ในรูปแบบปกติของจอร์แดน[ 30 ] นั่นคือ มีเมทริกซ์ผกผันได้อยู่เช่นนั้น[ 31 ] เมท ริกซ์ในกรณีนี้เรียกว่าเมทริกซ์โมดอลทั่วไป สำหรับ[ 32 ] ถ้าเป็นค่าลักษณะเฉพาะที่มีความซ้ำซ้อนเชิงพีชคณิตแล้วจะมี เวกเตอร์ลักษณะเฉพาะทั่วไปที่เป็น อิสระเชิงเส้นที่สอดคล้องกับ[ 33 ] ผลลัพธ์เหล่านี้ในทางกลับกัน ให้วิธีการที่ตรงไปตรงมาสำหรับ การ คำนวณฟังก์ชันเมทริกซ์บางอย่างของ[ 34 ] n {\displaystyle n} n {\displaystyle n} A {\displaystyle A} n {\displaystyle n} J {\displaystyle J} M {\displaystyle M} J = M − 1 A M {\displaystyle J=M^{-1}AM} M {\displaystyle M} A {\displaystyle A} λ {\displaystyle \lambda } μ {\displaystyle \mu } A {\displaystyle A} μ {\displaystyle \mu } λ {\displaystyle \lambda } A {\displaystyle A}
หมายเหตุ: สำหรับเมทริกซ์เหนือฟิลด์ ที่จะแสดงในรูปแบบปกติของจอร์แดน ค่าลักษณะเฉพาะทั้งหมดของเมทริกซ์จะต้องอยู่ในฟิลด์นั้น นั่นคือพหุนามลักษณะเฉพาะ จะต้องแยกตัวประกอบได้อย่างสมบูรณ์เป็นตัวประกอบเชิงเส้น ฟิลด์จะต้องเป็นฟิลด์ปิดเชิงพีชคณิต ตัวอย่างเช่น ถ้า เมทริกซ์ มี องค์ประกอบ ที่เป็นค่าจริง อาจจำเป็นต้องให้ค่าลักษณะเฉพาะและส่วนประกอบของเวกเตอร์ลักษณะเฉพาะมีค่าเป็น จำนวนเชิงซ้อน[ 35 ] [ 36 ] [ 37 ] n × n {\displaystyle n\times n} A {\displaystyle A} F {\displaystyle F} A {\displaystyle A} F {\displaystyle F} f ( x ) {\displaystyle f(x)} F {\displaystyle F} A {\displaystyle A}
เซตที่ครอบคลุม โดยเวกเตอร์ลักษณะทั่วไปทั้งหมดสำหรับรูปแบบ ที่กำหนดจะทำให้ เกิด ปริภูมิลักษณะทั่วไป สำหรับ[ 38 ] λ {\displaystyle \lambda } λ {\displaystyle \lambda }
ตัวอย่าง ต่อไปนี้เป็นตัวอย่างบางส่วนเพื่ออธิบายแนวคิดของเวกเตอร์ลักษณะเฉพาะทั่วไป รายละเอียดบางส่วนจะอธิบายในภายหลัง
ตัวอย่างที่ 1 ตัวอย่างนี้เรียบง่ายแต่แสดงให้เห็นประเด็นได้อย่างชัดเจน เมทริกซ์ประเภทนี้ถูกใช้บ่อยในตำราเรียน[ 39 ] [ 40 ] [ 41 ] สมมติว่า
A = ( 1 1 0 1 ) . {\displaystyle A={\begin{pmatrix}1&1\\0&1\end{pmatrix}}.} ดังนั้นจึงมีค่าลักษณะเฉพาะเพียงค่าเดียว คือและค่าความซ้ำซ้อนเชิงพีชคณิตของค่าลักษณะเฉพาะนี้คือ. λ = 1 {\displaystyle \lambda =1} m = 2 {\displaystyle m=2}
โปรดสังเกตว่าเมทริกซ์นี้อยู่ในรูปแบบปกติของจอร์แดน แต่ไม่ใช่เมทริกซ์ทแยงมุม ดังนั้น เมทริกซ์นี้จึงไม่สามารถทำให้เป็นเมทริกซ์ทแยงมุมได้ เนื่องจากมี องค์ประกอบ เหนือแนวทแยงมุม หนึ่ง ตัว จึงจะมีเวกเตอร์ลักษณะเฉพาะทั่วไปที่มีอันดับมากกว่า 1 เพียงหนึ่งตัว (หรืออาจสังเกตได้ว่าปริภูมิเวกเตอร์มีมิติ 2 ดังนั้นจึงมีเวกเตอร์ลักษณะเฉพาะทั่วไปที่มีอันดับมากกว่า 1 ได้มากที่สุดหนึ่งตัว) อีกทางเลือกหนึ่งคือ สามารถคำนวณมิติของปริภูมิว่าง ของ เมทริกซ์ ได้เป็นและดังนั้นจึงมีเวกเตอร์ลักษณะเฉพาะทั่วไปที่มีอันดับมากกว่า 1 V {\displaystyle V} A − λ I {\displaystyle A-\lambda I} p = 1 {\displaystyle p=1} m − p = 1 {\displaystyle m-p=1}
เวกเตอร์ลักษณะเฉพาะแบบธรรมดาจะคำนวณตามปกติ (ดู ตัวอย่างได้ในหน้า เวกเตอร์ลักษณะเฉพาะ ) โดยใช้เวกเตอร์ลักษณะเฉพาะนี้ เราจะคำนวณเวกเตอร์ลักษณะเฉพาะแบบทั่วไป โดยการแก้สมการ v 1 = ( 1 0 ) {\displaystyle \mathbf {v} _{1}={\begin{pmatrix}1\\0\end{pmatrix}}} v 2 {\displaystyle \mathbf {v} _{2}}
( A − λ I ) v 2 = v 1 . {\displaystyle (A-\lambda I)\mathbf {v} _{2}=\mathbf {v} _{1}.} การเขียนค่าต่างๆ ออกมา:
( ( 1 1 0 1 ) − 1 ( 1 0 0 1 ) ) ( v 21 v 22 ) = ( 0 1 0 0 ) ( v 21 v 22 ) = ( 1 0 ) . {\displaystyle \left({\begin{pmatrix}1&1\\0&1\end{pmatrix}}-1{\begin{pmatrix}1&0\\0&1\end{pmatrix}}\right){\begin{pmatrix}v_{21}\\v_{22}\end{pmatrix}}={\begin{pmatrix}0&1\\0&0\end{pmatrix}}{\begin{pmatrix}v_{21}\\v_{22}\end{pmatrix}}={\begin{pmatrix}1\\0\end{pmatrix}}.} สิ่งนี้ทำให้ง่ายขึ้นเป็น
v 22 = 1. {\displaystyle v_{22}=1.} องค์ประกอบนี้ไม่มีข้อจำกัดใดๆ เวกเตอร์ลักษณะเฉพาะทั่วไปอันดับ 2 คือโดยที่a สามารถมีค่าสเกลาร์ใดๆ ก็ได้ การเลือกa = 0 มักจะเป็นวิธีที่ง่ายที่สุด v 21 {\displaystyle v_{21}} v 2 = ( a 1 ) {\displaystyle \mathbf {v} _{2}={\begin{pmatrix}a\\1\end{pmatrix}}}
โปรดทราบว่า
( A − λ I ) v 2 = ( 0 1 0 0 ) ( a 1 ) = ( 1 0 ) = v 1 , {\displaystyle (A-\lambda I)\mathbf {v} _{2}={\begin{pmatrix}0&1\\0&0\end{pmatrix}}{\begin{pmatrix}a\\1\end{pmatrix}}={\begin{pmatrix}1\\0\end{pmatrix}}=\mathbf {v} _{1},} ดังนั้นนี่จึงเป็นเวกเตอร์ลักษณะเฉพาะทั่วไป เพราะว่า v 2 {\displaystyle \mathbf {v} _{2}}
( A − λ I ) 2 v 2 = ( A − λ I ) [ ( A − λ I ) v 2 ] = ( A − λ I ) v 1 = ( 0 1 0 0 ) ( 1 0 ) = ( 0 0 ) = 0 , {\displaystyle (A-\lambda I)^{2}\mathbf {v} _{2}=(A-\lambda I)[(A-\lambda I)\mathbf {v} _{2}]=(A-\lambda I)\mathbf {v} _{1}={\begin{pmatrix}0&1\\0&0\end{pmatrix}}{\begin{pmatrix}1\\0\end{pmatrix}}={\begin{pmatrix}0\\0\end{pmatrix}}=\mathbf {0} ,} ดังนั้น จึงเป็นเวกเตอร์ลักษณะเฉพาะธรรมดา และและเป็นอิสระเชิงเส้น และด้วยเหตุนี้จึงเป็นฐานสำหรับปริภูมิเวกเตอร์ v 1 {\displaystyle \mathbf {v} _{1}} v 1 {\displaystyle \mathbf {v} _{1}} v 2 {\displaystyle \mathbf {v} _{2}} V {\displaystyle V}
ตัวอย่างที่ 2 ตัวอย่างนี้ซับซ้อนกว่าตัวอย่างที่ 1 น่าเสียดายที่การสร้างตัวอย่างที่น่าสนใจที่มีลำดับต่ำนั้นค่อนข้างยาก[ 42 ] เมทริกซ์
A = ( 1 0 0 0 0 3 1 0 0 0 6 3 2 0 0 10 6 3 2 0 15 10 6 3 2 ) {\displaystyle A={\begin{pmatrix}1&0&0&0&0\\3&1&0&0&0\\6&3&2&0&0\\10&6&3&2&0\\15&10&6&3&2\end{pmatrix}}} มีค่าลักษณะเฉพาะ และที่มีความซ้ำซ้อนเชิงพีชคณิต และแต่มี ความซ้ำซ้อนเชิงเรขาคณิต และλ 1 = 1 {\displaystyle \lambda _{1}=1} λ 2 = 2 {\displaystyle \lambda _{2}=2} μ 1 = 2 {\displaystyle \mu _{1}=2} μ 2 = 3 {\displaystyle \mu _{2}=3} γ 1 = 1 {\displaystyle \gamma _{1}=1} γ 2 = 1 {\displaystyle \gamma _{2}=1}
ปริภูมิ ไอเก น ทั่วไป ของคำนวณได้ดังนี้ คือเวกเตอร์ไอเก น ปกติที่เกี่ยวข้องกับคือเวกเตอร์ไอเกนทั่วไปที่เกี่ยวข้องกับคือเวกเตอร์ไอเกนปกติที่เกี่ยวข้องกับและ คือเวกเตอร์ไอเกนทั่วไปที่เกี่ยวข้องกับ A {\displaystyle A} x 1 {\displaystyle \mathbf {x} _{1}} λ 1 {\displaystyle \lambda _{1}} x 2 {\displaystyle \mathbf {x} _{2}} λ 1 {\displaystyle \lambda _{1}} y 1 {\displaystyle \mathbf {y} _{1}} λ 2 {\displaystyle \lambda _{2}} y 2 {\displaystyle \mathbf {y} _{2}} y 3 {\displaystyle \mathbf {y} _{3}} λ 2 {\displaystyle \lambda _{2}}
( A − 1 I ) x 1 = ( 0 0 0 0 0 3 0 0 0 0 6 3 1 0 0 10 6 3 1 0 15 10 6 3 1 ) ( 0 3 − 9 9 − 3 ) = ( 0 0 0 0 0 ) = 0 , {\displaystyle (A-1I)\mathbf {x} _{1}={\begin{pmatrix}0&0&0&0&0\\3&0&0&0&0\\6&3&1&0&0\\10&6&3&1&0\\15&10&6&3&1\end{pmatrix}}{\begin{pmatrix}0\\3\\-9\\9\\-3\end{pmatrix}}={\begin{pmatrix}0\\0\\0\\0\\0\end{pmatrix}}=\mathbf {0} ,} ( A − 1 I ) x 2 = ( 0 0 0 0 0 3 0 0 0 0 6 3 1 0 0 10 6 3 1 0 15 10 6 3 1 ) ( 1 − 15 30 − 1 − 45 ) = ( 0 3 − 9 9 − 3 ) = x 1 , {\displaystyle (A-1I)\mathbf {x} _{2}={\begin{pmatrix}0&0&0&0&0\\3&0&0&0&0\\6&3&1&0&0\\10&6&3&1&0\\15&10&6&3&1\end{pmatrix}}{\begin{pmatrix}1\\-15\\30\\-1\\-45\end{pmatrix}}={\begin{pmatrix}0\\3\\-9\\9\\-3\end{pmatrix}}=\mathbf {x} _{1},} ( A − 2 I ) y 1 = ( − 1 0 0 0 0 3 − 1 0 0 0 6 3 0 0 0 10 6 3 0 0 15 10 6 3 0 ) ( 0 0 0 0 9 ) = ( 0 0 0 0 0 ) = 0 , {\displaystyle (A-2I)\mathbf {y} _{1}={\begin{pmatrix}-1&0&0&0&0\\3&-1&0&0&0\\6&3&0&0&0\\10&6&3&0&0\\15&10&6&3&0\end{pmatrix}}{\begin{pmatrix}0\\0\\0\\0\\9\end{pmatrix}}={\begin{pmatrix}0\\0\\0\\0\\0\end{pmatrix}}=\mathbf {0} ,} ( A − 2 I ) y 2 = ( − 1 0 0 0 0 3 − 1 0 0 0 6 3 0 0 0 10 6 3 0 0 15 10 6 3 0 ) ( 0 0 0 3 0 ) = ( 0 0 0 0 9 ) = y 1 , {\displaystyle (A-2I)\mathbf {y} _{2}={\begin{pmatrix}-1&0&0&0&0\\3&-1&0&0&0\\6&3&0&0&0\\10&6&3&0&0\\15&10&6&3&0\end{pmatrix}}{\begin{pmatrix}0\\0\\0\\3\\0\end{pmatrix}}={\begin{pmatrix}0\\0\\0\\0\\9\end{pmatrix}}=\mathbf {y} _{1},} ( A − 2 I ) y 3 = ( − 1 0 0 0 0 3 − 1 0 0 0 6 3 0 0 0 10 6 3 0 0 15 10 6 3 0 ) ( 0 0 1 − 2 0 ) = ( 0 0 0 3 0 ) = y 2 . {\displaystyle (A-2I)\mathbf {y} _{3}={\begin{pmatrix}-1&0&0&0&0\\3&-1&0&0&0\\6&3&0&0&0\\10&6&3&0&0\\15&10&6&3&0\end{pmatrix}}{\begin{pmatrix}0\\0\\1\\-2\\0\end{pmatrix}}={\begin{pmatrix}0\\0\\0\\3\\0\end{pmatrix}}=\mathbf {y} _{2}.} ผลลัพธ์ที่ได้คือฐานสำหรับปริภูมิเวกเตอร์ลักษณะเฉพาะทั่วไป ของ แต่ละตัว เมื่อรวมกันแล้ว โซ่ ของเวกเตอร์ลักษณะเฉพาะทั่วไป ทั้งสองจะครอบคลุมปริภูมิของเวกเตอร์คอลัมน์ 5 มิติทั้งหมด A {\displaystyle A}
{ x 1 , x 2 } = { ( 0 3 − 9 9 − 3 ) , ( 1 − 15 30 − 1 − 45 ) } , { y 1 , y 2 , y 3 } = { ( 0 0 0 0 9 ) , ( 0 0 0 3 0 ) , ( 0 0 1 − 2 0 ) } . {\displaystyle \left\{\mathbf {x} _{1},\mathbf {x} _{2}\right\}=\left\{{\begin{pmatrix}0\\3\\-9\\9\\-3\end{pmatrix}},{\begin{pmatrix}1\\-15\\30\\-1\\-45\end{pmatrix}}\right\},\left\{\mathbf {y} _{1},\mathbf {y} _{2},\mathbf {y} _{3}\right\}=\left\{{\begin{pmatrix}0\\0\\0\\0\\9\end{pmatrix}},{\begin{pmatrix}0\\0\\0\\3\\0\end{pmatrix}},{\begin{pmatrix}0\\0\\1\\-2\\0\end{pmatrix}}\right\}.} เมทริกซ์ "เกือบทแยงมุม" ในรูปแบบปกติของจอร์แดน ซึ่งคล้ายคลึงกับได้มาดังนี้: J {\displaystyle J} A {\displaystyle A}
M = ( x 1 x 2 y 1 y 2 y 3 ) = ( 0 1 0 0 0 3 − 15 0 0 0 − 9 30 0 0 1 9 − 1 0 3 − 2 − 3 − 45 9 0 0 ) , {\displaystyle M={\begin{pmatrix}\mathbf {x} _{1}&\mathbf {x} _{2}&\mathbf {y} _{1}&\mathbf {y} _{2}&\mathbf {y} _{3}\end{pmatrix}}={\begin{pmatrix}0&1&0&0&0\\3&-15&0&0&0\\-9&30&0&0&1\\9&-1&0&3&-2\\-3&-45&9&0&0\end{pmatrix}},} J = ( 1 1 0 0 0 0 1 0 0 0 0 0 2 1 0 0 0 0 2 1 0 0 0 0 2 ) , {\displaystyle J={\begin{pmatrix}1&1&0&0&0\\0&1&0&0&0\\0&0&2&1&0\\0&0&0&2&1\\0&0&0&0&2\end{pmatrix}},} โดยที่เมท ริก ซ์ โมดอลทั่วไป สำหรับคอลัมน์ของเป็นฐานแคนอนิก สำหรับและ[ 43 ] M {\displaystyle M} A {\displaystyle A} M {\displaystyle M} A {\displaystyle A} A M = M J {\displaystyle AM=MJ}
โซ่จอร์แดน นิยาม: ให้เป็นเวกเตอร์ลักษณะเฉพาะทั่วไปอันดับm ที่สอดคล้องกับเมทริกซ์และค่าลักษณะเฉพาะโซ่ที่สร้างโดย คือเซตของเวกเตอร์ที่กำหนดโดย x m {\displaystyle \mathbf {x} _{m}} A {\displaystyle A} λ {\displaystyle \lambda } x m {\displaystyle \mathbf {x} _{m}} { x m , x m − 1 , … , x 1 } {\displaystyle \left\{\mathbf {x} _{m},\mathbf {x} _{m-1},\dots ,\mathbf {x} _{1}\right\}}
โดยที่เป็นเวกเตอร์ลักษณะเฉพาะธรรมดาที่มีค่าลักษณะเฉพาะที่กำหนดให้เสมอดังนั้น โดยทั่วไปแล้ว x 1 {\displaystyle \mathbf {x} _{1}} λ {\displaystyle \lambda }
x j = ( A − λ I ) m − j x m = ( A − λ I ) x j + 1 ( j = 1 , 2 , … , m − 1 ) . {\displaystyle \mathbf {x} _{j}=(A-\lambda I)^{m-j}\mathbf {x} _{m}=(A-\lambda I)\mathbf {x} _{j+1}\qquad (j=1,2,\dots ,m-1).} 2
เวกเตอร์ที่กำหนดโดย ( 2 ) เป็นเวกเตอร์ลักษณะเฉพาะทั่วไปของอันดับj ที่สอดคล้องกับค่าลักษณะเฉพาะโซ่คือเซตของเวกเตอร์ที่เป็นอิสระเชิงเส้น[ 44 ] x j {\displaystyle \mathbf {x} _{j}} λ {\displaystyle \lambda }
ฐานหลัก นิยาม: เซตของเวกเตอร์ลักษณะเฉพาะทั่วไปที่เป็นอิสระเชิงเส้นจำนวนn ตัว จะเป็น ฐานหลักมาตรฐาน (canonical basis) ก็ต่อเมื่อเซตนั้นประกอบขึ้นจากโซ่จอร์แดน (Jordan chains) ทั้งหมด
ดังนั้น เมื่อเรากำหนดแล้วว่าเวกเตอร์ลักษณะเฉพาะทั่วไปที่มีอันดับm อยู่ในฐานมาตรฐาน ก็จะสรุปได้ว่า เวกเตอร์ m − 1 ที่อยู่ในโซ่จอร์แดนที่สร้างขึ้นโดยก็อยู่ในฐานมาตรฐานเช่นกัน[ 45 ] x m − 1 , x m − 2 , … , x 1 {\displaystyle \mathbf {x} _{m-1},\mathbf {x} _{m-2},\ldots ,\mathbf {x} _{1}} x m {\displaystyle \mathbf {x} _{m}}
ให้ λ เป็นค่าลักษณะเฉพาะของเมทริกซ์ λ ที่มีความซ้ำซ้อนเชิงพีชคณิตเท่ากับ n ก่อนอื่น ให้หาอันดับ (อันดับเมทริกซ์) ของเมทริกซ์ λ โดยกำหนดให้ λ เป็นจำนวนเต็มตัวแรก ที่ λ มีอันดับเท่ากับ n ( โดยที่ n คือจำนวนแถวหรือคอลัมน์ของλ นั่นคือ λ มีขนาดn × n ) λ i {\displaystyle \lambda _{i}} A {\displaystyle A} μ i {\displaystyle \mu _{i}} ( A − λ i I ) , ( A − λ i I ) 2 , … , ( A − λ i I ) m i {\displaystyle (A-\lambda _{i}I),(A-\lambda _{i}I)^{2},\ldots ,(A-\lambda _{i}I)^{m_{i}}} m i {\displaystyle m_{i}} ( A − λ i I ) m i {\displaystyle (A-\lambda _{i}I)^{m_{i}}} n − μ i {\displaystyle n-\mu _{i}} A {\displaystyle A} A {\displaystyle A}
ตอนนี้กำหนดคำจำกัดความ
ρ k = rank ( A − λ i I ) k − 1 − rank ( A − λ i I ) k ( k = 1 , 2 , … , m i ) . {\displaystyle \rho _{k}=\operatorname {rank} (A-\lambda _{i}I)^{k-1}-\operatorname {rank} (A-\lambda _{i}I)^{k}\qquad (k=1,2,\ldots ,m_{i}).} ตัวแปรนี้ระบุจำนวนเวกเตอร์ลักษณะเฉพาะทั่วไปที่เป็นอิสระเชิงเส้นที่มีอันดับk ซึ่งสอดคล้องกับค่าลักษณะเฉพาะที่จะปรากฏในฐานมาตรฐานสำหรับโปรดทราบว่า ρ k {\displaystyle \rho _{k}} λ i {\displaystyle \lambda _{i}} A {\displaystyle A}
rank ( A − λ i I ) 0 = rank ( I ) = n {\displaystyle \operatorname {rank} (A-\lambda _{i}I)^{0}=\operatorname {rank} (I)=n} [ 46 ]
การคำนวณเวกเตอร์ลักษณะเฉพาะทั่วไป ในส่วนก่อนหน้านี้ เราได้เห็นเทคนิคในการหาเวกเตอร์ลักษณะเฉพาะทั่วไปที่เป็นอิสระเชิงเส้นของฐานมาตรฐานสำหรับปริภูมิเวกเตอร์ที่เกี่ยวข้องกับเมทริกซ์เทคนิคเหล่านี้สามารถนำมารวมกันเป็นขั้นตอนได้ดังนี้: n {\displaystyle n} V {\displaystyle V} n × n {\displaystyle n\times n} A {\displaystyle A}
จงแก้สมการลักษณะเฉพาะ เพื่อหาค่าไอเกนและจำนวนเท่าเชิงพีชคณิตของ ค่าไอเก น เหล่านั้นA {\displaystyle A} λ i {\displaystyle \lambda _{i}} μ i {\displaystyle \mu _{i}} สำหรับแต่ละรายการλ i : {\displaystyle \lambda _{i}:} กำหนด;n − μ i {\displaystyle n-\mu _{i}} กำหนด;m i {\displaystyle m_{i}} กำหนดค่าสำหรับ;ρ k {\displaystyle \rho _{k}} ( k = 1 , … , m i ) {\displaystyle (k=1,\ldots ,m_{i})} กำหนดโซ่จอร์แดนแต่ละอันสำหรับ;λ i {\displaystyle \lambda _{i}}
ตัวอย่างที่ 3 เมทริกซ์
A = ( 5 1 − 2 4 0 5 2 2 0 0 5 3 0 0 0 4 ) {\displaystyle A={\begin{pmatrix}5&1&-2&4\\0&5&2&2\\0&0&5&3\\0&0&0&4\end{pmatrix}}} มีค่าไอเกนที่มีความซ้ำซ้อนเชิงพีชคณิตและค่าไอเกนที่มีความซ้ำซ้อนเชิงพีชคณิตนอกจากนี้เรายังมีเนื่องจากเรามี λ 1 = 5 {\displaystyle \lambda _{1}=5} μ 1 = 3 {\displaystyle \mu _{1}=3} λ 2 = 4 {\displaystyle \lambda _{2}=4} μ 2 = 1 {\displaystyle \mu _{2}=1} n = 4 {\displaystyle n=4} λ 1 {\displaystyle \lambda _{1}} n − μ 1 = 4 − 3 = 1 {\displaystyle n-\mu _{1}=4-3=1}
( A − 5 I ) = ( 0 1 − 2 4 0 0 2 2 0 0 0 3 0 0 0 − 1 ) , rank ( A − 5 I ) = 3. {\displaystyle (A-5I)={\begin{pmatrix}0&1&-2&4\\0&0&2&2\\0&0&0&3\\0&0&0&-1\end{pmatrix}},\qquad \operatorname {rank} (A-5I)=3.} ( A − 5 I ) 2 = ( 0 0 2 − 8 0 0 0 4 0 0 0 − 3 0 0 0 1 ) , rank ( A − 5 I ) 2 = 2. {\displaystyle (A-5I)^{2}={\begin{pmatrix}0&0&2&-8\\0&0&0&4\\0&0&0&-3\\0&0&0&1\end{pmatrix}},\qquad \operatorname {rank} (A-5I)^{2}=2.} ( A − 5 I ) 3 = ( 0 0 0 14 0 0 0 − 4 0 0 0 3 0 0 0 − 1 ) , rank ( A − 5 I ) 3 = 1. {\displaystyle (A-5I)^{3}={\begin{pmatrix}0&0&0&14\\0&0&0&-4\\0&0&0&3\\0&0&0&-1\end{pmatrix}},\qquad \operatorname {rank} (A-5I)^{3}=1.} จำนวนเต็มตัวแรกที่มี อันดับคือm 1 {\displaystyle m_{1}} ( A − 5 I ) m 1 {\displaystyle (A-5I)^{m_{1}}} n − μ 1 = 1 {\displaystyle n-\mu _{1}=1} m 1 = 3 {\displaystyle m_{1}=3}
ตอนนี้เราจะกำหนด
ρ 3 = rank ( A − 5 I ) 2 − rank ( A − 5 I ) 3 = 2 − 1 = 1 , {\displaystyle \rho _{3}=\operatorname {rank} (A-5I)^{2}-\operatorname {rank} (A-5I)^{3}=2-1=1,} ρ 2 = rank ( A − 5 I ) 1 − rank ( A − 5 I ) 2 = 3 − 2 = 1 , {\displaystyle \rho _{2}=\operatorname {rank} (A-5I)^{1}-\operatorname {rank} (A-5I)^{2}=3-2=1,} ρ 1 = rank ( A − 5 I ) 0 − rank ( A − 5 I ) 1 = 4 − 3 = 1. {\displaystyle \rho _{1}=\operatorname {rank} (A-5I)^{0}-\operatorname {rank} (A-5I)^{1}=4-3=1.} ดังนั้น จะมีเวกเตอร์ลักษณะเฉพาะทั่วไปที่เป็นอิสระเชิงเส้นสามตัว โดยแต่ละตัวมีอันดับ 3, 2 และ 1 เนื่องจากสอดคล้องกับสายโซ่เดียวของเวกเตอร์ลักษณะเฉพาะทั่วไปที่เป็นอิสระเชิงเส้นสามตัว เราจึงทราบว่ามีเวกเตอร์ลักษณะเฉพาะทั่วไปอันดับ 3 ที่สอดคล้องกับเช่นนั้น λ 1 {\displaystyle \lambda _{1}} x 3 {\displaystyle \mathbf {x} _{3}} λ 1 {\displaystyle \lambda _{1}}
( A − 5 I ) 3 x 3 = 0 {\displaystyle (A-5I)^{3}\mathbf {x} _{3}=\mathbf {0} } 3
แต่
( A − 5 I ) 2 x 3 ≠ 0 . {\displaystyle (A-5I)^{2}\mathbf {x} _{3}\neq \mathbf {0} .} 4
สมการ ( 3 ) และ ( 4 ) แสดงถึงระบบเชิงเส้น ที่สามารถแก้หาค่าได้ให้ x 3 {\displaystyle \mathbf {x} _{3}}
x 3 = ( x 31 x 32 x 33 x 34 ) . {\displaystyle \mathbf {x} _{3}={\begin{pmatrix}x_{31}\\x_{32}\\x_{33}\\x_{34}\end{pmatrix}}.} แล้ว
( A − 5 I ) 3 x 3 = ( 0 0 0 14 0 0 0 − 4 0 0 0 3 0 0 0 − 1 ) ( x 31 x 32 x 33 x 34 ) = ( 14 x 34 − 4 x 34 3 x 34 − x 34 ) = ( 0 0 0 0 ) {\displaystyle (A-5I)^{3}\mathbf {x} _{3}={\begin{pmatrix}0&0&0&14\\0&0&0&-4\\0&0&0&3\\0&0&0&-1\end{pmatrix}}{\begin{pmatrix}x_{31}\\x_{32}\\x_{33}\\x_{34}\end{pmatrix}}={\begin{pmatrix}14x_{34}\\-4x_{34}\\3x_{34}\\-x_{34}\end{pmatrix}}={\begin{pmatrix}0\\0\\0\\0\end{pmatrix}}} และ
( A − 5 I ) 2 x 3 = ( 0 0 2 − 8 0 0 0 4 0 0 0 − 3 0 0 0 1 ) ( x 31 x 32 x 33 x 34 ) = ( 2 x 33 − 8 x 34 4 x 34 − 3 x 34 x 34 ) ≠ ( 0 0 0 0 ) . {\displaystyle (A-5I)^{2}\mathbf {x} _{3}={\begin{pmatrix}0&0&2&-8\\0&0&0&4\\0&0&0&-3\\0&0&0&1\end{pmatrix}}{\begin{pmatrix}x_{31}\\x_{32}\\x_{33}\\x_{34}\end{pmatrix}}={\begin{pmatrix}2x_{33}-8x_{34}\\4x_{34}\\-3x_{34}\\x_{34}\end{pmatrix}}\neq {\begin{pmatrix}0\\0\\0\\0\end{pmatrix}}.} ดังนั้น เพื่อให้เป็นไปตามเงื่อนไข ( 3 ) และ ( 4 ) เราต้องมีและไม่มีข้อจำกัดใดๆ สำหรับและโดยการเลือกเราจะได้ x 34 = 0 {\displaystyle x_{34}=0} x 33 ≠ 0 {\displaystyle x_{33}\neq 0} x 31 {\displaystyle x_{31}} x 32 {\displaystyle x_{32}} x 31 = x 32 = x 34 = 0 , x 33 = 1 {\displaystyle x_{31}=x_{32}=x_{34}=0,x_{33}=1}
x 3 = ( 0 0 1 0 ) {\displaystyle \mathbf {x} _{3}={\begin{pmatrix}0\\0\\1\\0\end{pmatrix}}} ในฐานะเวกเตอร์ลักษณะทั่วไปอันดับ 3 ที่สอดคล้องกับโปรดทราบว่าเป็นไปได้ที่จะได้รับเวกเตอร์ลักษณะทั่วไปอันดับ 3 อื่นๆ อีกมากมายโดยการเลือกค่าที่แตกต่างกันของ , และโดยที่อย่างไรก็ตาม ตัวเลือกแรกของเราเป็นตัวเลือกที่ง่ายที่สุด[ 47 ] λ 1 = 5 {\displaystyle \lambda _{1}=5} x 31 {\displaystyle x_{31}} x 32 {\displaystyle x_{32}} x 33 {\displaystyle x_{33}} x 33 ≠ 0 {\displaystyle x_{33}\neq 0}
เมื่อใช้สมการ ( 1 ) เราจะได้และเป็นเวกเตอร์ลักษณะเฉพาะทั่วไปที่มีอันดับ 2 และ 1 ตามลำดับ โดยที่ x 2 {\displaystyle \mathbf {x} _{2}} x 1 {\displaystyle \mathbf {x} _{1}}
x 2 = ( A − 5 I ) x 3 = ( − 2 2 0 0 ) , {\displaystyle \mathbf {x} _{2}=(A-5I)\mathbf {x} _{3}={\begin{pmatrix}-2\\2\\0\\0\end{pmatrix}},} และ
x 1 = ( A − 5 I ) x 2 = ( 2 0 0 0 ) . {\displaystyle \mathbf {x} _{1}=(A-5I)\mathbf {x} _{2}={\begin{pmatrix}2\\0\\0\\0\end{pmatrix}}.} ค่าลักษณะเฉพาะอย่างง่าย สามารถจัดการได้โดยใช้เทคนิคมาตรฐาน และมีเวกเตอร์ลักษณะเฉพาะแบบปกติ λ 2 = 4 {\displaystyle \lambda _{2}=4}
y 1 = ( − 14 4 − 3 1 ) . {\displaystyle \mathbf {y} _{1}={\begin{pmatrix}-14\\4\\-3\\1\end{pmatrix}}.} พื้นฐานที่เป็นหลักการสำหรับคือ A {\displaystyle A}
{ x 3 , x 2 , x 1 , y 1 } = { ( 0 0 1 0 ) ( − 2 2 0 0 ) ( 2 0 0 0 ) ( − 14 4 − 3 1 ) } . {\displaystyle \left\{\mathbf {x} _{3},\mathbf {x} _{2},\mathbf {x} _{1},\mathbf {y} _{1}\right\}=\left\{{\begin{pmatrix}0\\0\\1\\0\end{pmatrix}}{\begin{pmatrix}-2\\2\\0\\0\end{pmatrix}}{\begin{pmatrix}2\\0\\0\\0\end{pmatrix}}{\begin{pmatrix}-14\\4\\-3\\1\end{pmatrix}}\right\}.} x 1 , x 2 {\displaystyle \mathbf {x} _{1},\mathbf {x} _{2}} และคือเวกเตอร์ลักษณะเฉพาะทั่วไปที่เกี่ยวข้องกับในขณะที่คือเวกเตอร์ลักษณะเฉพาะธรรมดาที่เกี่ยวข้องกับ x 3 {\displaystyle \mathbf {x} _{3}} λ 1 {\displaystyle \lambda _{1}} y 1 {\displaystyle \mathbf {y} _{1}} λ 2 {\displaystyle \lambda _{2}}
นี่เป็นตัวอย่างที่ค่อนข้างง่าย โดยทั่วไปแล้ว จำนวนเวกเตอร์ลักษณะเฉพาะทั่วไปที่เป็นอิสระเชิงเส้นที่มีอันดับจะไม่เท่ากันเสมอไป กล่าวคือ อาจมีโซ่หลายเส้นที่มีความยาวต่างกันซึ่งสอดคล้องกับค่าลักษณะเฉพาะเฉพาะค่าหนึ่ง[ 48 ] ρ k {\displaystyle \rho _{k}} k {\displaystyle k}
เมทริกซ์โมดอลทั่วไป ให้เป็นเมทริกซ์ขนาดn × n เมทริกซ์โมดอลทั่วไป สำหรับคือ เมทริกซ์ขนาด n × n ที่คอลัมน์ของมัน เมื่อพิจารณาว่าเป็นเวกเตอร์ จะก่อให้เกิดฐานมาตรฐานสำหรับและปรากฏในตามกฎต่อไปนี้: A {\displaystyle A} M {\displaystyle M} A {\displaystyle A} A {\displaystyle A} M {\displaystyle M}
โซ่จอร์แดนทั้งหมดที่ประกอบด้วยเวกเตอร์เดียว (กล่าวคือ เวกเตอร์ที่มีความยาวเดียว) จะปรากฏในคอลัมน์แรกๆ ของตารางM {\displaystyle M} เวกเตอร์ทั้งหมดของสายโซ่หนึ่งๆ จะปรากฏอยู่ด้วยกันในคอลัมน์ที่อยู่ติดกันของ.M {\displaystyle M} แต่ละสายโซ่จะปรากฏตามลำดับที่เพิ่มขึ้น (นั่นคือ เวกเตอร์ลักษณะทั่วไปของลำดับที่ 1 จะปรากฏก่อนเวกเตอร์ลักษณะทั่วไปของลำดับที่ 2 ของสายโซ่เดียวกัน ซึ่งจะปรากฏก่อนเวกเตอร์ลักษณะทั่วไปของลำดับที่ 3 ของสายโซ่เดียวกัน เป็นต้น) [ 49 ] M {\displaystyle M}
[ ⌜ λ 1 1 λ 1 λ 1 ⌝ ⌜ λ 2 1 λ 2 ⌝ [ λ 3 ] ⋱ ⌜ λ n 1 λ n ⌝ ⌜ λ 1 1 λ 1 1 λ 1 ⌝ ⌜ λ 2 1 λ 2 ⌝ [ λ 3 ] ⋱ ⌜ λ n 1 λ n ⌝ ⌞ λ 1 1 λ 1 λ 1 ⌟ ⌜ λ 2 1 λ 2 ⌝ [ λ 3 ] ⋱ ⌜ λ n 1 λ n ⌝ ⌜ λ 1 1 λ 1 1 λ 1 ⌝ ⌜ λ 2 1 n ⌝ [ λ 3 ] ⋱ ⌜ λ n 1 λ n ⌝ ⌜ λ 1 1 λ 1 1 λ 1 ⌟ ⌞ λ 2 λ 2 ⌟ [ λ 3 ] ⋱ ⌜ λ n 1 λ n ⌝ ⌜ λ 1 1 λ 1 1 λ 1 ⌝ ⌜ λ 2 1 λ 2 ⌝ [ λ 3 ] ⋱ ⌜ λ n 1 λ n ⌝ ⌜ λ 1 1 λ 1 1 λ 1 ⌝ ⌜ λ 2 1 λ 2 ⌝ [ λ 3 ] ⋱ ⌜ λ n 1 λ n ⌝ ⌜ λ 1 1 λ 1 1 λ 1 ⌝ ⌜ λ 2 1 λ 2 ⌝ [ λ 3 ] ⋱ ⌜ λ n 1 n ⌝ ⌞ λ 1 1 λ 1 1 λ 1 ⌝ ⌜ λ 2 1 λ 2 ⌝ [ λ 3 ] ⋱ ⌞ λ n λ n ⌟ ] {\displaystyle {\begin{bmatrix}{\color {red}\ulcorner }\lambda _{1}1{\hphantom {\lambda _{1}\lambda _{1}}}{\color {red}\urcorner }{\hphantom {\ulcorner \lambda _{2}1\lambda _{2}\urcorner [\lambda _{3}]\ddots \ulcorner \lambda _{n}1\lambda _{n}\urcorner }}\\{\hphantom {\ulcorner \lambda _{1}1}}\lambda _{1}1{\hphantom {\lambda _{1}\urcorner \ulcorner \lambda _{2}1\lambda _{2}\urcorner [\lambda _{3}]\ddots \ulcorner \lambda _{n}1\lambda _{n}\urcorner }}\\{\color {red}\llcorner }{\hphantom {\lambda _{1}1\lambda _{1}}}\lambda _{1}{\color {red}\lrcorner }{\hphantom {\ulcorner \lambda _{2}1\lambda _{2}\urcorner [\lambda _{3}]\ddots \ulcorner \lambda _{n}1\lambda _{n}\urcorner }}\\{\hphantom {\ulcorner \lambda _{1}1\lambda _{1}1\lambda _{1}\urcorner }}{\color {red}\ulcorner }\lambda _{2}1{\hphantom {n}}{\color {red}\urcorner }{\hphantom {[\lambda _{3}]\ddots \ulcorner \lambda _{n}1\lambda _{n}\urcorner }}\\{\hphantom {\ulcorner \lambda _{1}1\lambda _{1}1\lambda _{1}\lrcorner }}{\color {red}\llcorner }{\hphantom {\lambda _{2}}}\lambda _{2}{\color {red}\lrcorner }{\hphantom {[\lambda _{3}]\ddots \ulcorner \lambda _{n}1\lambda _{n}\urcorner }}\\{\hphantom {\ulcorner \lambda _{1}1\lambda _{1}1\lambda _{1}\urcorner \ulcorner \lambda _{2}1\lambda _{2}\urcorner }}{\color {red}[}\lambda _{3}{\color {red}]}{\hphantom {\ddots \ulcorner \lambda _{n}1\lambda _{n}\urcorner }}\\{\hphantom {\ulcorner \lambda _{1}1\lambda _{1}1\lambda _{1}\urcorner \ulcorner \lambda _{2}1\lambda _{2}\urcorner [\lambda _{3}]}}\ddots {\hphantom {\ulcorner \lambda _{n}1\lambda _{n}\urcorner }}\\{\hphantom {\ulcorner \lambda _{1}1\lambda _{1}1\lambda _{1}\urcorner \ulcorner \lambda _{2}1\lambda _{2}\urcorner [\lambda _{3}]\ddots }}{\color {red}\ulcorner }\lambda _{n}1{\hphantom {n}}{\color {red}\urcorner }\\{\hphantom {\llcorner \lambda _{1}1\lambda _{1}1\lambda _{1}\urcorner \ulcorner \lambda _{2}1\lambda _{2}\urcorner [\lambda _{3}]\ddots }}{\color {red}\llcorner }{\hphantom {\lambda _{n}}}\lambda _{n}{\color {red}\lrcorner }\end{bmatrix}}}
ตัวอย่างของเมทริกซ์ในรูปแบบปกติของจอร์แดนบล็อกสีแดงเรียกว่าบล็อกจอร์แดน
ให้ V เป็น ปริภูมิเวกเตอร์ n มิติ ให้ V เป็นแผนที่เชิงเส้นในL ( V ) ซึ่งเป็นเซตของแผนที่เชิงเส้นทั้งหมดจาก V ไปยังตัวมันเอง และให้ V เป็นเมทริกซ์แทน V เทียบกับฐานเรียงลำดับบางฐาน สามารถแสดงได้ว่า ถ้าพหุนามลักษณะเฉพาะ ของV แยกตัวประกอบได้เป็นตัวประกอบเชิงเส้น ดังนั้น V จะมีรูปแบบ V {\displaystyle V} ϕ {\displaystyle \phi } V {\displaystyle V} A {\displaystyle A} ϕ {\displaystyle \phi } f ( λ ) {\displaystyle f(\lambda )} A {\displaystyle A} f ( λ ) {\displaystyle f(\lambda )}
f ( λ ) = ± ( λ − λ 1 ) μ 1 ( λ − λ 2 ) μ 2 ⋯ ( λ − λ r ) μ r , {\displaystyle f(\lambda )=\pm (\lambda -\lambda _{1})^{\mu _{1}}(\lambda -\lambda _{2})^{\mu _{2}}\cdots (\lambda -\lambda _{r})^{\mu _{r}},} โดยที่ค่าไอเกนที่แตกต่างกันของแต่ละค่าคือความซ้ำซ้อนเชิงพีชคณิตของค่าไอเกนที่สอดคล้องกันและคล้ายกับเมทริกซ์ในรูปแบบปกติของจอร์แดน โดยที่แต่ละค่าปรากฏต่อเนื่องกันบนแนวทแยง และค่าที่อยู่เหนือแต่ละค่าโดยตรง(นั่นคือ บนซูเปอร์ไดอะโกนัล ) จะเป็น 0 หรือ 1: ในแต่ละบล็อก ค่าที่อยู่เหนือการปรากฏครั้งแรกของแต่ละค่าจะเป็น 0 เสมอ (ยกเว้นในบล็อกแรก) ค่าอื่นๆ ทั้งหมดบนซูเปอร์ไดอะโกนัลจะเป็น 1 ค่าอื่นๆ ทั้งหมด (นั่นคือ นอกแนวทแยงและซูเปอร์ไดอะโกนัล) จะเป็น 0 (แต่ไม่มีการกำหนดลำดับใดๆ ระหว่างค่าไอเกน หรือระหว่างบล็อกสำหรับค่าไอเกนที่กำหนด) เมทริกซ์นี้ใกล้เคียงที่สุดเท่าที่จะเป็นไปได้กับเมทริกซ์ทแยงมุมของถ้าสามารถทำให้เป็นเมทริกซ์ทแยงมุมได้ ค่าทั้งหมดที่อยู่เหนือแนวทแยงมุมจะเป็นศูนย์[ 50 ] โปรดทราบว่าตำราเรียนบางเล่มมีค่าที่อยู่บนซับไดอะโกนัล นั่นคือ อยู่ใต้แนวทแยงมุมหลักโดยตรง แทนที่จะอยู่บนซูเปอร์ไดอะโกนัล ค่าไอเกนยังคงอยู่บนแนวทแยงหลัก[ 51 ] [ 52 ] λ 1 , λ 2 , … , λ r {\displaystyle \lambda _{1},\lambda _{2},\ldots ,\lambda _{r}} A {\displaystyle A} μ i {\displaystyle \mu _{i}} λ i {\displaystyle \lambda _{i}} A {\displaystyle A} J {\displaystyle J} λ i {\displaystyle \lambda _{i}} μ i {\displaystyle \mu _{i}} λ i {\displaystyle \lambda _{i}} λ i {\displaystyle \lambda _{i}} J {\displaystyle J} A {\displaystyle A} A {\displaystyle A}
เมทริกซ์n × n ทุก เมทริกซ์ จะคล้ายกับเมทริกซ์ในรูปแบบปกติของจอร์แดน ซึ่งได้มาจากการแปลงความคล้ายคลึงโดยที่เป็นเมทริกซ์โมดอลทั่วไปสำหรับ[ 53 ] ( ดูหมายเหตุ ข้างต้น) A {\displaystyle A} J {\displaystyle J} J = M − 1 A M {\displaystyle J=M^{-1}AM} M {\displaystyle M} A {\displaystyle A}
ตัวอย่างที่ 4 จงหาเมทริกซ์ในรูปแบบปกติของจอร์แดนที่คล้ายกับ
A = ( 0 4 2 − 3 8 3 4 − 8 − 2 ) . {\displaystyle A={\begin{pmatrix}0&4&2\\-3&8&3\\4&-8&-2\end{pmatrix}}.} คำตอบ: สมการลักษณะเฉพาะของคือดังนั้น จึงเป็นค่าลักษณะเฉพาะที่มีความซ้ำซ้อนเชิงพีชคณิตเท่ากับสาม เมื่อปฏิบัติตามขั้นตอนในส่วนก่อนหน้า เราพบว่า A {\displaystyle A} ( λ − 2 ) 3 = 0 {\displaystyle (\lambda -2)^{3}=0} λ = 2 {\displaystyle \lambda =2}
rank ( A − 2 I ) = 1 {\displaystyle \operatorname {rank} (A-2I)=1} และ
rank ( A − 2 I ) 2 = 0 = n − μ . {\displaystyle \operatorname {rank} (A-2I)^{2}=0=n-\mu .} ดังนั้นและซึ่งหมายความว่าฐานมาตรฐานสำหรับจะประกอบด้วยเวกเตอร์ลักษณะเฉพาะทั่วไปที่เป็นอิสระเชิงเส้นหนึ่งตัวที่มีอันดับ 2 และเวกเตอร์ลักษณะเฉพาะทั่วไปที่เป็นอิสระเชิงเส้นสองตัวที่มีอันดับ 1 หรือเทียบเท่ากับโซ่ของเวกเตอร์สองตัวและโซ่ของเวกเตอร์หนึ่งตัวเมื่อกำหนดให้เราพบว่า ρ 2 = 1 {\displaystyle \rho _{2}=1} ρ 1 = 2 {\displaystyle \rho _{1}=2} A {\displaystyle A} { x 2 , x 1 } {\displaystyle \left\{\mathbf {x} _{2},\mathbf {x} _{1}\right\}} { y 1 } {\displaystyle \left\{\mathbf {y} _{1}\right\}} M = ( y 1 x 1 x 2 ) {\displaystyle M={\begin{pmatrix}\mathbf {y} _{1}&\mathbf {x} _{1}&\mathbf {x} _{2}\end{pmatrix}}}
M = ( 2 2 0 1 3 0 0 − 4 1 ) , {\displaystyle M={\begin{pmatrix}2&2&0\\1&3&0\\0&-4&1\end{pmatrix}},} และ
J = ( 2 0 0 0 2 1 0 0 2 ) , {\displaystyle J={\begin{pmatrix}2&0&0\\0&2&1\\0&0&2\end{pmatrix}},} โดยที่เป็นเมทริกซ์โมดอลทั่วไปสำหรับคอลัมน์ของเป็นฐานมาตรฐานสำหรับและ[ 54 ] โปรดทราบว่าเนื่องจากเวกเตอร์ลักษณะเฉพาะทั่วไปเองก็ไม่เป็นเอกลักษณ์ และเนื่องจากบางคอลัมน์ของทั้งและ อาจสลับกันได้ จึงสรุป ได้ ว่าทั้งและไม่เป็นเอกลักษณ์[ 55 ] M {\displaystyle M} A {\displaystyle A} M {\displaystyle M} A {\displaystyle A} A M = M J {\displaystyle AM=MJ} M {\displaystyle M} J {\displaystyle J} M {\displaystyle M} J {\displaystyle J}
ตัวอย่างที่ 5 ในตัวอย่างที่ 3 เราพบฐานมาตรฐานของเวกเตอร์ลักษณะเฉพาะทั่วไปที่เป็นอิสระเชิงเส้นสำหรับเมทริกซ์เมทริกซ์โมดอลทั่วไปสำหรับคือ A {\displaystyle A} A {\displaystyle A}
M = ( y 1 x 1 x 2 x 3 ) = ( − 14 2 − 2 0 4 0 2 0 − 3 0 0 1 1 0 0 0 ) . {\displaystyle M={\begin{pmatrix}\mathbf {y} _{1}&\mathbf {x} _{1}&\mathbf {x} _{2}&\mathbf {x} _{3}\end{pmatrix}}={\begin{pmatrix}-14&2&-2&0\\4&0&2&0\\-3&0&0&1\\1&0&0&0\end{pmatrix}}.} เมทริกซ์ในรูปแบบปกติของจอร์แดน คล้ายกับคือ A {\displaystyle A}
J = ( 4 0 0 0 0 5 1 0 0 0 5 1 0 0 0 5 ) , {\displaystyle J={\begin{pmatrix}4&0&0&0\\0&5&1&0\\0&0&5&1\\0&0&0&5\end{pmatrix}},} ดังนั้น. A M = M J {\displaystyle AM=MJ}
แอปพลิเคชัน
ฟังก์ชันเมทริกซ์ สามการดำเนินการพื้นฐานที่สุดที่สามารถทำได้กับเมทริกซ์จัตุรัส ได้แก่ การบวกเมทริกซ์ การคูณด้วยสเกลาร์ และการคูณเมทริกซ์[ 56 ] การดำเนินการเหล่านี้จำเป็นอย่างยิ่งสำหรับการกำหนด ฟังก์ชัน พหุนาม ของเมทริกซ์n × n [ 57 ] หาก เราจำได้จากแคลคูลัส พื้นฐาน ว่าฟังก์ชันหลายฟังก์ชันสามารถเขียนได้ในรูปอนุกรมแมคลาลิน เราก็สามารถกำหนดฟังก์ชันทั่วไปของเมทริกซ์ได้ง่ายขึ้น[ 58 ] ถ้าเมทริกซ์สามารถทำให้เป็นแนวทแยงได้ นั่นคือ A {\displaystyle A} A {\displaystyle A}
D = M − 1 A M , {\displaystyle D=M^{-1}AM,} กับ
D = ( λ 1 0 ⋯ 0 0 λ 2 ⋯ 0 ⋮ ⋮ ⋱ ⋮ 0 0 ⋯ λ n ) , {\displaystyle D={\begin{pmatrix}\lambda _{1}&0&\cdots &0\\0&\lambda _{2}&\cdots &0\\\vdots &\vdots &\ddots &\vdots \\0&0&\cdots &\lambda _{n}\end{pmatrix}},} แล้ว
D k = ( λ 1 k 0 ⋯ 0 0 λ 2 k ⋯ 0 ⋮ ⋮ ⋱ ⋮ 0 0 ⋯ λ n k ) {\displaystyle D^{k}={\begin{pmatrix}\lambda _{1}^{k}&0&\cdots &0\\0&\lambda _{2}^{k}&\cdots &0\\\vdots &\vdots &\ddots &\vdots \\0&0&\cdots &\lambda _{n}^{k}\end{pmatrix}}} และการประเมินอนุกรม Maclaurin สำหรับฟังก์ชันของนั้นง่ายขึ้นมาก[ 59 ] ตัวอย่างเช่น เพื่อให้ได้กำลังk ใดๆ ของเราจำเป็นต้องคำนวณคูณด้วย แล้วคูณผลลัพธ์ด้วย เท่านั้น [ 60 ] A {\displaystyle A} A {\displaystyle A} D k {\displaystyle D^{k}} D k {\displaystyle D^{k}} M {\displaystyle M} M − 1 {\displaystyle M^{-1}}
โดยใช้เวกเตอร์ลักษณะเฉพาะทั่วไป เราสามารถหาฟอร์มปกติของจอร์แดนสำหรับและผลลัพธ์เหล่านี้สามารถขยายไปสู่วิธีการที่ตรงไปตรงมาสำหรับการคำนวณฟังก์ชันของเมทริกซ์ที่ไม่สามารถทำให้เป็นแนวทแยงได้[ 61 ] (ดูการแยกส่วนฟังก์ชันเมทริกซ์#จอร์แดน ) A {\displaystyle A}
สมการเชิงอนุพันธ์ พิจารณาปัญหาการแก้ระบบสมการเชิงอนุพันธ์สามัญเชิงเส้น
x ′ = A x , {\displaystyle \mathbf {x} '=A\mathbf {x} ,} 5
ที่ไหน
x = ( x 1 ( t ) x 2 ( t ) ⋮ x n ( t ) ) , x ′ = ( x 1 ′ ( t ) x 2 ′ ( t ) ⋮ x n ′ ( t ) ) , {\displaystyle \mathbf {x} ={\begin{pmatrix}x_{1}(t)\\x_{2}(t)\\\vdots \\x_{n}(t)\end{pmatrix}},\quad \mathbf {x} '={\begin{pmatrix}x_{1}'(t)\\x_{2}'(t)\\\vdots \\x_{n}'(t)\end{pmatrix}},} และ A = ( a i j ) . {\displaystyle A=(a_{ij}).} ถ้าเมทริกซ์เป็นเมทริกซ์แนวทแยงมุมเพื่อให้สำหรับแล้วระบบ ( 5 ) จะลดลงเหลือระบบสม การ n สมการซึ่งมีรูปแบบดังนี้ A {\displaystyle A} a i j = 0 {\displaystyle a_{ij}=0} i ≠ j {\displaystyle i\neq j}
ในกรณีนี้ คำตอบทั่วไปคือ
x 1 = k 1 e a 11 t {\displaystyle x_{1}=k_{1}e^{a_{11}t}} x 2 = k 2 e a 22 t {\displaystyle x_{2}=k_{2}e^{a_{22}t}} ⋮ {\displaystyle \vdots } x n = k n e a n n t . {\displaystyle x_{n}=k_{n}e^{a_{nn}t}.} โดยทั่วไปแล้ว เราพยายามทำให้ระบบ ( 5 ) เป็นเมทริกซ์ทแยงมุมและลดรูปให้เหลือระบบเช่น ( 6 ) ดังต่อไปนี้ ถ้าสามารถทำให้เป็นเมทริกซ์ทแยงมุมได้ เราจะได้โดยที่เป็นเมทริกซ์โมดอลสำหรับแทนค่า สมการ ( 5 ) จะอยู่ในรูปแบบหรือ A {\displaystyle A} A {\displaystyle A} D = M − 1 A M {\displaystyle D=M^{-1}AM} M {\displaystyle M} A {\displaystyle A} A = M D M − 1 {\displaystyle A=MDM^{-1}} M − 1 x ′ = D ( M − 1 x ) {\displaystyle M^{-1}\mathbf {x} '=D(M^{-1}\mathbf {x} )}
y ′ = D y , {\displaystyle \mathbf {y} '=D\mathbf {y} ,} 7
ที่ไหน
x = M y . {\displaystyle \mathbf {x} =M\mathbf {y} .} 8
คำตอบของ ( 7 ) คือ
y 1 = k 1 e λ 1 t {\displaystyle y_{1}=k_{1}e^{\lambda _{1}t}} y 2 = k 2 e λ 2 t {\displaystyle y_{2}=k_{2}e^{\lambda _{2}t}} ⋮ {\displaystyle \vdots } y n = k n e λ n t . {\displaystyle y_{n}=k_{n}e^{\lambda _{n}t}.} จากนั้นจะได้ คำตอบของ ( 5 ) โดยใช้ความสัมพันธ์ ( 8 ) [ 62 ] x {\displaystyle \mathbf {x} }
ในทางกลับกัน หากเมทริกซ์ไม่สามารถทำให้เป็นเมทริกซ์ทแยงมุมได้ เราจะเลือกให้เป็นเมทริกซ์โมดอลทั่วไปสำหรับโดยที่เป็นรูปแบบปกติของจอร์แดนของ ระบบจะมีรูปแบบดังนี้ A {\displaystyle A} M {\displaystyle M} A {\displaystyle A} J = M − 1 A M {\displaystyle J=M^{-1}AM} A {\displaystyle A} y ′ = J y {\displaystyle \mathbf {y} '=J\mathbf {y} }
y 1 ′ = λ 1 y 1 + ϵ 1 y 2 ⋮ y n − 1 ′ = λ n − 1 y n − 1 + ϵ n − 1 y n y n ′ = λ n y n , {\displaystyle {\begin{aligned}y_{1}'&=\lambda _{1}y_{1}+\epsilon _{1}y_{2}\\&\vdots \\y_{n-1}'&=\lambda _{n-1}y_{n-1}+\epsilon _{n-1}y_{n}\\y_{n}'&=\lambda _{n}y_{n},\end{aligned}}}
9
โดยที่คือค่าลักษณะเฉพาะจากแนวทแยงหลักของและคือค่าหนึ่งและศูนย์จากแนวทแยงเหนือของระบบ ( 9 ) มักจะแก้ได้ง่ายกว่า ( 5 ) เราอาจแก้สมการสุดท้ายใน ( 9 ) สำหรับโดยจะได้ จากนั้นเราแทนค่าคำตอบนี้สำหรับลงในสมการรองสุดท้ายใน ( 9 ) และแก้หา ดำเนินการตามขั้นตอนนี้ต่อไป เราทำงานผ่าน ( 9 ) จากสมการสุดท้ายไปจนถึงสมการแรก โดยแก้ระบบทั้งหมดสำหรับ จากนั้นจะได้ คำตอบโดยใช้ความสัมพันธ์ ( 8 ) [ 63 ] λ i {\displaystyle \lambda _{i}} J {\displaystyle J} ϵ i {\displaystyle \epsilon _{i}} J {\displaystyle J} y n {\displaystyle y_{n}} y n = k n e λ n t {\displaystyle y_{n}=k_{n}e^{\lambda _{n}t}} y n {\displaystyle y_{n}} y n − 1 {\displaystyle y_{n-1}} y {\displaystyle \mathbf {y} } x {\displaystyle \mathbf {x} }
บทพิสูจน์ย่อย:
กำหนดให้ลำดับของเวกเตอร์ลักษณะเฉพาะทั่วไปที่มีความยาวดังต่อไปนี้r , {\displaystyle r,}
X 1 = v 1 e λ t {\displaystyle X_{1}=v_{1}e^{\lambda t}} X 2 = ( t v 1 + v 2 ) e λ t {\displaystyle X_{2}=(tv_{1}+v_{2})e^{\lambda t}} X 3 = ( t 2 2 v 1 + t v 2 + v 3 ) e λ t {\displaystyle X_{3}=\left({\frac {t^{2}}{2}}v_{1}+tv_{2}+v_{3}\right)e^{\lambda t}} ⋮ {\displaystyle \vdots } X r = ( t r − 1 ( r − 1 ) ! v 1 + . . . + t 2 2 v r − 2 + t v r − 1 + v r ) e λ t {\displaystyle X_{r}=\left({\frac {t^{r-1}}{(r-1)!}}v_{1}+...+{\frac {t^{2}}{2}}v_{r-2}+tv_{r-1}+v_{r}\right)e^{\lambda t}} ,ฟังก์ชันเหล่านี้แก้ระบบสมการ
X ′ = A X . {\displaystyle X'=AX.} การพิสูจน์:
กำหนด
v 0 = 0 {\displaystyle v_{0}=0} X j ( t ) = e λ t ∑ i = 1 j t j − i ( j − i ) ! v i . {\displaystyle X_{j}(t)=e^{\lambda t}\sum _{i=1}^{j}{\frac {t^{j-i}}{(j-i)!}}v_{i}.} จากนั้น เช่นเดียวกับและ, t 0 = 1 {\displaystyle {t^{0}}=1} 1 ′ = 0 {\displaystyle 1'=0}
X j ′ ( t ) = e λ t ∑ i = 1 j − 1 t j − i − 1 ( j − i − 1 ) ! v i + e λ t λ ∑ i = 1 j t j − i ( j − i ) ! v i {\displaystyle X'_{j}(t)=e^{\lambda t}\sum _{i=1}^{j-1}{\frac {t^{j-i-1}}{(j-i-1)!}}v_{i}+e^{\lambda t}\lambda \sum _{i=1}^{j}{\frac {t^{j-i}}{(j-i)!}}v_{i}} .ในทางกลับกัน เราก็มีและอื่นๆ v 0 = 0 {\displaystyle v_{0}=0}
A X j ( t ) = e λ t ∑ i = 1 j t j − i ( j − i ) ! A v i {\displaystyle AX_{j}(t)=e^{\lambda t}\sum _{i=1}^{j}{\frac {t^{j-i}}{(j-i)!}}Av_{i}} = e λ t ∑ i = 1 j t j − i ( j − i ) ! ( v i − 1 + λ v i ) {\displaystyle =e^{\lambda t}\sum _{i=1}^{j}{\frac {t^{j-i}}{(j-i)!}}(v_{i-1}+\lambda v_{i})} = e λ t ∑ i = 2 j t j − i ( j − i ) ! v i − 1 + e λ t λ ∑ i = 1 j t j − i ( j − i ) ! v i {\displaystyle =e^{\lambda t}\sum _{i=2}^{j}{\frac {t^{j-i}}{(j-i)!}}v_{i-1}+e^{\lambda t}\lambda \sum _{i=1}^{j}{\frac {t^{j-i}}{(j-i)!}}v_{i}} = e λ t ∑ i = 1 j − 1 t j − i − 1 ( j − i − 1 ) ! v i + e λ t λ ∑ i = 1 j t j − i ( j − i ) ! v i {\displaystyle =e^{\lambda t}\sum _{i=1}^{j-1}{\frac {t^{j-i-1}}{(j-i-1)!}}v_{i}+e^{\lambda t}\lambda \sum _{i=1}^{j}{\frac {t^{j-i}}{(j-i)!}}v_{i}} = X j ′ ( t ) {\displaystyle =X'_{j}(t)} ตามความจำเป็น
หมายเหตุ ^ บรอนสัน (1970 , หน้า 189)↑ Beauregard & Fraleigh (1973 , หน้า 310)^ เนอริง (1970 , หน้า 118)↑ Golub & Van Loan (1996 , หน้า 316)↑ Beauregard & Fraleigh (1973 , หน้า 319)^ บรอนสัน (1970 , หน้า 194–195)↑ Golub & Van Loan (1996 , หน้า 311)^ บรอนสัน (1970 , หน้า 196)^ บรอนสัน (1970 , หน้า 189)↑ โบเรอการ์ด แอนด์ ฟราเลห์ (1973 , หน้า 316–318)^ เนอริง (1970 , หน้า 118)^ บรอนสัน (1970 , หน้า 196)↑ แอนตัน (1987 , หน้า 301–302)↑ Beauregard & Fraleigh (1973 , หน้า 266)↑ ภาระและงานแสดง (1993 , หน้า 401)↑ Golub & Van Loan (1996 , หน้า 310–311)^ ฮาร์เปอร์ (1976 , หน้า 58)^ เฮอร์สไตน์ (1964 , หน้า 225)↑ ไครส์ซิก (1972 , หน้า 273, 684)^ เนอริง (1970 , หน้า 104)↑ ภาระและงานแสดง (1993 , หน้า 401)↑ โบเรอการ์ด แอนด์ ฟราเลห์ (1973 , หน้า 270–274)^ บรอนสัน (1970 , หน้า 179–183)^ บรอนสัน (1970 , หน้า 181)^ บรอนสัน (1970 , หน้า 179)↑ โบเรอการ์ด แอนด์ ฟราเลห์ (1973 , หน้า 270–274)^ บรอนสัน (1970 , หน้า 179–183)^ บรอนสัน (1970 , หน้า 189)^ บรอนสัน (1970 , หน้า 190, 202)^ บรอนสัน (1970 , หน้า 189, 203)^ บรอนสัน (1970 , หน้า 206–207)^ บรอนสัน (1970 , หน้า 205)^ บรอนสัน (1970 , หน้า 196)^ บรอนสัน (1970 , หน้า 189, 209–215)↑ Golub & Van Loan (1996 , หน้า 316)^ เฮอร์สไตน์ (1964 , หน้า 259)^ เนอริง (1970 , หน้า 118)^ เนอริง (1970 , หน้า 118)^ เนอริง (1970 , หน้า 118)^ เฮอร์สไตน์ (1964 , หน้า 261)↑ Beauregard & Fraleigh (1973 , หน้า 310)↑ เนิร์ง (1970 , หน้า 122, 123)^ บรอนสัน (1970 , หน้า 189–209)^ บรอนสัน (1970 , หน้า 194–195)^ บรอนสัน (1970 , หน้า 196, 197)^ บรอนสัน (1970 , หน้า 197, 198)^ บรอนสัน (1970 , หน้า 190–191)^ บรอนสัน (1970 , หน้า 197–198)^ บรอนสัน (1970 , หน้า 205)↑ Beauregard & Fraleigh (1973 , หน้า 311)^ คัลเลน (1966 , หน้า 114)^ แฟรงคลิน (1968 , หน้า 122)^ บรอนสัน (1970 , หน้า 207)^ บรอนสัน (1970 , หน้า 208)^ บรอนสัน (1970 , หน้า 206)↑ โบเรอการ์ด แอนด์ ฟราเลห์ (1973 , หน้า 57–61)^ บรอนสัน (1970 , หน้า 104)^ บรอนสัน (1970 , หน้า 105)^ บรอนสัน (1970 , หน้า 184)^ บรอนสัน (1970 , หน้า 185)^ บรอนสัน (1970 , หน้า 209–218)↑ โบเรอการ์ด แอนด์ ฟราเลห์ (1973 , หน้า 274–275)↑ Beauregard & Fraleigh (1973 , หน้า 317)