Symetrická matice

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 \boldsymbol A=\boldsymbol A^{\mathrm 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á. +more 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 \boldsymbol A řádu n nad tělesem T, se nazývá symetrická, pokud pro všechna i, j \in\{ 1, \ldots, n\} platí:

: a_{ij} = a_{ji}.

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

Příklady

Symetrickými maticemi jsou například:

: \begin{pmatrix} 2 \end{pmatrix} ,\quad \begin{pmatrix} 1 & 5 \\ 5 & 7 \end{pmatrix} ,\quad \begin{pmatrix} 1 & 1 \\ 1 & 0 \end{pmatrix} ,\quad \begin{pmatrix} 1 & 3 & 0 \\ 3 & 2 & 6 \\ 0 & 6 & 5 \end{pmatrix} .

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

: \begin{pmatrix} a & b \\ b & c \end{pmatrix} ,\quad \begin{pmatrix} a & b & c \\ b & d & e \\ c & e & f \end{pmatrix} ,\quad \begin{pmatrix} a & b & c & d \\ b & e & f & g \\ c & f & h & i \\ d & g & i & j \end{pmatrix} .

Speciální případy

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

* čtvercová nulová matice obsahuje jen nulové prvky * čtvercová jedničková matice obsahuje jen jednotkové prvky * diagonální matice jsou takové matice, které mají mimo diagonálu jen nulové prvky * Hankelova matice je matice, která je konstantní v rámci diagonál vedoucích seshora zprava doleva dolů * bisymetrická matice je matice osově symetrická podle hlavní diagonály i vedlejší diagonály

Vlastnosti

Popis

U symetrické matice stačí znát pouze prvky na diagonále a pod ní U symetrické matice \boldsymbol A \in T^{n \times n} stačí znát n prvků na diagonále a \tfrac{n(n-1)}{2} prvků na jedné straně diagonály (nad nebo pod). +more Hodnoty prvků na opačné straně diagonály lze odvodit ze symetrie matice. Symetrická matice může mít nejvýše : n + \frac{n(n-1)}{2} = \frac{n(n+1)}{2}.

různých prvků. Ve srovnání s nesymetrickými maticemi řádu n, které mohou mít až n^2 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.

Součet

Součet \boldsymbol A + \boldsymbol B dvou symetrických matic \boldsymbol A, \boldsymbol B \in T^{n \times n} je vždy symetrická matice, protože

: (\boldsymbol A + \boldsymbol B)^{\mathrm T} = \boldsymbol A^{\mathrm T} +\boldsymbol B^{\mathrm T} =\boldsymbol A +\boldsymbol B .

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

\operatorname{Symm}_n = \{\boldsymbol A \in T^{n \times n} \colon \boldsymbol A^{\mathrm T} = \boldsymbol A \}

prostoru čtvercových matic T^{n \times n}. Tento podprostor má dimenzi \tfrac{n^2+n}{2}. +more Jeho bázi lze vytvořit z matic \mathbf E_{ii} pro i\in\{1,\dots,n\}, a součtů \mathbf E_{ij}+\mathbf E_{ji} pro 1 \leq i . Uvedené matice \mathbf E_{ij} tvoří standardní bázi prostoru T^{n \times n} , čili mají jediný nenulový prvek \mathbf{e}_{ij}=1.

Rozklad

Pokud je charakteristika tělesa T různá od 2, lze libovolnou čtvercovou matici \boldsymbol M \in T^{n \times n} zapsat jednoznačně jako součet \boldsymbol M = \boldsymbol A + \boldsymbol B, kde matice \boldsymbol A je symetrická a matice \boldsymbol B je antisymetrická:

:\boldsymbol A = \frac{1}{2}(\boldsymbol M + \boldsymbol M^{\mathrm T}) a \boldsymbol B = \frac{1}{2}(\boldsymbol M - \boldsymbol M^{\mathrm T})

Antisymetrické matice tvoří vektorový podprostor prostoru čtvercových matic. Značí se \operatorname{Skew}_na má dimenzi \tfrac{n^2-n}{2}. +more Prostor čtvercových matic T^{n \times n} dimenze n^2 lze vyjádřit jako direktní součet : T^{n \times n} = \operatorname{Symm}_n \oplus \operatorname{Skew}_n.

prostorů symetrických a antisymetrických matic.

Součin

Součin \boldsymbol{AB} dvou symetrických matic \boldsymbol A,\boldsymbol B \in T^{n\times n} nemusí být opět symetrická matice. Součin symetrických matic je symetrický, právě když je součin \boldsymbol A a \boldsymbol B komutativní. +more Jinými slovy, pokud součin splňuje \boldsymbol{AB} = \boldsymbol{BA}, pak také platí: : (\boldsymbol{AB})^{\mathrm T} =\boldsymbol B^{\mathrm T}\boldsymbol A^{\mathrm T} = \boldsymbol{B A} = \boldsymbol{A B} .

Pro symetrickou matici \boldsymbol A proto platí, že symetrické jsou všechny její mocniny \boldsymbol A^k, kde k \in \N, i její maticová exponenciála e^{\boldsymbol A}.

Pro každou matici \boldsymbol M \in T^{m \times n} jsou matice \boldsymbol{M M}^{\mathrm T} typu m \times m, i matice \boldsymbol M^{\mathrm T} \boldsymbol M typu n \times n symetrické.

Kongruence a podobnost

Každá matice \boldsymbol B \in T^{n \times n}, která je kongruentní symetrické matici \boldsymbol A \in T^{n\times n}, je také symetrická, protože platí : \boldsymbol B^{\mathrm T} = (\boldsymbol S^{\mathrm T} \boldsymbol A \boldsymbol S)^{\mathrm T} = \boldsymbol S^{\mathrm T}\boldsymbol A^{\mathrm T} \boldsymbol S = \boldsymbol S^{\mathrm T}\boldsymbol A \boldsymbol S =\boldsymbol B,

přičemž \boldsymbol S \in T^{n \times 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 \boldsymbol A \in T^{n \times n} regulární, potom matice k ní inverzní \boldsymbol A^{-1} je symetrická, protože pro ni platí: : (\boldsymbol A^{-1})^{\mathrm T} = (\boldsymbol A^{\mathrm T})^{-1} = \boldsymbol A^{-1} .

V tomto případě jsou symetrické všechny mocniny \boldsymbol A^{-k} pro k \in \N.

Reálné symetrické matice

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

Normální matice

Reálná symetrická matice \boldsymbol A \in \R^{n \times n} je normální, protože platí : \boldsymbol A^{\mathrm T} \boldsymbol A = \boldsymbol A\boldsymbol A = \boldsymbol A \boldsymbol A^{\mathrm 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 \R každé číslo shoduje se svým komplexně sdruženým protějškem, neboli z=\overline z, splývají reálné symetrické matice s reálnými hermitovskými. Formálně: : \boldsymbol A^{\mathrm H} = \overline{\boldsymbol A}^{\mathrm T} =\boldsymbol A^{\mathrm T} = \boldsymbol A,

přičemž \boldsymbol A^{\mathrm H} je hermitovská transpozice matice \boldsymbol A a \overline{\boldsymbol A} je komplexně sdružená matice k \boldsymbol A.

Reálná symetrická matice \boldsymbol A \in \R^{n \times n} je vždy hermitovská mimo jiné i proto, že vzhledem k standardnímu skalárnímu součinu \langle \cdot, \cdot \rangle na \mathbb R^n splňuje: : \langle \boldsymbol{Ax}, \boldsymbol y \rangle = (\boldsymbol{Ax})^{\mathrm T}\boldsymbol y =\boldsymbol x^{\mathrm T} \boldsymbol A^{\mathrm T} \boldsymbol y = \boldsymbol x^{\mathrm T} \boldsymbol{ A y} = \langle \boldsymbol x, \boldsymbol{ A y} \rangle

pro všechny vektory \boldsymbol x,\boldsymbol y \in \R^n. Reálné symetrické matice jsou hermitovské i s ohledem na standardní skalární součin nad \C.

Vlastní čísla

Jednotková kružnice je transformována na elipsu pomocí reálné symetrické matice řádu 2. +more Poloosy elipsy odpovídají vlastním vektorům matice a jejich délky vlastním číslům. Vlastní čísla reálné symetrické matice \boldsymbol A \in \R^{n \times n}, tedy řešení rovnice \boldsymbol{ A x} = \lambda \boldsymbol x, jsou vždy reálná. Kdyby \lambda \in \Complex bylo komplexní vlastní číslo matice \boldsymbol A příslušné netriviálnímu vlastnímu vektoru \boldsymbol x \in \C^n, \boldsymbol x \neq \boldsymbol 0, pak z toho, že \boldsymbol A je hermitovská plyne: : \lambda \langle \boldsymbol x, \boldsymbol x \rangle = \langle \boldsymbol x, \lambda \boldsymbol x \rangle = \langle \boldsymbol x, \boldsymbol{Ax} \rangle = \langle \boldsymbol{Ax}, \boldsymbol x \rangle = \langle \lambda \boldsymbol x, \boldsymbol x \rangle = \overline\lambda \langle \boldsymbol x, \boldsymbol x \rangle.

Protože pro každé \boldsymbol x \neq \boldsymbol 0 platí \langle \boldsymbol x, \boldsymbol x \rangle \neq 0, musí vlastní číslo \lambda splňovat \lambda = \overline\lambda, a proto je reálné. V důsledku lze i příslušný vlastní vektor \boldsymbol x zvolit reálný.

Násobnosti vlastních čísel

Pro každou reálnou symetrickou matici \boldsymbol A \in \R^{n \times n} se algebraické a geometrické násobnosti všech vlastních čísel shodují. Důvod je následující. +more Pro vlastní číslo \lambda matice \boldsymbol A s geometrickou násobností k existuje ortonormální báze \{ \boldsymbol x_1, \ldots, \boldsymbol x_k \} prostoru vlastních vektorů příslušných k \lambda. Tuto bázi lze rozšířit pomocí vektorů\{ \boldsymbol x_{k+1}, \ldots,\boldsymbol x_n \} na ortonormální bázi celého prostoru \R^{n}. S pomocí ortogonální matice \boldsymbol S = (\boldsymbol x_1 \mid \cdots \mid \boldsymbol x_n ) je matice \boldsymbol A převedena na podobnou : \boldsymbol C = \boldsymbol S^{-1} \boldsymbol A \boldsymbol S = \boldsymbol S^{\mathrm T} \boldsymbol A \boldsymbol S = \left( \begin{array}{c|c} \lambda \mathbf I & \boldsymbol 0 \\ \hline \boldsymbol 0 & \boldsymbol X \end{array} \right).

což je bloková diagonální matice s bloky \lambda \mathbf I \in \R^{k \times k} a \boldsymbol X \in \R^{(n-k) \times (n-k)}. Vzhledem k tomu, že matice \boldsymbol A je hermitovská a vektory \boldsymbol x_1, \ldots, \boldsymbol x_ntvoří ortonormální bázi, platí pro prvky c_{ij} matice C s indexy \min \{ i,j \} \leq k, že: : c_{ij} = \langle \boldsymbol x_i, \boldsymbol A \boldsymbol x_j \rangle = \langle \boldsymbol A \boldsymbol x_i, \boldsymbol x_j \rangle = \lambda \langle \boldsymbol x_i, \boldsymbol x_j \rangle = \lambda \delta_{ij} ,

kde \delta_{ij} je Kroneckerovo delta. Vektory \boldsymbol x_{k+1}, \ldots, \boldsymbol x_n nejsou podle předpokladu vlastní vektory matice \boldsymbol A příslušné vlastnímu číslu \lambda, proto \lambda není žádným vlastním číslem matice \boldsymbol X. +more Vlastní číslo \lambda matice \boldsymbol C má podle vzorce pro determinant blokových matic shodnou algebraickou i geometrickou násobnost k. Totéž platí i pro matici \boldsymbol A díky vzájemné podobnosti s maticí \boldsymbol C.

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 \boldsymbol A bázi prostoru \R^n. Reálná symetrická matice je tedy vždy diagonalizovatelná, to znamená, že existuje regulární matice \boldsymbol S \in \R^{n \times n} a diagonální matice \boldsymbol D \in \R^{n \times n}splňující: : \boldsymbol S^{-1} \boldsymbol A \boldsymbol S = \boldsymbol D

Matice \boldsymbol S = (\boldsymbol x_1 \mid \cdots \mid \boldsymbol x_n) je sestavena z vlastních vektorů \boldsymbol x_1, \ldots, \boldsymbol x_n po sloupcích a matice \boldsymbol D = \operatorname{diag}(\lambda_1, \ldots, \lambda_n) má vlastní čísla \lambda_1, \ldots, \lambda_n příslušná těmto vlastním vektorům na diagonále. Vzhledem k tomu, že sloupce matice \boldsymbol S, neboli vlastní vektory lze libovolně přerovnat, může být odpovídající pořadí prvků na diagonále \boldsymbol D libovolné. +more 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é. +more 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 \boldsymbol x_i, \boldsymbol x_j příslušné dvěma různým vlastním číslům \lambda_i \neq \lambda_j reálné symetrické matici \boldsymbol A \in \R^{n \times n} jsou vzájemně kolmé. Uvedený vztah opět z následující vlastnosti hermitovských matic \boldsymbol A: : \lambda_i \langle \boldsymbol x_i, \boldsymbol x_j \rangle = \langle \lambda_i \boldsymbol x_i, \boldsymbol x_j \rangle = \langle \boldsymbol A \boldsymbol x_i, \boldsymbol x_j \rangle = \langle \boldsymbol x_i, \boldsymbol A \boldsymbol x_j \rangle = \langle \boldsymbol x_i, \lambda_j \boldsymbol x_j \rangle = \lambda_j \langle \boldsymbol x_i, \boldsymbol x_j \rangle .

Z předpokladu, že \lambda_i a \lambda_j jsou různá, pak plyne \langle \boldsymbol x_i, \boldsymbol x_j \rangle = 0. Vlastní vektory \boldsymbol A tvoří ortonormální bázi prostoru \R^n. +more Každou reálnou symetrickou matici lze proto ortogonálně diagonalizovat, neboli existuje ortogonální matice \boldsymbol S splňující: : \boldsymbol S^{\mathrm T} \boldsymbol A \boldsymbol S = \boldsymbol D.

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

Parametry

Každá reálná symetrická matice \boldsymbol A \in \R^{n \times n} diagonalizovatelná, a proto pro její stopu platí: : \operatorname{tr} \boldsymbol A = \lambda_1 + \ldots + \lambda_n

Její determinant tudíž splňuje:

: \det \boldsymbol A = \lambda_1 \cdot \ldots \cdot \lambda_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 : \operatorname{rank} \boldsymbol A = n - \left( \delta_{\lambda_1,0} + \ldots + \delta_{\lambda_n,0} \right) .

Reálná symetrická matice je regulární, právě když má všechna vlastní čísla nenulová. Spektrální norma reálné symetrické matice je : \| \boldsymbol A \|_2 = \max \{ | \lambda_1 |, \ldots, | \lambda_n | \}

a tedy rovna spektrálnímu poloměru matice. Frobeniova norma vyplývá z normality : \| \boldsymbol A \|_{\mathrm F} = \sqrt{ \lambda_1^2 + \ldots + \lambda_n^2}.

Definitnost

Pro reálnou symetrickou matici \boldsymbol A \in \R^{n \times n} a vektor \boldsymbol x \in \R^n se výraz : Q_{\boldsymbol A}(\boldsymbol x) = \boldsymbol x^{\mathrm T} \boldsymbol{Ax} = \langle \boldsymbol x, \boldsymbol{Ax} \rangle

nazývá kvadratická forma určená maticí \boldsymbol A. Podle toho, jestli je Q_{\boldsymbol A}(\boldsymbol x) pro všechna \boldsymbol x \neq 0 kladná, resp. +more nezáporná, záporná či nekladná, nazývá se matrice \boldsymbol A pozitivně definitní, resp. pozitivně semidefinitní, negativně definitní nebo negativně semidefinitní. Pokud Q_{\boldsymbol A}(\boldsymbol x) nabývá kladných i záporných hodnot, nazývá se matice \boldsymbol A 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é \boldsymbol A \in \R^{n \times n} odhadnout pomocí Rayleighova kvocientu. Konkrétně, pro všechna netriviální \boldsymbol x \in \R^n platí: : \min\{ \lambda_1, \ldots, \lambda_n\} \leq \frac{\langle \boldsymbol x, \boldsymbol {Ax} \rangle}{\langle \boldsymbol x, \boldsymbol x \rangle} \leq \max\{ \lambda_1, \ldots, \lambda_n\}

Rovnost platí, právě když je \boldsymbol x 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 \boldsymbol A, \boldsymbol B \in \R^{n \times n} s vlastními čísly seřazenými sestupně \lambda_1 \geq \ldots \geq \lambda_n a \mu_1 \geq \ldots \geq \mu_n platí odhad : \operatorname{tr}(\boldsymbol { AB}) \leq \lambda_1 \mu_1 + \ldots + \lambda_n \mu_n.

Rovnost je splněna, právě když matice \boldsymbol A a \boldsymbol B jsou současně diagonalizovatelné vzhledem k uspořádání vlastních čísel, neboli když existuje ortogonální matice \boldsymbol S \in \R^{n \times n} taková, že platí \boldsymbol A = \boldsymbol S \operatorname{diag}(\lambda_1, \ldots, \lambda_n) \boldsymbol S^{\mathrm T} a \boldsymbol B = \boldsymbol S \operatorname{diag}(\mu_1, \ldots, \mu_n) \boldsymbol S^{\mathrm T}. Uvedená nerovnost zobecňuje Cauchy-Schwarzovu nerovnost pro Frobeniův skalární součin a permutační nerovnost pro vektory.

Komplexní symetrické matice

Rozklad

Podobně jako u reálných matic lze prostor komplexních čtvercových matic {\mathbb C}^{n \times n} zapsat jako direktní součet prostorů symetrických a antisymetrických matic: : {\mathbb C}^{n \times n} = \operatorname{Symm}_n \oplus \operatorname{Skew}_n

Jde zároveň o ortogonální součet vzhledem k Frobeniově skalárnímu součinu, protože pro všechny matice \boldsymbol A \in \operatorname{Symm}_n a \boldsymbol B \in \operatorname{Skew}_n platí: : \langle \boldsymbol A, \boldsymbol B \rangle_{\mathrm F} = \operatorname{tr}(\boldsymbol A^{\mathrm H}\boldsymbol B) = \operatorname{tr}(\overline{\boldsymbol A}\boldsymbol B) = \operatorname{tr}(\boldsymbol B\overline{\boldsymbol A}) = \operatorname{tr}((\boldsymbol B\overline{\boldsymbol A})^{\mathrm T}) = -\operatorname{tr}(\boldsymbol A^{\mathrm H}\boldsymbol B) = -\langle \boldsymbol A, \boldsymbol B \rangle_{\mathrm F}

z čehož vyplývá \langle \boldsymbol A, \boldsymbol B \rangle_{\mathrm F} = 0. Ortogonalita rozkladu platí i pro reálný maticový prostor {\mathbb R}^{n \times n}.

Spektrum

Pro komplexní matice \boldsymbol A \in \Complex^{n \times n} nemá symetrie žádný zvláštní vliv na spektrum matice. Komplexní symetrická matice může mít nereálná vlastní čísla. +more Například komplexní symetrická matice \boldsymbol A = \begin{pmatrix} 1 & \mathrm i \\ \mathrm i & 1 \end{pmatrix} \in \Complex^{2 \times 2}má dvě vlastní čísla \lambda_{1,2} = 1 \pm \mathrm i .

Existují komplexní symetrické matice, které nelze diagonalizovat. Například matice \boldsymbol A = \begin{pmatrix} 1 & \mathrm i \\ \mathrm i & -1 \end{pmatrix} \in \Complex^{2 \times 2} má jediné vlastní číslo \lambda = 0 s algebraické násobnosti dvě a geometrické násobnosti jedna. +more Obecně platí, že jakákoli komplexní čtvercová matice je podobná komplexní symetrické matici. Spektrum komplexní symetrické matice proto nevykazuje žádné zvláštnosti.

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

Rozklad

Libovolnou komplexní symetrickou matici \boldsymbol A \in \Complex^{n \times n} lze pomocí Autonne-Takagiho faktorizace rozložit na součin

: \boldsymbol A = \boldsymbol U^{\mathrm T} \boldsymbol D \boldsymbol U ,

kde matice \boldsymbol U \in \Complex^{n \times n} je unitární, \boldsymbol D = \operatorname{diag}(\sigma_1, \ldots, \sigma_n) \in \R^{n \times n}je reálná diagonální. Prvky diagonální matice jsou singulární hodnoty \boldsymbol A, neboli odmocniny vlastních čísel matice \boldsymbol A^{\mathrm H}\boldsymbol A.

Použití

Symetrické bilineární formy

Každá bilineární forma b \colon V \times V \to T na vektorovém prostoru V dimenze n nad tělesem T může být vzhledem k bázi \{ \boldsymbol v_1, \ldots, \boldsymbol v_n \} prostoru V reprezentována čtvercovou maticí \boldsymbol A_b \in T^{n \times n} danou vztahem: : (\boldsymbol A_b)_{ij} = b(\boldsymbol v_i, \boldsymbol v_j)

Pokud je bilineární forma symetrická, pak platí b(\boldsymbol v,\boldsymbol w)=b(\boldsymbol w,\boldsymbol v) pro všechny \boldsymbol v,\boldsymbol w \in V, a matice \boldsymbol A_b je symetrická. Naopak každá symetrická matice \boldsymbol A \in T^{n \times n} definuje symetrickou bilineární formu b_{\boldsymbol A} \colon T^n \times T^n \to T vztahem: : b_{\boldsymbol A}(\boldsymbol x,\boldsymbol y) = \boldsymbol{x^{\mathrm T} A y}

Je-li matice \boldsymbol A \in \R^{n \times n} navíc pozitivně definitní, pak b_\boldsymbol{A} představuje skalární součin na euklidovském prostoru \R^n.

Samoadjungované zobrazení

Je-li (V, \langle \cdot, \cdot \rangle) reálný prostor se skalárním součinem dimenze n, pak může být každé lineární zobrazení f \colon V \to V vzhledem k ortonormální bázi \{\boldsymbol e_1, \ldots, \boldsymbol e_n \} prostoru V reprezentováno maticí zobrazení : \boldsymbol A_f = ( a_{ij} ) \in \R^{n \times n},

kde f(\boldsymbol e_j) = a_{1j}\boldsymbol e_1 + \ldots + a_{nj}\boldsymbol e_n pro j=1, \ldots, n. Matice zobrazení \boldsymbol A_f je symetrická, právě když je zobrazení f samoadjungované. +more To vyplývá ze vztahu : \langle f(\boldsymbol v), \boldsymbol w \rangle = (\boldsymbol A_f\boldsymbol x)^{\mathrm T}\boldsymbol y = \boldsymbol x^{\mathrm T}\boldsymbol A_f^{\mathrm T}\boldsymbol y = \boldsymbol x^{\mathrm T}\boldsymbol A_f\boldsymbol y = \boldsymbol x^{\mathrm T}(\boldsymbol A_f\boldsymbol y) = \langle \boldsymbol v, f(\boldsymbol w) \rangle,.

kde \boldsymbol x a \boldsymbol y jsou vektory souřadnic vektorů \boldsymbol v=x_1 \boldsymbol e_1+ \ldots + x_n \boldsymbol e_n a \boldsymbol w=y_1 \boldsymbol e_1 + \ldots + y_n \boldsymbol e_n.

Projekce a souměrnost

Ortogonální rozklady jsou popsány symetrickými maticemi Je-li opět ( V, \langle \cdot, \cdot \rangle ) reálný prostor se skalárním součinem dimenze n a U je jeho k-dimenzionální podprostor, přičemž \boldsymbol x_1, \ldots, \boldsymbol x_k jsou vektory ortonormální báze prostoru U, potom matice kolmé projekce na podprostor U je : \boldsymbol A_U = \boldsymbol x_1 \boldsymbol x_1^{\mathrm T} + \ldots + \boldsymbol x_k \boldsymbol x_k^{\mathrm T} \in \R^{n \times 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^\bot je díky reprezentaci \boldsymbol A_{U^\bot}=\mathbf I-\boldsymbol A_U vždy symetrická. +more S pomocí matic projekcí \boldsymbol A_U a \boldsymbol A_{U^\perp} může být libovolný vektor \boldsymbol v \in V rozložen na součet vzájemně kolmých vektorů \boldsymbol u \in U a \boldsymbol u^\perp \in U^\perp. Geometrická transformace souměrnosti podle podprostoru U má symetrickou matici \mathbf I-2\boldsymbol A_U.

Soustavy lineárních rovnic

Řešení soustavy lineárních rovnic \boldsymbol{Ax}=\boldsymbol b se symetrickou maticí soustavy \boldsymbol A může být zjednodušeno, pokud se využije symetrie matice \boldsymbol A, konkrétně jejího rozkladu: : \boldsymbol A = \boldsymbol{LDL}^{\mathrm T}

s dolní trojúhelníkovou matricí \boldsymbol L s jedničkami na diagonále a diagonální maticí \boldsymbol D. Tento rozklad se používá např. +more 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 \boldsymbol A \in \R^{n \times n} má polární rozklad

: \boldsymbol A = \boldsymbol{QP}

s ortogonální maticí \boldsymbol Q \in \R^{n \times n} a pozitivní semidefinitní symetrickou maticí \boldsymbol P \in \R^{n \times n}. Matice \boldsymbol P je druhá odmocnina z \boldsymbol{A^{\mathrm T}A} . +more Pokud je \boldsymbol A regulární, je \boldsymbol P pozitivně definitní a polární rozklad je pak dán \boldsymbol Q = \boldsymbol{A P}^{-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. +more Každou kvadriku lze definovat pomocí nenulové symetrické matice \boldsymbol A \in \R^{n \times n}, vektoru \boldsymbol b \in \R^n a absolutního členu c \in \R jako množinu bodů.

Q = \left\{ \boldsymbol x \in \R^n \mid \boldsymbol{x^{\mathrm T} A x} + 2 \boldsymbol{b^{\mathrm T} x} + c = 0 \right \}.

Analýza

Charakterizaci extrémů dvakrát spojitě derivovatelných funkcí f \colon D \subset \R^n \to \R lze provést pomocí Hessovy matice

:H_f(x) = \left(\frac{\partial^2f}{\partial x_i\partial x_j}(x)\right) \in \R^{n \times n}

Hessova matice je podle Schwarzovy věty symetrická. Podle toho, jestli je H_f(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. +more Matice sousednosti \boldsymbol A_G neorientovaného hranově váženého grafu G=(V,E,d) s množinou vrcholů V = \{ v_1, \ldots, v_n \} je z definice : \boldsymbol A_G = (a_{ij}) \in(\R\cup\infty)^{n \times n}, kde \boldsymbol A_{ij} = \begin{cases} d(e) & \text{pro} ~ e=\{ v_i,v_j \} \in E \\ \infty & \text{jinak} \end{cases}.

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é. +more Analýza těchto matic je předmětem spektrální teorie grafů.

Stochastika

Je-li X = ( X_1, \ldots, X_n ) náhodný vektor sestávající z n reálných náhodných veličin X_1, \ldots, X_n s konečným rozptylem, pak přidružená kovarianční matice : \Sigma_X = \left( \operatorname{Cov}(X_i, X_j) \right) \in \R^{n \times n}

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

Symetrizovatelná matice

Čtvercová matice \boldsymbol A se nazývá symetrizovatelná, pokud existuje regulární diagonální matice \boldsymbol D a symetrická matice \boldsymbol S takové, že \boldsymbol A = \boldsymbol{DS}.

Transpozice symetrizovatelné matice je symetrická, protože \boldsymbol A^{\mathrm T} =(\boldsymbol{DS})^{\mathrm T} =\boldsymbol{SD}= \boldsymbol D^{-1}(\boldsymbol{DSD}) a \boldsymbol {DSD} je symetrická.

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

# a_{ij} = 0 implikuje a_{ji} = 0 pro všechna i, j \in \{1,\dots, n\} a # a_{i_1 i_2} a_{i_2 i_3} \dots a_{i_k i_1} = a_{i_2 i_1} a_{i_3 i_2} \dots a_{i_1 i_k} pro jakoukoli konečnou posloupnost \left(i_1, i_2, \dots, i_k\right).

Odkazy

Reference

Literatura

Související články

Hermitovská matice * Lineární zobrazení * Podobnost matic * Stopa matice * Vlastní vektory a vlastní čísla

Externí odkazy

Kategorie:Matice