Symetrická matice

Z testwiki
Verze z 2. 1. 2025, 20:44, kterou vytvořil imported>Vít Karásek (-zavádející odkaz)
(rozdíl) ← Starší verze | zobrazit aktuální verzi (rozdíl) | Novější verze → (rozdíl)
Skočit na navigaci Skočit na vyhledávání
Schematické naznačení symetrie v čtvercové matici stupně pět

Symetrická matice je v lineární algebře každá čtvercová matice, která je osově souměrná podle své hlavní diagonály. Jedná o čtvercovou matici, která se shoduje se svou transponovanou maticí, neboli 𝑨=𝑨T.

Symetrické matice se v lineární algebře používají k popisu symetrických bilineárních forem. Matice samoadjungovaného lineárního zobrazení vzhledem k ortonormální bázi je vždy symetrická. Soustavy lineárních rovnic se symetrickými maticemi soustavy lze řešit efektivně a numericky stabilně. Dále se symetrické matice používají v ortogonálních projekcích a při polárním rozkladu matic.

Symetrické matice mají aplikace také v geometrii, analýze, teorii grafů a stochastice.

Definice

Čtvercová matice 𝑨 řádu n nad tělesem T, se nazývá symetrická, pokud pro všechna i,j{1,,n} platí:

aij=aji.

Matice, která není symetrická se nazývá asymetrická, neplést s antisymetrickou maticí.

Příklady

Symetrickými maticemi jsou například:

(2),(1557),(1110),(130326065).

Obecně mají symetrické matice o rozměrech 2×2, 3×3 a 4×4 následující podobu:

(abbc),(abcbdecef),(abcdbefgcfhidgij).

Speciální případy

Některé symetrické matice se zvláštními vlastnostmi mají vlastní název:

Vlastnosti

Popis

U symetrické matice stačí znát pouze prvky na diagonále a pod ní

U symetrické matice 𝑨Tn×n stačí znát n prvků na diagonále a n(n1)2 prvků na jedné straně diagonály (nad nebo pod). Hodnoty prvků na opačné straně diagonály lze odvodit ze symetrie matice. Symetrická matice může mít nejvýše

n+n(n1)2=n(n+1)2

různých prvků. Ve srovnání s nesymetrickými maticemi řádu n, které mohou mít až n2 různých prvků, jde zhruba o poloviční množství dat, a proto byly pro symetrické matice navrženy speciální formáty pro ukládání v počítači. [1]

Vektorový prostor symetrických matic

Součet 𝑨+𝑩 dvou symetrických matic 𝑨,𝑩Tn×n je vždy symetrická matice, protože

(𝑨+𝑩)T=𝑨T+𝑩T=𝑨+𝑩

Stejně tak i skalární násobek c𝑨 symetrické matice skalárem cT je opět symetrická matice. Protože je nulová matice také symetrická, tvoří množina symetrických matic řádu n vektorový podprostor

Symmn={𝑨Tn×n:𝑨T=𝑨}

prostoru čtvercových matic Tn×n. Tento podprostor má dimenzi n2+n2. Jeho bázi lze vytvořit z matic 𝐄ii pro i{1,,n}, a součtů 𝐄ij+𝐄ji pro 1i<jn. Uvedené matice 𝐄ij tvoří standardní bázi prostoru Tn×n , čili mají jediný nenulový prvek 𝐞ij=1.

Rozklad

Pokud je charakteristika tělesa T různá od 2, lze libovolnou čtvercovou matici 𝑴Tn×n zapsat jednoznačně jako součet 𝑴=𝑨+𝑩, kde matice 𝑨 je symetrická a matice 𝑩 je antisymetrická:

𝑨=12(𝑴+𝑴T)  a  𝑩=12(𝑴𝑴T)

Antisymetrické matice tvoří vektorový podprostor prostoru čtvercových matic. Značí se Skewna má dimenzi n2n2. Prostor čtvercových matic Tn×n dimenze n2 lze vyjádřit jako direktní součet

Tn×n=SymmnSkewn

prostorů symetrických a antisymetrických matic.

Součin

