10.11.2014 Views

Uge 9.pdf

Uge 9.pdf

Uge 9.pdf

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

SEKTION 9.1<br />

EGENVÆRDIER OG EGENVEKTORER<br />

9.1 Egenværdier og egenvektorer<br />

Definition 9.1.1<br />

1. Lad V være et F-vektorrum; og lad T : V → V være en lineær transformation. λ ∈ F<br />

er en egenværdi for T , hvis der findes v ∈ V \{0}, så<br />

T (v) =λv.<br />

v kaldes da en egenvektor for T associeret til, eller svarende til, λ.<br />

2. Lad A ∈ Mat n,n (F). λ ∈ F er en egenværdi for A hvis der findes z ∈ F n \{0}, så<br />

Az = λz.<br />

z kaldes da en egenvektor for A associeret til, eller svarende til, λ.<br />

Eksempler 9.1.2<br />

1. Ethvert tal λ ∈ R er en egenværdi for differentiation D : C ∞ (R) → C ∞ (R); f(t) =e λt<br />

er en tilsvarende egenvektor.<br />

(En egenvektor i et funktionsrum kaldes ofte en egenfunktion.)<br />

2. √ 2, − √ 2 er egenværdier for<br />

[ ] 1 1<br />

∈ Mat<br />

1 −1<br />

2,2 (R),<br />

idet<br />

[ ][ √ ]<br />

1 1 1+ 2<br />

1 −1 1<br />

= √ [ √<br />

1+ 2<br />

2<br />

1<br />

]<br />

,<br />

[ ][ √ ]<br />

1 1 1 − 2<br />

1 −1 1<br />

= − √ [ √<br />

1 − 2<br />

2<br />

1<br />

]<br />

.<br />

Lemma 9.1.3<br />

Lad V være et F-vektorrum med ordnet basis V = {v 1 , . . . , v n }, og lad T : V → V være en<br />

lineær transformation.<br />

Så er λ en egenværdi for T ⇔ λ er en egenværdi for M(T ) V,V ,<br />

og v ∈ V er en egenvektor for T associeret til λ ⇔ [v] V er en egenvektor for M(T ) V,V<br />

associeret til λ.<br />

Bevis<br />

Da koordinatisering mht. V giver en isomorfi V → F n , gælder<br />

T (v) =λv ⇔ [T (v)] V =[λv] V ⇔ M V,V (T )[v] V = λ[v] V .<br />

Da v = 0 ⇔ [v] V = 0 følger udsagnet umiddelbart.<br />

155


SEKTION 9.1<br />

EGENVÆRDIER OG EGENVEKTORER<br />

Pga. lemma 9.1.3 vil vores udregning af egenværdier og egenvektorer næsten altid foretages<br />

for matricer – men det er oftest egenskaber af lineære transformationer, som vi dermed prøver at<br />

afdække. Se f.eks. Application 1 og Application 2 i [L], 6.1.<br />

Det er nemlig sådan, at vi kan udregne egenværdier og egenvektorer for en matrix på algoritmisk<br />

vis:<br />

Proposition 9.1.4 ([L], s. 302)<br />

Lad A ∈ Mat n,n (F), λ ∈ F. Følgende er ækvivalente:<br />

(a) λ er en egenværdi for A.<br />

(b) Ligningssystemet (A − λI)x = 0 har en ikke-triviel løsning.<br />

(c) N(A − λI) ≠ {0}.<br />

(d) A − λI er singulær.<br />

(e) det(A − λI) = 0.<br />

Bevis<br />

(a)⇒(b): Der findes z ∈ F n \{0} så Az = λz, så 0 = Az − λz =(A − λI)z, og z er en<br />

ikke-triviel løsning til (A − λI)x = 0.<br />

(b)⇒(c): En ikke-triviel løsning z til (A − λI)x = 0 ligger i N(A − λI), og er ikke 0.<br />

(b)⇔(d): Sætning 1.4.8.<br />

(d)⇔(e): Sætning 8.1.15.<br />

(c)⇒(a): Lad z ∈ N(A − λI) \{0}, så 0 =(A − λI)z = Az − λz og Az = λz. Da z ≠ 0 er λ<br />

en egenværdi for A.<br />

Proposition 9.1.5<br />

Hvis det(A − λI) beregnes for et ubestemt λ, så fås et polynomium af grad n i λ,<br />

p A (λ) = det(A − λI);<br />

vi har<br />

p A (λ) =(−1) n λ n +(−1) n tr(A)λ n−1 + · · · + det(A)<br />

hvor tr(A) =a 11 + · · · + a nn er summen af diagonalindgangene i A.<br />

156


SEKTION 9.1<br />

EGENVÆRDIER OG EGENVEKTORER<br />

Bevis<br />

