Písomka - ortogonálna podobnosť

Moderator: Martin Sleziak

Post Reply
Martin Sleziak
Posts: 5687
Joined: Mon Jan 02, 2012 5:25 pm

Písomka - ortogonálna podobnosť

Post by Martin Sleziak »

Pre danú symetrickú maticu $A$ nájdite diagonálnu maticu $D$ a ortogonálnu maticu $P$, tak aby platilo $PAP^T=D$.
Líšila sa len zadaná matica:
$$A=
\begin{pmatrix}
2 &-4 & 2 \\
-4 & 2 & 2 \\
2 & 2 & 5
\end{pmatrix}$$
$$B=
\begin{pmatrix}
-1 &-4 & 2 \\
-4 &-1 & 2 \\
2 & 2 & 2
\end{pmatrix}$$
$$C=
\begin{pmatrix}
-1 &-4 & 2 \\
-4 &-1 & 2 \\
2 & 2 & 2
\end{pmatrix}$$
Martin Sleziak
Posts: 5687
Joined: Mon Jan 02, 2012 5:25 pm

Re: Písomka - ortogonálna podobnosť

Post by Martin Sleziak »

Riešenie

Napíšem sem riešenie pre skupinu $A$, ostatné skupiny by sa riešili podobne.

Najprv vypočítajme charakteristický polynóm:
$ch_A(x)=
\begin{vmatrix}
2-x&-4 & 2 \\
-4 &2-x& 2 \\
2 & 2 & 5-x
\end{vmatrix}=$ $
\begin{vmatrix}
6-x&x-6& 0 \\
-4 &2-x& 2 \\
2 & 2 & 5-x
\end{vmatrix}=$ $
(6-x)\begin{vmatrix}
1 &-1 & 0 \\
-4 &2-x& 2 \\
2 & 2 & 5-x
\end{vmatrix}=$ $
(6-x)\begin{vmatrix}
1 &-1 & 0 \\
0 &-2-x& 2 \\
0 & 4 & 5-x
\end{vmatrix}=$ $
(6-x)\begin{vmatrix}
-2-x& 2 \\
4 & 5-x
\end{vmatrix}=$ $
(6-x)[(x+2)(x-5)-8]=$ $
(6-x)(x^2-3x-18)=$ $
(6-x)(x-6)(x+3)=$ $
-(x-6)^2(x+3)$

Zistili sme, že vlastné čísla sú $6$ a $-3$. Už na tomto mieste vieme povedať, že $D=
\begin{pmatrix}
-3 & 0 & 0 \\
0 & 6 & 0 \\
0 & 0 & 6
\end{pmatrix}$.

Ak chceme nájsť maticu takú, že $PAP^{-1}=D$, tak $P$ má ako riadky vlastné vektory. Poďme teda hľadať vlastné vektory. (Potom sa vrátime k tomu, či vieme dostať ortonormálnu bázu z vlastných vektorov - to potrebujeme ak chceme dostať ortogonálnu maticu.)

Vlastný vektor k $-3$ je $(2,2,-1)$ (a všetky jeho nenulové násobky).
Spoiler:
$(A+3I)^T=A+3I=
\begin{pmatrix}
5 &-4 & 2 \\
-4 & 5 & 2 \\
2 & 2 & 8
\end{pmatrix}\sim$ $
\begin{pmatrix}
1 & 1 & 4 \\
-4 & 5 & 2 \\
1 & 1 & 4
\end{pmatrix}\sim$ $
\begin{pmatrix}
1 & 1 & 4 \\
0 & 9 &18 \\
0 & 0 & 0
\end{pmatrix}\sim$ $
\begin{pmatrix}
1 & 1 & 4 \\
0 & 1 & 2 \\
0 & 0 & 0
\end{pmatrix}\sim$ $
\begin{pmatrix}
1 & 0 & 2 \\
0 & 1 & 2 \\
0 & 0 & 0
\end{pmatrix}$
Ešte chcem nájsť vlastné vektory k $6$, čo môžeme urobiť riešením sústavy


$(A-6I)^T=A-6I=
\begin{pmatrix}
-4 &-4 & 2 \\
-4 &-4 & 2 \\
2 & 2 &-1
\end{pmatrix}\sim$ $
\begin{pmatrix}
2 & 2 &-1 \\
0 & 0 & 0 \\
0 & 0 & 0
\end{pmatrix}$

