Module 10 Lineaire Algebra - Wisnet

L
Module 10 Lineaire Algebra
Afbeeldingen (vervolg (b))
In deze les worden de eigenwaarden en eigenvectoren van lineaire afbeeldingen behandeld.
Inhoud van de leskern
1 Basistransformatie
Vak 57.52 Les 36.0
*0000==:;000=*

Leskern 57.52-36.0 2
1 Basistransformatie
1.1 Kentallen (coördinaten) ten opzichte van een basis
Zoals u weet vormt het stelsel fe1;e2g een basis voor de R2 : Als er
verder niets vermeld wordt, bedoelt men ook deze basis.
Iedere vector in de R2 kan uitgedrukt worden in deze basisvectoren.
Bijvoorbeeld: de vector x = ¡ 3¢
2 =3e1 +2e2: De getallen 3 en 2 zijn
de kentallen van de vector x ten opzichte van de basis fe1;e2g :
Het is echter ook mogelijk om in de R2 twee andere onafhankelijke
vectoren te kiezen die de R2 opspannen. Deze onafhankelijke vectoren
hoeven in het geheel niet loodrecht op elkaar te staan en het hoeft ook
niet zo te zijn dat ze lengte 1 hebben. Men kan immers elke vector van
de R2 schrijven als een combinatie van zo’n onafhankelijk stelsel
vectoren.
Voorbeeld 1
Neem de twee onafhankelijke vectoren b 1 en b 2 als nieuwe basis voor
de R 2 .
b 1 =
µ 1
1
en b 2 =
µ ¡1
1
Het is mogelijk om iedere vector in deze basisvectoren uit te drukken.
Druk bijvoorbeeld de vector x = ¡ ¢ 3
2 uit in deze basisvectoren, dat wil
zeggen
= x1 e1 + x2 e2 = y1 b1 + y2 b2 µ x1
x2
en bereken de kentallen y1 en y2:
- Men noemt x1 en x2 de kentallen van x ten opzichte van fe 1 ;e 2 g :
- Men noemt y1 en y2 de kentallen van dezelfde vector maar dan ten
opzichte van de nieuwe basis fb 1;b 2g :
Zie ook afbeelding 1.
Afbeelding 1
Antwoord:
b2
x2-as
-1/2 b2
b1
5/2 b1
x1-as
x

µ y1
y2
57.52-36.0 3
De volgende vectorvergelijking moet opgelost worden waarin y1 en y2
de onbekenden zijn.
µ
µ µ
3
1 ¡1 y1
= y1 b
2
1 + y2 b2 =
=)
1 1 y2
µ ¡1 µ µ µ
12 1
1 ¡1 3
=
=
2 3
=
1 1 2
2
De matrix waarvan de kolommen gevormd worden door de nieuwe
basisvectoren wordt B genoemd. Voor de berekening van de inverse
matrix B ¡1 kunt u Maple gebruiken
> with(linalg):
> B:=matrix([[1,-1],[1,1]])
" #
1 ¡1
B :=
1 1
> inverse(B)
2
6
4
1
2
¡1
2
> [y1,y2]=evalm(inverse(B) &* vector([3,2]))
· ¸
5 ¡1
[y1;y2] = ;
2 2
1
2
1
2
3
7
5
¡ 1 2
1
2
µ 52
Het komt er dus op neer dat u de vectoren van de nieuwe basis b 1 en b 2
in de kolommen van een matrix zet die u B noemt. Elke vector x 2 R 2
kanuitgedruktwordenindebasisvectorenb 1 en b 2: Uhoeftslechtsde
inverse van B te kennen en de nieuwe kentallen y1 en y2 van de
desbetreffende vector x zijn:
µ y1
y2
= B ¡1
µ
x1
x2
Voorbeeld 2
Bereken de kentallen van e 1 en e 2 ten opzichte van de basis fb 1 ;b 2 g
van voorbeeld 1.
Antwoord:
De nieuwe kentallen worden
B ¡1
µ
1
0
=
µ
12
¡
1
2
1
2
1
2
e1 = 1
2 b1 ¡ 1
2 b2 B ¡1
µ
0
1
=
µ
1
2
1
2
¡ 1
2
1
2
e 2 = 1
2 b 1 + 1
2 b 2
µ 1
0
µ 0
1
=
=
µ 12
µ 1
2
1
2
¡ 1
2
=)
=)
¡ 1 2

