Peatüki algus: Maatriksid
Eelmine: Cayley-Hamiltoni teoreem
Järgmine: LINEAARALGEBRA ARVUTUSMEETODID


Maatriksargumendiga funktsioonid

Olgu antud maatriks tex2html_wrap_inline2348 ja kompleksmuutuja funktsioon
displaymath2174
On palju võimalusi maatriksargumendiga funktsiooni f(A) defineerimiseks, lähtudes kompleksmuutuja funktsioonist f(z). Lihtsaim neist võimalustest tundub olema muutuja "z" vahetu asendamine muutujaga "A". Näiteks,
displaymath2175
ja
displaymath2176
ning
displaymath2177

displaymath2178

displaymath2179

displaymath2180
Osutub, et paljude probleemide lahendamisel ei ole selline lähenemine otstarbekas.

Definitsioon 2.9.1. Kui tex2html_wrap_inline2348 ja f(z) on analüütiline lahtises piirkonnas tex2html_wrap_inline2362 ning tex2html_wrap_inline2364 on kinnine lihtne joon (ei lõika iseennast) piirkonnas tex2html_wrap_inline2362 ja maatriksi A spekter tex2html_wrap_inline2370 sisaldub joone poolt hõlmatavas piirkonnas tex2html_wrap_inline2374 siis
 equation859
kusjuures integraali maatriksile rakendatakse elementide kaupa.

Märkus 2.9.1. Valem (28) on kompleksmuutuja funktsioonide korral tõestatud Cauchy' integraalvalemi analoog.

Näide 2.9.1. Olgu f(z)=z ja tex2html_wrap_inline2378 Uurime, kuidas toimub arvutamine eeskirja (28) põhjal. Kuna maatriksi A omaväärtusteks on tex2html_wrap_inline2382 ja tex2html_wrap_inline2384 siis valime joone tex2html_wrap_inline2386 kus tex2html_wrap_inline2388 Funktsioon f(z)=z on analüütiline piirkonnas tex2html_wrap_inline2392 Leiame esiteks,
displaymath2181

displaymath2182
ja siis,

displaymath2183

displaymath2184

displaymath2185

Lause 2.9.1. Kui tex2html_wrap_inline2394 siis
displaymath2186
kusjuures tex2html_wrap_inline2396 on ruumi tex2html_wrap_inline2398 vektor, mille k-s komponent on üks ja ülejäänud on nullid.

Tõestus. Olgu tex2html_wrap_inline2402 siis
displaymath2187
Kuna maatriksit integreeritakse elementide kaupa, siis saame, et

displaymath2188

Lause 2.9.2. Kui tex2html_wrap_inline2348 ja tex2html_wrap_inline2406 st on rahuldatud definitsioonis 2.9.1 esitatud tingimused, ning
 equation944
siis
 
equation951

Tõestus. Seoste (28) , (29) ja XX-1=I põhjal leiame, et
displaymath2189

displaymath2190

displaymath2191

displaymath2192
Kuna
displaymath2193
ja
displaymath2194

displaymath2195

displaymath2196
siis
displaymath2197
mida oligi vaja tõestada.

Lause 2.9.3. Kui tex2html_wrap_inline2348 ja tex2html_wrap_inline2414 on maatriksi A Jordani normaalkuju, kusjuures
displaymath2198
on tex2html_wrap_inline2086 Jordani blokk ning tex2html_wrap_inline2420 ja f(z) on analüütiline lahtisel hulgal, mis sisaldab maatriksi A spektrit tex2html_wrap_inline2426 siis
 equation993
kusjuures
 
equation997

Tõestus. Lause 2.9.2 põhjal piisab uurida vaid väärtust F(G), kus tex2html_wrap_inline2430 on tex2html_wrap_inline2432 Jordani blokk ja tex2html_wrap_inline2434 Olgu maatriks zI-G regulaarne. Kuna
displaymath2199
siis
displaymath2200
ja
displaymath2201