Poznámka: Nemalo by nás prekvapiť, že všetky riešenia tejto sústavy (všetky vlastné vektory k $6$) sú kolmé na vektor $(2,2,-1)$ (vlastný vektor k $-3$). Z prednášky poznáme vetu, že pre symetrickú maticu sú vlastné vektory k rôznym vlastným číslam na seba kolmé.

Odtiaľto vieme (rôznymi spôsobmi) vyjadriť bázu priestoru riešení sústavy, napríklad ako $(1,-1,0)$, $(1,0,2)$.
Chceli by sme dostať takú bázu, kde budú tieto vektory na seba kolmé. Mohli by sme napríklad použiť Gram-Schmidtov ortogonalizačný proces.

Pre dva vektory to asi môžeme spraviť aj jednoduchšie. Ak sme zobrali ako prvý vektor do bázy $(1,-1,0)$, tak ešte chceme nájsť druhý vektor, ktorý je kolmý na $(1,-1,0)$ a súčasne je riešením danej sústavy. Inak povedané, vektor ktorý hľadáme má vyhovovať sústave
$\begin{pmatrix}
2 & 2 &-1 \\
1 &-1 & 0 \\
0 & 0 & 0
\end{pmatrix}$
Túto sústavu ľahko vyriešime, nájdeme, že riešením je napríklad $(1,1,4)$.

Vlastné vektory k $6$ sú $(1,-1,0)$ a $(1,1,4)$; takto zvolené vektory sú navyše na seba kolmé. (A oba sú kolmé na vlastný vektor k $-3$)

Ak ešte každý z vlastných vektorov predelíme veľkosťou, dostaneme tak ortonormálnu bázu pozostávajúcu z vlastných vektorov.
Potom z nich môžeme poskladať ortogonálnu maticu $P$
$P=
\begin{pmatrix}
\frac23 & \frac23 & -\frac13 \\
\frac1{\sqrt2} & -\frac1{\sqrt2} & 0 \\
\frac1{3\sqrt2} & \frac1{3\sqrt2} & \frac4{3\sqrt2}
\end{pmatrix}
$
Pre túto maticu platí $PAP^T=D$.

(Azda by som mal poznamenať, že matica $P$ nie je jednoznačne určená - mohli ste dostať inú ortogonálnu maticu a je to správne riešenie. Matica $D$ je až na poradie prvkov na diagonále určená jednoznačne.)

Tu je kontrola, že to naozaj vychádza. (Aj keď toto by som asi neočakával, že niekto bude rátať naozaj ručne na papier. Na druhú stranu, skontrolovať ručne, či riadky sú naozaj vlastné vektory nie je veľmi ťažké - a tam môžeme použiť vektory pred normalizáciou.)

http://www.wolframalpha.com/input/?i=%5 ... 9%5D%5D%29
http://www.wolframalpha.com/input/?i=%5 ... 9%5D%5D%29

Ostatné skupiny

V ostatných skupinách sa zadané matice líšia len o násobok jednej matice.

Ak $PAP^T=D$, kde $P$ je ortogonálna, tak dostaneme
$$P(A+cI)P^T=PAP^T+cPP^T=D+cI,$$
čiže výsledok je presne rovnaký, treba len k $D$ pripočítať vhodný násobok jednotkovej matice.

Komentáre k riešeniu

Jedno vlastné číslo sa dalo pomerne ľahko uhádnuť.
Vlastné čísla sú tie čísla, pre ktoré $A-\lambda I$ má hodnosť menšiu ako 3. Ak si všimneme, že sa dá nájsť také $\lambda$, kde matica $A-\lambda I$ bude mať dva rovnaké riadky, tak to je určite vlastné číslo.
Konkrétne v skupine A máme
$A-6I=
\begin{pmatrix}
-4 &-4 & 2 \\
-4 &-4 & 2 \\
2 & 2 &-1
\end{pmatrix}$
Táto matica nie je regulárna, teda $6$ je vlastné číslo.
Ak si navyše všimneme, že táto matica má hodnosť $h(A-6I)=1$, tak vieme, že $6$ je dvojnásobné vlastné číslo. Teda diagonálna matica bude vyzerať takto:
$D=
\begin{pmatrix}
d_1 & 0 & 0 \\
0 & 6 & 0 \\
0 & 0 & 6
\end{pmatrix}$
Z toho, že podobné matice majú rovnakú stopu, dostaneme $6+6+d=9$ a $d=-3$.
Čiže v podstate sa dali nájsť vlastné čísla a matica $D$ aj bez rátania charakteristického polynómu, aj keď s trochou hádania.
(Ale určite je prirodzenejší postup cez charakteristický polynóm ako takéto "tipovanie".)

