I Algebra liniowa z geometrią 1

Rozdział 11 Do czego zmierzamy?

SŁOWA KLUCZOWE ROZDZIAŁU wektor własny ◇ wartość własna ◇ diagonalizacja ◇ baza wektorów własnych ◇ diagonalizacja ortogonalna ◇ macierz stochastyczna ◇ dominująca wartość własna ◇ wielomian charakterystyczny ◇ zespolona postać Jordana ◇ rzeczywista postać Jordana ◇ krotność geometryczna i algebraiczna ◇ twierdzenie Cayleya–Hamiltona ◇ diagonalizacja ortogonalna macierzy symetrycznej

Przedstawimy teraz motywację dla naszych dalszych rozważań. Jak już pewnie zauważyliście szczególnie łatwe do analizy są macierze diagonalne

D=diag⁢(λ1,…,λn),λ1,…,λn∈ℝ.

Przykładowo,

Dk=diag⁢(λ1k,…,λnk),det⁢D=λ1⋅…⋅λn.

Jeśli D jest nieosobliwa, to D-1=diag⁢(1/λ1,…,1/λn).

Załóżmy, że macierz A∈Mn×n⁢(ℝ) jest podobna do macierzy diagonalnej D, czyli istnieje taka macierz nieosobliwa

S=[v1…vn]∈Mn×n⁢(ℝ),

że

S-1⁢A⁢S=D.

Wtedy A=S⁢D⁢S-1, więc

Ak=S⁢Dk⁢S-1,det⁢A=det⁢D,

czyli A jest równie prosta do analizy jak macierz diagonalna D.

PROBLEM 1 Scharakteryzować macierze A∈Mn×n⁢(ℝ), które są diagonalizowalne tzn. istnieje taka macierz nieosobliwa S=[v1…vn]∈Mn×n⁢(ℝ), że S-1⁢A⁢S=D=diag⁢(λ1,…,λn)=[λ1⋅e1…λn⋅en].

Załóżmy, że S-1⁢A⁢S=D. Wtedy A⁢S=S⁢D, czyli

[A⁢v1…A⁢vn] =A⁢S=S⁢D
=S⁢[λ1⋅e1…λn⋅en]
=[λ1⋅S⁢e1…λn⋅S⁢en]
=[λ1⋅v1…λn⋅vn].
Stąd
A⁢v1=λ1⋅v1,…,A⁢vn=λn⋅vn.

Z powyższych rozważań wynika następujące:

Twierdzenie 11.1 (Diagonalizacja).

Dla macierzy A∈Mn×n⁢(R) poniższe warunki są równoważne

  • (i)

    istnieje taka macierz nieosobliwa S∈Mn×n⁢(ℝ), że S-1⁢A⁢S jest diagonalna,

  • (ii)

    istnieje baza v1,…,vn dla ℝn oraz takie skalary λ1,…,λn∈ℝ, że

    A⁢vi=λi⋅vi,i=1,…,n.

Warto jakoś nazwać własność jaką spełniają wektory vi w punkcie (ii). Będą one pełniły podstawową rolę w naszych dalszych rozważaniach.

Definicja 11.1.

Wektor własny i wartość własna

Powiemy, że niezerowy wektor 0≠v∈ℝn jest wektorem własnym macierzy rzeczywistej A∈Mn×n⁢(ℝ), jeśli istnieje taka liczba rzeczywista λ∈ℝ, że

A⁢v=λ⋅v.

Liczbę λ∈ℝ nazywamy wartością własną macierzy A∈Mn×n⁢(ℝ) dla wektora własnego v∈ℝn.

Twierdzenie 11.1 możemy sformułować następująco:

Twierdzenie 11.2 (Diagonalizacja=baza wektorów własnych).

Dla macierzy A∈Mn×n⁢(R) poniższe warunki są równoważne

  • (i)

    istnieje taka macierz nieosobliwa S∈Mn×n⁢(ℝ), że S-1⁢A⁢S jest diagonalna,

  • (ii)

    istnieje baza v1,…,vn dla ℝn złożona z wektorów własnych A.

Podamy geometryczną interpretację wektora własnego.

Wniosek 11.1.

Dla niezerowego wektora v∈Rn i macierzy A∈Mn×n⁢(R) poniższe warunki są równoważne

  • (1)

    v jest wektorem własnym A,

  • (2)

    prosta generowana przez v jest niezmiennicza dla A tzn.

    A⁢(span⁢{v})⊂span⁢{v}.
  • U

    Nie każda macierz rzeczywista ma jakiś wektor własny. Wystarczy znaleźć macierz rzeczywistą A, która nie ma prostych niezmienniczych. Przykładowo, macierz obrotu o kąt π2, czyli

    A=[0-110],

    nie ma żadnej prostej niezmienniczej.

    Możemy na to spojrzeć jeszcze trochę inaczej. Wektor 0≠v∈ℝ2 jest wektorem własnym A, gdy dla pewnej liczby rzeczywistej λ∈ℝ zachodzi warunek

    (A-λ⁢I)⁢v=0,

    czyli ker⁢(A-λ⁢I)≠0, więc det⁢(A-λ⁢I)=0 dla pewnej liczby rzeczywistej λ∈ℝ. Zauważmy, że w naszym przykładzie

    det⁢(A-λ⁢I)=λ2+1≠0,

    dla λ∈ℝ.

    Macierz rzeczywistą A możemy traktować w naturalny sposób jako macierz zespoloną. Możemy, więc spytać o istnienie zespolonych wektorów własnych v∈ℂ2. Zauważmy, że det⁢(A-λ⁢I)=0 dla λ=±i. Rozważmy macierz

    A-i⁢I=[-i-11-i]

    i znajdźmy taki niezerowy wektor v=[a+i⁢b,c+i⁢d]∈ℂ2, że (A-i⁢I)⁢v=0. Mamy

    (A-i⁢I)⁢v=[-i⁢a+b-c-i⁢d,a+i⁢b-i⁢c+d]T=[0,0]T,

    czyli

    a+d=0,b-c=0.

    Wektor v ma więc postać

    v=[a+i⁢b,b-i⁢a]T=[a+i⁢b,-i⁢(a+i⁢b)]T=(a+i⁢b)⁢[1,-i]T.

    Wynika stąd, że

    A⁢[1,-i]T=i⁢[1,-i]T,

    czyli [1,-i]T jest zespolonym wektorem własnym A z wartością własną i. Analogicznie sprawdzamy, że [1,i]T jest zespolonym wektorem własnym A dla wartości własnej -i. Macierz A jest diagonalizowalna jako macierz zespolona. Dla

    S=[11i-i]

    mamy

    S-1⁢A⁢S=[i00-i].
  • U

    Pojęcie wektora własnego możemy też zdefiniować w bardziej ogólnym kontekście. Niech L:V⟶V będzie odwzorowaniem liniowym. Powiemy, że niezerowy wektor v∈V jest wektorem własnym L jeśli

    L⁢(v)=λ⋅v.

    Przestrzeń V nie musi być skończenie wymiarowa. Przykładowo, niech C∞⁢(ℝ,ℝ) będzie przestrzenią wszystkich funkcji f:ℝ⟶ℝ mających ciągłą pochodną dowolnego rzędu. Wtedy odwzorowanie (operator różniczkowania)

    L:C∞⁢(ℝ,ℝ)∋f⟼f′∈C∞⁢(ℝ,ℝ)

    jest liniowe. Ponadto, dla f=eλ⁢x∈C∞⁢(ℝ,ℝ) mamy

    L⁢(f) =(eλ⁢x)′
    =λ⁢eλ⁢x
    =λ⋅f,

    czyli f=eλ⁢x jest wektorem własnym L dla wartości własnej λ.

    Uzasadnimy teraz, że jeśli L⁢(f)=λ⋅f, to f=C⁢eλ⁢x dla pewnego C∈ℝ. Rzeczywiście, równość f′=λ⁢f oznacza, że f′-λ⁢f=0, czyli również

    f′⁢e-λ⁢x-λ⁢e-λ⁢x⁢f⏟=(f⁢e-λ⁢x)′=0,

    czyli f⁢e-λ⁢x jest funkcją stałą, więc f⁢e-λ⁢x=C.

    Wynika stąd, że jądro ker⁢(L-λ⁢I) ma wymiar 1 oraz

    ker⁢(L-λ⁢I)=span⁢{eλ⁢x}.
  • U

    Niech L:V⟶V będzie odwzorowaniem liniowym i niech v1,…,vn∈V będzie bazą V. Rozważmy macierz A∈Mn×n⁢(𝔽) odwzorowania L w tej bazie. Wektor v=x1⋅v1+…+xn⋅vn jest wektorem własnym L dla wartości własnej λ∈𝔽 wtedy i tylko wtedy, gdy wektor współrzędnych [x1,…,xn]T∈𝔽n jest wektorem własnym macierzy A dla wartości własnej λ.

