msleziak.com

Tvrdenie.$\newcommand{\Tra}{\operatorname{Tr}}
\newcommand{\inv}[1]{{#1}^{-1}}$ Ak $P$ je matica ortogonálnej projekcie na podpriestor $S$ (v $\mathbb R^n$ so štandardným skalárnym súčinom), tak platí
$$\Tra(P)=\dim(S).$$

Pripomeňme, že $\Tra(P)$ označuje stopu matice $P$, t.j. $\Tra(P)=p_{11}+p_{22}+\dots+p_{nn}$.

Takéto tvrdenie sme spomenuli na cvičení - pri výpočte matice projekcie to je vec, ktorú vieme skontrolovať ľahko. (A ak to nesedí, vieme, že treba hľadať chybu.)

V skutočnosti platí dokonca aj silnejšie tvrdenie - pre každú reálnu (dokonca aj komplexnú) maticu $n\times n$ takú, že $P^2=P$ platí $\Tra(P)=h(P)$.
Niečo o tomto tvrdení sa dá nájsť tu: viewtopic.php?t=1799

Na tomto mieste sa ale zatiaľ pozrime iba na takéto menej všeobecné tvrdenie - že to platí pre maticu ortogonálnej projekcie.

Otázka:$\newcommand{\Tra}{\operatorname{Tr}}
\newcommand{\inv}[1]{{#1}^{-1}}$ Ak $|\vec u|=1$, vedeli by ste zistiť, čomu sa rovná $\Tra(\vec u^T\vec u)$?
Vieme z toho dostať, že stopa matice ortogonálnej projekcie je rovná $\dim(S)$?

Ešte si môžeme ukázať iné zdôvodnenie.$\newcommand{\Tra}{\operatorname{Tr}}\newcommand{\inv}[1]{{#1}^{-1}}$

Vec, ktorá je užitočná aj inde a môže pomôcť aj tu, je to, že pre ľubovoľné štvorcové matice platí
$$\Tra(XY)=\Tra(YX).$$
(Môžete sa zamyslieť nad tým, prečo to je tak - vlastne to je jednoduché úloha na precvičenie definície súčinu matíc.)

Máme $$\vec xP=
\begin{cases}
\vec x & \text{ak }x\in S, \\
\vec 0 & \text{ak }x\in S^\bot.
\end{cases}
$$
Vytvorme bázu $\vec b_1,\dots,\vec b_n$ priestoru $\mathbb R^n$ takú, že $\vec b_1,\dots,\vec b_k$ je ortonormálna báza podpriestoru $S$ a $\vec b_{k+1},\dots,\vec b_n$ je ortonormálna báza podpriestoru $S^\bot$. (Označili sme $k=\dim(S)$.)

Ak týchto $n$ vektorov poukladám ako riadky do matice $B$.
$$BPB^T=D,$$
kde $D=\operatorname{diag}(\underset{\text{$k$-krát}}{\underbrace{1,\dots,1}},\underset{\text{$(n-k)$-krát}}{\underbrace{0,\dots,0}})$, t.j. je to diagonálna matica, kde na diagonále máme $k$ jednotiek a zvyšok núl.
(Skúste si rozmyslieť prečo.) Teraz už zostáva len rozmyslieť si, že platí
\begin{align*}
h(BPB^T)&=h(P)\\
\Tra(BPB^T)&=\Tra(P)
\end{align*}
Z týchto rovností už totiž vidíme, že $\Tra(P)=\Tra(A)=k=h(P)=\dim(S)$.

Dokázať rovnosť pre hodnosti by malo byť ľahké na základe vecí, ktoré vieme z minulého semestra. Ak sa chceme pozrieť na stopu, tak nám tam nejako pomôže fakt, že $\Tra(XY)=\Tra(YX)$.

Pripomeňme, že $BB^T=B^TB=I$. (Z prednášky už vieme, že matice, kde riadky tvoria ortonormálnu bázu, sú presne ortogonálne matice.)

Vedeli by ste s použitím týchto vecí už dokázať rovnosť stôp? Priznám, že takéto riešenie je trochu podvod - v tom zmysle, že keď už budeme mať prebraté nejaké veci o podobnosti matíc (v druhej polovici LS), tak takéto riešenie už bude vyzerať veľmi prirodzene.
Ale v tomto okamihu môže vyzerať veľmi zvláštne, že prečo sme si zobrali akurát maticu $BPB^T$ alebo prečo by sme tam čakali, že vyjde rovnaká stopa.
Snáď ale aj s tým, čo vieme zatiaľ, sa na to dá pozerať tak, že nám niekto prezradil, že máme skúsiť takúto maticu - a potom už overiť tieto veci nie je až také veľmi ťažké.