Součin 𝑨𝑩 dvou symetrických matic 𝑨,𝑩Tn×n nemusí být opět symetrická matice. Součin symetrických matic je symetrický, právě když je součin 𝑨 a 𝑩 komutativní. Jinými slovy, pokud součin splňuje 𝑨𝑩=𝑩𝑨, pak také platí:

(𝑨𝑩)T=𝑩T𝑨T=𝑩𝑨=𝑨𝑩 .

Pro symetrickou matici 𝑨 proto platí, že symetrické jsou všechny její mocniny 𝑨k, kde k, i její maticová exponenciála e𝑨.

Pro každou matici 𝑴Tm×n jsou matice 𝑴𝑴T typu m×m, i matice 𝑴T𝑴 typu n×n symetrické.

Kongruence a podobnost

Každá matice 𝑩Tn×n, která je kongruentní symetrické matici 𝑨Tn×n, je také symetrická, protože platí

𝑩T=(𝑺T𝑨𝑺)T=𝑺T𝑨T𝑺=𝑺T𝑨𝑺=𝑩,

přičemž 𝑺Tn×n je odpovídající regulární matice.

Na druhou stranu existují i nesymetrické matice, které jsou podobné symetrické matici.

Inverze

Pokud je symetrická matice 𝑨Tn×n regulární, potom matice k ní inverzní 𝑨1 je symetrická, protože pro ni platí:

(𝑨1)T=(𝑨T)1=𝑨1 .

V tomto případě jsou symetrické všechny mocniny 𝑨k pro k.

Reálné symetrické matice

Symetrické matice s reálnými prvky mají řadu dalších vlastností.

Normální matice

Reálná symetrická matice 𝑨n×n je normální, protože platí

𝑨T𝑨=𝑨𝑨=𝑨𝑨T.

Každá reálná symetrická matice komutuje se svou transpozicí. Existují však i normální matice, které nejsou symetrické, například antisymetrické matice.

Hermitovské matice

Protože se na každé číslo shoduje se svým komplexně sdruženým protějškem, neboli z=z, splývají reálné symetrické matice s reálnými hermitovskými. Formálně:

𝑨H=𝑨T=𝑨T=𝑨,

přičemž 𝑨H je hermitovská transpozice matice 𝑨 a 𝑨 je komplexně sdružená matice k 𝑨.

Reálná symetrická matice 𝑨n×n je vždy hermitovská mimo jiné i proto, že vzhledem k standardnímu skalárnímu součinu , na n splňuje:

𝑨𝒙,𝒚=(𝑨𝒙)T𝒚=𝒙T𝑨T𝒚=𝒙T𝑨𝒚=𝒙,𝑨𝒚

pro všechny vektory 𝒙,𝒚n. Reálné symetrické matice jsou hermitovské i s ohledem na standardní skalární součin nad .

Vlastní čísla

Jednotková kružnice je transformována na elipsu pomocí reálné symetrické matice řádu 2. Poloosy elipsy odpovídají vlastním vektorům matice a jejich délky vlastním číslům.

Vlastní čísla reálné symetrické matice 𝑨n×n, tedy řešení rovnice 𝑨𝒙=λ𝒙, jsou vždy reálná. Kdyby λ bylo komplexní vlastní číslo matice 𝑨 příslušné netriviálnímu vlastnímu vektoru 𝒙n, 𝒙0, pak z toho, že 𝑨 je hermitovská plyne:

λ𝒙,𝒙=𝒙,λ𝒙=𝒙,𝑨𝒙=𝑨𝒙,𝒙=λ𝒙,𝒙=λ𝒙,𝒙.

Protože pro každé 𝒙0 platí 𝒙,𝒙0, musí vlastní číslo λ splňovat λ=λ, a proto je reálné. V důsledku lze i příslušný vlastní vektor 𝒙 zvolit reálný.

Násobnosti vlastních čísel