Możemy pójść jeszcze krok dalej. Wektory własne macierzy diagonalnej D to wektory bazy standardowej e1,…,en. Ma ona własność

(ei|ej)=δi⁢j.

Takie bazy ℝn będziemy nazywać ortonormalnymi. Wektory ei oraz ej są ortogonalne i mają normę 1.

PROBLEM 2 Scharakteryzować macierze A∈Mn×n⁢(ℝ), które są ortogonalnie diagonalizowalne tzn. istnieje taka macierz nieosobliwa S=[v1…vn]∈Mn×n⁢(ℝ), że S-1⁢A⁢S=D=diag⁢(λ1,…,λn)=[λ1⋅e1…λn⋅en] oraz (vi|vj)=δi⁢j.

Zaczniemy od przyglądnięcia się takim macierzom

S=[v1…vn]∈Mn×n⁢(ℝ),

że (vi|vj)=δi⁢j. Bezpośrednio z definicji iloczynu macierzy i definicji macierzy transponowanej wynika, że wtedy

ST⁢S=I,czyli S-1=ST.

Takie macierze będziemy nazywać ortogonalnymi. W Problemie 2 pytamy więc o istnienie takiej macierzy ortogonalnej S∈Mn×n⁢(ℝ), że

S-1⁢A⁢S=ST⁢A⁢S=D

jest diagonalna. Zauważmy, że wtedy A=S⁢D⁢ST oraz

AT =(S⁢D⁢ST)T
=(ST)T⁢DT⁢ST
=S⁢D⁢ST
=A,

czyli A musi być symetryczna. Jest to warunek konieczny, aby zaszła sytuacja opisana w Problemie 2. W następnych rozdziałach pokażemy, że macierze symetryczne to jedyne macierze rzeczywiste rozwiązujące Problem 2.

11.1 Rozkład Jordana w wymiarze 2

Spróbujmy nabrać intuicji w przypadku wymiaru 2. Zakładamy, że ciało F jest równe ℝ lub ℂ. Rozważmy macierz A∈M2×2⁢(𝔽).

Definicja 11.2.

Wektor własny i wartość własna

Liczbę λ∈𝔽 nazywamy wartością własną macierzy A∈M2×2⁢(𝔽), jeśli istnieje taki niezerowy wektor v∈𝔽2, że

A⁢v=λ⋅v.

Każdy taki wektor v nazywamy wektorem własnym A odpowiadającym wartości własnej λ.

  • U

    Niech 𝔽=ℝ. Interpretacja geometryczna jest następująca. Jeśli v∈ℝ2 jest wektorem własnym dla wartości własnej λ∈ℝ, to prosta l=span⁢{v} jest niezmiennicza dla A tzn. jeśli w∈l, to A⁢w∈l.

Jeśli v∈𝔽2 jest wektorem własnym dla wartości własnej λ∈𝔽, to (A-λ⁢I)⁢v=0. Oznacza to, że ker⁢(A-λ⁢I)≠{0}, czyli det⁢(A-λ⁢I)=0. Odwrotnie, jeśli det⁢(A-λ⁢I)=0 dla pewnego skalara λ∈𝔽, to macierz A-λ⁢I jest osobliwa. W szczególności, ker⁢(A-λ⁢I)≠{0}, czyli istnieje wektor własny v∈ker⁢(A-λ⁢I).

Wniosek 11.2.

Dla A∈M2×2⁢(F) i λ∈F następujące warunki są równoważne

  • (1)

    λ jest wartością własną A;

  • (2)

    ker⁢(A-λ⁢I)≠{0};

  • (3)

    A-λ⁢I jest osobliwa;

  • (4)

    det⁢(A-λ⁢I)=0.

Warunek (4) oznacza, że wartość własna λ∈𝔽 jest pierwiastkiem wielomianu

pA⁢(x)=det⁢(A-x⁢I)

Wielomian pA nazywamy wielomianem charakterystycznym macierzy A. Bezpośredni rachunek pokazuje, że

pA⁢(x) =det⁢[a11-xa12a21a22-x]
=x2-(a11+a22)⁢x+(a11⁢a22-a12⁢a22)
=x2-tr⁢(A)⁢x+det⁢(A).
  • U

    Macierze podobne mają takie same wielomiany charakterystyczne. Rzeczywiście, jeśli B=S-1⁢A⁢S, to

    det⁢(B-x⁢I) =det⁢(S-1⁢A⁢S-x⁢I)
    =det⁢(S-1⁢(A-x⁢I)⁢S)
    =det⁢(A-x⁢I).

    Jeśli macierze A i B mają takie same wielomiany charakterystyczne, to nie muszą być podobne. Przykładowo, A=0 oraz B=[0100] mają taki sam wielomian charakterystyczny p⁢(x)=x2, ale nie są podobne, bo macierz zerowa jest podobna tylko do siebie samej.

Dla 𝔽=ℂ zasadnicze twierdzenie algebry gwarantuje, że wielomian pA ma dwa (niekoniecznie różne) pierwiastki zespolone λ1, λ2∈ℂ. Są one wtedy wartościami własnymi macierzy A∈M2×2⁢(ℂ). Macierz zespolona ma więc zawsze dwie wartości własne. Nie oznacza to, że ma ona dwa liniowo niezależne wektory własne. Przykładowo, wartościami własnymi macierzy

A=[λ10λ],λ∈ℂ,

są λ1=λ2=λ. Mamy więc dwukrotną wartość własną λ. Zobaczmy jak wyglądają wektory własne. Niezerowy wektor v=[v1,v2]T∈ℂ2 jest wektorem własnym dla λ jeśli (A-λ⁢I)⁢v=0, czyli

[00]=[0100]⁢[v1v2]=[v20].

Stąd v2=0, czyli wektor własny v ma postać

[v10]=v1⋅[10]=v1⋅e1.

W przypadku 𝔽=ℝ, z algebraicznego punktu widzenia dla wielomianu charakterystycznego mamy trzy możliwości:

  • (i)

    pA ma dwa różne pierwiastki rzeczywiste;

  • (ii)

    pA ma dwukrotny pierwiastek rzeczywisty;

  • (iii)

    pA ma dwa sprzężone (nierzeczywiste) pierwiastki zespolone.

W przypadku (iii) macierz rzeczywista A nie ma rzeczywistych wartości własnych.

Wniosek 11.3.

Niech A∈M2×2⁢(C). Jeśli λ1,λ2∈C są pierwiastkami równania charakterystycznego

x2-tr⁢(A)⁢x+det⁢(A)=0,

to zachodzą wzory Viete’a

tr⁢(A)=λ1+λ2,det⁢(A)=λ1⋅λ2.

Dowód.

Z założenia mamy

(x-λ1)⁢(x-λ2)=pA⁢(x)=x2-tr⁢(A)⁢x+det⁢(A),