Komentáre k odovzdaným riešeniam, chyby, ktoré sa objavovali

V podstate asi najpodstatnejšie bolo to, že ste niektorí zrátali takú maticu, že $PAP^T=D$ alebo takú, že $PAP^{-1}=D$, ale pre vašu maticu neplatilo $P^T=P^{-1}$. (Nebola ortogonálna.)
Dával som za takéto riešenie (v podstate riešenie inej úlohy) 5 bodov. Niekedy aj viac, ak ste povedzme aspoň napísali, čo by ste robili ďalej, aby ste dostali ortogonálnu maticu, ale už to nezrátali.

Ešte tiež poznamenám, že pri hľadaní koreňov a delení koreňovým činiteľom sa azda oplatilo používať Hornerovu schému.

Objavila sa aj písomka, kde ste zobrali riadky zadanej matice a na ne pustil Gram-Schmidtov proces. Týmto spôsobom rátate niečo úplne iné.

Charakteristický polynóm matice rozmerov $n\times n$ má stupeň $n$. Ak vám pre maticu $3\times3$ vyšiel polynóm stupňa nižšieho ako $3$, treba niekde hľadať chybu.

Ak mi vyšlo, že $\lambda$ je charakteristický polynóm a súčasne som pri riešení sústavy nenašiel žiadne vlastné vektory, treba niekde hľadať chybu. (Buď som sa pomýlili pri hľadaní vlastných vektorov/riešení sústavy, alebo $\lambda$ nie je vlastné číslo a chyba je niekde pri výpočte charakteristického polynómu alebo pri hľadaní jeho koreňov.)

Riadky ortogonálnej matice sú navzájom kolmé vektory dĺžky jedna. Niektorí z vás našli vlastné vektory, ktoré boli na seba kolmé - ale zabudli ste ich predeliť veľkosťou.

Obtiažnosť úlohy

Určite som ochotný uznať, že ak ste rátali charakteristický polynóm pomocou Sarrusovho pravidla, tak tam vychádzali pomerne veľké čísla (koeficienty niekde v rozsahu medzi 50 a 100). Takže sa tam dalo pomerne ľahko pomýliť.
A aj ak ho zrátam a nepomýlim sa, tak stále zostáva problém nájsť korene.
Aj tak si však myslím, že táto úloha bola zrátateľná v pohode pod 30 minút - ak ste zrátali charakteristický polynóm napríklad tak, ako som uviedol vyššie, tak už sa výpočet charakteristického polynómu výrazne zjednoduší a navyše hneď máte jeden koreň. (Stačí hľadať už len korene kvadratického polynómu.)
Viackrát, keď sme na cviku rátali charakteristický polynóm, ste mohli vidieť, že práve takýto postup nám mohol zjednodušiť život.

Ak ste rátali charakteristický polynóm cez Sarrusa, tak je to ťažšie. Ale snáď aj tu je istá výhoda (a pomôže nájsť chybu) to, že jeden koreň hneď vidíte. (Počas písomky som práve kvôli tomuto hovoril hint o tom, že sa treba pozrieť, či pre nejaké $\lambda$ nebude $A-\lambda I$ matica, ktorá má dva riadky rovnaké.)

Klobúk dolu, že ste mnohí zvládli zrátať charakteristický polynóm aj takýmto spôsobom a nepomýliť sa. Nemyslím si však, že by som zadal neprimerane ťažkú písomku - ako som ukázal vyššie, vlastné čísla sa dali nájsť aj oveľa rýchlejším spôsobom.
Martin Sleziak
Posts: 5687
Joined: Mon Jan 02, 2012 5:25 pm

