הבדלים בין גרסאות בדף "משפט ז'ורדן"

מתוך Math-Wiki
קפיצה אל: ניווט, חיפוש
(אלגוריתם לז'ירדון מטריצה)
(אלגוריתם לז'ירדון מטריצה)
שורה 25: שורה 25:
  
  
*מצא את הפולינום המינימלי של המטריצה A. נסמן את הערכים העצמיים של המטריצה ב <math>\lambda_1,...,\lambda_n</math>
+
*נמצא את הפולינום המינימלי של המטריצה A. נסמן את הערכים העצמיים של המטריצה ב <math>\lambda_1,...,\lambda_n</math>
  
  
*עבור כל ע"ע <math>\lambda</math> מצא בסיס מז'רדן עבור המרחב העצמי המוכלל <math>K_\lambda</math> באופן הבא:
+
*עבור כל ע"ע <math>\lambda</math> נמצא בסיס מז'רדן עבור המרחב העצמי המוכלל <math>K_\lambda</math> באופן הבא:
  
  
שורה 34: שורה 34:
  
  
:*נביט במטריצה <math>(A-\lambda I)^{k-1}</math> ונבחר את עמודות מספר <math>i_1,...,i_p</math> המהוות בסיס למרחב העמודות <math>C([A-\lambda I]^{k-1})</math>
+
:*נמצא בסיס ל <math>V_\lambda\cap C([A-\lambda I]^{k-1})</math>
  
  
:*נפתור את מערכת המשוואות <math>x_1AC_{i_1}([A-\lambda I]^{k-1}) + ... + x_pC_{i_p}([A-\lambda I]^{k-1})</math>
+
::*נביט במטריצה <math>(A-\lambda I)^{k-1}</math> ונבחר עמודות <math>C_{i_1}([A-\lambda I]^{k-1}),...,C_{i_p}([A-\lambda I]^{k-1})</math> המהוות בסיס למרחב העמודות <math>C([A-\lambda I]^{k-1})</math>
 +
 
 +
 
 +
::*נפתור את מערכת המשוואות <math>x_1(A-\lambda I)C_{i_1}([A-\lambda I]^{k-1}) + ... + x_p(A-\lambda I)C_{i_p}([A-\lambda I]^{k-1})=0</math>  
 +
 
 +
 
 +
::*לכל וקטור <math>x=(x_1,...,x_n)</math> בבסיס למרחב הפתרונות למערכת נסמן <math>u_x=x_1e_{i_1}+...+x_pe_{i_p}</math>. '''הערה''': שימו לב כי תמיד מתקיים <math>C_i(A)=Ae_i</math> כאשר <math>e_i</math> הוקטור ה-i בבסיס הסטנדרטי.
 +
 
 +
 
 +
::*עבור כל וקטור x בבסיס למרחב הפתרונות נוסיף את כל הוקטורים במסלול <math>(A-\lambda I)^{k-1}u_x, (A-\lambda I)^{k-2}u_x,...,(A-\lambda I)u_x,u_x</math> לבסיס המז'רדן בסדר משמאל לימין.
  
  

גרסה מ־15:10, 4 בדצמבר 2012

בלוק ז'ורדן

בלוק ז'ורדן הינו מטריצה ריבועית מהצורה

J_n(\lambda):=\begin{pmatrix} \lambda & 1 & 0 & \cdots & 0 \\ 0 & \lambda & 1 & \ddots & \vdots \\ \vdots & \ddots & \ddots & \ddots &\vdots \\\vdots & \ddots & \ddots &\lambda & 1\\ 0 & \cdots & \cdots & 0 & \lambda \end{pmatrix}


לדוגמא,

J_3(0)=\begin{pmatrix}0 & 1 & 0 \\ 0 & 0 & 1 \\ 0 & 0 & 0\end{pmatrix}, J_3(2)=\begin{pmatrix}2 & 1 & 0 \\ 0 & 2 & 1 \\ 0 & 0 & 2\end{pmatrix}

נזכר בסימון של סכום ישר של מטריצות, לדוגמא: J_1(2)\oplus J_2(0)\oplus J_2(2)= \begin{pmatrix} 2 & 0 & 0 & 0 & 0 \\ 0 & 0 & 1 & 0 & 0 \\  0 & 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 2 & 1 \\ 0 & 0 & 0 & 0 & 2\end{pmatrix}

משפט ז'ורדן

תהי A מטריצה ריבועית, כך שהפולינום האופייני שלה מתפרק לגורמים לינאריים. אזי A דומה למטריצה אלכסונית בלוקים, כאשר כל בלוקיה הם בצורת ג'ורדן. בנוסף, צורה זו יחידה עד כדי סדר הבלוקים.

הוכחה ומציאת מטריצה מז'רדנת

סיכום בנושא משפט ז'ורדן על ידי דר' בועז צבאן


אלגוריתם לז'ירדון מטריצה

תהי A מטריצה כך שהפולינום האופייני שלה מתפרק לגורמים לינאריים.


  • נמצא את הפולינום המינימלי של המטריצה A. נסמן את הערכים העצמיים של המטריצה ב \lambda_1,...,\lambda_n


  • עבור כל ע"ע \lambda נמצא בסיס מז'רדן עבור המרחב העצמי המוכלל K_\lambda באופן הבא:


  • נסמן ב k את החזקה של הגורם האי פריק (x-\lambda) בפולינום המינימלי


  • נמצא בסיס ל V_\lambda\cap C([A-\lambda I]^{k-1})


  • נביט במטריצה (A-\lambda I)^{k-1} ונבחר עמודות C_{i_1}([A-\lambda I]^{k-1}),...,C_{i_p}([A-\lambda I]^{k-1}) המהוות בסיס למרחב העמודות C([A-\lambda I]^{k-1})


  • נפתור את מערכת המשוואות x_1(A-\lambda I)C_{i_1}([A-\lambda I]^{k-1}) + ... + x_p(A-\lambda I)C_{i_p}([A-\lambda I]^{k-1})=0


  • לכל וקטור x=(x_1,...,x_n) בבסיס למרחב הפתרונות למערכת נסמן u_x=x_1e_{i_1}+...+x_pe_{i_p}. הערה: שימו לב כי תמיד מתקיים C_i(A)=Ae_i כאשר e_i הוקטור ה-i בבסיס הסטנדרטי.


  • עבור כל וקטור x בבסיס למרחב הפתרונות נוסיף את כל הוקטורים במסלול (A-\lambda I)^{k-1}u_x, (A-\lambda I)^{k-2}u_x,...,(A-\lambda I)u_x,u_x לבסיס המז'רדן בסדר משמאל לימין.


  • אחד את הבסיסים המז'רדנים למרחבים המוכללים לכדי בסיס B למרחב כולו, זהו הבסיס המז'רדן של המטריצה


  • שים את איברי הבסיס B בעמודות מטריצה P. מתקיים כי J=P^{-1}AP הינה צורת הז'ורדן של המטריצה A.

דוגמאות

ז'ירדון של מטריצה ניליפוטנטית

מצאו בסיס מז'רדן למטריצה הבאה:

A=\begin{pmatrix}
0 & 1 & 0 & 1 & -1 \\
0 & 0 & 0 & 0 & 1 \\
0 & -1 & 0 & -1 & 1 \\
1 & 0 & 1 & 0 & 0 \\
0 & 0 & 0 & 0 & 0 \\

\end{pmatrix}


  • ראשית, נחשב את הפולינום האופייני p_A(x)=x^5, כלומר זוהי מטריצה ניליפוטנטית
  • שנית, נמצא את הפולינום המינימלי m_A(x)=x^3, בפרט המטריצה ניליפוטנטית מסדר 3
  • כעת נמצא בסיס ל C(A^{3-1}) מהצורה A^2v_1,A^2v_2,...,A^2v_k באופן הבא:
    • נבחר עמודות של המטריצה A^2 המהוות בסיס ל- C(A^2)
    • כל עמודה i שבחרנו ניתן להציג כ- A^2e_i


A^2=\begin{pmatrix}
1 & 0 & 1 & 0 & 1 \\
0 & 0 & 0 & 0 & 0 \\
-1 & 0 & -1 & 0 & -1 \\
0 & 0 & 0 & 0 & 0 \\
0 & 0 & 0 & 0 & 0 \\

\end{pmatrix}


לכן בסיס למרחב העמודות הינו A^2e_1

  • כעת המסלול A^2e_1,Ae_1,e_1 הוא חלק של הבסיס המז'רדן משמאל לימין. שימו לב שסדר הוקטורים בבסיס המז'רדן חשוב מאד.


  • השלב הבא הוא להשלים את הבסיס שמצאנו (A^2e_1) לבסיס למרחב N(A)\cap C(A^{3-2})=N(A)\cap C(A) מהצורה Av_1,Av_2,...,Av_p באופן הבא:
    • נבחר בסיס u_1,...,u_r למרחב העמודות C(A)
    • נפתור את המערכת A(a_1u_1+...+a_ru_r) על מנת למצוא בסיס ל N(A)\cap C(A)
    • נשמיט וקטורים על מנת שלא תהא תלות לינארית בבסיס שבחרנו עד כה


בדוגמא שלנו, העמודה הראשונה, השנייה והחמישית מהוות בסיס למרחב העמודות של A:

u_1= (0,0,0,1,0)
u_2= (1,0,-1,0,0)
u_3= (-1,1,1,0,0)


כעת נפתור את המערכת a_1Au_1+a_2Au_2+a_3Au_3=0, זהו בדיוק מרחב האפס של המטריצה שעמודותיה הן Au_1,Au_2,Au_3:


N  \begin{pmatrix}
1 & 0 & 1  \\
0 & 0 & 0  \\
-1 & 0 & -1  \\
0 & 0 & 0  \\
0 & 0 & 0  \\

\end{pmatrix}  = span\{(0,1,0),(-1,0,1)\}


כיוון שאלו המקדמים a_1,a_2,a_3 אנו מקבלים את בסיס ל N(A)\cap C(A):

\{u_2,u_3-u_1\}=\{Ae_2,A(e_5-e_1)\}


הערה: שימו לב שu3=Ae_5 כיוון שזו העמודה החמישית


כיוון ש Ae_2=A^2e_1 אנו משמטים איבר זה ונשארים עם A(e_5-e_1)


  • המסלול A(e_5-e_1),e_5-e_1 משלים לנו את הבסיס המז'רדן.


סיכום

הבסיס המז'רדן הינו

A^2e_1,Ae_1,e_1,A(e_5-e_1),e_5-e_1

נסמן בP את המטריצה שעמודותיה הן איברי הבסיס

P=\begin{pmatrix}
1 & 0 & 1 & -1 & -1 \\
0 & 0 & 0 & 1 & 0 \\
-1 & 0 & 0 & 1 & 0 \\
0 & 1 & 0 & -1 & 0 \\
0 & 0 & 0 & 0 & 1 \\

\end{pmatrix}


אזי (לא במפתיע) מתקיימת המשוואה הבאה:


P^{-1}AP=\begin{pmatrix}
0 & 1 & 0 & 0 & 0 \\
0 & 0 & 1 & 0 & 0 \\
0 & 0 & 0 & 0 & 0 \\
0 & 0 & 0 & 0 & 1 \\
0 & 0 & 0 & 0 & 0 \\

\end{pmatrix}

כלומר זו צורת הז'ורדן של המטריצה A.

ז'ירדון של מטריצה עם ע"ע יחיד

מצאו בסיס מז'רדן למטריצה הבאה:

A=\begin{pmatrix}
2 & -1 & 0 & 0 & 0 & 0 \\
0 & 2 & 0 & 0 & 0 & 0 \\
0 & 0 & 2 & 0 & 1 & 0 \\
0 & -1 & 0 & 2 & -1 & 0 \\
0 & 0 & 0 & 0 & 2 & 0 \\
\frac{1}{2} & 1 & -\frac{1}{2} & -\frac{1}{2} & 2.5 & 2 \\

\end{pmatrix}
  • ראשית נמצא את הפולינום האופייני p_A(x)=(x-2)^6, כלומר 2 הינו הערך העצמי היחיד
  • לפי משפט קיילי המילטון (A-2I)^6=0 ולכן A-2I ניליפוטנטית.
  • נמצא לה צורת ז'ורדן J=P^{-1}(A-2I)P = P^{-1}AP - P^{-1}2IP = P^{-1}AP-2I
  • לכן צורת הז'ורדן של המטריצה A הינה J+2I, כאשר הבסיס המז'רדן הוא אותו בסיס המז'רדן את A-2I.


A-2I=\begin{pmatrix}
0 & -1 & 0 & 0 & 0 & 0 \\
0 & 0 & 0 & 0 & 0 & 0 \\
0 & 0 & 0 & 0 & 1 & 0 \\
0 & -1 & 0 & 0 & -1 & 0 \\
0 & 0 & 0 & 0 & 0 & 0 \\
\frac{1}{2} & 1 & -\frac{1}{2} & -\frac{1}{2} & 2.5 & 0 \\

\end{pmatrix}


  • כעת (A-2I)^2=0, לכן נמצא בסיס לC(A-2I)
  • העמודה הראשונה, השנייה והחמישית פורסות את מרחב העמודות של המטריצה ולכן הבסיס הינו (A-2I)e_1,(A-2I)e_2,(A-2I)e_5
  • בסיס זה מייצר שלושה מסלולים מאורך שתים, ולכן מצאנו מיד בסיס מז'רדן:
(A-2I)e_1,e_1,(A-2I)e_2,e_2,(A-2I)e_5,e_5


P=\begin{pmatrix}
0 & 1 & -1 & 0 & 0 & 0 \\
0 & 0 & 0 & 1 & 0 & 0 \\
0 & 0 & 0 & 0 & 1 & 0 \\
0 & 0 & -1 & 0 & -1 & 0 \\
0 & 0 & 0 & 0 & 0 & 1 \\
\frac{1}{2} & 0 & 1 & 0 &  2.5 & 0 \\

\end{pmatrix}


ושוב, הפלא ופלא, מתקיים:


P^{-1}AP=\begin{pmatrix}
2 & 1 & 0 & 0 & 0 & 0 \\
0 & 2 & 0 & 0 & 0 & 0 \\
0 & 0 & 2 & 1 & 0 & 0 \\
0 & 0 & 0 & 2 & 0 & 0 \\
0 & 0 & 0 & 0 & 2 & 1 \\
0 & 0 & 0 & 0 & 0 & 2 \\

\end{pmatrix}