wystarczy więc porównać współczynniki wielomianów po obu stronach równości. ∎

  • U

    Dla dowolnych macierzy X,Y∈M2×2⁢(ℂ) macierze X⁢Y i Y⁢X mają takie same wielomiany charakterystyczne. Wynika, to z faktu, że

    tr⁢X⁢Y=tr⁢Y⁢X,det⁢X⁢Y=det⁢Y⁢X.
Lemat 11.1.

Jeśli A∈M2×2⁢(F) ma dwie różne rzeczywiste wartości własne λ1≠λ2 oraz vi∈F2 są wektorami własnymi odpowiadającymi λi∈F, to v1, v2 są liniowo niezależne. W szczególności, v1,v2 tworzą bazę F2.

Dowód.

W przeciwnym razie v2=p⋅v1 dla pewnego skalara p∈𝔽. Wtedy

λ2⋅v2 =A⁢(v2)
=A⁢(p⋅v1)
=p⋅A⁢(v1)
=p⁢λ1⋅v1
=λ1⋅v2,

czyli λ1=λ2, co prowadzi do sprzeczności. Skorzystaliśmy tu z faktu, że v2≠0. ∎

Twierdzenie 11.3.

Jeśli A∈M2×2⁢(F) ma dwa liniowo niezależne wektory własne v1,v2∈F2 odpowiadające wartościom własnym λ1, λ2∈F (niekoniecznie różnym), to dla macierzy B=[v1v2] mamy

Λ:=B-1⁢A⁢B=[λ100λ2]

Dowód.

Z definicji B mamy, że B⁢(ei)=vi dla i=1,2. Oczywiście B jest odwracalna oraz dla i=1,2 mamy

B-1⁢A⁢B⁢ei =B-1⁢A⁢vi
=B-1⁢(λi⋅vi)
=λi⋅ei.

∎

  • U

    Twierdzenie 11.3 orzeka, że jeśli 𝔽2 ma bazę złożoną z wektorów własnych, to A∈M2×2⁢(𝔽) ma w tej bazie postać diagonalną Λ. Macierz Λ nazywamy postacią Jordana macierzy A.

Rysunek 11.1: Działanie macierzy A na wektorach własnych v1, v2 o wartościach własnych λ1=4, λ2=-3.
Przykład 11.1.

Znajdziemy postać Jordana macierzy

A=[323-2].

Szukamy najpierw wartości własnych czyli pierwiastków równania charakterystycznego

pA⁢(x)=x2-x-12=0.

Są nimi liczby λ1=4 oraz λ2=-3. Aby znaleźć wektor własny dla λ1=4 znajdujemy jądro macierzy

A-4⁢I=[-123-6].

Rozwiązując równanie (A-4⁢I)⁢x=0 otrzymujemy, że x=[2⁢x2,x2]T. Jednym z wektorów własnych jest v1=[2,1]T. Analogicznie, v2=[-1,3]T jest wektorem własnym dla λ2. W bazie v1, v2 macierz A ma postać diagonalną

Λ=[400-3].
  • !

    Macierz zespolona A∈M2×2⁢(ℂ) może mieć rzeczywistą wartość własną, ale może nie mieć odpowiadającego jej rzeczywistego wektora własnego w ℝ2. Dla przykładu rozważmy macierz

    A=[0-ii0].

    Wielomian charakterystyczny ma postać

    pA⁢(x)=x2-1,

    więc A ma wartości własne λ1=1 oraz λ2=-1. Wektor własny [x,y]T∈ℂ2 dla λ1=1 jest niezerowym rozwiązaniem równania

    [-1-ii-1]⁢[xy]=[00],

    czyli

    -x-i⁢y=0,i⁢x-y=0.

    Stąd, y=i⁢x oraz każdy wektor własny jest postaci

    c⋅[1,i]T,c∈ℂ∖{0}.

    Nie może on być wektorem z ℝ2.

Zajmiemy się teraz przypadkiem, gdy λ∈𝔽 jest dwukrotnym pierwiastkiem wielomianu charakterystycznego pA. Mogą zachodzić dwie możliwości:

  • •

    dim⁢ker⁢(A-λ⁢I)=2. Mówimy wtedy, że krotność geometryczna wartości własnej λ jest równa jej krotności algebraicznej. Wówczas A ma dwa liniowo niezależne wektory własne odpowiadające λ i jesteśmy w sytuacji opisanej Twierdzeniem 11.3. Wtedy postać Jordana Λ macierzy A jest diagonalna:

    Λ:=B-1⁢A⁢B=[λ00λ]=λ⁢I.

    Z powyższej równości wynika, że wtedy A=λ⁢I.

  • •

    dim⁢ker⁢(A-λ⁢I)=1. Istnieje wtedy tylko jeden liniowo niezależny wektor własny v∈𝔽2 dla λ∈𝔽.

Wniosek 11.4.

Jeśli λ∈F jest dwukrotnym pierwiastkiem charakterystycznym macierzy A∈M2×2⁢(F) i A≠λ⁢I, to dim⁢ker⁢(A-λ⁢I)=1 i A ma tylko jeden liniowo niezależny wektor własny.

Lemat 11.2 (Twierdzenie Cayleya–Hamiltona).

Dla macierzy A∈M2×2⁢(F) zachodzi równość

pA⁢(A)=A2-tr⁢(A)⁢A+det⁢(A)⁢I=[0000].
Dowód.

Zastosujemy brutalną siłę, czyli bezpośredni rachunek. Niech A=[abcd]. Ponieważ pA⁢(x)=(x-a)⁢(x-d)-b⁢c, więc

pA⁢(A) =(A-a⁢I)⁢(A-d⁢I)-b⁢c⁢I
=[0bcd-a]⁢[a-dbc0]-[b⁢c00b⁢c]
=[b⁢c0c⁢(a-d)+(d-a)⁢cb⁢c]-[b⁢c00b⁢c]
=[0000].

∎

Przykład 11.2.

Niech A∈M2×2⁢(ℂ). Twierdzenie Cayleya–Hamiltona pozwala łatwo wyznaczyć An dla n∈ℕ. Niech λ1, λ2∈ℂ będą wartościami własnymi A. Wtedy

A2-(λ1+λ2)⁢A+λ1⁢λ2⁢I=0,

czyli

(A-λ1⁢I)⁢(A-λ2⁢I)=0.

Jeśli λ1≠λ2, to rozważamy macierze

X=A-λ2⁢Iλ1-λ2,Y=A-λ1⁢Iλ2-λ1.

Wtedy

X2=X,X⁢Y=Y⁢X=0,Y2=Y

Środkowa równość wynika bezpośrednio z Twierdzenia Cayleya–Hamiltona. Sprawdzimy równość X2=X. Mamy

X2 =(A-λ2⁢I)2(λ1-λ2)2
=(A-λ2⁢I)⁢(A-λ1⁢I+(λ1-λ2)⁢I)(λ1-λ2)2
=(A-λ2⁢I)⁢(A-λ1⁢I)(λ1-λ2)2+(A-λ2)⁢(λ1-λ2)⁢I(λ1-λ2)2
=A-λ2⁢I(λ1-λ2)
=X.

Stąd dla k≥2 mamy

Xk=X,Yk=Y,λ1≠λ2.

Dla λ1=λ2 rozważamy

Z=A-λ1⁢I.

Wtedy Zk=0 dla k≥2. Mamy więc

