$\newcommand{\metric}[2]{\langle #1,#2 \rangle}$$ \newcommand{\R}{\mathbb{R}}$ He aquí otra prueba. El argumento es un ajuste de la prueba de la descomposición del valor singular de una matriz. No pretendo que sea la prueba más corta de la historia, pero creo que es bastante informativa.
El punto de vista geométrico. Una simétrica $n\times n$ -matriz $A$ da lugar a una forma bilineal en $\mathbb{R}^n$ : $$ \metric{v}{w} = w^T A v. $$ Asociada a esta forma está la forma cuadrática $Q(v)=\metric{v}{v}=v^T A v$ .
Consideremos esta forma bilineal como un producto interior generalizado sobre $\R^n$ . En este sentido, piense en $\metric{v}{v}$ como la "longitud al cuadrado" de un vector $v\in\R^n$ . Si todos los valores propios de $A$ son positivos, entonces $\metric{v}{v}$ es siempre un número positivo; en este caso $\metric{v}{w}$ es un producto interno. Si $A$ tiene valores propios negativos de cero, este punto de vista es menos intuitivo. No obstante, estas formas bilineales se dan de forma natural (por ejemplo, el espaciotiempo de Minkowski en $\R^4$ ).
Surge una pregunta natural. ¿Cómo se comportan las "distancias" en $\R^n$ equipado con esta forma bilineal? Para ello, investigamos el mapa $$ f: S^{n-1}\subset\R^n\to \R: v \mapsto v^T A v. $$ Vamos a proceder como sigue. Primero determinamos el vector $v_1 \in S^{n-1}$ tal que $f(v_1)$ es máxima. Entonces vamos a encontrar un vector $v_2 \in S^{n-1} \cap \mathrm{span}\{v_1\}^\perp$ tal que $f(v_2)$ es maximal (en el complemento ortogonal de $v_1$ ). Repitiendo este argumento obtenemos una base ortonormal $\{v_1,\ldots, v_n\}$ con valores propios $f(v_1),\ldots, f(v_n)$ . Como la función tiene valor real, los valores propios son reales.
Teorema. La matriz $A$ admite una base ortonormal $\{v_1,\ldots, v_n\}$ de vectores propios con valores propios reales.
Prueba. Paso básico. Desde $S^{n-1}$ es compacto, existe un $v_1\in S^{n-1}$ tal que $f(v_1)$ es máxima. Ahora tomemos a $w \in S^{n-1}\cap \mathrm{span}\{v_1\}^\perp$ . Consideremos la curva $$ \alpha\colon (-\epsilon,\epsilon)\to\R^n : t \mapsto \cos t\, v_1 + \sin t\, w. $$ Tenga en cuenta que $\alpha(0)=v_1$ y $\alpha'(0)=w$ . Consideremos ahora la composición $g(t)=f(\alpha(t))$ . Su derivada es $$ \begin{align*} g'(t) &= (-\sin t\, v_1 + \cos t \, w)^T A (\cos t\, v_1 + \sin t\,w) \\ & \qquad + (\cos t\, v_1 + \sin t\,w)^T A (-\sin t\, v_1 + \cos t \, w) \\ &= 2 (-\sin t\, v_1 + \cos t \, w)^T A (\cos t\, v_1 + \sin t\,w) . \end{align*} $$ Aquí utilizamos la simetría de $A$ . Tenga en cuenta que $g$ es máxima en $t=0$ desde $\alpha(0)=v_1$ . Por lo tanto $g'(0)=2 w^T A v_1$ debe ser cero. Dado que $w$ es un vector unitario arbitrario perpendicular a $v_1$ , $Av_1$ debe ser múltiplo de $v_1$ . Concluimos que $v_1$ es un valor propio con valor propio real $f(v_1)$ .
Paso de inducción. Supongamos que ya hemos encontrado un conjunto ortonormal $\{v_1,\ldots, v_k\}$ de valores propios con valores propios reales. Entonces simplemente hacemos el mismo argumento que arriba, pero lo aplicamos a la función $f$ restringida a la subesfera $S^{n-1}\cap \mathrm{span}\{v_1,\ldots,v_k\}^\perp$ .