4.4 Koordinatsystemer

4.4 Koordinatsystemer
Definisjon: B = b1 , b2 , . . . , bn ⊂ V kalles en basis for et
vektorrom V dersom B er lineært uavhengig og utspenner V .
TEOREM 7: (Teoremet om entydig representasjon) Anta at
B = {b1 , b2 , . . . , bn } er en basis for vektorrommet V . For hver
x ∈ V fins da entydig bestemte skalarer c1 , c2 , . . . , cn slik at
x = c1 b1 + c2 b2 + · · · + cn bn .
I
I
I
Vektene c1 , c2 , . . . , cn kalles koordinatene til x med hensyn på
(eller, relativt til) B.
Vektoren [x]B = (c1 , c2 , . . . , cn ) i Rn kalles koordinatvektoren
til x m.h.p. B.
Avbildningen x → [x]B fra V til Rn kalles koordinatavbildningen m.h.p. B.
MERK: Vektorene i B betraktes som ordnet (”ordnet basis”);
bytter vi rekkefølgen vil også koordinatvektoren ombyttes
tilsvarende.
1 / 15
Eksempler.
1) Anta x ∈ Rn og B = {e1 , . . . , en }. Da er [x]B = x .
Så vektorer i Rn er i utgangspunktet koordinatisert m.h.p.
standardbasisen
2) Betrakt b1 = (1, −1), b2 = (2, 1), B = {b1 , b2 } og x = (3, 3).
Da er
[x]B = (−1 , 2).
3) Betrakt P2 og B = {1, t, t 2 }.
Hvis p(t) = a0 + a1 t + a2 t 2 , så er [p]B = (a0 , a1 , a2 ) .
F.eks. er [3 − t + 2 t 2 ]B = (3, −1, 2).
2 / 15
Spesielt om koordinater i Rn
La B = {b1 , b2 , . . . , bn }h være en basisi for Rn og la x ∈ Rn . Da er
x = PB [x]B der PB = b1 b2 . . . bn .
I
Matrisen PB = b1 b2 . . . bn kalles koordinatskiftematrisen
(eller basisskiftematrisen) fra B til standardbasisen i Rn . Den
er invertibel, og [x]B = (PB )−1 x for alle x ∈ Rn .
I
Gitt x ∈ Rn kan vi bestemme koordinatvektoren c = [x]B ved
å løse likningen
PB c = x
med hensyn på c.
I
Dette viser at koordinatavbildningen (m.h.p. B) er lineær, at
dens standardmatrise er (PB )−1 , og at den er 1 − 1 og på Rn .
3 / 15
TEOREM 8: Anta at B = {b1 , b2 , . . . , bn } er en (ordnet) basis for
et vektorrom V . Da er x → [x]B en lineær avbildning fra V til Rn
som er 1-1 og på.
Definisjon: En isomorfi fra et vektorrom V til et vektorrom W er
en lineær avbildning fra V til W som 1 − 1 og på W .
Teorem 8 sier altså at koordinatavbildningen m.h.p. en (ordnet)
basis for V er en isomorfi fra V til Rn .
Korollar til Teorem 8 [fra Notat 1]:
Anta at B = b1 , b2 , . . . , bn er en (ordnet) basis for V .
Betrakt en delmengde S = {u1 , u2 , . . . , up } av V .
Definer
o
n
SB = [u1 ]B , [u2 ]B , . . . , [up ]B
(som er en delmengde av Rn ) og [SB ] ∈ Mn×p ved
h
i
[SB ] = [u1 ]B [u2 ]B · · · [up ]B .
4 / 15
Da har vi at:
I
S utspenner V ⇔ SB utspenner Rn .
I
S er lineært uavhengig i V ⇔ SB er lineært uavhengig i Rn .
I
S er en basis for V ⇔ SB er en basis for Rn
⇔ p = n og matrisen [SB ] er invertibel.
Kommentarer:
I
Viktig konsekvens: Dersom et vektorrom V har en basis med
n elementer, så er antall vektorer i en hvilken som helst annen
basis for V også lik n.
Tallet n kalles dimensjonen til V , se neste avsnitt.
I
Når S er en basis for V , kalles den invertible matrisen [SB ] for
koordinatskiftematrisen fra S til B.
Slike matriser skal vi studere nærmere i avsnitt 4.7.
5 / 15
Nok en kommentar
La R = rref([SB ]). Husk at
I
SB utspenner Rn ⇔ R har pivoter i alle sine rader.
I
SB er lineært uavhengig ⇔ R har pivoter i alle sine kolonner.
I
SB er en basis for Rn ⇔ p = n og R = In .
I praksis bruker vi korollaret til Teorem 8 ved å bestemme SB , lage
matrisen [SB ] og deretter finne R = rref([SB ]).
6 / 15
Eksempel. Vi illustrerer dette med P2 , B = {1, t, t 2 } ,
S = {1 − t 2 , t − t 2 , 2 − t − t 2 } Da er
SB = [1 − t 2 ]B , [t − t 2 ]B , [2 − t − t 2 ]B

 
 