Re: Písomka - ortogonálna podobnosť

Post by Martin Sleziak »

Pridám sem nejaké linky týkajúce sa hľadania koreňov charakteristického polynómu: viewtopic.php?t=890

Ak hľadáte korene polynómu vyššieho stupňa, oplatí sa pripomenúť si, čo viete o racionálnych koreňoch polynómov s celočíselnými koeficientami: viewtopic.php?t=1091
Martin Sleziak
Posts: 5687
Joined: Mon Jan 02, 2012 5:25 pm

Re: Písomka - ortogonálna podobnosť

Post by Martin Sleziak »

Aj tento semester bolo presne rovnaké zadanie (s inými maticami).

Napíšem sem len to, aké vyšli vlastné čísla a ako vyzerajú vlastné vektory (už kolmé). Na zostavenie matice $P$ ich treba ešte vynormovať (predeliť veľkosťou)

$$A=
\begin{pmatrix}
0 & 1 & 2 \\
1 & 0 & 2 \\
2 & 2 & 3 \\
\end{pmatrix}
$$
Vlastné hodnoty sú $-1$, $-1$, $-5$. Vlastné vektory k $-1$: $[(1,1,-1),(1,-1,0)]$. Vlastné vektory k $5$: $(1,1,2)$.

$$A=
\begin{pmatrix}
1 & 2 & 2 \\
2 & 1 & 2 \\
2 & 2 & 1 \\
\end{pmatrix}
$$
Vlastné hodnoty $-1$, $-1$, $5$. Vlastné vektory k $-1$: $[(1,1,-2),(1,-1,0)]$. Vlastné vektory k $5$: $[(1,1,1)]$.

$$A=
\begin{pmatrix}
3 & 1 & 1 \\
1 & 3 & 1 \\
1 & 1 & 3
\end{pmatrix}$$
Vlastné čísla $2$, $2$, $5$. Vlastné vektory k $2$: $[(1,-1,0),(1,1,-2)]$. Vlastné vektory k $5$: $[(1,1,1)]$.

$$A=
\begin{pmatrix}
1 & 2 & 2 \\
2 & 1 & 2 \\
2 & 2 & 1
\end{pmatrix}$$
Vlastné čísla: $-1$, $-1$, $5$. Vlastné vektory k $-1$: $[(1,-1,0),(1,1,-2)]$; vlastné vektory k $5$: $[(1,1,1)]$.
Martin Sleziak
Posts: 5687
Joined: Mon Jan 02, 2012 5:25 pm

Re: Písomka - ortogonálna podobnosť

Post by Martin Sleziak »

Pre danú maticu $A$ nájdite ortogonálnu maticu $P$ a diagonálnu maticu $D$ tak, aby platilo $PAP^T=D$. (Ako výsledok chcem obe matice -- ortogonálnu maticu $P$ aj diagonálnu maticu $D$.)
$$A=
\begin{pmatrix}
2 & 1 &-1 \\
1 & 2 &-1 \\
-1 &-1 & 2 \\
\end{pmatrix}
$$
Teda vlastne to čo chceme urobiť je:
  • Nájsť vlastné hodnoty.
  • Nájsť k nim vlastné vektory.
  • Ak je niektorá vlastná hodnota viacnásobná, tak ešte potrebujeme ortogonalizovať.
  • Ak už máme ortonormálnu bázu pozostávajúcu z vlastných vektorov, tak z nich môžeme zostaviť maticu $P$.
Vlastné hodnoty.

