Přeskočit na obsah

Wikipedista:Jirka Fiala/Pískoviště

Z Wikipedie, otevřené encyklopedie

V lineární algebře, spektrální rozklad je faktorizace matice do kanonické formy, whereby matice je reprezentována v podmínkách jeho eigenvalues a eigenvectors . Tímto způsobem lze faktorizovat pouze diagonalizovatelné matice . Když je matice faktorizovaná normální nebo skutečná symetrická matice, rozklad se nazývá „spektrální rozklad“, odvozený ze spektrální věty.

Základní teorie maticových vlastních vektorů a vlastních čísel[editovat | editovat zdroj]

(Nenulový) vektor dimenze je vlastním vektorem čtvercové matice , pokud splňuje lineární rovnici tvaru pro nějaké skalární . Potom se nazývá vlastní hodnota odpovídající . Geometricky řečeno, vlastní vektory jsou vektory, které pouze prodlužuje nebo zmenšuje, a velikost, o kterou se prodlužují/zmenšují, je vlastní hodnota. Výše uvedená rovnice se nazývá rovnice vlastních čísel nebo problém vlastních čísel.

To dává rovnici pro vlastní čísla Charakteristickým polynomem nazýváme (\lambda )</math> a rovnice, zvaná charakteristická rovnice, je rovnice polynomu -tého řádu v neznámé . Tato rovnice bude mít odlišných řešení, kde . Množina řešení, tedy vlastní čísla, se nazývá spektrum </ref> [1]

Pokud je pole skalárů algebraicky uzavřené, pak můžeme faktor jako Celé číslo se nazývá algebraická násobnost vlastní hodnoty . Součet algebraických násobků k  :

Pro každé vlastní číslo máme specifickou rovnici vlastního čísla Pro každou rovnici vlastních čísel bude existovat _{} \le _{}</math> lineárně nezávislých řešení. Lineární kombinace _{}</math> (kromě toho, které dává nulový vektor) jsou vlastní vektory spojené s vlastní hodnotou . Celé číslo _{}</math> se nazývá geometrická násobnost . Je důležité mít na paměti, že algebraická násobnost _{}</math> a geometrická násobnost _{}</math> se mohou a nemusí rovnat, ale vždy platí _{} \le _{}</math> . Nejjednodušší případ je samozřejmě, když _{} = _{} = 1</math> . Celkový počet lineárně nezávislých vlastních vektorů, _{''}</math>, lze vypočítat sečtením geometrických násobků Vlastní vektory lze indexovat vlastními čísly pomocí dvojitého indexu, přičemž '_{ij}</math> je tým vlastním vektorem té vlastní hodnoty. Vlastní vektory lze také indexovat pomocí jednoduššího zápisu jediného indexu '_{}</math>, kde = 1, 2, ..., _{''}1= = 1, 2, ..., _{''}</math> .

spektrální rozklad matice[editovat | editovat zdroj]

Nechť je čtvercová matice s lineárně nezávislými vlastními vektory (kde = 1, ..., ). Potom lze faktorizovat jako kde je čtvercová matice , jejíž tý sloupec je vlastním vektorem , a je diagonální matice, jejíž diagonální prvky jsou odpovídající vlastní hodnoty, . Všimněte si, že tímto způsobem lze faktorizovat pouze diagonalizovatelné matice . Například vadná matrice (což je smyková matice ) nelze diagonalizovat.

vlastních vektorů je obvykle normalizováno, ale nemusí být. Nenormalizovaná množina vlastních vektorů, může být také použita jako sloupce . To lze pochopit, když si všimneme, že velikost vlastních vektorů v se při rozkladu ruší přítomností '^{−1}</math> . Jestliže jedno z vlastních čísel má více lineárně nezávislých vlastních vektorů (to znamená, že geometrická násobnost je větší než 1), pak tyto vlastní vektory pro toto vlastní číslo mohou být vybrány jako vzájemně ortogonální ; pokud však dva vlastní vektory patří dvěma různým vlastním číslům, může být nemožné, aby byly navzájem ortogonální (viz příklad níže). Jedním speciálním případem je, že pokud je normální matice, pak podle spektrální věty je vždy možné diagonalizovat v ortonormální bázi {q_{i}} .