1
0
2
o
n
=  0  ,  1  ,  −1  . Så
−1
−1
−1


1
0
2
1 −1 .
[SB ] =  0
−1 −1 −1

Utregning gir

1 0 2
R = rref([SB ]) =  0 1 −1  . Vi ser at SB
0 0 0
ikke utspenner R3 og at SB er lineært avhengig. Dermed kan vi
konkludere med at S ikke utspenner P2 og at S er lineært avhengig.
7 / 15
4.5 Dimensjon
Vi lar V betegne et vektorrom.
TEOREM 9 og 10. Anta at V har en basis B med n vektorer.
La S være en endelig delmengde av V .
1. Anta S har flere enn n vektorer. Da er S lineært avhengig.
2. Anta S har færre enn n vektorer. Da er V ikke utspent av S.
3. Enhver basis for V har også n elementer.
Definisjon: Et vektorrom V kalles endeligdimensjonalt dersom det
er utspent av en endelig mengde. Hvis ikke, sier vi at V er
uendeligdimensjonalt.
I
I
I
Anta V 6= {0} er endeligdimensjonalt. Da har V en basis (ved
Teorem 5). Antall elementer i en hvilken som helst basis for V
kalles dimensjonen til V , og betegnes med dim V .
Vi setter dim{0} = 0.
Eksempler: dim Pn = n + 1, dim Mm×n = m n
8 / 15
Om dim Col A og dim Nul A
La A være en m × n matrise. Sett R = rref(A). Da gjelder:
I
Kolonnevektorene i A som svarer til pivotene i R danner en
basis for Col A.
Så dim Col A = antall pivoter i R.
I
Vektorene vi får ved å ”separere” de frie variablene når vi skal
angi Nul A på parameterform danner en basis for Nul A.
Så dim Nul A = antall frie variable i systemet A x = 0
= antall kolonner i R som ikke er pivotkolonner.
Merk: dim Col A kalles rangen til A og er tema for avsn. 4.6.
9 / 15
To nyttige teoremer, der V er endeligdimensjonalt, med
dim V = n ≥ 1 .
TEOREM 11: La H være et underrom av V .
Da kan enhver lineært uavhengig delmengde av H utvides (om
nødvendig) til en basis for H.
Videre er H endeligdimensjonalt og dim H ≤ dim V .
TEOREM 12 (”Basis teoremet”). La S være en endelig delmengde
av V , og anta at S består av n vektorer.
I
Hvis S er lineært uavhengig, så er S en basis for V .
I
Hvis S utspenner V , så er S en basis for V .
Eksempel: Betrakt V = Span{sin2 t, cos2 t} og S = {1, cos 2t}.
Siden dimV = 2 og S består av 2 lineært uavhengige vektorer i V ,
så er S er en basis for V .
10 / 15
4.9 Anvendelser: Markovkjeder
Markov kjeder er en spesiell type diskret dynamisk system. Brukes
bl.a. i finans, fysikk, biologi, økonomi, demografiske modeller
(populasjonsmodeller). Er også grunnlaget for rangeringen i
Google, algoritmen kalles Pagerank. Man kan lære mer om
Markovkjeder i STK2130 - Modellering av stokastiske prosesser
Eksempel: Populasjonsmodell.
I
I
I
I
I
Land inndelt i n områder.
Tilstandsvektor x ∈ IRn der xi angir andelen somPbor i område
i (i ≤ n). Da er komponentene ikkenegative og j xj = 1. En
slik vektor kalles ofte en stokastisk vektor.
Antar at hvert år (eller ved angitt tidspunkt) vil en viss andel
flytte fra en region til en annen region. Dette beskrives ved
andeler (sannsynligheter).
pij er andelen av befolkningen i region j som flytter til region i
i løpet av ett år.
Vi antar at denne andelen er den samme for hvert år.
11 / 15
En matrise P = [pij ] kalles (kolonne)stokastisk dersom den er
komponentvis ikkenegativ og hver kolonnesum er lik 1, dvs. hver
kolonne er en stokastisk vektor. Brukes i sannsynlighetsregning:
Markov kjede, pij er da overgangssannsynlighet fra tilstand j til
tilstand i. Slike matriser kalles også Markovmatriser. Eksempel:


1 0.2
0
0 0
 0 0.5 0.2
0 0 



P2 =  0 0.3 0.5 0.2 0 

 0
0 0.3 0.5 0 
0
0
0 0.3 1
Tilbake til populasjonsmodell:
Hvis xk er tilstandsvektor ved starten av tidsperiode k, så blir
tilstandsvektor ved neste periode k + 1 gitt ved
xk+1 = Pxk
12 / 15
Hva skjer langt fram i tid?
xk = Pxk−1 = P(Pxk−2 ) = P 2 xk−2 = · · · = P k x0
Så den asymptotiske fordelingen blir bestemt av matrisen P k .
Elementet i posisjon (i, j) i denne matrisen P k kan tolkes som
sannsynligheten for at en person flytter fra region j til region i i
løpet av k tidsperioder. P k er også en stokastisk matrise.
Oppgave: vis at produktet av to stokastiske matriser er en
stokastisk matrise.
I Matlab skriver man bare x=P∗x i en for-løkke. I eksemplet foran
”får vi konvergens” etter 5 iterasjoner og x = (0.5714, 0.4286).
13 / 15
Definisjon. Hvis P er en stokastisk matrise og q er en stokastisk
vektor slik at
q = Pq
så kalles q en likevektsvektor (equilibrium, steady-state, stasjonær
fordeling).
Teorem: Enhver stokastisk matrise har en likevektsvektor.
Skal senere knytte dette til egenvektorer (som behandles i neste
kapittel).
Vi kan finne en likevektsvektor ved å løse det lineære
likningssystemet.
(P − I )x = O
14 / 15
Numerisk (for store matriser) finner man gjerne likevektsvektor ved
iterasjonen x=P∗x. Men det er da viktig å utnytte sparsity i
matrisen. Samme type beregning ligger bak Google og PageRank
algoritmen!
Definisjon. En stokastisk matrise P kalles regulær dersom, for en
viss k, er alle elementene i P k positive.
I så fall vil alle elementene være positive også i P k+1 , P k+2 , . . ..
Dette svarer til at man kan komme seg fra enhver tilstand til
enhver annen tilstand ved passende mange tilstandsoverganger.
F.eks. hvis pij > 0 for alle i, j, så er P regulær.
TEOREM 18. Hvis P er en regulær stokastisk matrise, så har P en
entydig likevektsvektor q. Videre, hvis x0 er en vilkårlig startvektor
og xk+1 = Pxk for k = 0, 1, . . ., så vil Markovkjeden {xk }
konvergere mot q når k → ∞.
15 / 15