Pro každou reálnou symetrickou matici 𝑨n×n se algebraické a geometrické násobnosti všech vlastních čísel shodují. Důvod je následující. Pro vlastní číslo λ matice 𝑨 s geometrickou násobností k existuje ortonormální báze {𝒙1,,𝒙k} prostoru vlastních vektorů příslušných k λ. Tuto bázi lze rozšířit pomocí vektorů{𝒙k+1,,𝒙n} na ortonormální bázi celého prostoru n. S pomocí ortogonální matice 𝑺=(𝒙1𝒙n) je matice 𝑨 převedena na podobnou

𝑪=𝑺1𝑨𝑺=𝑺T𝑨𝑺=(λ𝐈00𝑿)

což je bloková diagonální matice s bloky λ𝐈k×k a 𝑿(nk)×(nk). Vzhledem k tomu, že matice 𝑨 je hermitovská a vektory 𝒙1,,𝒙ntvoří ortonormální bázi, platí pro prvky cij matice C s indexy min{i,j}k, že:

cij=𝒙i,𝑨𝒙j=𝑨𝒙i,𝒙j=λ𝒙i,𝒙j=λδij,

kde δij je Kroneckerovo delta. Vektory 𝒙k+1,,𝒙n nejsou podle předpokladu vlastní vektory matice 𝑨 příslušné vlastnímu číslu λ, proto λ není žádným vlastním číslem matice 𝑿. Vlastní číslo λ matice 𝑪 má podle vzorce pro determinant blokových matic shodnou algebraickou i geometrickou násobnost k. Totéž platí i pro matici 𝑨 díky vzájemné podobnosti s maticí 𝑪. [2]

Diagonalizovatelnost

Vzhledem k tomu, že se algebraické a geometrické násobnosti všech vlastních čísel shodují, a protože vlastní vektory příslušné různým vlastním číslům jsou lineárně nezávislé, tvoří vlastní vektory reálné symetrické matice 𝑨 bázi prostoru n. Reálná symetrická matice je tedy vždy diagonalizovatelná, to znamená, že existuje regulární matice 𝑺n×n a diagonální matice 𝑫n×nsplňující:

𝑺1𝑨𝑺=𝑫

Matice 𝑺=(𝒙1𝒙n) je sestavena z vlastních vektorů 𝒙1,,𝒙n po sloupcích a matice 𝑫=diag(λ1,,λn) má vlastní čísla λ1,,λn příslušná těmto vlastním vektorům na diagonále. Vzhledem k tomu, že sloupce matice 𝑺, neboli vlastní vektory lze libovolně přerovnat, může být odpovídající pořadí prvků na diagonále 𝑫 libovolné. V důsledku si dvě reálné symetrické matice jsou podobné, právě když mají stejná vlastní čísla. Kromě toho jsou dvě reálné symetrické matice současně diagonalizovatelné, právě když spolu komutují.

Ortogonální diagonalizace

U symetrických matic platí, že vlastní vektory (modrý a fialový) příslušné různým vlastním číslům (zde 3 a 1) jsou na sebe kolmé. Při provedení transformace odpovídající matici se modré vektory třikrát prodlouží, zatímco fialové vektory si svou délku zachovají.

Vlastní vektory 𝒙i,𝒙j příslušné dvěma různým vlastním číslům λiλj reálné symetrické matici 𝑨n×n jsou vzájemně kolmé. Uvedený vztah opět z následující vlastnosti hermitovských matic 𝑨:

λi𝒙i,𝒙j=λi𝒙i,𝒙j=𝑨𝒙i,𝒙j=𝒙i,𝑨𝒙j=𝒙i,λj𝒙j=λj𝒙i,𝒙j .

Z předpokladu, že λi a λj jsou různá, pak plyne 𝒙i,𝒙j=0. Vlastní vektory 𝑨 tvoří ortonormální bázi prostoru n. Každou reálnou symetrickou matici lze proto ortogonálně diagonalizovat, neboli existuje ortogonální matice 𝑺 splňující:

𝑺T𝑨𝑺=𝑫

Tato reprezentace tvoří základ pro transformaci hlavní osy a je nejjednodušší verzí spektrální věty.

Parametry

Každá reálná symetrická matice 𝑨n×n diagonalizovatelná, a proto pro její stopu platí:

tr𝑨=λ1++λn

Její determinant tudíž splňuje:

det𝑨=λ1λn

Hodnost reálné symetrické matice je rovna počtu nenulových vlastních čísel. Za pomoci Kroneckerovy delty ji lze vyjádřit výrazem

rank𝑨=n(δλ1,0++δλn,0) .

Reálná symetrická matice je regulární, právě když má všechna vlastní čísla nenulová. Spektrální norma reálné symetrické matice je

𝑨2=max{|λ1|,,|λn|}

a tedy rovna spektrálnímu poloměru matice. Frobeniova norma vyplývá z normality

𝑨F=λ12++λn2.

Definitnost

Šablona:Podrobně Pro reálnou symetrickou matici 𝑨n×n a vektor 𝒙n se výraz

Q𝑨(𝒙)=𝒙T𝑨𝒙=𝒙,𝑨𝒙

nazývá kvadratická forma určená maticí 𝑨. Podle toho, jestli je Q𝑨(𝒙) pro všechna 𝒙0 kladná, resp. nezáporná, záporná či nekladná, nazývá se matrice 𝑨 pozitivně definitní, resp. pozitivně semidefinitní, negativně definitní nebo negativně semidefinitní. Pokud Q𝑨(𝒙) nabývá kladných i záporných hodnot, nazývá se matice 𝑨 indefinitní. Definitnost reálné symetrické matice závisí na znaméncích jejích vlastních čísel. Pokud jsou všechna vlastní čísla kladná, je matice pozitivně definitní, pokud jsou všechna záporná, je matice negativně definitní atd. Trojice čísel daná počtem kladných, záporných a nulových vlastních čísel se nazývá signatura matice. Podle Sylvesterova zákona setrvačnosti je signatura zachována u kongruentních reálných symetrických matic.

Odhady vlastních čísel

Podle Courant-Fischerovy věty lze nejmenší a největší vlastní číslo symetrické 𝑨n×n odhadnout pomocí Rayleighova kvocientu. Konkrétně, pro všechna netriviální 𝒙n platí:

min{λ1,,λn}𝒙,𝑨𝒙𝒙,𝒙max{λ1,,λn}

Rovnost platí, právě když je 𝒙 je vlastní vektor příslušný k danému vlastnímu číslu. V důsledku lze nejmenší a největší vlastní číslo reálné symetrické matice určit minimalizací nebo maximalizací Rayleighova kvocientu.

Další možnost pro odhad vlastních čísel nabízejí Geršgorinovy kruhy, které u reálných symetrických matic mají tvar intervalů.

Pro dvě reálné symetrické matice 𝑨,𝑩n×n s vlastními čísly seřazenými sestupně λ1λn a μ1μn platí odhad

tr(𝑨𝑩)λ1μ1++λnμn.

Rovnost je splněna, právě když matice 𝑨 a 𝑩 jsou současně diagonalizovatelné vzhledem k uspořádání vlastních čísel, neboli když existuje ortogonální matice 𝑺n×n taková, že platí 𝑨=𝑺diag(λ1,,λn)𝑺T a 𝑩=𝑺diag(μ1,,μn)𝑺T. Uvedená nerovnost zobecňuje Cauchy-Schwarzovu nerovnost pro Frobeniův skalární součin a permutační nerovnost pro vektory. [3]

Komplexní symetrické matice

Rozklad

Podobně jako u reálných matic lze prostor komplexních čtvercových matic n×n zapsat jako direktní součet prostorů symetrických a antisymetrických matic:

n×n=SymmnSkewn

Jde zároveň o ortogonální součet vzhledem k Frobeniově skalárnímu součinu, protože pro všechny matice 𝑨Symmn a 𝑩Skewn platí:

𝑨,𝑩F=tr(𝑨H𝑩)=tr(𝑨𝑩)=tr(𝑩𝑨)=tr((𝑩𝑨)T)=tr(𝑨H𝑩)=𝑨,𝑩F

z čehož vyplývá 𝑨,𝑩F=0. Ortogonalita rozkladu platí i pro reálný maticový prostor n×n.

