Orthonormalizing a Given Basis
Introduction
Abstract We introduce the Gram-Schmidt procedure which turns any basis of a subspace \(V\) of \(\RNrSpc{n}\) into a orthonormal basis of \(V\).
Outline In the previous sections we saw how convenient it can be to work with a basis of orthonormal vectors vs. an arbitrary basis of a vector space \(V\). It is therefore useful to know that an arbitrary ordered basis can always be transformed into an orthonormal one. The procedure which accomplishes this feat is called Gram-Schmidt orthonormalization, named after Jorgen Pedersen Gram and Eduard Schmidt. It had, however, been used earlier in specialized contexts. So, here goes:
If \(\mathcal{A}=(\Vect{a}_1,\dots ,\Vect{a}_n)\) is an arbitrary ordered set of linearly independent vectors in \(\RNrSpc{n}\), let \(V\) denote its span. The first step toward transforming \(\mathcal{A}\) into an ONB of \(V\) is just normalizing \(\Vect{a}_1\):
\[\Vect{v}_1 \DefEq \dfrac{\Vect{a}_1}{ \Norm{ \Vect{a}_1 } }\]
The second step is to split \(\Vect{a}_2 = \Vect{b}_2 + \Vect{u}_2\), with \(\Vect{b}_2\in \SpanOfSet{ \Vect{v}_1 }\) and \(\Vect{u}_2 \in \SpanOfSet{ \Vect{v}_1 }^{\bot}\). Then \(\Vect{b}_2\) contributes nothing new to \(\SpanOfSet{\Vect{v}_1}\); so forget about it. However, \(\Vect{v}_1\), and \(\Vect{v}_2\DefEq \frac{\Vect{u}_2}{ \Norm{\Vect{u}_2} }\) are orthonormal vectors by design, and they satisfy
\[\SpanOfSet{\Vect{a}_1,\Vect{a}_2} = \SpanOfSet{\Vect{v}_1,\Vect{v}_2}\]
Next we split \(\Vect{a}_3 = \Vect{b}_3 + \Vect{u}_3\), with \(\Vect{b}_3\in \SpanOfSet{ \Vect{v}_1,\Vect{v}_2 }\) and \(\Vect{u}_3 \in \SpanOfSet{ \Vect{v}_1,\Vect{v}_2 }^{\bot}\). Then \(\Vect{b}_3\) contributes nothing new to \(\SpanOfSet{\Vect{v}_1, \Vect{v}_2}\); so forget about it. However, \(\Vect{v}_1,\Vect{v}_2\), and \(\Vect{v}_3\DefEq \frac{\Vect{u}_3}{ \Norm{\Vect{u}_3} }\) are orthonormal vectors by design, and they satisfy
\[\SpanOfSet{\Vect{a}_1,\Vect{a}_2,\Vect{a}_3} = \SpanOfSet{\Vect{v}_1,\Vect{v}_2,\Vect{v}_3}\]
We continue in this manner until we have processed all the vectors in \(\mathcal{A}\), and the result is an ONB \(\Vect{v}_1,\dots ,\Vect{v}_r\) for \(V\). Of course, at some point we must
check that this procedure actually accomplishes what it claims to do.
As a consequence of the Gram-Schmidt orthonormalization procedure, we see that every orthonormal set in a subspace \(V\) or \(\RNrSpc{n}\)
can be complemented to an ONB of \(V\).
TheoremGram-Schmidt orthonormalization
Given an ordered linearly independent set \(\EuScript{A}\DefEq (\Vect{a}_1,\dots ,\Vect{a}_r)\) of vectors in \(\RNrSpc{n}\), the ordered set of vectors \(\EuScript{B}\DefEq (\Vect{v}_1,\dots ,\Vect{v}_r)\) defined below is an ordered ONB of \(\span(\EuScript{A})\).
\(\Vect{v}_1\) | \(\DefEq \) | \(\dfrac{ \Vect{a}_1 }{ \Norm{ \Vect{a}_1} }\) |
\(\Vect{v}_2\) | \(\DefEq \) | \(\dfrac{ \Vect{a}_2 - (\DotPr{ \Vect{a}_2 }{ \Vect{v}_1 })\Vect{v}_1 }{ \Norm{ \Vect{a}_2 - (\DotPr{ \Vect{a}_2 }{ \Vect{v}_1 })\Vect{v}_1} }\) |
\(\vdots\) | \(\) | \(\qquad \vdots \qquad\qquad \vdots\) |
\(\Vect{v}_r\) | \(\DefEq \) | \(\dfrac{ \Vect{a}_r - (\DotPr{ \Vect{a}_r }{ \Vect{v}_1 })\Vect{v}_1 - (\DotPr{ \Vect{a}_r }{ \Vect{v}_2 })\Vect{v}_2 - \cdots - (\DotPr{ \Vect{a}_r }{ \Vect{v}_{r-1} })\Vect{v}_{r-1} }{ \Norm{ \Vect{a}_r - (\DotPr{ \Vect{a}_r }{ \Vect{v}_1 })\Vect{v}_1 - (\DotPr{ \Vect{a}_r }{ \Vect{v}_2 })\Vect{v}_2 - \cdots - (\DotPr{ \Vect{a}_r }{ \Vect{v}_{r-1} })\Vect{v}_{r-1} } }\) |
Moreover, \(\span\Set{ \Vect{a}_1,\dots ,\Vect{a}_j } = \span\Set{ \Vect{v}_1,\dots ,\Vect{v}_j }\) for each \(1\leq j\leq r\) and, if \(\Vect{a}_1,\dots ,\Vect{a}_j\) are already orthonormal, then \(\Vect{a}_k=\Vect{v}_k\) for each \(1\leq k\leq j\).
The Gram-Schmidt orthonormalization theorem has the following immediate consequence
CorollaryExistence of an ONB
In a nonzero subvector space \(V\) of \(\RNrSpc{n}\) every orthonormal subset \(S\) of \(V\) can be complemented to an orthonormal basis of \(V\).
http://emath.ualberta.ca/LinrAlgbra/LinAlgInRn.xml/pages/Sec_LeastSquaresOrthonormalizing a Given Basishttp://emath.ualberta.ca/LinrAlgbra/LinAlgInRn.xml/pages/Sec_ChangeCoordinates