Rozklad lze odvodit ze základní vlastnosti vlastních vektorů: Lineárně nezávislé vlastní vektory s nenulovými vlastními hodnotami tvoří základ (ne nutně ortonormální) pro všechny možné součiny ', pro ' ∈ ''^{}</math>, který je stejný jako obraz (nebo rozsah ) příslušné transformace matice, a také sloupcový prostor matice . Počet lineárně nezávislých vlastních vektorů s nenulovými vlastními čísly je roven hodnosti matice a také rozměru obrazu (nebo rozsahu) příslušné transformace matice a také jejímu sloupcovému prostoru.

Lineárně nezávislé vlastní vektory s vlastní hodnotou nula tvoří základ (který lze zvolit jako ortonormální) pro nulový prostor (také známý jako jádro) maticové transformace .

Příklad[editovat | editovat zdroj]

Reálná matice 2 × 2 lze rozložit na diagonální matici násobením nesingulární matice Pak pro nějakou skutečnou diagonální matici .

Vynásobením obou stran rovnice vlevo  : Výše uvedená rovnice se dá rozložit na dvě simultánní rovnice : Rozložení vlastních hodnot a  : Pronájem to nám dává dvě vektorové rovnice: A může být reprezentován jedinou vektorovou rovnicí zahrnující dvě řešení jako vlastní čísla: kde představuje dvě vlastní čísla a a představuje vektory a .