$\chi_A(t)=\det(A-tI)=
\begin{vmatrix}
2-t& 1 &-1 \\
1 &2-t&-1 \\
-1 &-1 &2-t\\
\end{vmatrix}=-(t-1)^2(t-4)$
Vlastné hodnoty sú $1$, $1$, $4$.
Spoiler:
$\det(A-tI)=
\begin{vmatrix}
2-t& 1 &-1 \\
1 &2-t&-1 \\
-1 &-1 &2-t\\
\end{vmatrix}=$ $
\begin{vmatrix}
1-t&t-1& 0 \\
1 &2-t&-1 \\
-1 &-1 &2-t\\
\end{vmatrix}=$ $
(1-t)\begin{vmatrix}
1 &-1 & 0 \\
1 &2-t&-1 \\
-1 &-1 &2-t\\
\end{vmatrix}=$ $
(1-t)\begin{vmatrix}
1 &-1 & 0 \\
0 &3-t&-1 \\
0 &-2 &2-t\\
\end{vmatrix}=$ $
(1-t)\begin{vmatrix}
3-t&-1 \\
-2 &2-t\\
\end{vmatrix}=$ $
(1-t)[(t-3)(t-2)-2]=$ $
(1-t)(t^2-5t+4)=$ $
(1-t)(t-1)(t-4)
$
Vlastné vektory k $4$: Vlastné vektory k vlastnej hodnote $\lambda$ dostaneme riešením homogénnej sústavy s maticou $(A-\lambda I)^T$. (Pretože $A$ je symetrická, máme $(A-\lambda I)^T=(A-\lambda I)$.)

$A-4I=\begin{pmatrix}
-2 & 1 &-1 \\
1 &-2 &-1 \\
-1 &-1 &-2 \\
\end{pmatrix}\sim\dots\sim
\begin{pmatrix}
1 & 0 & 1 \\
0 & 1 & 1 \\
0 & 0 & 0 \\
\end{pmatrix}$

Množina riešení tejto sústavy je $[(1,1,-1)]$, každý nenulový vektor z tejto množiny je vlastný vektor k vlastnému číslu $4$.
Spoiler:
$A-4I=\begin{pmatrix}
-2 & 1 &-1 \\
1 &-2 &-1 \\
-1 &-1 &-2 \\
\end{pmatrix}\sim$ $
\begin{pmatrix}
-3 & 3 & 0 \\
1 &-2 &-1 \\
-1 &-1 &-2 \\
\end{pmatrix}\sim$ $
\begin{pmatrix}
1 &-1 & 0 \\
1 &-2 &-1 \\
-1 &-1 &-2 \\
\end{pmatrix}\sim$ $
\begin{pmatrix}
1 &-1 & 0 \\
0 &-1 &-1 \\
0 &-2 &-2 \\
\end{pmatrix}\sim$ $
\begin{pmatrix}
1 &-1 & 0 \\
0 & 1 & 1 \\
0 & 0 & 0 \\
\end{pmatrix}\sim$ $
\begin{pmatrix}
1 & 0 & 1 \\
0 & 1 & 1 \\
0 & 0 & 0 \\
\end{pmatrix}
$
Ortogonálne vlastné vektory k $1$: Vieme, že vlastné vektory k rôznym vlastným hodnotám sú na seba kolmé. Tu máme iba dve vlastné hodnoty, už sme vypočítali vlastné vektory k jeden z nich.
Stačí teda zobrať vektory kolmé na $[(1,1,-1)]$, t.j. $[(1,0,1),(0,1,1)]$. (Môžete si vyskúšať, že presne to isté dostaneme riešením homogénnej sústavy s maticou $A-I$)
Z nich môžeme vytvoriť rôzne ortogonálne bázy, napríklad riadky týchto matíc dajú ortogonálnu bázu tohoto podpriestoru:
$\begin{pmatrix}
1 & 1 & 2 \\
1 &-1 & 0 \\
\end{pmatrix}\sim$ $
\begin{pmatrix}
1 & 0 & 1 \\
1 &-2 &-1 \\
\end{pmatrix}\sim$ $
\begin{pmatrix}
2 &-1 & 1 \\
0 & 1 & 1 \\
\end{pmatrix}$.
(Na hľadanie ortogonálnej resp. ortonormálnej bázy sme videli viacero postupov, sem som nepísal celý výpočet.)

Matice $P$ a $D$.
Matica $D$ obsahuje vlastné čísla (toľkokrát, koľko je ich násobnosť):
$$D=\operatorname{diag}(1,1,4)=
\begin{pmatrix}
1 & 0 & 0 \\
0 & 1 & 0 \\
0 & 0 & 4 \\
\end{pmatrix}
$$
Ortogonálnu maticu $P$ zostavíme tak, že to nej poukladáme ako riadky ortonormálne vlastné vektory. (V takom poradí, v akom sme dali do diagonálnej matice $D$ vlastné hodnoty.)
$$P=
\begin{pmatrix}
\frac2{\sqrt6} &-\frac1{\sqrt6} & \frac1{\sqrt6} \\
0 & \frac1{\sqrt2} & \frac1{\sqrt2} \\
\frac1{\sqrt3} & \frac1{\sqrt3} &-\frac1{\sqrt3} \\
\end{pmatrix}
$$