Spektrum

Pro komplexní matice 𝑨n×n nemá symetrie žádný zvláštní vliv na spektrum matice. Komplexní symetrická matice může mít nereálná vlastní čísla. Například komplexní symetrická matice 𝑨=(1ii1)2×2má dvě vlastní čísla λ1,2=1±i .

Existují komplexní symetrické matice, které nelze diagonalizovat. Například matice 𝑨=(1ii1)2×2 má jediné vlastní číslo λ=0 s algebraické násobnosti dvě a geometrické násobnosti jedna. Obecně platí, že jakákoli komplexní čtvercová matice je podobná komplexní symetrické matici. Spektrum komplexní symetrické matice proto nevykazuje žádné zvláštnosti. [4]

Komplexním rozšířením reálných symetrických matic, pokud jde o matematické vlastnosti, jsou hermitovské matice.

Rozklad

Libovolnou komplexní symetrickou matici 𝑨n×n lze pomocí Autonne-Takagiho faktorizace rozložit na součin

𝑨=𝑼T𝑫𝑼,

kde matice 𝑼n×n je unitární, 𝑫=diag(σ1,,σn)n×nje reálná diagonální. Prvky diagonální matice jsou singulární hodnoty 𝑨, neboli odmocniny vlastních čísel matice 𝑨H𝑨. [5]

Použití

Symetrické bilineární formy

Každá bilineární forma b:V×VT na vektorovém prostoru V dimenze n nad tělesem T může být vzhledem k bázi {𝒗1,,𝒗n} prostoru V reprezentována čtvercovou maticí 𝑨bTn×n danou vztahem:

(𝑨b)ij=b(𝒗i,𝒗j)

Pokud je bilineární forma symetrická, pak platí b(𝒗,𝒘)=b(𝒘,𝒗) pro všechny 𝒗,𝒘V, a matice 𝑨b je symetrická. Naopak každá symetrická matice 𝑨Tn×n definuje symetrickou bilineární formu b𝑨:Tn×TnT vztahem:

b𝑨(𝒙,𝒚)=𝒙T𝑨𝒚

Je-li matice 𝑨n×n navíc pozitivně definitní, pak b𝑨 představuje skalární součin na euklidovském prostoru n.

Samoadjungované zobrazení

Je-li (V,,) reálný prostor se skalárním součinem dimenze n, pak může být každé lineární zobrazení f:VV vzhledem k ortonormální bázi {𝒆1,,𝒆n} prostoru V reprezentováno maticí zobrazení

𝑨f=(aij)n×n,

kde f(𝒆j)=a1j𝒆1++anj𝒆n pro j=1,,n. Matice zobrazení 𝑨f je symetrická, právě když je zobrazení f samoadjungované. To vyplývá ze vztahu

f(𝒗),𝒘=(𝑨f𝒙)T𝒚=𝒙T𝑨fT𝒚=𝒙T𝑨f𝒚=𝒙T(𝑨f𝒚)=𝒗,f(𝒘),

kde 𝒙 a 𝒚 jsou vektory souřadnic vektorů 𝒗=x1𝒆1++xn𝒆n a 𝒘=y1𝒆1++yn𝒆n.

Projekce a souměrnost

Ortogonální rozklady jsou popsány symetrickými maticemi

Je-li opět (V,,) reálný prostor se skalárním součinem dimenze n a U je jeho k-dimenzionální podprostor, přičemž 𝒙1,,𝒙k jsou vektory ortonormální báze prostoru U, potom matice kolmé projekce na podprostor U je

𝑨U=𝒙1𝒙1T++𝒙k𝒙kTn×n.

Tato matice je symetrická, neboť je dána součtem symetrických matic. Také matice kolmé projekce do ortogonálního doplňku U je díky reprezentaci 𝑨U=𝐈𝑨U vždy symetrická. S pomocí matic projekcí 𝑨U a 𝑨U může být libovolný vektor 𝒗V rozložen na součet vzájemně kolmých vektorů 𝒖U a 𝒖U. Geometrická transformace souměrnosti podle podprostoru U má symetrickou matici 𝐈2𝑨U.