displaymath2202
ning arvestades tingimust tex2html_wrap_inline2438 peab paika lause väide.

Näide 2.9.2. Leida tex2html_wrap_inline2442 kui
displaymath2203

Kuna tex2html_wrap_inline2444 ja funktsioon tex2html_wrap_inline2446 on analüütiline punkti 0 ümbruses, siis on väärtuse tex2html_wrap_inline2450 leidmiseks rakendatav lauses 2.9.3 esitatud algoritm. Kasutame maatriksi A Jordani lahutuse leidmiseks paketti ''Maple'':
displaymath2204
Järelikult, J=diag(J1,J2), kusjuures
displaymath2205
Leiame, esiteks, valemi (3) abil maatriksid tex2html_wrap_inline2456 ja tex2html_wrap_inline2458 :
displaymath2206
ja
displaymath2207
Seejärel leiame valemi (2) abil meid huvitava maatriksi:

displaymath2208

displaymath2209

Järeldus 2.9.1. Kui tex2html_wrap_inline2348 ja tex2html_wrap_inline2462 ning siis
displaymath2210

Tõestus. Tegemist on lause 2.9.3 erijuhuga, kus kõik Jordani blokid on tex2html_wrap_inline2466

Näide 2.9.3. Kui maatriksi tex2html_wrap_inline2348 omaväärtused on tex2html_wrap_inline2470 tex2html_wrap_inline2472 ja tex2html_wrap_inline2474 on neile vastavad lineaarselt sõltumatud omavektorid, st tex2html_wrap_inline2474 moodustavad baasi (nn maatriksi A omabaasi) ruumis tex2html_wrap_inline2480 siis tex2html_wrap_inline2482 ja funktsioonide tex2html_wrap_inline2484 tex2html_wrap_inline2486 tex2html_wrap_inline2488 analüütilisusest kogu komplekstasandi lõplikus osas järeldub, et
displaymath2211
kusjuures tex2html_wrap_inline2490 ja

displaymath2212

displaymath2213

Vaatleme järgnevalt probleemi, mis tekib funktsiooni f(A) lähendamisel funktsiooniga g(A). Sellist laadi probleem tekib näiteks f(A) asendamisel tema q-astme Taylori polünoomiga.

Lause 2.9.4. Kui tex2html_wrap_inline2348 ja tex2html_wrap_inline2502 kusjuures
displaymath2214
on tex2html_wrap_inline2086 Jordani blokk ning tex2html_wrap_inline2420 ja funktsioonid f(z) ja g(z) on analüütilised lahtisel hulgal, mis sisaldab maatriksi A spektrit tex2html_wrap_inline2426 siis
 
equation1149

Tõestus. Kui valida h(z)=f(z)-g(z), siis
displaymath2215

displaymath2216
Lause 2.9.3 ja võrratuse tex2html_wrap_inline2518 abil leiame, et
displaymath2217
ning seega peab paika lause väide.

Näide 2.9.4. Olgu
displaymath2218
Hindame vahet tex2html_wrap_inline2522

Kuna tex2html_wrap_inline2524 ja funktsioonid tex2html_wrap_inline2526 ning g(z)=z on analüütilised punkti 0.1 ümbruses, siis võime kasutada lauses 2.9.4 saadud hinnangut (33). Esiteks, kasutame paketti ''Maple'' maatriksi A Jordani lahutuse leidmiseks:
displaymath2219
Järelikult, maatriksi A Jordani lahutuses on vaid üks Jordani blokk, st
displaymath2220
Teiseks, leiame paketi ''Maple'' abil maatriksi X konditsiooniarvu:
displaymath2221
Kuna


displaymath2222


displaymath2223
ja