Môžeme overiť (či už manuálnym výpočtom alebo pomocou nejakého softvéru), že skutočne platí $PP^T=I$ a $PAP^T=D$.
Spoiler:
$PP^T=
\begin{pmatrix}
\frac2{\sqrt6} &-\frac1{\sqrt6} & \frac1{\sqrt6} \\
0 & \frac1{\sqrt2} & \frac1{\sqrt2} \\
\frac1{\sqrt3} & \frac1{\sqrt3} &-\frac1{\sqrt3} \\
\end{pmatrix}
\begin{pmatrix}
\frac2{\sqrt6} & 0 & \frac1{\sqrt3} \\
-\frac1{\sqrt6} & \frac1{\sqrt2} & \frac1{\sqrt3} \\
\frac1{\sqrt6} & \frac1{\sqrt2} &-\frac1{\sqrt3} \\
\end{pmatrix}=$ $
\begin{pmatrix}
1 & 0 & 0 \\
0 & 1 & 0 \\
0 & 0 & 1 \\
\end{pmatrix}$

$PAP^T=
\begin{pmatrix}
\frac2{\sqrt6} &-\frac1{\sqrt6} & \frac1{\sqrt6} \\
0 & \frac1{\sqrt2} & \frac1{\sqrt2} \\
\frac1{\sqrt3} & \frac1{\sqrt3} &-\frac1{\sqrt3} \\
\end{pmatrix}
\begin{pmatrix}
2 & 1 &-1 \\
1 & 2 &-1 \\
-1 &-1 & 2 \\
\end{pmatrix}
\begin{pmatrix}
\frac2{\sqrt6} & 0 & \frac1{\sqrt3} \\
-\frac1{\sqrt6} & \frac1{\sqrt2} & \frac1{\sqrt3} \\
\frac1{\sqrt6} & \frac1{\sqrt2} &-\frac1{\sqrt3} \\
\end{pmatrix}=$ $
\begin{pmatrix}
\frac2{\sqrt6} &-\frac1{\sqrt6} & \frac1{\sqrt6} \\
0 & \frac1{\sqrt2} & \frac1{\sqrt2} \\
\frac1{\sqrt3} & \frac1{\sqrt3} &-\frac1{\sqrt3} \\
\end{pmatrix}
\begin{pmatrix}
\frac2{\sqrt6} & 0 & \frac1{\sqrt3} \\
-\frac1{\sqrt6} & \frac1{\sqrt2} & \frac1{\sqrt3} \\
\frac1{\sqrt6} & \frac1{\sqrt2} &-\frac1{\sqrt3} \\
\end{pmatrix}
\begin{pmatrix}
\frac2{\sqrt6} &-\frac1{\sqrt6} & \frac1{\sqrt6} \\
0 & \frac1{\sqrt2} & \frac1{\sqrt2} \\
\frac4{\sqrt3} & \frac4{\sqrt3} &-\frac4{\sqrt3} \\
\end{pmatrix}=$ $
\begin{pmatrix}
\frac2{\sqrt6} & 0 &-\frac4{\sqrt3} \\
-\frac1{\sqrt6} & \frac1{\sqrt2} & \frac4{\sqrt3} \\
\frac1{\sqrt6} & \frac1{\sqrt2} &-\frac4{\sqrt3} \\
\end{pmatrix}=$ $
\begin{pmatrix}
1 & 0 & 0 \\
0 & 1 & 0 \\
0 & 0 & 4 \\
\end{pmatrix}$
Chyby, ktoré sa vyskytli v odovzdaných riešeniach

Niektorí ste postupovali ako pri úprave na kanonický tvar. Tak skutočne viete dostať nejaké matice $P$ a $D$, pre ktoré platí $PAP^T=D$. Nemáte však zaručené, že takýmto postupom dostanete ortogonálnu maticu.
Post Reply