Posunutí na levou stranu a vyloučení Protože je nesingulární, je důležité, aby bylo nenulové. Proto, Tím pádem což nám dává řešení vlastních hodnot pro matici jako nebo a výsledná diagonální matice z vlastního rozkladu je tedy {{Nowrap|

Uvedení řešení zpět do výše uvedených simultánních rovnic Řešení rovnic máme Matice potřebná pro spektrální rozklad je tedy to je:

Matice inverzní prostřednictvím vlastního rozkladu[editovat | editovat zdroj]

Pokud lze matici rozložit a pokud žádné z jejích vlastních hodnot není nulové, pak je invertibilní a její inverzní je dána vztahem Li je symetrická matice, protože je tvořena z vlastních vektorů , je tedy zaručeno, že jde o ortogonální matici . Navíc, protože je diagonální matice, její inverzní hodnotu lze snadno vypočítat:

Praktické důsledky[editovat | editovat zdroj]

Když se na matici naměřených reálných dat použije spektrální rozklad, může být inverzní hodnota méně platná, pokud jsou všechna vlastní čísla použita bez úprav ve výše uvedeném tvaru. Je to proto, že jak jsou vlastní čísla relativně malá, jejich příspěvek k inverzi je velký. Ty, které jsou blízko nuly nebo jsou v "šumu" měřicího systému, budou mít nepřiměřený vliv a mohly by bránit řešení (detekci) pomocí inverzní.

Byla navržena dvě zmírnění: zkrácení malých nebo nulových vlastních hodnot a rozšíření nejnižší spolehlivé vlastní hodnoty na hodnoty pod nimi. Viz také Tichonovova regularizace jako statisticky motivovaná, ale neobjektivní metoda pro odvalování vlastních hodnot, když v nich dominuje šum.

První metoda zmírnění je podobná řídkému vzorku původní matrice a odstraňuje složky, které nejsou považovány za cenné. Pokud se však řešení nebo proces detekce blíží úrovni šumu, může zkrácení odstranit komponenty, které ovlivňují požadované řešení.

Druhé zmírnění rozšiřuje vlastní hodnotu, takže nižší hodnoty mají mnohem menší vliv na inverzi, ale stále přispívají, takže řešení blízko šumu budou stále nalezena.

Spolehlivou vlastní hodnotu lze nalézt za předpokladu, že vlastní hodnoty extrémně podobné a nízké hodnoty jsou dobrou reprezentací šumu měření (který se u většiny systémů předpokládá jako nízký).

Pokud jsou vlastní čísla seřazena podle hodnoty, pak lze spolehlivou vlastní hodnotu najít minimalizací Laplaciánu seřazených vlastních hodnot: kde vlastní hodnoty jsou indexovány s k označení třídění. Poloha minimalizace je nejnižší spolehlivá vlastní hodnota. V měřicích systémech je druhou odmocninou této spolehlivé vlastní hodnoty průměrný šum na součástech systému.

Funkční kalkul[editovat | editovat zdroj]

spektrální rozklad umožňuje mnohem snadnější výpočet mocninných řad matic. Pokud  ()</math> je dáno pak to víme Protože je diagonální matice, lze funkce velmi snadno vypočítat: Mimodiagonální prvky  (\lambda )</math> jsou nula; tedy  (\lambda )</math> je také diagonální matice. Proto výpočet  ('')</math> redukuje na pouhý výpočet funkce na každém z vlastních čísel.

Podobná technika funguje obecněji s holomorfním funkčním kalkulem, using shora . Opět to zjišťujeme

Příklady[editovat | editovat zdroj]

což jsou příklady funkcí  . dále  je matice exponenciální .

Rozklad pro spektrální matice[editovat | editovat zdroj]

Spektrální matice jsou matice, které mají odlišné vlastní hodnoty a kompletní sadu vlastních vektorů. Tato charakteristika umožňuje, aby spektrální matice byly plně diagonální, což znamená, že je lze rozložit na jednodušší formy pomocí vlastního rozkladu. Tento proces rozkladu odhaluje základní poznatky o struktuře a chování matice, zejména v oblastech, jako je kvantová mechanika, zpracování signálu a numerická analýza. [2]

Normální matice[editovat | editovat zdroj]

Čtvercová matice s komplexní hodnotou je normální (tzn. , kde je konjugovaná transpozice ) právě tehdy, když ji lze rozložit jako , kde je unitární matice (tj ) a diag( </ref> Sloupce z tvoří ortonormální základ a jsou vlastními vektory s odpovídajícími vlastními hodnotami . [3]

Uvažujme například normální matici 2 x 2 .

Vlastní hodnoty jsou a .

(Normalizované) vlastní vektory odpovídající těmto vlastním číslům jsou a .

Diagonalizace je , kde , a .

Ověření je .

Tento příklad ilustruje proces diagonalizace normální matice nalezením vlastních čísel a vlastních vektorů, které tvoří unitární matici , diagonální matice a ověření rozkladu.

Podmnožiny důležitých tříd matic

Reálné symetrické matice[editovat | editovat zdroj]

Speciálním případem je, že pro každých reálných symetrických matic jsou vlastní čísla reálná a vlastní vektory lze vybrat reálné a ortonormální . Reálnou symetrickou matici lze tedy rozložit jako , kde je ortogonální matice, jejíž sloupce jsou skutečné, ortonormální vlastní vektory , a je diagonální matice, jejíž vstupy jsou vlastní čísla </ref>

Diagonalizovatelné matice[editovat | editovat zdroj]

Diagonalizovatelné matice lze rozložit pomocí vlastního rozkladu za předpokladu, že mají úplnou sadu lineárně nezávislých vlastních vektorů. Mohou být vyjádřeny jako , kde je matice, jejíž sloupce jsou vlastními vektory a je diagonální matice skládající se z odpovídajících vlastních hodnot . [3]

Pozitivní určité matice[editovat | editovat zdroj]

Kladně definitní matice jsou matice, pro které jsou všechna vlastní čísla kladná. Lze je rozložit jako pomocí Choleského rozkladu, kde je nižší trojúhelníková matice. [4]

Unitární a hermitovské matice[editovat | editovat zdroj]

Unitární matice vyhovují (skutečný případ) popř (složitý případ), kde označuje konjugovanou transpozici a označuje konjugovanou transpozici. Diagonalizují pomocí unitárních transformací . [3]

Hermitovské matrice vyhovují , kde označuje konjugovanou transpozici. Mohou být diagonalizovány pomocí unitárních nebo ortogonálních matic . [3]

Užitečná fakta[editovat | editovat zdroj]

Užitečná fakta týkající se vlastních čísel[editovat | editovat zdroj]

  • Součin vlastních čísel se rovná determinantu
  • Součet vlastních čísel se rovná stopě
  • Pokud jsou vlastní hodnoty '\lambda _{}</math> a je invertibilní, pak jsou vlastní hodnoty '^{−1}</math> jednoduše Nelze pochopit (syntaktická chyba): {\displaystyle ''\lambda ''{{su|b=<math>i} |p=−1</math></br> Nelze pochopit (syntaktická chyba): {\displaystyle ''\lambda ''{{su|b=<math>i} |p=−1</math>
  • Jestliže vlastní čísla jsou , pak vlastní čísla  ('')</math> jsou prostě  (\lambda _{})</math> pro jakoukoli holomorfní funkci .

Užitečná fakta o vlastních vektorech[editovat | editovat zdroj]

  • Je-li hermitovské a plné, může být základ vlastních vektorů vybrán tak, aby byl vzájemně ortogonální . Vlastní hodnoty jsou skutečné.
  • Vlastní vektory '^{−1}</math> jsou stejné jako vlastní vektory .
  • Vlastní vektory jsou definovány pouze do multiplikativní konstanty. To znamená, že pokud '' je také vlastní vektor pro jakékoli skalární . Zejména Nelze pochopit (syntaktická chyba): {\displaystyle −'<math>v} a ^{}' (pro libovolné θ ) jsou také vlastní vektory.
  • V případě degenerovaných vlastních čísel (vlastní hodnota mající více než jeden vlastní vektor) mají vlastní vektory další volnost lineární transformace, to znamená, že jakákoli lineární (ortonormální) kombinace vlastních vektorů sdílejících vlastní hodnotu (v degenerovaném podprostoru) je sama sebou. vlastní vektor (v podprostoru).

Užitečná fakta týkající se vlastního rozkladu[editovat | editovat zdroj]

  • can be eigendecomposed if and only if the number of linearly independent eigenvectors, _{''}</math>, equals the dimension of an eigenvector: _{''} =
  • If the field of scalars is algebraically closed and if (\lambda )</math> has no repeated roots, that is, if then can be eigendecomposed.
  • The statement " can be eigendecomposed" does not imply that has an inverse as some eigenvalues may be zero, which is not invertible.
  • The statement " has an inverse" does not imply that can be eigendecomposed. A counterexample is , which is an invertible defective matrix.

Užitečná fakta o inverzní matici[editovat | editovat zdroj]

  • can be inverted if and only if all eigenvalues are nonzero:
  • If and _{''} = , the inverse is given by

Numerické výpočty[editovat | editovat zdroj]

Numerický výpočet vlastních čísel[editovat | editovat zdroj]

Předpokládejme, že chceme vypočítat vlastní čísla dané matice. Pokud je matice malá, můžeme je spočítat symbolicky pomocí charakteristického polynomu . To je však u větších matic často nemožné, v takovém případě musíme použít numerickou metodu .

V praxi se vlastní čísla velkých matic nepočítají pomocí charakteristického polynomu. Výpočet polynomu je sám o sobě nákladný a přesné (symbolické) kořeny polynomu vysokého stupně může být obtížné vypočítat a vyjádřit: Abel-Ruffiniho teorém naznačuje, že kořeny polynomů vysokého stupně (5 nebo více) obecně nemohou být vyjádřen jednoduše pomocí -tých kořenů. Proto jsou obecné algoritmy k nalezení vlastních vektorů a vlastních hodnot iterativní .

Iterativní numerické algoritmy pro aproximaci kořenů polynomů existují, jako je Newtonova metoda, ale obecně je nepraktické vypočítat charakteristický polynom a poté tyto metody aplikovat. Jedním z důvodů je, že malé zaokrouhlovací chyby v koeficientech charakteristického polynomu mohou vést k velkým chybám ve vlastních číslech a vlastních vektorech: kořeny jsou extrémně špatně podmíněnou funkcí koeficientů. [5]

Jednoduchou a přesnou iterační metodou je mocninná metoda : vybere se náhodný vektor a posloupnost jednotkových vektorů se vypočítá jako Tato posloupnost bude téměř vždy konvergovat k vlastnímu vektoru odpovídajícímu vlastní hodnotě největší velikosti za předpokladu, že má nenulovou složku tohoto vlastního vektoru v bázi vlastního vektoru (a také za předpokladu, že existuje pouze jedna vlastní hodnota největší velikosti). Tento jednoduchý algoritmus je užitečný v některých praktických aplikacích; například Google jej používá k výpočtu hodnocení stránek dokumentů v jejich vyhledávači. Výkonová metoda je také výchozím bodem pro mnoho sofistikovanějších algoritmů. Například tím, že ponecháme nejen poslední vektor v sekvenci, ale místo toho se podíváme na rozsah všech vektorů v sekvenci, můžeme získat lepší (rychlejší konvergující) aproximaci vlastního vektoru, a tato myšlenka je základem Arnoldiho iterace . [5] Alternativně je důležitý QR algoritmus také založen na jemné transformaci mocninné metody. [5]

Numerický výpočet vlastních vektorů[editovat | editovat zdroj]

Jakmile jsou vlastní čísla vypočtena, lze vlastní vektory vypočítat řešením rovnice pomocí Gaussovy eliminace nebo jakékoli jiné metody pro řešení maticových rovnic .

V praktických metodách vlastních čísel ve velkém měřítku se však vlastní vektory obvykle počítají jinými způsoby, jako vedlejší produkt výpočtu vlastních čísel. Například při mocninné iteraci je vlastní vektor ve skutečnosti vypočítán před vlastní hodnotou (která se obvykle vypočítává pomocí Rayleighova podílu vlastního vektoru). [5] V QR algoritmu pro hermitovskou matici (nebo jakoukoli normální matici) jsou ortonormální vlastní vektory získány jako součin matic z kroků v algoritmu. [5] (U obecnějších matic poskytuje algoritmus QR nejprve Schurův rozklad, ze kterého lze získat vlastní vektory zpětným substitučním postupem. [6] ) Pro hermitovské matice je algoritmus vlastních hodnot Divide-and-conquer efektivnější než algoritmus QR. pokud jsou požadovány jak vlastní vektory, tak vlastní hodnoty. [5]

Viz také[editovat | editovat zdroj]

[[Kategorie:Maticové rozklady]] [[Kategorie:Teorie matic]] [[Kategorie:Údržba:Články s nekontrolovanými překlady]]

  1. [s.l.]: [s.n.] Dostupné online. ISBN 978-0-471-50728-4. 
  2. Dostupné online. 
  3. a b c d Dostupné online.  Chybná citace: Neplatná značka <ref>; název „:0“ použit vícekrát s různým obsahem
  4. [s.l.]: [s.n.] ISBN 9781611977431. 
  5. a b c d e f [s.l.]: [s.n.] ISBN 978-0-89871-361-9.  Chybná citace: Neplatná značka <ref>; název „Trefethen“ použit vícekrát s různým obsahem
  6. [s.l.]: [s.n.] Dostupné online. ISBN 978-0-387-98959-4.