Soustavy lineárních rovnic

Řešení soustavy lineárních rovnic 𝑨𝒙=𝒃 se symetrickou maticí soustavy 𝑨 může být zjednodušeno, pokud se využije symetrie matice 𝑨, konkrétně jejího rozkladu:

𝑨=𝑳𝑫𝑳T

s dolní trojúhelníkovou matricí 𝑳 s jedničkami na diagonále a diagonální maticí 𝑫. Tento rozklad se používá např. při Choleského rozkladu pozitivně-definitivních symetrických matic.

Metody CG a MINRES jsou příklady moderních přístupů pro numerické řešení rozsáhlých soustav lineárních rovnic s řídkou symetrickou maticí soustavy.

Polární rozklad

Každá čtvercová matice 𝑨n×npolární rozklad

𝑨=𝑸𝑷

s ortogonální maticí 𝑸n×n a pozitivní semidefinitní symetrickou maticí 𝑷n×n. Matice 𝑷 je druhá odmocnina z 𝑨T𝑨 . Pokud je 𝑨 regulární, je 𝑷 pozitivně definitní a polární rozklad je pak dán 𝑸=𝑨𝑷1.

Aplikace

Geometrie

Kvadriky lze popsat symetrickými maticemi

Kvadrika v n -rozměrném euklidovském prostoru je množina kořenů kvadratického polynomu v n neznámých. Každou kvadriku lze definovat pomocí nenulové symetrické matice 𝑨n×n, vektoru 𝒃n a absolutního členu c jako množinu bodů

Q={𝒙n𝒙T𝑨𝒙+2𝒃T𝒙+c=0}.

Analýza

Charakterizaci extrémů dvakrát spojitě derivovatelných funkcí f:Dn lze provést pomocí Hessovy matice

Hf(x)=(2fxixj(x))n×n

Hessova matice je podle Schwarzovy věty symetrická. Podle toho, jestli je Hf(x) je pozitivně definitní, negativně definitní nebo indefinitní leží v bodě x lokální minimum, lokální maximum nebo sedlový bod.

Teorie grafů

Neorientovaný hranově­ vážený graf má symetrickou matici sousednosti.

Matice sousednosti 𝑨G neorientovaného hranově váženého grafu G=(V,E,d) s množinou vrcholů V={v1,,vn} je z definice

𝑨G=(aij)()n×n, kde 𝑨ij={d(e)proe={vi,vj}Ejinak

vždy symetrická. Matice odvozené z matice sousednosti součty nebo mocninami, jako například Laplaceova matice, matice sousednosti nebo matice vzdálenosti jsou také symetrické. Analýza těchto matic je předmětem spektrální teorie grafů.

Stochastika

Je-li X=(X1,,Xn) náhodný vektor sestávající z n reálných náhodných veličin X1,,Xn s konečným rozptylem, pak přidružená kovarianční matice

ΣX=(Cov(Xi,Xj))n×n

je matice všech párových kovariancí těchto náhodných veličin. Protože pro všechna i,j{1,,n} platí: Cov(Xi,Xj)=Cov(Xj,Xi), je kovarianční matice symetrická.

Symetrizovatelná matice

Čtvercová matice 𝑨 se nazývá symetrizovatelná, pokud existuje regulární diagonální matice 𝑫 a symetrická matice 𝑺 takové, že 𝑨=𝑫𝑺.

Transpozice symetrizovatelné matice je symetrická, protože 𝑨T=(𝑫𝑺)T=𝑺𝑫=𝑫1(𝑫𝑺𝑫) a 𝑫𝑺𝑫 je symetrická.

Matice 𝑨 je symetrizovatelná, právě když jsou splněny následující podmínky:

  1. aij=0 implikuje aji=0 pro všechna i,j{1,,n} a
  2. ai1i2ai2i3aiki1=ai2i1ai3i2ai1ik pro jakoukoli konečnou posloupnost (i1,i2,,ik).

Odkazy

Reference

Šablona:Překlad

Literatura

Související články

Šablona:Autoritní data