Den (i,j)’te indgang i A − λI er {<br />

a ij<br />

a ii − λ<br />

Korollar 8.2.5 foretæller, at<br />

det(A) = ∑<br />

hvis i ≠ j<br />

hvis i = j<br />

σ∈S n<br />

a σ(1),1 . . . a σ(n),1 sgn(σ)<br />

anvendt på A − λI ser vi, at det(A − λI) er en sum af produkter af n elementer, som er enten<br />

en a ij ,i≠ j eller en a ii − λ.<br />

p A (λ) = det(A − λI) er således et polynomium af grad højst n; det er dog af grad n, idet<br />

leddet (a 11 − λ)(a 22 − λ) . . . (a nn − λ) er det eneste af de indgående produkter, hvor λ n<br />

optræder. Vi ser, at den optræder med koefficient (−1) n .<br />

Hvis σ ∈ S n ikke er identiteten, så er σ(i) ≠ i for mindst to i ∈{1, . . . , n}, (hvis σ afbilder<br />

n − 1 af 1, . . . , n til sig selv, så må den også afbilde det sidste element til sig selv, idet den er<br />

surjektiv), så λ n−1 optræder i p A (λ) også kun i produktet (a 11 − λ)(a 22 − λ) . . . (a nn − λ): vi<br />

ser, at den optræder med koefficient<br />

(−1) n−1 (a 11 + · · · + a nn ) = (−1) n−1 tr(A).<br />

Endelig er konstant-leddet i p A (λ) det(A), idet dette konstante led er<br />

p A (0) = det(A − 0I) = det(A).<br />

Notation 9.1.6<br />

p A kaldes det karakteristiske polynomium for A.<br />

Korollar 9.1.7<br />

Lad A ∈ Mat n,n (F). λ 0 er en rod for p A (dvs. p A (λ 0 )=0) ⇔ λ 0 er en egenværdi for A, og<br />

hvis λ 0 er en egenværdi for A så er v en tilsvarende egenvektor ⇔ v ∈ N(A − λ 0 I) \{0}.<br />

Vi kan således finde samtlige egenværdier og egenvektorer for A ∈ Mat n,n (F):<br />

1. Find alle rødderne λ 1 , . . . , λ k for p A .<br />

2. For i =1, . . . , k, find en basis {z i1 , . . . , z i,ni } for N(A − λ i I) (her n i = dim N(A − λ i I)).<br />

Så er λ 1 , . . . , λ k de eneste egenværdier for A og for i =1, . . . , k er egenvektorerne for λ i givet som<br />

ikke-trivielle lineære kombinationer af z i1 , . . . , z i,ni .<br />

157


SEKTION 9.1<br />

EGENVÆRDIER OG EGENVEKTORER<br />

Notation 9.1.8<br />

N(A − λ i I) kaldes egenrummet E A (λ i ) svarende til λ i : og dets dimension kaldes den geometriske<br />

multiplicitet Geo A (λ i ) af λ i .<br />

Den algebraiske multiplicitet Alg A (λ 0 ) af λ 0 er derimod det antal gange, λ−λ 0 går op i p A (λ):<br />

hvis p A (λ) =(λ − λ 0 ) m0 q(λ), hvor q(λ 0 ) ≠0, så er m 0 = Alg A (λ 0 ).<br />

Eksempel 9.1.9<br />

Lad<br />

⎡<br />

9 0<br />

⎤<br />

7<br />

A = ⎣−3 2 −3⎦ ∈ Mat 3,3 (R).<br />

−8 0 −6<br />

⎡<br />

⎤<br />

9 − λ 0 7<br />

p A (λ) = det ⎣ −3 2− λ −3 ⎦<br />

−8 0 −6 − λ<br />

[ ]<br />

=(−1) 2+2 9 − λ 7<br />

(2 − λ)<br />

−8 −6 − λ<br />

= (2 − λ)((9 − λ)(−6 − λ) + 56) = (2 − λ)(λ 2 − 3λ + 2)<br />

= (2 − λ) 2 (1 − λ).<br />

Egenværdierne er derfor λ 1 =1, med algebraisk multiplicitet 1, og λ 2 =2med algebraisk<br />

multiplicitet 2.<br />

⎡<br />

⎤ ⎡<br />

⎤<br />

8 0 7 R 3 →R 3 +R 1 8 0 7<br />

R<br />

A − = ⎣−3 1 −3⎦<br />

2→8R<br />

∼<br />

2<br />

⎣−24 8 −24⎦<br />

−8 0 −7<br />

0 0 0<br />

R 2→R<br />

∼<br />

2+3R 1<br />

⎡<br />

⎣ 8 0 7<br />

⎤<br />

0 8 −3⎦<br />

0 0 0<br />

R 1 → 1 8 R 1<br />

R 2→ 1 8 R2<br />

∼<br />

⎡<br />

1 0<br />

⎤<br />

7/8<br />

⎣0 1 −3/8⎦ ,<br />

0 0 0<br />

så<br />

Geo(λ 1 ) = 1.<br />

og<br />

Geo(λ 2 )=2.<br />

E A (λ 1 )=N(A −<br />

⎛⎡<br />

) = Span ⎝⎣ −7 ⎤⎞<br />

3 ⎦⎠ ;<br />

8<br />

⎡<br />

⎤ ⎡ ⎤<br />

7 0 7 1 0 1<br />

A − 2 = ⎣−3 0 −3⎦ ∼ ⎣0 0 0⎦<br />

−8 0 −8 0 0 0<br />

⎛⎡<br />

⎤ ⎡ ⎤⎞<br />

0 −1<br />

E A (λ 2 )=N(A − 2 ) = Span ⎝⎣1⎦ , ⎣ 0 ⎦⎠ ;<br />

0 1<br />

158


SEKTION 9.1<br />

EGENVÆRDIER OG EGENVEKTORER<br />

Eksempel 9.1.10<br />

Lad<br />

[ ]<br />

cos θ − sin θ<br />

A =<br />

∈ Mat<br />

sin θ cos θ<br />

2,2 (R).<br />

A er SMR for rotationen R θ af planen gennem en vinkel θ, hvis θ ≠0,π, så er Av, v ikke<br />

parallelle for v ∈ R 2 \{0}, så A har ingen egenværdi.<br />

Hvis vi udregner<br />

ser vi, at p A (λ) har rødder<br />

p A (λ) = (cos θ − λ) 2 + sin 2 θ = λ 2 − 2λ cos θ +1<br />

2 cos θ ± √ 4 cos 2 θ − 4<br />

2<br />

så rødderne er ikke reelle når sin θ ≠0, dvs. når θ ≠0,π.<br />

= cos θ ± i sin θ,<br />

Sætning 9.1.11 (Algebraens Fundamentalsætning)<br />

Lad p være et ikke-konstant komplekst polynomium. Da har p en rod. Faktisk kan p skrives<br />

på formen<br />

p(x) =a(x − x 1 ) n1 . . . (x − x k ) n k<br />

hvor a ∈ C \{0} og x 1 , . . . , x r ∈ C er p’s rødder.<br />

Vi vil ikke vise denne sætning i dette kursus, der henvises i stedet for til bøger eller kurser i<br />

kompleks analyse.<br />

Notation 9.1.12<br />

n i er multipliciteten af x i som rod af p; p har n = n 1 + · · · + n k rødder talt med multiplicitet.<br />

Korollar 9.1.13<br />

Lad A ∈ Mat n,n (C). Da har Anegenværdier, talt med algebraisk multiplicitet.<br />

Når en matrix har reelle indgange, så kan den også opfattes som en kompleks matrix, fordi de<br />

reelle tal er indlejrede i de komplekse. Så vi kan, og ofte vil, betragte komplekse egenværdier for en<br />

reel kvadratisk matrix. Mange udredninger og udregninger bliver meget nemmere, når komplekse<br />

egenværdier tages i betragtning, også selv om det er reel information, der søges.<br />

159


SEKTION 9.1<br />

EGENVÆRDIER OG EGENVEKTORER<br />

Lemma 9.1.14<br />

Lad A ∈ Mat n,n (R); lad λ ∈ C \ R være en egenværdi for A, med tilsvarende egenvektor v.<br />

1. ¯λ er en egenværdi for A, med tilsvarende egenvektor ¯v.<br />

2. λ, ¯λ har den samme geometriske multiplicitet som egenværdier for A.<br />

Bevis<br />

1. Vi har Av = λv. Konjuger: Av = λv; så A¯v = ¯λ¯v, idet A har reelle indgange.<br />

2. Antag, at {v 1 , . . . , v k } er en basis for E λ (A). Så er ¯v 1 , . . . , ¯v k ∈ E¯λ(A). De er uafhængige,<br />

for<br />

Så dim E¯λ(A) ≥ dim E λ (A).<br />

c 1 ¯v 1 + · · · + c k ¯v k =0⇒ ¯c 1 v 1 + · · · +¯c k v k =0<br />

⇒ ¯c 1 =0, . . . , ¯c k =0<br />

⇒ c 1 =0, . . . , c k =0,<br />

Med λ erstattet med ¯λ i det ovenstående fås dim E¯λ(A) ≥ dim E¯λ(A);<br />

da ¯λ = λ har vi vist<br />

dim E λ (A) = dim E¯λ(A).<br />

Lemma 9.1.15<br />

Lad A ∈ Mat n,n (R); lad λ ∈ R være en egenværdi for A. Så har λ den samme geometrisk<br />

multiplicitet, uanset om A betragtes som en reel eller en kompleks matrix.<br />

Bevis<br />

Lad<br />

E R A(λ) ={x ∈ R | (A − λI)x = 0}, E C A(λ) ={z ∈ C | (A − λI)z = 0};<br />

så dim EA R (λ) er den geometriske multiplicitet af λ som egenværdi af A betragtet som reel<br />

matrix, mens dim EA C (λ) er den geometriske multiplicitet af λ som egenværdi af A betragtet<br />

som kompleks matrix.<br />

Lad A−λ ∼ H i RREF. Da A−λ har reelle indgange kan denne rækkeækvivalens fås ved<br />

reelle rækkeoperationer, og i så fald gælder A − λ ∼ H uanset om A − λ og H betragtes<br />

som reelle eller komplekse matricer.<br />

Vi har da, ifølge Proposition 3.2.8, 1, at<br />

som påstået.<br />

dim E R A(λ) =Antal søjler uden pivot i H = dim E C A(λ),<br />

160


SEKTION 9.1<br />

EGENVÆRDIER OG EGENVEKTORER<br />

Lemma 9.1.16 ([L], 6.1.1)<br />

Lad A, B ∈ Mat n,n (F) være similære.<br />

1. p A = p B .<br />

2. Hvis λ 0 er en egenværdi for A og B, så er Geo A (λ 0 ) = Geo B (λ 0 ).<br />

Bevis<br />

1. Der gælder, at B = S −1 AS, hvor S ∈ Mat n,n (F) er invertibel.<br />

Vi har for et ubestemt λ<br />

Så<br />

S −1 (A − λ )S = S −1 AS − S −1 (λ )S = B − λ .<br />

p B (λ) = det(B − λ )<br />

= det(S −1 (A − λ )S)<br />

= det(S −1 ) det(A − λ ) det(S)<br />

= det(A − λ )<br />

= p A (λ).<br />

(Vi har brugt, at det(S −1 ) det(S) = det(S −1 S) = det( ) = 1).<br />

2. Lad {v 1 , . . . , v k } være en basis for N(B − λ 0 ). Vi har da for i =1, . . . , n<br />

så Sv 1 , . . . , Sv k ∈ N(A − λ 0 ).<br />

Sv 1 , . . . , Sv k er uafhængige: for antag, at<br />

Så er<br />

(A − λ 0 )Sv i = S(B − λ 0 )v i =0,<br />

c 1 Sv 1 + · · · + c k Sv k = 0, så S(c 1 v 1 , . . . , c k v k )=0.<br />

c 1 v 1 + · · · + c k v k = S −1 S(c 1 v 1 + · · · + c k v k )=S −1 0 = 0;<br />

og c 1 =0, . . . , c k =0fordi v 1 , . . . , v k er uafhængige.<br />

Så<br />

dim N(A − λ 0 ) ≥ dim N(B − λ 0 ).<br />

På samme måde ses, at hvis {w 1 , . . . , w l } er en basis for N(A − λ ), så er<br />

S −1 w 1 , . . . , S −1 w l ∈ N(B − λ ), og uafhængige, så<br />

dim N(A − λ 0 ) ≤ dim(B − λ 0 ).<br />

Så<br />

Geo A (λ 0 ) = dim N(A − λ 0 ) = dim N(B − λ 0 ) = Geo B (λ 0 ).<br />

161


SEKTION 9.2<br />

DIAGONALISERING<br />

9.2 Diagonalisering<br />

Definition 9.2.1<br />

1. Lad L : V → V være en lineær transformation.<br />

L er diagonaliserbar hvis der findes en basis V = {v 1 , . . . , v n } for V bestående af<br />

egenvektorer for L.<br />

Hvis egenværdien svarende til v i er λ i , for i =1, . . . , n, så er<br />

M V,V (L) = [[L(v 1 )] V , . . . , [L(v n )] V ]<br />

= [[λ 1 v 1 ] V , . . . , [λ n v n ] V ]<br />

=[λ 1 e 1 , . . . , λ n e n ]<br />

⎡<br />

⎤<br />

λ 1 0<br />

⎢<br />

= ⎣<br />

. ..<br />

⎥<br />

⎦ ,<br />

0 λ n<br />

en diagonalmatrix.<br />

Definition 9.2.2, fortsat<br />

2. Lad A ∈ Mat n,n (F). A er diagonaliserbar, hvis der findes en basis V = {v 1 , . . . , v n } for<br />

F n bestående af egenvektorer for A.<br />

L A er således diagonaliserbar, og, som ovenfor, hvis egenværdien svarende til v i er λ i for<br />

i =1, . . . , n, så er<br />

⎡<br />

⎤<br />

λ 1<br />

⎢<br />

M V,V (L A )= ⎣<br />

. ..<br />

⎥<br />

⎦ .<br />

λ n<br />

Ifølge Korollar 4.2.9 er<br />

M V,V (L A )=K V,E M E,E (L A )K E,V ,<br />

hvor K E,V er koordinatskiftematricen til E-koordinater (dvs. standard-koordinater) fra V-koordinater<br />

(se slides 3.3).<br />

Skriv V =[v 1 , . . . , v n ] i søjleform.<br />

Ifølge Lemma 3.3.5 (se også [L], s. 158-9) er K E,V = V, K V,E = V −1 ; da M E,E(LA ) = A har vi<br />

⎡<br />

⎤<br />

λ 1<br />

. ..<br />

λ n<br />

⎢<br />

⎣<br />

⎥<br />

⎦ = V −1 AV.<br />

Leon ([L], s. 326) definerer A til at være diagonaliserbar vha. denne formel. Definitionerne er ensbetydende:<br />

162


SEKTION 9.2<br />

DIAGONALISERING<br />

Lemma 9.2.3 ([L], 6.3.2)<br />

Lad A ∈ Mat n,n (F). Følgende er ækvivalente:<br />

(1) Der findes en basis for F n bestående af egenvektorer for A.<br />

(2) Der findes n lineært uafhængige egenvektorer for A.<br />

(3) Der findes en invertibel matrix V ∈ Mat n,n (F) så V −1 AV er en diagonalmatrix.<br />

Bevis<br />

(1)⇒(2): er oplagt.<br />

(2)⇒(1): 3.1.4, (eller [L], 3.4.3,I).<br />

(2)⇒(3): Lad v 1 , . . . , v n være lineært uafhængige egenvektorer for A; lad λ 1 , . . . , λ n være<br />

de tilsvarende egenværdier. Skriv V =[v 1 , . . . , v n ] i søjleform. Vi har<br />

hvor<br />

er diagonal.<br />

AV = A[v 1 , . . . , v n ]=[Av 1 , . . . , Av n ]<br />

⎡<br />

⎤<br />

λ 1 0<br />

⎢<br />

=[λ 1 v 1 , . . . , λ n v n ] = [v 1 , . . . , v n ] ⎣<br />

. ..<br />

⎥<br />

⎦<br />

0 λ n<br />

= VD<br />

⎡<br />

⎤<br />

λ 1 0<br />

⎢<br />

D = ⎣<br />

. ..<br />

⎥<br />

⎦<br />

0 λ n<br />

Da V har uafhængige søjler er den invertibel. Vi har derfor V −1 AV = V −1 VD= D.<br />

(3)⇒(2): Antag, at der findes invertibel X ∈ Mat n,n (F) så X −1 AX = D, hvor D er diagonal,<br />

⎡<br />

⎤<br />

λ 1 0<br />

⎢<br />

D = ⎣<br />

. ..<br />

⎥<br />

⎦ .<br />

0 λ n<br />

Skriv X =[x 1 , . . . , x n ] i søjleform. Da X er invertibel, er x 1 , . . . , x n uafhængige.<br />

Vi har AX = X(X −1 AX) = XD, dvs. [Ax 1 , . . . , Ax n ] = [λ 1 x 1 , . . . , λ n x n ]. Men så er<br />

Ax i = λ i x i for i =1, . . . , n; x 1 , . . . , x n er n uafhængige egenvektorer for A.<br />

Lemmaet undgår behændigt at nævne koordinatskift; alligevel er det nyttigt at huske, specielt i<br />

anvendelser, at et koordinatskift er involveret.<br />

163


SEKTION 9.2<br />

DIAGONALISERING<br />

Det er nemt at regne med potenser af en diagonaliserbar matrix:<br />

Lemma 9.2.4<br />

Lad A ∈ Mat n,n (F) være diagonaliserbar, og lad v 1 , . . . , v n være en basis af F n bestående<br />

af egenvektorer for A med tilsvarende egenværdier λ 1 , . . . , λ n .<br />

1. Der gælder, at<br />

for alle c 1 , . . . , c n ∈ F og k ∈ N.<br />

A k (c 1 v 1 + · · · + c n v n )=c 1 λ k 1v 1 + · · · + c n λ k nv n<br />

(∗)<br />

2. Lad V =[v 1 , . . . , v n ]. Der gælder, at<br />

⎡<br />

A k = V −1 ⎢<br />

⎣<br />

λ k 1<br />

⎤<br />

⎥<br />

⎦ V<br />

. ..<br />

λ k n<br />

for k ∈ N.<br />

Bevis<br />

1. For k ∈ N er<br />

A k (c 1 v 1 + · · · + c n v n )=c 1 A k v 1 + · · · + c n A k b n<br />

= c 1 λ k 1v 1 + · · · + c n λ k nv n<br />

2. (∗) kan omskrives som<br />

⎡ ⎤<br />

1<br />

A k ⎢<br />

V ⎣c<br />

.<br />

c n<br />

⎥<br />

⎦ = V<br />

⎡<br />

⎢<br />

⎣<br />

λ k 1<br />

⎤ ⎡ ⎤<br />

1<br />

. ..<br />

⎥ ⎢<br />

⎦ ⎣c<br />

.<br />

λ k n c n<br />

⎥<br />

⎦ for alle<br />

⎡ ⎤<br />

1<br />

⎢<br />

⎣c<br />

.<br />

c n<br />

⎥<br />

⎦ ∈ F n ,<br />

så<br />

A k V = V<br />

⎡<br />

⎢<br />

⎣<br />

λ k 1<br />

. ..<br />

λ k n<br />

⎤<br />

⎡<br />

⎥<br />

⎦ og A k ⎢<br />

= V ⎣<br />

λ k 1<br />

. ..<br />

λ k n<br />

⎤<br />

⎥<br />

⎦ V −1 .<br />

164


SEKTION 9.2<br />

DIAGONALISERING<br />

Eksempel 9.2.5 ([L] 3.5, Application 1)<br />

Vi har set på dette eksempel tidligere (i begyndelsen af sektion 3.3). Det handler om befolkningsændringer<br />

i en storby: hvert år flytter 6% af befolkningen fra midtbyen til omegnskommuner,<br />

mens 2% flytter den anden vej. Hvis 30% af befolkningen bor i midtbyen, 70% i omegnskommunerne<br />

nu, hvad bliver fordelingen ad årene?<br />

Vi skrev<br />

[ ]<br />

0, 94 0, 02<br />

A =<br />

, x<br />

0, 06 0, 98 0 =<br />

[ 0, 30<br />

0, 70<br />

så giver x n = Ax n−1 ,n=1, 2, . . . denne fordeling efter n år.<br />

Vores løsningsmetode var egentlig at identificere egenvektorer<br />

[ [ ]<br />

1 −1<br />

u 1 = , u<br />

3]<br />

2 =<br />

1<br />

for A; vi så, at Au 1 = u 1 , Au 2 =0, 92u 2 . (Vi kan naturligvis også finde egenværdierne fra<br />

det karakteristiske polynomium<br />

» –<br />

0, 94 − λ 0, 02<br />

p A(λ) =<br />

= λ 2 − 1, 92λ +0, 92 = (λ − 1)(λ − 0, 92),<br />

0, 06 0, 98 − λ<br />

og egenvektorer ved en nulrums- (egenrums-)udregning<br />

f.eks.<br />

Vi har da<br />

Da<br />

E 1 = N A−I = N 2 −0.06<br />

4<br />

0.02<br />

0.06 −0.02<br />

3<br />

5<br />

]<br />

;<br />

[ 1<br />

= R · )<br />

3]<br />

A n (au 1 + bu 2 )=aA n u 1 + bA n u 2 = au 1 + b(0, 92) n u 2 .<br />

(det er her, der skiftes koordinater!) er<br />

x 0 =0, 25u 1 − 0, 05u 2 (1)<br />

x n = A n x 0 =0, 25u 1 − 0, 05(0, 92) n u 2 → 0, 25u 1 når n →∞. (2)<br />

Da u 1 , u 2 er uafhængige, er A diagonaliserbar, og, hvis U =[u 1 , u 2 ], så er<br />

[ ]<br />

[ ]<br />

A n 1 0<br />

= U<br />

0 (0, 92) n U −1 , så A n 1 0<br />

x 0 = U<br />

0 (0, 92) n U −1 x 0 .<br />

og<br />

Da U −1 = K U,E er U −1 x 0 =[x 0 ] U =<br />

Da U =<br />

[ ] 0, 25<br />

−0, 05<br />

[<br />

]<br />

A n 0, 25<br />

x 0 = U<br />

−(0, 92) n .<br />

0, 05<br />

[ ] 1 −1<br />

udregnes A n x<br />

3 1<br />

0 =<br />

men formuleringen i (2) er måske nemmere at overskue.<br />

(fra (1))<br />

[ 0, 25 + (0, 92) n 0, 05<br />

0, 75 − (0, 92) n 0, 05<br />

]<br />

165


SEKTION 9.2<br />

DIAGONALISERING<br />

Uafhængigheden af egenvektorer svarende til forskellige egenværdier gælder generelt:<br />

Lemma 9.2.6 ([L], 6.3.1)<br />

Lad A ∈ Mat n,n (F); antag, at λ 1 , . . . , λ k er forskellige egenværdier for A.<br />

Lad v 1 , . . . , v k være tilsvarende egenvektorer. Så er v 1 , . . . , v k uafhængige.<br />

Bevis<br />

Lad S = Span(v 1 , . . . , v k ); og lad r = dim S. Vi må vise, at r = k. Antag, modsætningsvis, at<br />

r < k. Efter omnummerering kan vi antage, at v 1 , . . . , v r udgør en basis for S, så der findes<br />

Så er<br />

dvs.<br />

c 1 , . . . , c r ∈ F med v r+1 = c 1 v 1 + · · · + c r v r (1)<br />

Av r+1 = A(c 1 v 1 + · · · + c r v r )=c 1 Av 1 + · · · + c r Av r ,<br />

Træk λ r+1 gange ligning (1) fra ligning (2):<br />

λ r+1 v r+1 = c 1 λ 1 v 1 + · · · + c r λ r v r (2)<br />

0=c 1 (λ 1 − λ r+1 )v 1 + · · · + c r (λ r − λr + 1)v r .<br />

Da v 1 , . . . , v r er uafhængige er koefficienterne c 1 (λ 1 − λ r+1 ), . . . , c r (λ r − λ r+1 ) i denne lineære<br />

relation alle 0. Da λ i ≠ λ j for i ≠ j må dette betyde, at c 1 , . . . , c r alle er 0. Ligning (1)<br />

giver så, at v r+1 =0. Modstrid – idet v r+1 er en egenvektor. Så vores antagelse var forkert,<br />

og r = k; v 1 , . . . , v k er uafhængige.<br />

Vi konkluderer, at mange matricer er diagonaliserbare:<br />

Korollar 9.2.7<br />

Lad A ∈ Mat n,n (F); antag, at A har n forskellige egenværdier.<br />

Så er A diagonaliserbar.<br />

166


SEKTION 9.2<br />

DIAGONALISERING<br />

Vi kan vise mere omkring uafhængigheden af egenvektorer:<br />

Korollar 9.2.8<br />

Lad A ∈ Mat n,n (F); antag, at λ 1 , . . . , λ k er forskellige egenværdier for A. For i =1, . . . , k,<br />

lad v ij ,j=1, . . . , m i være uafhængige egenvektorer for A svarende til λ i .<br />

Så er {v ij :1≤ i ≤ k, 1 ≤ j ≤ m i } uafhængige.<br />

Bevis<br />

Betragt en lineær relation<br />

k∑ ∑m i<br />

c ij v ij =0 (1)<br />

i=1 j=1<br />

med c ij ∈ F. Skriv x i = ∑ m i<br />

j=1 c ijv ij ; så relationen bliver til<br />

x 1 + · · · + x k =0. (2)<br />

Vi påstår, at x i =0for i =1, . . . , k. Antag, modsætningsvis, at nogle af x i ’erne ikke er 0;<br />

efter omnummerering kan vi antage, at x 1 , . . . , x r ≠0, x r+1 , . . . , x k =0.<br />

Så (2) bliver til<br />

x 1 + · · · + x r =0. (3)<br />

For 1 ≤ i ≤ r er x i en egenvektor for A svarende til λ i , da den er ikke 0 og en lineær<br />

kombination af elementer af E A (λ i ), så selv et element af E A (λ i ). Ifølge Lemma 9.2.6 er<br />

x 1 , . . . , x r uafhængige. (3) giver en modstrid. Så vores antagelse var forkert, og x i = 0 for<br />

i =1, . . . , k.<br />

Altså<br />

0 = x i = c i1 v i1 + · · · + c ini v ini for i =1, . . . , k.<br />

Da v i1 , . . . , v ini er uafhængige fås, at<br />

c i1 =0, . . . , c ini =0for i =1, . . . , k;<br />

så<br />

er uafhængige.<br />

v ij , 1 ≤ i ≤ k, 1 ≤ j ≤ n i ,<br />

Korollar 9.2.9<br />

Lad A ∈ Mat n,n (F); antag, at λ 1 , . . . , λ k er forskellige egenværdier for A. Så er<br />

Geo(λ 1 )+· · · + Geo(λ k ) ≤ n.<br />

Bevis<br />

Ifølge Korollar 9.2.8 er foreningen af baser for E A (λ 1 ), . . . , E A (λ k ) en uafhængig mængde i<br />

F n , så har højst n elementer.<br />

167


SEKTION 9.2<br />

DIAGONALISERING<br />

Vi får nu en karakterisering af diagonaliserbare matricer:<br />

Proposition 9.2.10<br />

Lad A ∈ Mat n,n (F), og lad λ 1 , . . . , λ k være de forskellige egenværdier for A.<br />

A er diagonaliserbar ⇔ n = Geo(λ 1 )+· · · + Geo(λ k ).<br />

Bevis<br />

⇒: Lad {v 1 , . . . , v n } være en basis for F n bestående af egenvektorer. Antag, for i =1, . . . , k,<br />

at n i af dem er indeholdt i egenrummet E A (λ i ). Da disse n i elementer er uafhængige,<br />

er Geo(λ i ) = dim E A (λ i ) ≥ n i . Da alle egenvektorer må have én af λ 1 , . . . , λ k som<br />

tilsvarende egenværdi, har vi<br />

n = n 1 + · · · + n k ;<br />

så n ≤ Geo(λ 1 )+· · · + Geo(λ k ). Vi har faktisk lighed, ifølge Korollar 9.2.9.<br />

⇐: Lad V i være en basis for E A (λ i ) for i =1, . . . , k; lad V = V 1 ∪ · · · ∪ V k . Ifølge Korollar<br />

9.2.8 består V af uafhængige elementer.<br />

Da V i har Geo(λ i ) elementer har V Geo(λ i )+· · · + Geo(λ k )=n elementer. Så V er en<br />

basis for F n , en basis bestående af egenvektorer for A.<br />

Med tanke på beregningsteknikken for egenværdier og egenvektorer er det nærliggende at inddrage<br />

de algebraiske multipliciteter.<br />

Lemma 9.2.11<br />

Lad A ∈ Mat n,n (F); og lad λ 0 være en egenværdi for A. Så er<br />

Geo A (λ 0 ) ≤ Alg A (λ 0 ).<br />

Bevis<br />

Lad v 1 , . . . , v k være en basis for E A (λ 0 ); så Geo(λ 0 )=k. Udvid til en basis v 1 , . . . , v n for<br />

F n . Lad V =[v 1 , . . . , v n ] i søjleform. V er invertibel. Vi har, for i =1, . . . , k.<br />

(V −1 AV )e i = V −1 Av i = V −1 λ 0 v i = λ 0 V −1 V e i = λ 0 e i .<br />

I blokform har vi da<br />

V −1 AV =<br />

k<br />

n − k<br />

k<br />

[<br />

λ0 I B<br />

0 C<br />

n− k<br />

]<br />

.<br />

168


SEKTION 9.2<br />

DIAGONALISERING<br />

Bevis, fortsat<br />

Da A, V −1 AV er similære er<br />

p A (λ) =p V −1 AV (λ)<br />

[ (λ0 − λ)I B<br />

= det<br />

0 C − λI<br />

Ved at udvikle determinanten og de efterfølgende minorer k gange efter første søjle fås<br />

så Alg(λ 0 ) ≥ k = Geo(λ 0 ).<br />

p A (λ) =(λ 0 − λ) k det(C − λI),<br />

]<br />

.<br />

Det giver en alternativ karakterisering af diagonaliserbare matricer:<br />

Sætning 9.2.12<br />

Lad A ∈ Mat n,n (F); lad λ 1 , . . . λ k ∈ F være de forskellige egenværdier for A. Så er A diagonaliserbar<br />

⇔<br />

(1) Alg(λ 1 )+· · · + Alg(λ k )=n<br />

(2) Alg(λ i ) = Geo(λ i ) for i =1, . . . , k.<br />

Bemærkning<br />

Hvis F = C gælder (1) altid.<br />

Bevis<br />

⇐: hvis (1) og (2) gælder, så er Geo(λ 1 )+· · · + Geo(λ n )=n og A er diagonaliserbar ifølge<br />

Proposition 9.2.10.<br />

⇒: hvis A er diagonaliserbar er Geo(λ 1 )+· · · + Geo(λ n )=n. Da p A er et polynomium af<br />

grad n har det højst n rødder, talt med multiplicitet, så Alg(λ 1 )+· · · + Alg(λ k ) ≤ n.<br />

Ifølge Lemma 9.2.12 er<br />

Alg(λ i ) ≥ Geo(λ i ) for i =1, . . . , k<br />

(∗)<br />

så n ≥ Alg(λ 1 )+· · · + Alg(λ k ) ≥ Geo(λ 1 )+· · · + Geo(λ k )=n. Vi må altså have ligheder<br />

hele vejen i denne kæde af uligheder; så n = Alg(λ 1 )+· · · + Alg(λ k ), mens ligheden<br />

Alg(λ 1 )+· · · + Alg(λ n ) = Geo(λ 1 )+· · · + Geo(λ k )<br />

sammen med (∗) implicerer, at Alg(λ i ) = Geo(λ i ) for i =1, . . . , k.<br />

169


SEKTION 9.2<br />

DIAGONALISERING<br />

Vi får en generel fremgangsmåde til at undersøge, om en kvadratisk matrix A ∈ Mat n,n (F) er<br />

diagonaliserbar, og til i givet fald at udføre diagonaliseringen:<br />

1. udregn det karakteristiske polynomium p A (λ) = det(A − λI) for A,<br />

2. find de forskellige rødder λ 1 , . . . , λ k ∈ F for p A (λ),<br />

3. Løs for hvert i ≤ k det homogene ligningssystem (A − λ i I)v = 0, og find derved en basis<br />

{v i1 , . . . , v i,di }, bestående af d i vektorer, for N(A − λ i I)=E A (λ i ).<br />

Lad n i være multipliciteten af λ i som rod af p A .<br />

Hvis n 1 + · · · + n k


SEKTION 9.2<br />

DIAGONALISERING<br />

Eksempel 9.2.14, fortsat<br />

Vi har<br />

⎡<br />

⎤<br />

16 − 20λ 2 1 1<br />

p A (λ) = 1<br />

20 4 det ⎢ 2 16 − 20λ 1 1<br />

⎥<br />

⎣ 1 1 16 − 20λ 2 ⎦ (R 4 → R 1 + R 2 + R 3 + R 4 )<br />

20 − 20λ 20 − 20λ 20 − 20λ 20 − 20λ<br />

⎡<br />

⎤<br />

16 − 20λ 2 1 1<br />

= 1 − λ<br />

20 3 det ⎢ 2 16 − 20λ 1 1<br />

⎥<br />

⎣ 1 1 16 − 20λ 2⎦<br />

1 1 1 1<br />

⎡<br />

⎤<br />

15 − 20λ 1 0 1<br />

= 1 − λ<br />

20 3 det ⎢ 1 15 − 20λ 0 1<br />

⎥<br />

⎣ −1 −1 14 − 20λ 2⎦ (S i → S i − S 4 ,i=1, 2, 3)<br />

0 0 0 1<br />

⎡<br />

⎤<br />

= 1 − λ 15 − 20λ 1 0<br />

[ ]<br />

20 3 det ⎣ 1 15 − 20λ 0 ⎦<br />

(1 − λ)(14 − 20λ) 15 − 20λ 1<br />

=<br />

20<br />

−1 −1 14 − 20λ<br />

3 det<br />

1 15 − 20λ<br />

= 1<br />

20 3 (1 − λ)(14 − 20λ)((15 − 20λ)2 − 1 2 ) = 1<br />

20 3 (1 − λ)(14 − 20λ)2 (16 − 20λ)<br />

= (1 − λ)( 7<br />

10 − λ)2 ( 8<br />

10 − λ).<br />

Vi beregner E A ( 7<br />

10 ): ⎡ ⎤ ⎡ ⎤<br />

2 2 1 1 1 1 2 2<br />

A − 7<br />

10 = 1 ⎢2 2 1 1<br />

⎥<br />

20 ⎣1 1 2 2⎦ ∼ ⎢2 2 1 1<br />

⎥<br />

⎣0 0 0 0⎦<br />

1 1 2 2 0 0 0 0<br />

⎡<br />

⎤ ⎡ ⎤<br />

1 1 2 2 1 1 0 0<br />

∼ ⎢0 0 −3 −3<br />

⎥<br />

⎣0 0 0 0 ⎦ ∼ ⎢0 0 1 1<br />

⎥<br />

⎣0 0 0 0⎦<br />

0 0 0 0 0 0 0 0<br />

Så<br />

På lignende vis fås<br />

( ) ( 7<br />

E A = N A − 7<br />

10<br />

10<br />

⎛⎡<br />

⎤ ⎡<br />

)<br />

−1<br />

= Span ⎜⎢<br />

1<br />

⎥<br />

⎝⎣<br />

0 ⎦ , ⎢<br />

⎣<br />

0<br />

0<br />

0<br />

−1<br />

1<br />

⎛⎡<br />

⎤⎞<br />

1<br />

E A (1) = N(A − ) = Span ⎜⎢1<br />

⎥⎟<br />

⎝⎣1⎦⎠ ,<br />

1<br />

⎛⎡<br />

⎤⎞<br />

−1<br />

E A ( 8<br />

10 )=N(A − 8 10 ) = Span ⎜⎢−1<br />

⎥⎟<br />

⎝⎣<br />

1 ⎦⎠ .<br />

1<br />

⎤⎞<br />

⎥⎟<br />

⎦⎠ .<br />

171


SEKTION 9.2<br />

DIAGONALISERING<br />

Eksempel 9.2.14, fortsat<br />

A er diagonaliserbar,<br />

hvor<br />

⎡<br />

1<br />

V −1 AV = ⎢<br />

⎣<br />

8/10<br />

7/10<br />

⎤<br />

⎥<br />

⎦<br />

7/10<br />

⎡<br />

⎤<br />

⎡<br />

⎤<br />

1 −1 −1 0<br />

1 1 1 1<br />

V = ⎢1 −1 1 0<br />

⎥<br />

⎣1 1 0 −1⎦ og V −1 = 1 ⎢−1 −1 1 1<br />

⎥<br />

4 ⎣−2 2 0 0⎦ .<br />

1 1 0 1<br />

0 0 −2 2<br />

172

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!