57.52-36.0 4
Ga op afbeelding 1 na dat de kentallen van de vectoren e 1 en e 2 ten
opzichte van de basis fb 1;b 2g kloppen.
- Een vector 2 R n met n kentallen ten opzichte van fe 1;e 2; :::::; e ng
kan overgevoerd worden naar dezelfde vector y 2 R n maar nu met
kentallen ten opzichte van een nieuwe basis fb 1;b 2; :::::; b ng door
middel van y = B ¡1 x waarbij B de matrix is waarvan de
kolommen gevormd worden door de vectoren van deze nieuwe
basis.
UzultbegrijpendateenbasisindeR n bestaat uit n vectoren met n
kentallen, dus de matrix B is vierkant. Deze n vectoren moeten
onafhankelijk zijn. Dit houdt in dat de matrix B altijd inverteerbaar is.
Opdracht 1
Bereken de kentallen van de vector x =(1; 2; 3) ten opzichte van de
basis:
(Ã !
1
à !
0
à !)
2
fb1 ;b2 ;b3g = 2
¡3
; 0
¡2
; 6
0
1.2 Afbeeldingen ten opzichte van een basis
U weet dat in de kolommen van een afbeeldingsmatrix van een
afbeelding van R 2 naar R 2 de beelden staan van de basisvectoren e 1 en
e 2: Dat wil zeggen dat de afbeelding beschreven is ten opzichte van
dezebasis.Meestalzegtmendaternietbijalshetomde
standaard-basis gaat.
Voorbeeld 3
Neemt u de afbeelding van R 2 naar R 2 vandespiegelingindelijn
x1 = x2: Zie ook afbeelding 2.
Afbeelding 2
b2
x2-as
b1
x1-as
De matrix (ten opzichte van fe1;e2g) van deze afbeelding S (spiegeling
in de lijn x1 = x2) is
µ
0 1
S =
1 0

57.52-36.0 5
en heeft als kolommen de beelden van e1 en e2: Het zou eigenlijk
mooier zijn om niet fe1;e2g als basis te nemen voor juist deze
afbeelding S maar over te gaan op een nieuwe basis fb1;b2g met
b 1 =
µ 1
1
en b 2 =
µ ¡1
1
Op afbeelding 2 zult u kunnen zien dat deze basis voor juist deze
afbeelding zéér geschikt is, want bij deze spiegeling worden de beelden
van deze basisvectoren op zichzelf of op hun eigen lijn afgebeeld!! De
matrix (ten opzichte van fb1;b2g) van déze spiegeling moet als
kolommen de beelden van b1 en b2 hebben! Bij deze afbeelding geldt
dus dat b1 op zijn plaats blijft en dat b2 alleen maar omgeklapt wordt.
De beelden van b1 en b2 zijn dus Sb1 = b1 en Sb2 = ¡b2: De
vectoren b1 en b2 zijn nu juist de eigenvectoren van deze spiegeling!
Als we deze beeldvectoren ook nog met de kentallen willen schrijven
ten opzichte van de nieuwe basis fb1;b2g dan worden de beelden van
deze basisvectoren na spiegelen in de lijn x1 = x2 dus
B ¡1 µ
1
Sb1 = en B
0
¡1 µ
0
Sb2 =
¡1
Let op!, de betekenis van ¡ 1¢
0 ten opzichte van de nieuwe basis fb1;b2g wil zeggen ¡ ¢ 1
0 =1¢ b1 +0¢ b2: Vormen we nu met de beelden van b1 en b2 de kolommen van een
matrix waarbij ook nog de kentallen ten opzichte van de nieuwe basis
fb1;b2g gegeven zijn, dan krijgen we de matrix die de genoemde
spiegeling beschrijft met als basis het stelsel fb1 ;b2g : Hetmoetu
opvallen dat het resultaat een diagonaal-matrix is.
B ¡1 µ
1 0
SB=
0 ¡1
- Algemeen geldt nu het volgende voor een afbeelding A (ten
opzichte van de normale basis). Zet men de beelden van een stelsel
nieuwe basisvectoren (met kentallen ten opzichte van deze
basisvectoren) in de kolommen van de afbeeldingsmatrix, dan
krijgt men de matrix B ¡1 AB:(De kolommen van matrix B
worden gevormd door de nieuwe basisvectoren.) Deze
samengestelde matrix stelt weer precies dezelfde afbeelding voor.
Men moet zich echter bewust zijn van het feit dat deze afbeelding
niet meer ten opzichte van de gewone basis beschreven wordt, maar
ten opzichte van een nieuwe basis die natuurlijk wel erbij vermeld
dient te worden.
Opdracht 2
Gegeven is een afbeeldingsmatrix A.
µ
5 ¡2
A =
6 ¡2
a. Bereken de matrix van deze afbeelding, maar nu ten opzichte van
de basis ©¡ 1¢
¡ ¡1¢ª
1 ; 1 waarvan de matrix B1 genoemd wordt.
b. Bereken de matrix van deze afbeelding, maar nu ten opzichte van
de basis ©¡ 1¢
¡ 2¢ª
2 ; 3 waarvan de matrix B2 genoemd wordt.

