Orthogonal Complement

Resources

Main Idea

Within Linear Algebra and Functional Analysis, given some subset W of an Inner Product space V, the orthogonal complement is the set of all vectors W⊥ that are orthogonal to every vector in W, as per some inner product on the parent space V. More formally, the orthogonal complement to W is

W⊥={h∈V:⟨h,g⟩=0,∀g∈U}.

The orthogonal decomposition writes a vector as a sum of a vector in W and a vector in the space orthogonal to W, W⊥.

Properties

Linear Algebra Case

Let W⊂Rn, then

W⊥:={v∈Rn:⟨u,v⟩=0 ∀u∈W}.

We can show that if W has an orthogonal basis {a1,a2,…,ap}, then it is equivalent for elements of W⊥ to be orthogonal to only these basis vectors, as opposed to everything in W. This can be shown by linearity of the inner product and the decomposition of any vector into a linear combination of its basis vectors.

A matrix A defines a subspace through it's Range / Column Space: W=range(A)=col(A). More formally,

W={Ax:x∈Rn}.

For u∈Rn, we wish to find its orthogonal decomposition as u=uW+uW⊥. By using an orthogonal basis, we can find the uW term,

uW=⟨a1,u⟩⟨a1,a1⟩a1+⟨a2,u⟩⟨a2,a2⟩a2+⋯+⟨ap,u⟩⟨ap,ap⟩ap.

This is a sum of projections along the basis vectors of the subspace of W. Then, uW⊥=u−uW.

Next, we wish to find W⊥. First, begin with its definition, and plug in u=Ax, which must be true for some x, as that is how u∈W is originally defined:

W⊥={v∈Rn:⟨Ax,v⟩=0 ∀x∈Rn}={v∈Rn:xTATv=0 ∀x∈Rn}.

Because this second condition holds for all x, we can "divide it out", giving

W⊥={v∈Rn:ATv=0}.

Thus, we have shown that the orthogonal complement of the range of a matrix is the Nullspace / Kernel of the transpose of the matrix!

Note, dim W≠dim W⊥, and W,W⊥⊂Rn. There are further properties about (W⊥)⊥ and the (set) addition of these spaces.

If we instead view A as an operator, A:Rn→Rn, then there exists (multiple?) A~:Rn→Rn such that the range of A is orthogonal to the range of A~ (over the domain Rn). More generally, the operators do not even need to share the same domain, just the same codomain. Here, A~ is the matrix associated with the nullspace of AT (i.e. proper dimension and such that ATA~=0).

For ease of discussion, let the columns of A, ai, be orthonormal. Then, repeating the formula from above:

uW=⟨a1,u⟩a1+⟨a2,u⟩a2+⋯+⟨ap,u⟩ap,=a1a1Tu+a2a2Tu+⋯+apapTu,=(a1a1T+a2a2T+⋯+apapT)u,=AATu.

Then,

uW⊥=(I−AAT)u.

From this, we can see that the operators that map from a vector to the subspace and it's orthogonal complement are AAT and I−AAT, respectively.