Gram-Schmidt Process

Dependencies:

  1. Inner product space
  2. Orthogonality and orthonormality
  3. Span
  4. Linear independence

The Gram-Schmidt process takes a set of n linearly independent vectors as input and outputs a set of n orthogonal vectors which have the same span.

Let W={w1,w2,,wn} be a linearly independent set of vectors from an inner product space. Then the Gram-Schmidt process outputs V={v1,v2,,vn} where

vi=wij=1i1ci,jvjwhereci,j=wi,vjvj2

Then:

Also, scaling each vi by a (possibly different) non-zero scalar maintains orthogonality and does not change the span.

Proof by induction

Let Vk={v1,v2,,vk} and Wk={w1,w2,,wk}.

Let the predicate P(k) be the 'logical and' of the following statements:

  1. span(Vk)=span(Wk).
  2. ik,vi0.
  3. i<jk,vi,vj=vj,vi=0.

Base case:

v1=w1.

  1. span(V1)=span({v1})=span({w1})=span(W1).
  2. Since W is linearly independent, w10v10.
  3. The third part of P(1) is trivially true.

Therefore, P(1) holds.

Inductive step:

Assume P(k) is true (inductive hypothesis).

Part 1

ik,wispan(Wk)=span(Vk)span(Vk+1)

wk+1=vk+1+j=1kck+1,jvjspan(Vk+1)

Since, each of w1,w2,,wk+1 is a linear combination of Vk+1. Their linear combination is also a linear combination of Vk+1. Therefore, span(Wk+1)span(Vk+1).

ik,vispan(Vk)=span(Wk)span(Wk+1)

j=1kck+1,jvjspan(Vk)=span(Wk)span(Wk+1)vk+1=wk+1j=1kck+1,jvjspan(Wk+1)

Since, each of v1,v2,,vk+1 is a linear combination of Wk+1. Their linear combination is also a linear combination of Wk+1. Therefore, span(Vk+1)span(Wk+1).

Therefore, span(Vk+1)=span(Wk+1), which is part 1 of P(k+1).

Part 2

vk+1=0wk+1=i=1kck+1,ivispan(Vk)=span(Wk)

wk+1span(Wk) means Wk+1 is linearly dependent, which is a contradiction. Therefore, vk+10. Combining this fact with part 2 of P(k), we get ik+1,vi0, which is the part 2 of P(k+1).

Part 3

Let ik.

vk+1,vi=wk+1j=1kwk+1,vjvj2vj,vi(linearity in first argument)=wk+1,vij=1kwk+1,vjvj2vj,vi(ijvj,vi=0)=wk+1,viwk+1,vivi2vi,vi=wk+1,viwk+1,vi=0

(by conjugate symmetry)vi,vk+1=vk+1,vi=0=0

Combining the above fact with part 3 of P(k), we get part 3 of P(k+1).

Therefore, P(n) is true by mathematical induction.

Dependency for:

  1. Standard normal random vector on vector space
  2. Symmetric operator iff hermitian
  3. Symmetric operator on V has a basis of orthonormal eigenvectors
  4. Joining orthogonal linindep sets

Info:

Transitive dependencies:

  1. Group
  2. Ring
  3. Field
  4. Vector Space
  5. Linear independence
  6. Span
  7. Inner product space
  8. Orthogonality and orthonormality