Hola
Como 3-4 veces me hablan de la matriz de Gram en el foro
no lo vi en el teórico, hay otra forma de probarlo? intente usando ese teorema que dice que una matriz es ortogonal si las filas forman una base ortogonal, pero no logre nada, lo que se podría usar es que en una c.lineal de un vector respecto a la base ortonormal puedo calcular las coordenadas fácilmente con el pi, pero no logro que eso sea una base ortonormal 
Independientemente de como te lo demuestren, que una matriz de cambio de base entre bases ortonormales es una matriz ortogonal, es decir, su inversa coincide con su traspuesta, es una propiedad muy importante y en ningún curso sobre espacios vectoriales con producto interno debería de pasarse por alto.
Una prueba sin usar la matriz de Gram en espacios vectoriales reales:
- Dada una base ortonormal \( \{u_1,u_2,\ldots,u_n\} \) supongo que sabes que un vector \( v \) respecto a esa base puede escribirse como:
\( v=<v,u_1>u_1+<v,u_2>u_2+\ldots+<v,u_n>u_n \)
Spoiler
Es bastante inmediato de probar.
Si \( v=a_1u_1+a_2u_2+\ldots+a_nu_n \) haciendo a ambos lados el producto escalar por \( u_i \):
\( <v,u_i>=a_1<u_1,u_i>+a_2<u_2,u_i>+\ldots+a_n<u_n,u_n> \)
Y por ser la base ortonormal esos productos \( <u_j,u_i> \) son nulos excepto si \( i=j \) en cuyo caso vale \( 1 \) y por eso \( <v-u_i>=a_i \).
- Dadas dos bases \( B=\{u_1,u_2,\ldots,u_n\} \) y \( B'=\{v_1,v_2,\ldots,v_n\} \) la matriz de cambio de base de la segunda a la primera es aquella que tiene como columnas las coordenadas de los vectores de la segunda base con respecto a la primera. Pero por ser ésta ortonormal y por lo que dije antes:
\( v_1=<v_1,u_1>u_1+<v_1,u_2>u_2+\ldots+<v_1,u_n>u_n \)
\( v_2=<v_2,u_1>u_1+<v_2,u_2>u_2+\ldots+<v_2,u_n>u_n \)
\( \ldots \)
\( v_n=<v_n,u_1>u_1+<v_n,u_2>u_2+\ldots+<v_n,u_n>u_n \)
y por tanto la matriz de cambio de base:
\( M=\begin{pmatrix}
<v_1,u_1>&<v_2,u_1>&\ldots&<v_n,u_1>\\
<v_1,u_2>&<v_2,u_2>&\ldots&<v_n,u_2>\\
\vdots&\vdots&\ddots&\vdots\\
<v_1,u_n>&<v_2,u_n>&\ldots&<v_n,u_n>\\
\end{pmatrix} \)
- Por otra parte es sabido que la matriz de cambio de base en sentido opuesto, es decir, de la primera base a la segunda es la inversa de la anterior. Pero también por definición es aquella que tiene como columnas las coordenadas de los vectores de la primera base con respecto a la segunda. Como la segunda también es ortogonal:
\( u_1=<u_1,v_1>v_1+<u_1,v_2>v_2+\ldots+<u_1,v_n>v_n \)
\( u_2=<u_2,v_1>v_1+<u_2,v_2>v_2+\ldots+<u_2,v_n>v_n \)
\( \ldots \)
\( u_n=<u_n,v_1>v_1+<u_n,v_2>v_2+\ldots+<u_n,v_n>v_n \)
y por tanto:
\( M^{-1}=\begin{pmatrix}
<u_1,v_1>&<u_2,v_1>&\ldots&<u_n,v_1>\\
<u_1,v_2>&<u_2,v_2>&\ldots&<u_n,v_2>\\
\vdots&\vdots&\ddots&\vdots\\
<u_1,v_n>&<u_2,v_n>&\ldots&<u_n,v_n>\\
\end{pmatrix}=\begin{pmatrix}
<v_1,u_1>&<v_1,u_2>&\ldots&<v_1,u_n>\\
<v_2,u_1>&<v_2,u_2>&\ldots&<v_2,u_n>\\
\vdots&\vdots&\ddots&\vdots\\
<v_n,u_1>&<v_n,u_2>&\ldots&<v_n,u_n>\\
\end{pmatrix}=M^t \)
donde en la igualdad intermedia hemos usado que el producto escalar real es conmutativo \( <u,v>=<v,u> \).
Saludos.