displaymath2224
siis hinnangu (30) abil leiame, et
displaymath2225
On teada, et maatriksi A Jordani lahutuses esinev maatriks X pole üheselt määratud. Üritame valida maatriksi X nii, et konditsiooniarv k2(X) oleks minimaalne. Kasutades maatriksi A Jordani lahutuse leidmiseks Filipovi algoritmi (vaadake lauset 2.5.2.1), saame, et
displaymath2226
kusjuures
displaymath2227
Osutub, et maatriksi A Jordani lahutuseks on ka
displaymath2228
kusjuures
displaymath2229
Seega parim hinnang, mida saame lause 2.9.4 abil, on
displaymath2230
Teisalt, selle näite korral on rakendatav lauses 2.9.3 esitatud algoritm väärtuse tex2html_wrap_inline2550 leidmiseks. Kasutades valemit (32), saame, et
displaymath2231

displaymath2232
Valemi (31) abil arvutame meid huvitava funktsiooni väärtuse:
displaymath2233

displaymath2234
Järelikult,
displaymath2235

displaymath2236
ja
displaymath2237
Antud näite tulemusena võib väita, et lauses 2.9.5 tõestatud hinnang (33) on antud näite korral suhteliselt jäme.

Lause 2.9.5. Kui funktsiooni f(z) Maclaurini arendus
displaymath2238
koondub ringis, mis sisaldab maatriksi tex2html_wrap_inline2348 spektrit tex2html_wrap_inline2426 siis

displaymath2239

Tõestame selle väite täiendaval lisatingimusel, et maatriksil A leidub omavektoreist koosnev baas. Siis järelduse 2.9.1 põhjal
displaymath2240

displaymath2241

displaymath2242

Lause 2.9.6. Kui funktsiooni f(z) Maclaurini arendus
displaymath2238
koondub ringis, mis sisaldab maatriksi tex2html_wrap_inline2348 spektrit tex2html_wrap_inline2426 siis

displaymath2244

Tõestus. Defineerime maatriksi E(s) seosega
 equation1326
Kuitex2html_wrap_inline2568 tex2html_wrap_inline2570, siis tex2html_wrap_inline2568 on analüütiline ja seega,
 equation1335
kusjuures tex2html_wrap_inline2574 Võrreldes muutuja s astmeid seostes (34) ja (35), leiame, et tex2html_wrap_inline2578 omab kuju
displaymath2245
Kui tex2html_wrap_inline2580 siis

displaymath2246

Ülesanne 2.9.1. Näidake, et suvalise maatriksi tex2html_wrap_inline2348 korral
displaymath2247
ja

displaymath2248

Ülesanne 2.9.2. Rakendage lauset 2.9.6 vea hindamisel ligikaudsete seoste
displaymath2249
ja
displaymath2250
korral.

Lause 2.9.7 (Sylvesteri teoreem). Kui maatriksi tex2html_wrap_inline2348 kõik omaväärtused tex2html_wrap_inline2586 on erinevad, siis
 equation1393
või


 equation1400
kus tex2html_wrap_inline2588 tex2html_wrap_inline2590 on determinant, mis on saadud Vandermonde'i determinandist
displaymath2251
asendades k-nda reavektori
displaymath2252
vektoriga

displaymath2253

Näide 2.9.4. Leida tex2html_wrap_inline2594 kuitex2html_wrap_inline2596
Leiame kõigepealt maatriksi A omaväärtused
displaymath2254
ja siis kasutame valemit (36)
eqnarray1427

displaymath2255
ning kasutates valemit (37)

eqnarray1450


displaymath2256


displaymath2257
Sama probleemi lahendame samuti valemi tex2html_wrap_inline2600 kus S on maatriksi A omavektoreist koostatud maatriks. Leiame A omavektorid
displaymath2258
ja
displaymath2259
ning maatriksi
displaymath2260
Seega

eqnarray1497


Peatüki algus: Maatriksid
Eelmine: Cayley-Hamiltoni teoreem
Järgmine: LINEAARALGEBRA ARVUTUSMEETODID