Veamos primero esto en dos dimensiones. Después de esto, deberías saber cómo hacerlo en tres.
Supongamos que se trabaja en el plano y se tienen dos vectores linealmente independientes $v$ y $w$ . Usted quiere hacer $v$ y $w$ ortogonales entre sí en términos del producto interior euclidiano estándar. ¿Cómo se puede hacer?
Pues fíjate que podemos restar un determinado múltiplo de $w$ de $v$ Llamémoslo $cw$ donde $c$ es alguna constante para que $v - cw$ será ortogonal a $w$ . En otras palabras, después de algunas manipulaciones algebraicas encontramos que $c$ debe ser igual a
$$\frac{\langle v,w \rangle }{\langle w,w \rangle}$$
simplemente resolviendo la ecuación $\langle v - cw, w \rangle = 0$ para $c$ . Entonces ahora tendrás una base ortogonal para $\mathbb{R}^2$ , es decir, los vectores
$$w \quad \text{and} \quad v - \frac{\langle v,w \rangle }{\langle w,w \rangle} w.$$
Para encontrar una base ortonormal, basta con dividir por la longitud de cada uno de los vectores.
Sur $\mathbb{R}^3$ sólo tienes que aplicar este proceso de forma recursiva como se muestra en el enlace de la wikipedia en los comentarios anteriores. Sin embargo, primero tienes que comprobar que tus vectores son linealmente independientes. Puedes comprobarlo calculando el determinante de la matriz cuyas columnas son los vectores que has indicado en tu pregunta.