Minimum Norm Problem

Choosing the best x\pmb x from a set of solutions.

The Problem

Consider the problem

Ax=b\begin{align*} A\pmb x = \pmb b \end{align*}

where A∈Cm×nA \in \mathbb C^{m \times n} is full rank and m<nm < n.

This is an under-determined system with an infinite number of solutions.

Suppose x0\pmb x_0 is a solution to the problem such that Ax0=bA \pmb x_0 = \pmb b.

Since AA has more columns than rows, the N(A)\mathcal N(A) is non-trivial

Therefore we now have an infinite set of solutions. For every v∈N(A)\pmb v \in \mathcal N(A) we have another solution A(x0+v)=bA(\pmb x_0 + \pmb v) = \pmb b.

These solutions form an affine set of solutions.

To pin down a specific x\pmb x we put a norm on x\pmb x that we want to minimize subject to satisfying this system of equations.

This new formulation is written as

arg min⁡x ∥x∥  s.t.  Ax=b\begin{align*} \argmin_{\pmb x}\ \lVert \pmb x \rVert\ \ s.t.\ \ A\pmb x = \pmb b \end{align*}

Solving

We know that all solutions to the problem must satisfy the constraint so let's investigate AA and its null space a bit further.

Let the rows of AA be written as:

A=[y1H⋮ymH]\begin{align*} A = \begin{bmatrix}\pmb y^H_1 \\ \vdots \\ \pmb y^H_m\end{bmatrix} \end{align*}

Then using knowledge that v∈N(A)\pmb v \in \mathcal N(A) we have

Av=[y1H⋮ymH]v=[y1Hv⋮ymHv]=[0⋮0]\begin{align*} A \pmb v = \begin{bmatrix}\pmb y^H_1 \\ \vdots \\ \pmb y^H_m\end{bmatrix} \pmb v &= \begin{bmatrix}\pmb y^H_1 \pmb v\\ \vdots \\ \pmb y^H_m \pmb v\end{bmatrix} = \begin{bmatrix}0\\ \vdots \\ 0\end{bmatrix} \end{align*}

Since yiHv=0\pmb y_i^H \pmb v = 0, each vector in the null space is orthogonal to yiH\pmb y_i^H.

Therefore N(A)=M⊥\mathcal N(A) = M^\perp where M=span(y1…,ym)M = \text{span}(\pmb y_1 \dots, \pmb y_m).

Thus we now now know that the candidate solutions come from x0+M⊥\pmb x_0 + M^\perp.

We don't have enough tools to prove this next part beyond proof by picture but the "shortest" solution is orthogonal to M⊥M^\perp. In other words x∈(M⊥⊥=M=span(y1…,ym))\pmb x \in (M^{\perp\perp}=M=\text{span}(\pmb y_1 \dots, \pmb y_m)). Therefore the best x\pmb x can be expressed as

x=∑i=1mciyi\begin{align*} \pmb x = \sum_{i=1}^m c_i \pmb y_i \end{align*}

Now returning to the constraint Ax=bA \pmb x = \pmb b and using our new found knowledge of x\pmb x we write

Ax=b[y1H⋮ymH]x=[b1⋮bm][y1Hx⋮ymHx]=[b1⋮bm][y1H∑i=1mciyi⋮ymH∑i=1mciyi]=[b1⋮bm][∑i=1mciy1Hyi⋮∑i=1mciymHyi]=[b1⋮bm][c1y1Hy1+⋯+cmy1Hym⋮c1ymHy1+⋯+cmymHym]=[b1⋮bm][y1Hy1…y1Hym⋮ymHy1…ymHym][c1⋮cm]=[b1⋮bm]AAHc=b  ⟹  c=(AAH)−1b\begin{align*} A\pmb x &= \pmb b \\[15pt] \begin{bmatrix}\pmb y^H_1 \\ \vdots \\ \pmb y^H_m\end{bmatrix} \pmb x &= \begin{bmatrix}b_1 \\ \vdots \\ b_m\end{bmatrix} \\[25pt] \begin{bmatrix}\pmb y^H_1 \pmb x \\ \vdots \\ \pmb y^H_m \pmb x\end{bmatrix} &= \begin{bmatrix}b_1 \\ \vdots \\ b_m\end{bmatrix} \\[25pt] \begin{bmatrix}\pmb y^H_1 \sum_{i=1}^m c_i \pmb y_i \\ \vdots \\ \pmb y^H_m \sum_{i=1}^m c_i \pmb y_i\end{bmatrix} &= \begin{bmatrix}b_1 \\ \vdots \\ b_m\end{bmatrix} \\[25pt] \begin{bmatrix} \sum_{i=1}^m c_i \pmb y^H_1 \pmb y_i \\ \vdots \\ \sum_{i=1}^m c_i\pmb y^H_m \pmb y_i\end{bmatrix} &= \begin{bmatrix}b_1 \\ \vdots \\ b_m\end{bmatrix} \\[25pt] \begin{bmatrix} c_1 \pmb y^H_1 \pmb y_1 + \dots + c_m \pmb y^H_1 \pmb y_m \\ \vdots \\ c_1 \pmb y^H_m \pmb y_1 + \dots + c_m \pmb y^H_m \pmb y_m\end{bmatrix} &= \begin{bmatrix}b_1 \\ \vdots \\ b_m\end{bmatrix} \\[25pt] \begin{bmatrix} \pmb y^H_1 \pmb y_1 & \dots & \pmb y^H_1 \pmb y_m \\ \vdots \\ \pmb y^H_m \pmb y_1 & \dots & \pmb y^H_m \pmb y_m\end{bmatrix} \begin{bmatrix} c_1 \\ \vdots \\ c_m\end{bmatrix} &= \begin{bmatrix}b_1 \\ \vdots \\ b_m\end{bmatrix} \\[25pt] AA^H \pmb c &= \pmb b \\ \implies \pmb c &= (AA^H)^{-1} \pmb b \end{align*}

Bringing it all home

x=∑i=1mciyi=AHc=AH(AAH)−1b\begin{align*} \pmb x &= \sum_{i=1}^m c_i \pmb y_i\\ &= A^H \pmb c \\ &= A^H(AA^H)^{-1} \pmb b \end{align*}