Käänteismatriisi
Lineaarinen yhtälöryhmä voidaan kirjoittaa kerroinmatriisin A ja
vakiotermivektorin b avulla muodossa
Ax=b. Jos kerroinmatriisi sattuisi olemaan
neliömatriisi ja A=In, niin ryhmän ratkaisu olisi
yksinkertaisesti x=Inx=b. Muussa
tapauksessa voidaan kysyä, löydetäänkö sellainen matriisi B,
että tulo BA=In. Tällöin yhtälön ratkaisuksi kävisi
x=Inx=(BA)x=B(Ax)=Bb.
Tässä päättelyssä on olennaista, että A on
n×n-neliömatriisi.
Huomautus.
Kaikki neliömatriisit eivät ole kääntyviä ja vain
neliömatriisit voivat olla kääntyviä annetun määritelmän mukaan.
Matriisin kääntyvyyden riittävät ja välttämättömät ehdot voidaan
esitellä vasta myöhemmin.
Lause.
Kääntyvän matriisin A käänteismatriisi
A−1 on yksikäsitteinen.
Oletetaan, että matriisilla A on käänteismatriisit
B1 ja B2, eli B1A=In=AB1 ja
B2A=In=AB2. Tällöin
B1=B1In=B1(AB2)=(B1A)B2=InB2=B2.
Siis B1=B2, joten käänteismatriisi on yksikäsitteinen.
◻
Lause.
Olkoot A ja B kääntyviä samankokoisia
neliömatriiseja, sekä r≠0 reaaliluku. Tällöin
- A−1 on kääntyvä ja (A−1)−1=A,
- rA on kääntyvä ja (rA)−1=1rA−1,
- AB on kääntyvä ja (AB)−1=B−1A−1,
- AT on kääntyvä ja (AT)−1=(A−1)T=A−T,
- An, missä n≥0 on kokonaisluku, on kääntyvä ja
(An)−1=(A−1)n=A−n.
Osoitetaan esimerkkinä kohdat 3 ja 5. Kohdan 3
todistamiseksi riittää osoittaa, että matriisi B−1A−1
(kooltaan n×n) toteuttaa käänteismatriisin ehdot, eli
(B−1A−1)(AB)=B−1(A−1A)B=B−1InB=B−1B=In(AB)(B−1A−1)=A(BB−1)A−1=AInA−1=AA−1=In.
Siis matriisi AB on kääntyvä ja
(AB)−1=B−1A−1.
Todistetaan kohta 5 induktiolla
hyödyntäen kohtaa 3.
Alkuaskel. Kun n=0, väitetään matriisin A0=In
olevan kääntyvä ja että käänteismatriisi
I−1n=(A−1)0=In. Koska InIn=In, niin
I−1n=In ja väite on tosi.
Induktioaskel. Tehdään induktio-oletus, jonka mukaan matriisi
Ak on kääntyvä ja (Ak)−1=(A−1)k, kun
k≥0 on kokonaisluku. Pyritään osoittamaan, että matriisi
Ak+1 on kääntyvä ja että
(Ak+1)−1=(A−1)k+1. Induktio-oletuksen ja
kohdan 3 nojalla matriisi AkA=Ak+1 on kääntyvä ja
(Ak+1)−1=(AkA)−1=A−1(Ak)−1io=A−1(A−1)k=(A−1)k+1,
eli induktioväite on tosi.
Induktioperiaatteen nojalla siis matriisi An on kääntyvä ja
(An)−1=(A−1)n aina, kun n≥0 on
kokonaisluku. Loput jätetään harjoitustehtäväksi. ◻
Edellisen lauseen viimeinen kohta motivoi laajentamaan matriisin
potenssin määritelmää myös negatiivisille kokonaisluvuille. Jos
A on kääntyvä matriisi ja k on positiivinen
kokonaisluku, niin asetetaan
A−k=(A−1)k=(Ak)−1.
Vastaavanlaista merkintää käytetään joskus transpoosin
käänteismatriisille
A−T=(AT)−1=(A−1)T.
Esimerkki.
Ratkaise X, kun
- AXA2=A−1,
- AXB=(BA)2,
ja kaikki matriisit ovat kääntyviä n×n-matriiseja.
Koska matriisitulo ei ole vaihdannainen, tulee mainita
kerrotaanko yhtälöä oikealta vai vasemmalta puolelta.
Kerrotaan yhtälöä puolittain vasemmalta matriisilla A−1,
jolloin
AXA2=A−1⇔A−1AXA2=A−1A−1⇔XA2=A−2.
Kerrotaan sitten yhtälöä puolittain oikealta matriisilla
A−2, jolloin
XA2=A−2⇔XA2A−2=A−2A−2⇔X=A−4.
Siis matriisi X=A−4.
Kerrotaan yhtälöä puolittain vasemmalta matriisilla A−1 ja
oikealta matriisilla B−1, jolloin
AXB=(BA)2⇔A−1AXBB−1=A−1(BA)2B−1⇔X=A−1(BA)2B−1.
Koska matriisitulo ei ole vaihdannainen, ratkaisu ei etene
pidemmälle.
Yhtälöryhmän ratkaisu saadaan käänteismatriisin avulla seuraavasti.
Lause.
Jos A on kääntyvä n×n-matriisi, niin
lineaarisella yhtälöryhmällä Ax=b, missä
b on avaruuden Rn vektori, on
yksikäsitteinen ratkaisu x=A−1b avaruudessa
Rn.
Oletetaan, että A on kääntyvä, jolloin
käänteismatriisi A−1 on olemassa. Vektori
x=A−1b on tällöin yhtälöryhmän
Ax=b ratkaisu, sillä
Ax=A(A−1b)=(AA−1)b=Inb=b.
Ratkaisu on myös yksikäsitteinen, sillä jos y on toinen
ratkaisu, niin
y=Iny=(A−1A)y=A−1(Ay)=A−1b=x,
eli ratkaisut ovat samat. ◻
2×2-matriisin käänteismatriisille saadaan seuraava
helpohko, jopa ulkoa muistettavissa oleva sääntö.
Suoralla laskulla voitaisiin osoittaa, että kyseessä on
todella käänteismatriisi. Tutkitaan asiaa kuitenkin hieman
konstruktiivisemmin. Oletetaan, että matriisi
B=[wxyz]
toteuttaa käänteismatriisin määritelmän AB=In=BA. Tällöin
siis
[abcd][wxyz]=[aw+byax+bzcw+dycx+dz]=[1001].
Matriisit ovat samat silloin, kun niiden sarakkeet ovat samat, eli kun
[aw+bycw+dy]=[abcd][wy]=[10]ja[ax+bzcx+dz]=[abcd][xz]=[01].
Saadaan siis kaksi yhtälöryhmää, joiden kokonaismatriisit ovat
[ab1cd0]ja[ab0cd1],
ja joiden on toteuduttava yhtä aikaa. Jaetaan tarkastelu kahteen osaan
sen mukaan, onko kerroin a=0 vai ei.
a≠0. Viedään ensimmäistä kokonaismatriisia kohti
redusoitua riviporrasmuotoa.
[ab1cd0]1aR1→[1ba1acd0]R2−cR1→[1ba1a0d−bca−ca]
Jos tässä d−bca=ad−bca=0, niin
yhtälöryhmällä ei ole ratkaisua. Oletetaan siis, että
ad−bc≠0, jolloin
[1ba1a0d−bca−ca]aad−bcR2→[1ba1a01−cad−bc]R1−baR2→[10dad−bc01−cad−bc].
Samalla prosessilla ja oletuksilla saadaan jälkimmäiselle
kokonaismatriisille
rref[ab0cd1]=[10−bad−bc01aad−bc].
a=0. Jos nyt olisi lisäksi c=0, niin pitäisi
toteutua
[bydy]=[10]ja[bzdz]=[01].
Ensimmäisen yhtälön ensimmäisestä komponentista päätellään, että
b≠0, jolloin yhtälö bz=0 tarkoittaa, että
z=0. Mutta tällöin dz=0=1, mikä on ristiriita.
On siis oltava c≠0. Myös b≠0, sillä
muussa tapauksessa olisi oltava
[0cw+dy]=[10].
Viedään ensimmäinen kokonaismatriisi redusoituun riviporrasmuotoon.
[0b1cd0]R1↔R2→[cd00b1]1cR1→[1dc00b1]1bR2→[1dc0011b]R1−dcR2→[10−dbc011b]
Tässä oletukset b≠0 ja c≠0 voidaan
korvata oletuksella ad−bc=−bc≠0. Samalla
prosessilla ja oletuksilla saadaan jälkimmäiselle kokonaismatriisille
rref[0b0cd1]=[101c010].
Molemmissa tapauksissa siis oletus ad−bc≠0 riittää
yksikäsitteisten ratkaisujen olemassaoloon. Jälkimmäisen tapauksen
ratkaisussa voidaan lisäksi luonnollisesti kirjoittaa
w=−dbc=dad−bc,y=1b=−cad−bc,x=1c=−bad−bcjaz=0=aad−bc,
jotta se saadaan yhdenmukaiseksi ensimmäisen tapauksen ratkaisun kanssa.
Tämän esityksen mukainen ratkaisu luvuille w, y,
x ja z on siis aina olemassa, kun
ad−bc≠0, ja tällöin matriisi
B=[wxyz]=1ad−bc[d−b−ca]
on matriisin A käänteismatriisi. ◻
Esimerkki.
- Olkoon :math:`A = begin{bmatrix}
2 & 4 \ 1 & 3
end{bmatrix}` ja b=[−12]. Etsi A−1 ja ratkaise yhtälö
Ax=b.
Ratkaise yhtälöryhmä
{2x1+x2=2x1−x2=1
käänteismatriisin avulla.
Koska 2⋅3−4⋅1=2≠0, matriisi A
on kääntyvä ja
A−1=12[3−4−12]=[32−2−121].
Koska matriisi A on kääntyvä, yhtälöllä
Ax=b on yksikäsitteinen ratkaisu
x=A−1b=12[3−4−12][−12]=[−11252].
Yhtälöryhmä voidaan kirjoittaa matriisiyhtälönä
[211−1][x1x2]=[21].
Koska 2⋅(−1)−1⋅1=−3≠0, kerroinmatriisi
on kääntyvä ja yhtälöryhmällä on yksikäsitteinen ratkaisu
[x1x2]=−13[−1−1−12][21]=[10].
Alkeismatriisit
Yhtälöryhmää Ax=b ratkaistaessa viedään
kokonaismatriisi [A∣b] redusoituun riviporrasmuotoon
rivimuunnosten (rivien vaihto ja skaalaus, sekä lisääminen vakiolla
kerrottuna) avulla. Tähän prosessiin liittyvä teoria voidaan esittää
myös alkeismatriisien ja matriisitulon avulla.
Esimerkki.
Eräitä 3×3-alkeismatriiseja ovat
E13=[001010100],E2(k)=[1000k0001]jaE23(k)=[10001k001].
Seuraavan varsin ilmeisen tuloksen todistus on merkinnöiltään tekninen,
ja siksi sivuutetaan.
Lause.
Olkoon E rivimuunnokseen Q liittyvä
n×n-alkeismatriisi. Tällöin
AQ→EA
aina, kun A on n×r-matriisi.
Matriisiin A saadaan siis sovellettua rivimuunnos Q
kertomalla se vasemmalta sopivalla alkeismatriisilla E. Koska
rivimuunnokset ovat kääntyviä, saadaan seuraava tulos.
Lause.
Olkoon E rivimuunnokseen Q liittyvä
n×n-alkeismatriisi. Tällöin E on kääntyvä ja
E−1 on rivimuunnokseen Q−1 liittyvä
alkeismatriisi.
Oletetaan, että rivimuunnokseen Q−1 liittyy
alkeismatriisi F ja sovelletaan muunnoksia Q ja
Q−1 yksikkömatriisiin In.
InQ→EInQ−1→FEIn=FE
Toisaalta tuloksen on oltava alkuperäinen matriisi, joten
FE=In. Vastaavasti soveltamalla muunnoksia käänteisessä
järjestyksessä päätellään, että EF=In, joten matriisi
E on kääntyvä ja F=E−1.
Varsinainen lause on jo todistettu, mutta kirjoitetaan vielä ylös eri
rivimuunnoksiin liittyvien alkeismatriisien käänteismatriisit.
E−1ij=Eji,Ej(k)−1=Ej(1k),Eij(k)−1=Eij(−k).
Perustelut jätetään harjoitustehtäväksi. ◻
Alkeismatriisien avulla voidaan myös ilmaista riittäviä ja
välttämättömiä ehtoja neliömatriisin kääntyvyydelle. Seuraavan
kääntyvien matriisien peruslauseen tärkeyttä on vaikea korostaa
liikaa.
Lause.
Olkoon A n×n-neliömatriisi. Tällöin
seuraavat väitteet ovat yhtäpitäviä.
- A on kääntyvä.
- Yhtälöllä Ax=b on yksikäsitteinen ratkaisu
aina, kun b on avaruuden Rn
vektori.
- Yhtälöllä Ax=0 on ainoastaan
triviaaliratkaisu.
- Matriisin A redusoitu riviporrasmuoto on In.
- A on alkeismatriisien tulo.
Tämän muotoisen lauseen todistamiseksi riittää osoittaa,
että kohdasta 1 seuraa kohta 2, kohdasta 2 kohta 3, ja niin edelleen.
Lopuksi todistetaan, että viimeisestä kohdasta seuraa ensimmäinen.
Seuraavassa numerointi viittaa osatodistuksen oletuksena olevaan
väitteeseen.
Aiempi lause käsittelee tämän
kohdan.
Oletetaan, että yhtälöllä Ax=b on
yksikäsitteinen ratkaisu aina, kun b on avaruuden
Rn vektori. Tällöin myös yhtälöllä
Ax=0 on yksikäsitteinen ratkaisu. Mutta
triviaaliratkaisu x=0 on aina homogeenisen
lineaarisen yhtälöryhmän ratkaisu, joten nyt se on ainoa ratkaisu.
Oletetaan, että yhtälöllä Ax=0 on vain
triviaaliratkaisu. Yksikäsitteinen ratkaisu syntyy vain silloin, kun
yhtälöryhmässä ei ole vapaita muuttujia, ja tällöin matriisin
rref[A∣0] jokaisessa sarakkeessa
on johtava ykkönen. On siis oltava
rref(A)=In.
Oletetaan, että rref(A)=In. Tällöin on
löydyttävä rivimuunnokset Q1,Q2,…,Qk ja niitä
vastaavat alkeismatriisit E1,E2,…,Ek, joille
AQ1,Q2,…,Qk→Ek⋯E2E1A=rref(A)=In.
Alkeismatriisit ovat kääntyviä ja niiden käänteismatriisit ovat
alkeismatriiseja, joten
A=(Ek⋯E2E1)−1In=E−11E−12⋯E−1k,
alkeismatriisien tulo.
Oletetaan, että A=E1E2⋯El, missä
E1,E2,…,El ovat alkeismatriiseja. Tällöin
A on kääntyvien matriisien tulona kääntyvä.
◻
Kääntyvien matriisien peruslauseen todistuksen välitön sovellus on
neliömatriisin A kääntyvyyden tutkiminen ja käänteismatriisin
selvittäminen Gaussin eliminoinnilla.
Seuraus.
Olkoon A n×n-neliömatriisi. Jos
AQ1,Q2,…,Qk→In,
niin matriisi A on kääntyvä ja
InQ1,Q2,…,Qk→A−1.
Oletuksista seuraa, että
rref(A)=In, joten kääntyvien matriisien
peruslauseen nojalla A on kääntyvä. Jos rivimuunnoksia
Q1,Q2,…,Qk vastaavat alkeismatriisit ovat
E1,E2,…,Ek, niin (Ek⋯E2E1)A=In,
eli
A−1=Ek⋯E2E1=(Ek⋯E2E1)In,
ja väite on todistettu. ◻
Huomautus.
Neliömatriisi on siis kääntyvä täsmälleen silloin, kun
se voidaan muuntaa Gaussin eliminoinnilla yksikkömatriisiksi.
Jos siis A on n×n-neliömatriisi, sen kääntyvyyttä
voidaan tutkia seuraavasti Gaussin eliminoinnilla.
- Kirjoita kokonaismatriisi [A∣In].
- Sovella Gaussin eliminointia
matriisiin [A∣In] viedäksesi matriisin A
redusoituun riviporrasmuotoon.
- Jos rref(A)=In, niin kokonaismatriisi on
muodossa [In∣A−1]. Muuten käänteismatriisia ei ole.
Huomautus.
Vertaa edellä esitettyä kokonaismatriisimenetelmää
tapaan, jolla etsittiin 2×2-matriisin
käänteismatriisi.
Esimerkki.
Etsi matriisin A=[201010102] käänteismatriisi.
Sovelletaan Gaussin eliminointia kokonaismatriisiin
[A∣I3].
[201100010010102001]12R1→[10121200010010102001]R3−R1→[101212000100100032−1201]23R3→[10121200010010001−13023]R1−12R3→[100230−13010010001−13023]
Koska rref(A)=I3, käänteismatriisi
A−1=[230−13010−13023] on olemassa.
Esimerkki.
Tutki matriisien A=[2−2423−1312] ja B=[1−2−123040−2] kääntyvyyttä.
Sovelletaan Gaussin eliminointia kokonaismatriisiin
[A∣I3].
[2−2410023−1010312001]R2−R1→[2−2410005−5−110312001]R3−32R1→[2−2410005−5−11004−4−3201]R3−45R3→[2−2410005−5−110000−710−451]
Koska matriisin A puolelle syntyi nollarivi, matriisi A
ei voi olla riviekvivalentti yksikkömatriisin I3 kanssa. Tämän
vuoksi A on singulaarinen. Sovelletaan sitten Gaussin
eliminointia kokonaismatriisiin [B∣I3].
[1−2−110023001040−2001]R1−12R3→[−1−2010−1223001040−2001]R1+23R3→[1300123−1223001040−2001]3R1→[10032−3223001040−2001]R2−2R1→[10032−32030−6−3340−2001]13R2→[10032−32010−2−1140−2001]R3−4R1→[10032−32010−2−1100−2−12−87]−12R3→[10032−32010−2−1100164−72]
Koska matriisin B paikalle syntyy yksikkömatriisi I3,
niin käänteismatriisi on olemassa ja
B−1=[32−32−2−1164−72].
Esimerkki.
Osoita, että jos kääntyvän matriisin A
käänteismatriisi A−1 on symmetrinen, niin myös A on
symmetrinen.
Oletetaan, että A on kääntyvä
n×n-neliömatriisi ja että A−1=(A−1)T,
eli käänteismatriisi on symmetrinen. Tällöin käänteismatriisin
ominaisuuksien nojalla
AT=InAT=(AA−1)AT=A(A−1)TAT=A(AT)−1AT=AIn=A,
eli myös matriisi A on symmetrinen. ◻