An={(λ1⁢X+λ2⁢Y)n=λ1n⁢X+λ2n⁢Y,gdy λ1≠λ2(λ1⁢I+Z)n=λ1n⁢I+n⁢λ1n-1⁢Z,gdy λ1=λ2
Wniosek 11.5.

Jeśli A∈M2×2⁢(C) jest nieosobliwa, to

A-1=1det⁢A⁢(tr⁢A⋅I-A),tr⁢A-1=tr⁢Adet⁢A.

Ponadto,

det⁢A=12⁢((tr⁢A)2-tr⁢A2).
Dowód.

Z twierdzenia Cayleya–Hamiltona mamy

A2-tr⁢(A)⁢A+det⁢(A)⁢I=[0000].

Mnożąc powyższą równość przez A-1 otrzymujemy, że

A-1=1det⁢A⁢(tr⁢A⋅I-A).

Biorąc ślady obu stron w tej równości dostajemy, że tr⁢A-1=tr⁢Adet⁢A. Wzór na det⁢A otrzymujemy obliczając ślady macierzy po obu stronach w równości Cayleya–Hamiltona. ∎

Wniosek 11.6.

Niech A∈M2×2⁢(C). Jeśli A∉{α⁢I:α∈C} oraz A2-a⁢A+b⁢I= dla pewnych a,b∈C, to

a=tr⁢A,b=det⁢A.
Dowód.

Z założenia i twierdzenia Cayleya–Hamiltona wynika, że

(a-tr⁢A)⁢A=(b-det⁢A)⁢I.

Zauważmy, że a-tr⁢A=0, bo inaczej A=b-det⁢Aa-tr⁢A⁢I, co jest sprzeczne z założeniem. Stąd również b-det⁢A=0. ∎

  • U

    Macierz A∈M2×2⁢(ℂ) jest nilpotentna, jeśli An=0 dla pewnego n∈ℕ. Pokażemy, że wtedy 0 jest jedyną wartością własną A oraz A2=0. Zauważmy, że jeśli λ∈ℂ jest wartością własną dla wektora własnego v≠0, to 0=An⁢v=λn⋅v, czyli λ=0. W szczególności, tr⁢A=0 i det⁢A=0, czyli pA⁢(x)=x2. Z twierdzenia Cayleya–Hamiltona wynika, że A2=0.

Lemat 11.3.

Jeśli λ∈F jest dwukrotnym pierwiastkiem wielomianu charakterystycznego pA, to

(A-λ⁢I)2=[0000].
Dowód.

Ze wzorów Viete’a mamy tr⁢(A)=2⁢λ i det⁢(A)=λ2. Z Lematu Cayleya–Hamiltona wynika, że

A2-2⁢λ⁢A+λ2⁢I=[0000].

Ponieważ macierze A i λ⁢I komutują ze sobą, więc

(A-λ⁢I)2=A2-2⁢λ⁢A+λ2⁢I.

∎

Twierdzenie 11.4.

Jeśli λ∈F jest dwukrotną wartością własną macierzy A∈M2×2⁢(F) i A≠λ⁢I, to istnieje taka baza v,w dla F2, że

A⁢v=λ⋅v,A⁢w=v+λ⋅w

W szczególności, dla macierzy S=[vw] mamy następującą postać Jordana dla A:

Λ=S-1⁢A⁢S=[λ10λ]

Dowód.

Niech 0≠v∈ker⁢(A-λ⁢I) będzie dowolnym wektorem własnym dla λ. Każdy wektor własny A jest postaci p⋅v, gdzie v jest pewnym wektorem własnym A dla λ oraz 0≠p∈𝔽. Niech 0≠u∈𝔽2 będzie dowolnym wektorem liniowo niezależnym z v. W szczególności, u nie jest wektorem własnym odpowiadającym λ, więc (A-λ⁢I)⁢u≠0. Ponieważ,

0=(A-λ⁢I)2⁢u=(A-λ⁢I)⁢(A-λ⁢I)⁢u,

więc (A-λ⁢I)⁢u∈ker⁢(A-λ⁢I) jest wektorem własnym odpowiadającym λ. Stąd

(A-λ⁢I)⁢u=p⋅v,

dla pewnego 0≠p∈ℝ. Dla wektora w=1p⋅u otrzymujemy, że (A-λ⁢I)⁢w=v, czyli A⁢w=v+λ⋅w. ∎

POSTAĆ JORDANA MACIERZY ZESPOLONEJ Macierz zespolona A∈M2×2⁢(ℂ) ma jedną z dwóch postaci Jordana [λ100λ2],[λ10λ].
Wniosek 11.7.
  • (i)

    Jeśli 0≠A∈M2×2⁢(ℂ) jest macierzą nilpotentną, to

    A=S⁢[0100]⁢S-1

    dla pewnej macierzy nieosobliwej S.

  • (ii)

    Jeśli A∈M2×2⁢(ℂ) jest taka, że A2=A, to A=0 lub A=I lub

    A=S⁢[1000]⁢S-1

    dla pewnej macierzy nieosobliwej S.

  • (iii)

    Jeśli A∈M2×2⁢(ℂ) jest taka, że A2=I, to A=±I lub

    A=S⁢[100-1]⁢S-1

    dla pewnej macierzy nieosobliwej S.

  • (iv)

    Jeśli A∈M2×2⁢(ℂ) jest taka, że A2=-I, to A=±i⋅I lub

    A=S⁢[0-110]⁢S-1

    dla pewnej macierzy nieosobliwej S.

Wniosek 11.8.

Dla dowolnej macierzy A∈M2×2⁢(C) mamy

  • (i)

    A jest podobna do AT.

  • (ii)

    A jest podobna do macierzy symetrycznej.

  • (iii)

    A jest iloczynem macierzy symetrycznych.

Dowód.

Punkt (i). Możemy założyć, że A jest w postaci Jordana. Jeśli jest ona diagonalna, to nie ma czego dowodzić. Załóżmy więc, że

A=[λ10λ].

Wtedy dla S=[0110] mamy S-1⁢A⁢S=AT.

Punkt (ii). Ponownie możemy założyć, że A=[λ10λ]. Szukamy takiej macierzy nieosobliwej S=[abcd]∈M2×2⁢(ℂ), że S-1⁢A⁢S jest symetryczna. Dla Δ=det⁢S≠0 mamy

S-1⁢A⁢S=[c⁢d+λ⁢Δd2-c2-c⁢d+λ⁢Δ].

Ponieważ szukamy macierzy symetrycznej, więc musi zachodzić warunek c2+d2=0. Wystarczy więc przyjąć takie a,b,c,d, że a⁢d-b⁢c≠0 i c2+d2=0. Zwróćmy uwagę, że S nie może być rzeczywista. Przykładowo, dla d=i=a, c=1 i b=0 mamy

S=[i01i],S-1⁢A⁢S=[λ-i11λ+i].

Punkt (iii). Zauważmy najpierw, że możemy założyć, że A jest w postaci Jordana Λ, bo jeśli Λ=B⁢C jest iloczynem macierzy symetrycznych, to

A=S⁢Λ⁢S-1=(S⁢B⁢ST)⁢((S-1)T⁢C⁢S-1)

jest również iloczynem macierzy symetrycznych. Dla Λ=[λ100λ2] wystarczy przyjąć

B=[λ1001],C=[100λ2].

Jeśli Λ=[λ10λ], to przyjmujemy

B=[1110],C=[0λλ1-λ].

∎

  • U

    Rozważmy macierz A w postaci bloku Jordana

    A=[λ10λ].

    Dla macierzy

    S=[0110]

    mamy S-1=S oraz

    S-1⁢A⁢S =[0110]⁢[λ10λ]⁢[0110]
    =[λ01λ]
    =AT.

Naturalne jest pytanie o postać macierzy A∈M2×2⁢(ℝ), gdy jej wielomian charakterystyczny pA ma dwie zespolone wartości własne λ1=α+i⁢β, λ2=α-i⁢β oraz β≠0. Macierz A∈M2×2⁢(ℝ) możemy traktować jako element przestrzeni M2×2⁢(ℂ). Istnieje więc wektor własny v∈ℂ2 odpowiadający wartości własnej λ1=α+i⁢β. Wektor v∈ℂ2 wygodnie jest zapisać w postaci

v=[a1+i⁢a2b1+i⁢b2]=[a1b1]+i⁢[a2b2]:=v1+i⁢v2,v1,v2∈ℝ2.

Wówczas

A⁢v=A⁢v1+i⁢A⁢v2.

Ponadto,

A⁢(v1-i⁢v2)=(α-i⁢β)⋅(v1-i⁢v2),

czyli wektor v1-i⁢v2 jest wektorem własnym odpowiadającym wartości własnej λ2=α-i⁢β.

Pokażemy, że wektory v1,v2∈ℝ2 są liniowo niezależne. Rzeczywiście, przypuśćmy, że v1=p⋅v2 dla pewnej liczby rzeczywistej p∈ℝ. Wtedy

A⁢(v1+i⁢v2)=(α+i⁢β)⋅(v1+i⁢v2)=(α+i⁢β)⁢(p+i)⋅v2

oraz

A⁢(v1+i⁢v2)=A⁢((p+i)⁢v2)=(p+i)⋅A⁢(v2),

więc A⁢(v2)=(α+i⁢β)⋅v2, co prowadzi do sprzeczności z warunkiem β≠0.

Zauważmy, że

A⁢(v1)+i⁢A⁢(v2)=A⁢(v1+i⁢v2)=(α+i⁢β)⋅(v1+i⁢v2),

więc

A⁢(v1)=α⋅v1-β⋅v2,A⁢(v2)=β⋅v1+α⋅v2.

Dla macierzy S=[v1v2] otrzymujemy równość

Λ:=S-1⁢A⁢S=[αβ-βα].

Dla S1=[v2v1] mamy

S1-1⁢A⁢S1=[α-ββα].
Wniosek 11.9 (Rzeczywista postać Jordana).

Jeśli λ1=α+i⁢β, λ2=α-i⁢β z β≠0 są pierwiastkami wielomianu charakterystycznego pA dla macierzy A∈M2×2⁢(R), to istnieje taka nieosobliwa macierz S∈M2×2⁢(R), że

Λ=S-1⁢A⁢S=[αβ-βα].

Ponadto, S=[v1v2] dla dowolnego takiego niezerowego wektora v1+i⁢v2∈C, że

A⁢(v1+i⁢v2)=(α+i⁢β)⋅(v1+i⁢v2).
Przykład 11.3 (Ciąg Fibonacciego).

Ciąg Fibonacciego un jest zdefiniowany rekurencyjnie wzorem

un+1=un+un-1,u1=1,u0=0.

Dla macierzy

A=[1110]

mamy więc

[1110]⁢[unun-1]=[un+1un].

Wynika stąd, że

[un+1un]=An⁢[u1u0]=An⁢e1.

Wielomian charakterystyczny A ma postać pA⁢(x)=x2-x-1, więc wartościami własnymi A są

λ+=5+12,λ-=1-52.

Wektory własne mają postać

v+=[(1+5)/2,1]T,v-=[(1-5)/2,1]T.

Ponieważ

e1=15⋅(v+-v-),

więc

[un+1un]=15⁢(λ+n⋅v+-λ-n⋅v-),

czyli otrzymujemy formułę Bineta

un=15⁢(1+52)n-15⁢(1-52)n.
Przykład 11.4.

Rozważmy macierz

A=[31-11].

Jej wielomian charakterystyczny ma postać

pA⁢(x)=x2-4⁢x+4=(x-2)2.

Stąd λ=2 jest dwukrotną wartością własną macierzy A. Ponieważ

A-2⁢I=[11-1-1],

więc [x,y]T∈ker⁢(A-2⁢I) jeśli x+y=0, czyli

ker⁢(A-2⁢I)=span⁢{[1,-1]T}.

Wektor v=[1,-1]T jest wektorem własnym i każdy inny wektor własny jest liniowo zależny z v.

Szukamy teraz takiego wektora w=[x,y]T∈ℝ2, że

(A-2⁢I)⁢w=v=[1,-1]T.

Musimy więc rozwiązać układ równań

x+y=1,-x-y=-1.

Rozwiązaniami są wektory postaci [x,1-x]T=x⁢[1,-1]T+[0,1]T. Możemy więc przyjąć, że

w=[1,0]T=e1.

Dla macierzy

S=[11-10],

mamy

S-1⁢A⁢S=[0-111]⁢[31-11]⁢[11-10]=[2102].
Przykład 11.5.

Macierz

A=[0-410]

ma zespolone wartości własne λ1=2⁢i oraz λ2=-2⁢i. Znajdziemy wektor własny v1=[a+b⁢i,c+d⁢i]T∈ℂ2 dla λ1=2⁢i. Szukamy niezerowych rozwiązań równania

0=(A-2⁢i⁢I)⁢v =[-2⁢i-41-2⁢i]⁢[a+b⁢ic+d⁢i]
=[-2⁢i⁢a+2⁢b-4⁢c-4⁢d⁢ia+b⁢i-2⁢i⁢c+2⁢d]
=[-2⁢i⁢a+2⁢(b-2⁢c)-2⁢i⁢(a+2⁢d)a+2⁢d+i⁢(b-2⁢c)].

Stąd a+2⁢d=0 oraz b-2⁢c=0. Możemy więc przyjąć, że

v1=[2,-i]T=[2,0]T+i⁢[0,-1]T.

Dla wartości własnej λ2=-2⁢i wektorem własnym jest więc v2=[2,i]T. Dla macierzy

S=[22-ii]

mamy

S-1⁢A⁢S=14⁢i⁢[i-2i2]⁢[0-410]⁢[22-ii]=[2⁢i00-2⁢i]

Aby znaleźć rzeczywistą postać Jordana rozważamy macierz

P=[200-1]

której kolumnami są część rzeczywista i urojona wektora v1=[2,-i]T=[2,0]T+i⁢[0,-1]T. Wtedy

P-1⁢A⁢P=[1/200-1]⁢[0-410]⁢[200-1]=[02-20].

11.2 Rozkład Jordana w wymiarze 3

Podamy dowód Twierdzenia Jordana dla macierzy zespolonej A∈M3×3⁢(ℂ). Zastosujemy podejście, które uogólnia się na wyższe wymiary.

Twierdzenie 11.5 (Rozkład Jordana macierzy zespolonej).

Niech A∈M3×3⁢(C). Istnieje taka macierz nieosobliwa S∈M3×3⁢(C), że S-1⁢A⁢S jest jednej z poniższych postaci:

[λ1000λ2000λ3],[λ1000λ2100λ2],[λ1100λ1100λ1].

Ponadto, λi∈C są wartościami własnymi macierzy A.

Lemat 11.4.

Załóżmy, że wartości własne λ1,λ2,λ3∈C macierzy A∈M3×3⁢(C) są różne. Wtedy odpowiadające im wektory własne v1,v2,v3 są liniowo niezależne oraz dla S=[v1v2v3] mamy

S-1⁢A⁢S=diag⁢(λ1,λ2,λ3).
Dowód.

Zauważmy najpierw, że v1 i v2 są liniowo niezależne. Przypuśćmy, że v2=x⋅v1 dla pewnego x∈ℂ. Wtedy x≠0, bo v2≠0. Mamy

λ2⋅v2 =A⁢v2
=A⁢(x⋅v1)
=x⁢λ1⋅v1
=λ1⋅v2,

czyli λ1=λ2, sprzeczność.

Przypuśćmy, że v3∈span⁢{v1,v2}. Wtedy v3=x1⋅v1+x2⋅v2 dla pewnych x1,x2∈ℂ. Przynajmniej jedna z liczb x1, x2 jest różna od zera, bo v3≠0. Otrzymujemy, że

λ3⋅v3 =A⁢v3
=A⁢(x1⋅v1+x2⋅v2)
=x1⁢λ1⋅v1+x2⁢λ2⋅v2,

czyli

λ3⁢(x1⋅v1+x2⋅v2)=x1⁢λ1⋅v1+x2⁢λ2⋅v2.

W efekcie

x1⁢(λ3-λ1)⋅v1+x2⁢(λ2-λ1)⋅v2=0,

co jest sprzeczne z liniową niezależnością v1 i v2. ∎

Lemat 11.5.

Niech A∈M3×3⁢(C). Istnieje taka macierz nieosobliwa S∈M3×3⁢(C), że macierz S-1⁢A⁢S jest górnie trójkątna.

Dowód.

Niech λ∈ℂ będzie wartością własną A i niech v∈ℂ3 będzie odpowiadającym jej wektorem własnym. Uzupełniamy v do bazy v, u, w dla ℂ3. Rozważmy macierz W=[vuw]. Macierz W-1⁢A⁢W ma postać blokową

W-1⁢A⁢W=[λ**0ab0cd].

Z rozkładu Jordana w wymiarze 2 istnieje taka macierz nieosobliwa V∈M2×2⁢(ℂ), że

V-1⁢[abcd]⁢V=T

jest górnie trójkątna. Definiujemy macierz nieosobliwą

U=[100V]∈M3×3⁢(ℂ).

Wtedy

U-1⁢W-1⁢A⁢W⁢U =U-1⁢[λ**0ab0cd]⁢U
=[100V-1]⁢[λ**0ab0cd]⁢[100V]
=[λ*0T],

jest górnie trójkątna i teza zachodzi dla S=W⁢U. ∎

Lemat 11.6.

Rozważmy macierz górnie trójkątną z 0 jako jedyną wartością własną

A=[0ab00c000].

Wtedy zachodzi jeden z warunków

  • (i)

    A=0,

  • (ii)

    dim⁢ker⁢A=2 i A jest podobna do macierzy

    [010000000],
  • (iii)

    dim⁢ker⁢A=1 oraz A jest podobna do macierzy

    [010001000].
Dowód.

Sprawdzamy, że

A2=[0ab00c000]⁢[0ab00c000]=[00a⁢c000000],
A3=[0ab00c000]⁢[00a⁢c000000]=0.

Mamy trzy możliwości

  • (i)

    A=0,

  • (ii)

    A≠0, A2=0,

  • (iii)

    A2≠0, A3=0.

W przypadku (ii) a=0 lub c=0, bo A2=0 i macierz A jest jednej z postaci

A=[00b00c000] c≠0,A=[0ab000000] a≠0,A=[00b000000] b≠0.

We wszystkich przypadkach rank⁢A=dim⁢im⁢A=1 oraz

im⁢A⊂ker⁢A,dim⁢ker⁢A=2.

Inkluzja im⁢A⊂ker⁢A wynika stąd, że A2=0. Niech v1∈im⁢A będzie niezerowym wektorem. Istnieje więc taki (niezerowy) wektor v2∈ℂ3, że A⁢v2=v1. Wybieramy teraz wektor v3∈ker⁢A liniowo niezależny z v1. Rozważmy macierz

S=[v1v2v3].

Ponieważ v2∉ker⁢A, więc v1,v2,v3 są liniowo niezależne, czyli S jest nieosobliwa. Ponadto,

S-1⁢A⁢S=[010000000]

W przypadku (iii) mamy a≠0 i c≠0, więc

rank⁢A=dim⁢imA=2.

Ponadto,

im⁢A2⊂ker⁢A,dim⁢ker⁢A=1,więc im⁢A2=ker⁢A.

Niech v∈ker⁢A będzie wektorem własnym dla wartości własnej 0. Ponieważ

im⁢A2=ker⁢A,

więc istnieje taki wektor u∈ℂ3, że A2⁢u=v. Wektory

v,A⁢u,u

są liniowo niezależne. Rzeczywiście, przypuśćmy, że x⋅v+y⋅A⁢u+z⋅u=0 dla pewnych x,y,z∈ℂ. Wtedy

0=A2⁢(x⋅v+y⋅A⁢u+z⋅u)=z⋅v,

czyli z=0, bo v≠0. W konsekwencji,

0=A⁢(x⋅v+y⋅A⁢u)=y⋅v,

czyli y=0. Z równości x⋅v=0 również x=0. Dla

S=[vA⁢uu]

mamy

S-1⁢A⁢S=[010001000]

∎

Wniosek 11.10.

Załóżmy, że λ∈C jest 3-krotną wartością własną macierzy A∈M3×3⁢(C). Wtedy

  • (i)

    jeśli dim⁢ker⁢(A-λ⁢I)=3, to A=λ⁢I,

  • (ii)

    jeśli dim⁢ker⁢(A-λ⁢I)=2, to A jest podobna do macierzy

    [λ100λ000λ].
  • (iii)

    jeśli dim⁢ker⁢(A-λ⁢I)=1, to A jest podobna do macierzy

    [λ100λ100λ].
Dowód.

Z Lematu 11.5 możemy założyć, że A jest górnie trójkątna. Stosujemy Lemat 11.6 do macierzy A-λ⁢I. ∎

Lemat 11.7.

Załóżmy, że A∈M3×3⁢(C) jest macierzą górnie trójkątną postaci

A=[λab00c000],λ≠0.

Wtedy

  • (i)

    dim⁢ker⁢A<3,

  • (ii)

    jeśli dim⁢ker⁢A=2, to S-1⁢A⁢S=[λ00000000] dla pewnej macierzy nieosobliwej S∈M3×3⁢(ℂ),

  • (iii)

    jeśli dim⁢ker⁢A=1, to S-1⁢A⁢S=[λ00001000] dla pewnej macierzy nieosobliwej S∈M3×3⁢(ℂ).

Dowód.

Zauważmy, że

dim⁢im⁢A={1,gdy c=0,2,gdy c≠0, 

czyli

dim⁢ker⁢A={2,gdy c=0,1,gdy c≠0, 

Załóżmy, że dim⁢ker⁢A=2. Niech v2,v3∈ker⁢A będą bazą dla jądra. Jeśli v1 jest wektorem własnym dla λ, to v1,v2,v3 są liniowo niezależne. Przypuśćmy bowiem, że v1=x2⋅v2+x3⋅v3. Wtedy

λ⋅v1=A⁢v1=x2⋅A⁢v2+x3⋅A⁢v3=0,

sprzeczność. Dla macierzy nieosobliwej S=[v1v2v3] mamy

S-1⁢A⁢S=[λ00000000].

Niech teraz dim⁢ker⁢A=1, czyli c≠0. Zauważmy, że teza zajdzie, gdy znajdziemy taką bazę v1,v2,v3, że

A⁢v1=λ⋅v1,A⁢v2=0,A⁢v3=v2.

Z wyborem v1 i v2 nie mamy problemu: bierzemy wektory własne dla λ i 0. Musimy dobrać wektor v3, aby A⁢v3=v2. W tym celu wystarczy zauważyć, że

ker⁢A⊂im⁢A.

Rzeczywiście, A⁢[x,y,z]T=[λ⁢x+a⁢y+b⁢z,c⁢z,0]T, czyli [x,y,z]T∈ker⁢A, gdy z=0 i x=aλ⁢y, czyli jądro jest generowane przez wektor [a/λ,1,0]T. Sprawdzamy łatwo, że jest on w obrazie A generowanym przez wektory [λ,0,0]T i [b,c,0]T.

Pozostaje sprawdzić, że v1,v2,v3 są liniowo niezależne. Wiemy, że v1 i v2 są liniowo niezależne. Przypuśćmy, że v3=x1⋅v1+x2⋅v2. Wtedy

v2 =A⁢v3
=x1⋅A⁢v1+x2⋅A⁢v2
=x1⁢λ⋅v1,

co jest sprzeczne z liniową niezależnością v1 i v2. ∎

Wniosek 11.11.

Załóżmy, że λ,λ1∈C są różnymi wartościami własnymi macierzy A∈M3×3⁢(C) oraz λ1 ma krotność 2. Wtedy

  • (i)

    dim⁢ker⁢(A-λ1⁢I)<3,

  • (ii)

    jeśli dim⁢ker⁢(A-λ1⁢I)=2, to A jest podobna do macierzy

    [λ000λ1000λ1].
  • (iii)

    jeśli dim⁢ker⁢(A-λ1⁢I)=1, to A jest podobna do macierzy

    [λ000λ1100λ1].
Dowód.

Z Lematu 11.5 możemy założyć, że A jest górnie trójkątna. Stosujemy rozumowanie jak w dowodzie Lematu 11.7. ∎

Zajmiemy się teraz macierzami rzeczywistymi A∈M3×3⁢(ℝ). Wtedy wielomian charakterystyczny pA jest stopnia 3 i ma współczynniki rzeczywiste. Jeśli liczba zespolona λ∈ℂ jest pierwiastkiem pA, to jest nim również λ¯. Wynika stąd, że istnieje co najmniej jeden pierwiastek rzeczywisty pA. Jest on oczywiście wartością własną pA. Jeśli A ma trzy (liczone z krotnościami) rzeczywiste wartości własne, to A ma taką samą postać Jordana jak w przypadku zespolonym.

Pozostaje nam rozważyć przypadek, gdy pA ma rzeczywisty pierwiastek λ1∈ℝ i parę pierwiastków λ2=α+i⁢β oraz λ2¯=α-i⁢β z β≠0.

Wniosek 11.12 (Rzeczywista postać Jordana).

Jeśli dla A∈M3×3⁢(R) ma rzeczywisty pierwiastek λ1=λ∈R i parę pierwiastków sprzężonych λ2=α+i⁢β oraz λ2¯=α-i⁢β z β≠0, to istnieje taka nieosobliwa macierz S∈M3×3⁢(R), że

S-1⁢A⁢S=[λ000αβ0-βα].
Dowód.

Niech v∈ℝ3 będzie wektorem własnym dla λ, a w=w1+i⁢w2 zespolonym wektorem własnym dla α+i⁢β. Wtedy w¯=w1-i⁢w2 jest wektorem własnym dla α-i⁢β. Z naszych rozważań w wymiarze 2 wystarczy pokazać, że v,w1,w2 są liniowo niezależne nad ℝ, bo wtedy teza zachodzi dla S=[vw1w2]. Wiemy już, że w1 i w2 są liniowo niezależne. Przypuśćmy, że

v=x1⋅w1+x2⋅w2,

dla pewnych x1,x2∈ℝ. Ponieważ

w1=12⋅(w+w¯),w2=12⁢i⋅(w-w¯),

więc

v=x1-i⁢x22⋅w+x1+i⁢x22⋅w¯,

czyli v,w,w¯ są liniowo zależne nad ℂ. Prowadzi to do sprzeczności, bo wektory te odpowiadają różnym wartościom własnym. ∎

Przykład 11.6.

Macierz

A=[-30003-2011]

ma wartości własne λ1=-3, λ2=2+i, λ3=λ2¯=2-i. Odpowiadającymi im wektorami własnymi są przykładowo

v1=[1,0,0]T=e1,v2=[0,1+i,1]T,v3=[0,1-i,1]T.

Dla

S=[10001+i1-i011]

mamy

S-1⁢A⁢S=[-30002+i002-i].

Aby otrzymać postać rzeczywistą rozważamy macierz

P=[100011010]

której druga i trzecia kolumna, to odpowiednio, część rzeczywista i urojona wektora własnego v2=[0,1,1]T+i⁢[0,1,0]T. Otrzymujemy, że

P-1⁢A⁢P=[-3000210-12].
Przykład 11.7.

Wartości własne macierzy

A=[100-120112]

to λ1=1 i λ2=λ3=2. Sprawdzamy łatwo, że odpowiadające im wektory własne, to

v1=[1,1,-2]T,v2=[0,0,1]T.

Krotność algebraiczna wartości własnej λ2=2 jest równa 2, ale jej krotność geometryczna jest równa 1. Macierz A nie ma więc bazy wektorów własnych, czyli nie jest diagonalizowalna. Aby wyznaczyć postać Jordana A szukamy takiego wektora w, że

(A-2⁢I)⁢w=v2.

Ponieważ

A-2⁢I=[-100-100110],

więc możemy przyjąć, że w=[0,1,0]T. Dla macierzy

S=[100101-210]

mamy

S-1⁢A⁢S=[100201-110]⁢[100-120112]⁢[100101-210]=[100011001].
Przykład 11.8.

Macierz

A=[2100200-12]

ma tylko jedną wartość własną λ=2 o krotności algebraicznej 3. Zauważmy, że

A-2⁢I=[0100000-10],

więc

ker⁢(A-2⁢I)=span⁢{e1,e3}.

Ponadto,

im⁢(A-2⁢I)=span⁢{[1,0,-1]T}⊂ker⁢(A-2⁢I).

Szukamy teraz takiego wektora w, że

(A-2⁢I)⁢w=[1,0,-1]T.

Łatwo sprawdzamy, że możemy przyjąć w=[0,1,0]T=e2. Dla bazy

v1=[1,0,-1]T,v2=[0,1,0]T,v3=[1,0,0]T

Dla

S=[101010-100]

mamy

S-1⁢A⁢S=[00-1010100]⁢[2100200-12]⁢[101010-100]=[210020002]
Przykład 11.9.

Przyglądniemy się teraz macierzom antysymetrycznym A∈M3×3⁢(ℝ). Ponieważ AT=-A, więc A ma postać

A=[0ab-a0c-b-c0],a,b,c∈ℝ.

Zakładamy, że A jest niezerowa, czyli a2+b2+c2≠0. Ponieważ A jest antysymetryczna i wymiar n=3 jest nieparzysty, więc det⁢A=0. W szczególności, 0 jest zawsze wartością własną. Ma ona krotność 1, bo jak łatwo sprawdzić pozostałe wartości własne, to ±i⁢a2+b2+c2. Wektor własny dla wartości własnej 0 jest równy

vA=[-c,b,-a]T.

Ponadto, jeśli A,B∈M3×3⁢(ℝ) są antysymetryczne, to

tr(AB)=-2(vA|vB).

Co więcej, macierz

[A,B]:=A⁢B-B⁢A

jest antysymetryczna oraz

v[A,B]=vA×vB.

11.3 Diagonalizacja w wymiarach 2 i 3

Niech

A=[accb]

będzie rzeczywistą macierzą symetryczną. Równanie charakterystyczne dla A ma postać

pA⁢(x)=x2-(a+b)⁢x+a⁢b-c2=0.

Ponieważ

Δ:=(a+b)2-4⁢(a⁢b-c2)=(a-b)2+4⁢c2≥0,

więc wartości własne λ1, λ2 macierzy A są rzeczywiste oraz

λ1=a+b+Δ2,λ2=a+b-Δ2.
Twierdzenie 11.6.

Dla macierzy rzeczywistej A∈M2×2⁢(R) następujące warunki są równoważne

  • (i)

    A=AT,

  • (ii)

    ℝ2 ma bazę ortonormalną złożoną z wektorów własnych A,

  • (iii)

    Istnieje taka macierz ortogonalna Q∈M2×2⁢(ℝ), że

    QT⁢A⁢Q=D,

    gdzie D=diag⁢(λ1,λ2) jest diagonalna z wartościami własnymi A na przekątnej.

Dowód.

(i)⇒(i⁢i)

Załóżmy, że A jest symetryczna. Jeśli c=0, to e1, e2 jest bazą ortonormalną złożoną z wektorów własnych macierzy A=diag⁢(a,b).

Jeśli c≠0, to Δ>0 i A ma dwie różne rzeczywiste wartości własne λ1, λ2. Niech v1, v2 będą odpowiadającymi im wektorami własnymi o normie 1. Pokażemy, że wektory v1, v2 są ortogonalne. Mamy

λ1(v1|v2) =(λ1v1|v2)
=(Av1|v2)
=(v1|Av2)
=(v1|λ2v2)
=λ2(v1|v2),

więc (v1|v2)=0, bo λ1≠λ2.

(i⁢i)⇒(i⁢i⁢i) Wystarczy przyjąć Q=[v1v2] dla bazy ortonormalnej v1, v2 złożonej z wektorów własnych.

(i⁢i⁢i)⇒(i)

Wiemy, że A=Q⁢D⁢QT. Wtedy

AT =(QT)T⁢DT⁢QT
=Q⁢D⁢QT
=A.

∎

Rozważmy teraz macierz symetryczną A∈M3×3⁢(ℝ). Istnieje λ1∈ℝ rzeczywista wartość własna dla A. Niech v1 będzie odpowiadającym jej wektorem własnym o normie 1. Niech V⊂ℝ3 będzie płaszczyzną prostopadłą do wektora v1. Wybierzmy bazę ortonormalną v2, v3∈V dla V. Wtedy v1, v2, v3 jest bazą ortonormalną dla ℝ3.

Pokażemy, że A⁢v2, A⁢v3∈V, czyli A⁢(V)⊂V. Ze względu na symetrię wystarczy pokazać, że A⁢v2∈V. Ponieważ A jest symetryczna, więc

(v1|Av2) =(Av1|v2)
=(λ1v1|v2)
=λ1(v1|v2)
=0.

Rozważmy macierz ortogonalną W=[v1v2v3]. Wtedy

W-1⁢A⁢W=WT⁢A⁢W=[λ1000ac0db].

Ponieważ WT⁢A⁢W jest symetryczna, więc d=c, czyli

WT⁢A⁢W=[λ1000ac0cb]

Wiemy już, że istnieje taka macierz ortogonalna B∈M2×2⁢(ℝ), że

BT⁢[accb]⁢B=diag⁢(λ2,λ3).

Dla macierzy ortogonalnej

U=[100B]∈M3×3⁢(ℝ)

mamy

UT⁢WT⁢A⁢W⁢U=diag⁢(λ1,λ2,λ3).

Macierz Q=W⁢U jest ortogonalna i jej kolumny są bazą ortonormalną złożoną z wektorów własnych.

Wniosek 11.13.

Dla macierzy rzeczywistej A∈M3×3⁢(R) następujące warunki są równoważne

  • (i)

    A=AT,

  • (ii)

    ℝ3 ma bazę ortonormalną złożoną z wektorów własnych A,

  • (iii)

    Istnieje taka macierz ortogonalna Q∈M3×3⁢(ℝ), że

    QT⁢A⁢Q=D,

    gdzie D=diag⁢(λ1,λ2,λ3) jest diagonalna z wartościami własnymi A na przekątnej.

11.4 Twierdzenie Cayleya–Hamiltona w wymiarze 3

Twierdzenie 11.7 (Twierdzenie Cayleya–Hamiltona).

Niech λ1,λ2,λ3∈C będą wartościami własnymi macierzy A∈M3×3⁢(C). Wtedy

(A-λ1⁢I)⁢(A-λ2⁢I)⁢(A-λ3⁢I)=0.
Dowód.

Z Lematu 11.5 istnieje taka macierz nieosobliwa S∈M3×3⁢(ℂ), że S-1⁢A⁢S jest górnie trójkątna. Zauważmy, że

(S-1⁢A⁢S-λ1⁢I) (S-1⁢A⁢S-λ2⁢I)⁢(S-1⁢A⁢S-λ3⁢I)=
(S-1⁢A⁢S-λ1⁢S-1⁢S)⁢(S-1⁢A⁢S-λ2⁢S-1⁢S)⁢(S-1⁢A⁢S-λ3⁢S-1⁢S)=
S-1⁢(A-λ1⁢I)⁢S⁢S-1⁢(A-λ2⁢I)⁢S⁢S-1⁢(A-λ3⁢I)⁢S=
S-1⁢(A-λ1⁢I)⁢(A-λ2⁢I)⁢(A-λ3⁢I)⁢S.

Wynika stąd, że teza zachodzi dla A wtedy i tylko wtedy, gdy zachodzi dla S-1⁢A⁢S. Możemy więc założyć, że A jest górnie trójkątna postaci

A=[λ1ab0λ2c00λ3]

Niech v=[x,y,z]T∈ℂ3 będzie dowolnym wektorem. Zauważmy, że

(A-λ3⁢I)⁢v=[x1,y1,0]T,

dla pewnych x1,y1∈ℝ. Następnie,

(A-λ2⁢I)⁢[x1,y1,0]T=[x2,0,0]T,

dla pewnego x2∈ℂ. Ostatecznie,

(A-λ1⁢I)⁢[x2,0,0]T=[0,0,0]T,

czyli

(A-λ1⁢I)⁢(A-λ2⁢I)⁢(A-λ3⁢I)⁢v=0.

∎

  • U

    Z pewnością zauważyliście, że dokładnie te same argumenty pokazują, że dla macierzy górnie trójkątnej A∈Mn×n⁢(ℂ) zachodzi

    pA⁢(A)=0.

    Dowód opiera się na obserwacji, że dla macierzy górnie trójkątnej A podprzestrzenie span⁢{e1,…,ek} są niezmiennicze dla A. Pokażemy później, że każda macierz A∈Mn×n⁢(ℂ) jest podobna do macierzy górnie trójkątnej, więc Twierdzenie Cayleya–Hamiltona zachodzi dla dowolnej macierzy A∈Mn×n⁢(ℂ).

  • 𝗧𝗘𝗦𝗧→

    Niech A=[abcd]∈M2×2⁢(ℝ). Oceń prawdziwość poniższych zdań:

    • –

      Jeśli a+b=c+d=1, to [1,1]T jest wektorem własnym A.

    • –

      Jeśli A ma prostą niezmienniczą, to A ma rzeczywiste wartości własne.

    • –

      A ma rzeczywisty wektor własny.

    • –

      Jeśli A jest górnie trójkątna tzn. c=0, to A jest diagonalizowalna.

    • –

      Jeśli det⁢A<0, to A jest diagonalizowalna.

    • –

      Jeśli det⁢A<0, to A jest ortogonalnie diagonalizowalna.

    • –

      Jeśli A3=0, to A2=0.

    • –

      Jeśli ℝ2 ma bazę złożoną z wektorów własnych A, to A jest nieosobliwa.

    • –

      0 jest jedyną wartością własną A wtedy i tylko wtedy, gdy A2=0.

    • –

      Jeśli tr⁢A=0, to A jest diagonalizowalna lub A2=0.

    • –

      Jeśli det⁢A=0, to A ma prostą niezmienniczą.

    • –

      A i A2 mają takie same wartości własne.

    • –

      A i A2 mają takie same wektory własne.

    • –

      Jeśli A jest nieosobliwa, to A i A-1 mają takie same wartości własne.

    • –

      Jeśli A jest nieosobliwa, to A i A-1 mają takie same wektory własne.

    • –

      A i AT mają takie same wartości własne.

    • –

      A i AT mają takie same wektory własne.

    • –

      Jeśli S∈M2×2⁢(ℝ) jest nieosobliwa, to A i S-1⁢A⁢S mają takie same wartości własne.

    • –

      Jeśli S∈M2×2⁢(ℝ) jest nieosobliwa, to A i S-1⁢A⁢S mają takie same wektory własne.

    • –

      Jeśli A jest nieosobliwa i diagonalizowalna, to A-1 jest diagonalizowalna.

    • –

      Jeśli macierz kwadratowa A jest diagonalizowalna, to AT jest również diagonalizowalna.

    • –

      Jeśli macierz kwadratowa ma dwa identyczne wiersze (kolumny), to 0 jest jej wartością własną.

⟵