Una solución un tanto escurridiza: Dejemos que los subespacios sean $U$ y $V$ dentro del espacio euclidiano $E$ . Sea $\hat{P}:E \rightarrow V$ denotan la proyección, y $P:U \rightarrow V$ denota la restricción de $\hat{P}$ a $U$ . Entonces $A = P^{T}P:U \rightarrow U$ es simétrica, por lo que se puede formar una base ortonormal $e_1,\ldots,e_k$ de $U$ formado por los vectores propios de $A$ . Ordenar la base de manera que los vectores en el espacio nulo de $A$ vienen en último lugar. (nota que $Nul(A) = Nul(P)$ )
Entonces defina $f_i = Pe_i$ para cada uno de los $e_i$ tal que $Pe_i$ son distintos de cero. Tenemos $\langle f_i, f_j \rangle = \langle Pe_i,Pe_j \rangle = \langle e_i, P^{T}Pe_j \rangle = \langle e_i, Ae_j \rangle = \lambda_j \langle e_i, e_j \rangle$ De ahí que el $f_i$ son ortogonales entre sí, y podemos completar esto con una base ortogonal de $V$ .
Queda por demostrar que $\langle e_i, f_i \rangle \geq 0$ y $\langle e_i, f_j \rangle = 0$ si $i \neq j$ .
Si $Pe_i = 0$ entonces $e_i$ es ortogonal a cada vector en $V$ y así $\langle e_i, f_j \rangle = 0$ para todos $j$ .
En caso contrario, tenemos $e_i = f_i + \tilde{f}_i$ , donde $\tilde{f}_i$ es ortogonal a $V$ . Entonces $\langle e_i, f_j \rangle = \langle f_i + \tilde{f}_i, f_j \rangle = \langle f_i, f_j \rangle$ y el resultado se deduce ya que el $f$ forman una base ortogonal.
Por último, normalice el $f$ 's.