jordan
57.52-36.0 6
c. Bereken de eigenwaarden van afbeeldingsmatrix A.
d. Bereken de eigenwaarden van de afbeeldingsmatrix B1 AB1 ten
opzichte van basis B1 en ook van de afbeeldingsmatrix B2 AB2
ten opzichte van basis B2:
e. Bereken de eigenvectoren van A uitgedrukt in de basis fe 1;e 2g :
f. Bereken de eigenvectoren van de afbeelding B1 AB1 ten opzichte
van basis B1 en van de afbeelding B2 AB2 ten opzichte van basis
B2: Let op!, deze eigenvectoren zijn dan ook uitgedrukt in de
vectoren van de nieuwe basis.
g. Controleer dat de kolommen van de matrix B ¡1
1 AB1 gevormd
worden door de beelden van de bijbehorende basisvectoren,
uitgedrukt in deze basisvectoren.
1.3 Diagonaalmatrix
Het is niet voor niets geweest dat we in voorbeeld 3 de nieuwe
basisvectoren b1 en b2 hadden genomen! De spiegelingsmatrix
B ¡1SB(ten opzichte van fb1;b2g) wordtnuheelmooieen
diagonaalmatrix.
Misschien is u al opgevallen dat b1 en b2 juist de eigenvectoren waren
van de spiegelingsmatrix S:
Er geldt immers Sb1 = b1 en dus is b1 een eigenvector met
eigenwaarde 1.
Er geldt immers Sb2 = ¡b2 en dus is b2 een eigenvector met
eigenwaarde ¡1.
In de diagonaalmatrix staan heel netjes de eigenwaarden op de
diagonaalplaatsen.
Ook in opdracht 2 ziet u dat de basis van eigenvectoren de
diagonaalmatrix (B2) oplevert bij basistransformatie naar de basis van
eigenvectoren.
- De matrix van S (van voorbeeld 3) ten opzichte van fe1 ;e2g stelt
dezelfde afbeelding voor als de matrix B ¡1SB(met dezelfde
eigenvectoren en dezelfde eigenwaarden), maar de laatste is ten
opzichte van de basis fb1;b2g.Hierinvormendevectorenb1en b2 de kolommen van matrix B (uitgedrukt in deze basisvectoren).
- Voor iedere afbeelding A van Rn naar Rn kan men in Rn een
nieuwe basis kiezen (van n onafhankelijke vectoren) die Rn opspant (het domein van de afbeelding). De vectoren van deze
nieuwe basis vormen de kolommen van de vierkante matrix B. De
afbeelding wordt beschreven door matrix A (ten opzichte van
fe1; ::::::eng),maarkanookbeschrevenwordendoordematrix B ¡1ABten opzichte van fb1; ::::::bng :
U kunt zelf kiezen voor een willekeurige basis, maar de basis die
gevormd wordt door de eigenvectoren vandeafbeeldingishetmooist.
- De afbeeldingsmatrix ten opzichte van de basis van eigenvectoren
wordt dan de diagonaalmatrix met op de diagonaalplaatsen de
eigenwaarden.
Maple heeft een gemakkelijk commando om van een afbeeldingsmatrix
A van de R3 naar de R3 de diagonaalmatrix B ¡1AB te maken en
metéén ook de matrix B van eigenvectoren te geven. Dit commando kan
echter alléén gebruikt worden als u er zeker van bent dat u met een
stelsel onafhankelijke eigenvectoren te maken hebt. Zo niet, dan krijgt u
met dit commando ook geen diagonaalmatrix. Het commando luidt
jordan en hier volgt metéén een voorbeeld.

57.52-36.0 7
Voorbeeld 4
Gegeven de afbeeldingsmatrix van de spiegeling van voorbeeld 3.
Bepaal de diagonaalmatrix B ¡1 ABen ga na dat in de kolommen van
de matrix B de eigenvectoren van A staan.
Antwoord:
Bij het opgeven van de matrix A kunt u als optie ’B’ meegeven die dan
staat voor de matrix van eigenvectoren. De letter B (deze letter kiest u
zelf) staat tussen quotes, want anders zou Maple deze kunnen opvatten
als een eerder gedenieerde B.
> with(linalg):
> A:=matrix([[0,1],[1,0]])
" #
0 1
A :=
1 0
> jordan(A,’B’)
> evalm(B)
" #
¡1 0
0 1
2
1
6 2
4 ¡1
2
Uzietdatdediagonaalmatrix links boven de eigenwaarde ¸ = ¡1
vermeldt. In de bijbehorende matrix B vindt u als eerste kolom de
eigenvector ( 1 1
2 ; ¡ 2 ): Deze eigenvector noteren wij meestal zonder
breuken als (1; ¡1): Het maakt niet uit welke lengte deze eigenvector
heeft, het moet een representant zijn. De tweede eigenwaarde ¸ =1
van matrix A staat op de tweede diagonaalplaats van de
diagonaalmatrix en de tweede eigenvector staat in de tweede kolom van
matrix B. Deze tweede eigenvector staat als ( 1 2 ; 1 2 ) genoteerd, maar wij
schrijven liever zonder breuken (1; 1): De controle voor de
eigenwaarden en eigenvectoren kunt u ook met eigenvects doen.
1
2
1
2
> eigenvects(A)
[¡1; 1; f[¡1; 1]g]; [1; 1; f[1; 1]g]
Opdracht 3
Onderzoek of er een basis van eigenvectoren te vinden is bij de
volgende matrices en bepaal zo mogelijk de diagonaalmatrix.
P =
à 3 ¡1 ¡1
0 1 0
1 ¡1 1
!
3
7
5
0
1
1 3 ¡ 2 2 0
en Q = @ 3 1
2 ¡ 2 0 A
0 0 3