Choosing the best estimate from a set of feasible solutions.
October 8, 2021 · 3 min read
The Problem
Consider the problem
Ax=b
where A∈Cm×n is full rank and m<n.
This is an under-determined system with an infinite number of solutions.
Suppose x0 is a solution to the problem such that Ax0=b.
Since A has more columns than rows, the N(A) is non-trivial
Therefore we now have an infinite set of solutions. For every v∈N(A) we have another solution A(x0+v)=b.
These solutions form an affine set of solutions.
To pin down a specific x we put a norm on x that we want to minimize subject to satisfying this system of equations.
This new formulation is written as
xargmin∥x∥s.t.Ax=b
Solving
We know that all solutions to the problem must satisfy the constraint so let’s investigate A and its null space a bit further.
Let the rows of A be written as:
A=y1H⋮ymH
Then using knowledge that v∈N(A) we have
Av=y1H⋮ymHv=y1Hv⋮ymHv=0⋮0
Since yiHv=0, each vector in the null space is orthogonal to yiH.
Therefore N(A)=M⊥ where M=span(y1…,ym).
Thus we now now know that the candidate solutions come from x0+M⊥.
We don’t have enough tools to prove this next part beyond proof by picture but the “shortest” solution is orthogonal to M⊥. In other words x∈(M⊥⊥=M=span(y1…,ym)). Therefore the best x can be expressed as
x=i=1∑mciyi
Now returning to the constraint Ax=b and using our new found knowledge of x we write