Suppose three measurements constrain two unknowns:
The linear system asks for , , and , which are incompatible. A least-squares solution makes a balanced compromise. The matrix has full column rank, so
The fitted output and residual are
Thus the two direct measurements are each overshot by one, and the measured sum is undershot by one. The residual has Euclidean norm and is orthogonal to the column space, the plane in output coordinates .
Figure 1.
The pseudoinverse projects inconsistent data onto the attainable output plane
. For , the closest attainable output is . The segment between them is the residual , perpendicular to the column space. The drawing is an orthonormal two-dimensional slice of output space; the attainable plane intersects this slice in the blue straight line. Details.
For any perturbation , the orthogonality of to gives
Since the kernel of is trivial, every nonzero increases the residual. The least-squares solution is unique in this example. The formula using is convenient for these small exact calculations; an actual numerical computation can use a QR decomposition or singular value decomposition to avoid the extra conditioning cost of the least-squares normal equations. For full column rank, the squared singular values of are the eigenvalues of , so its condition number in the Euclidean norm is . [ref]

Articles by others on the same topic (0)

There are currently no matching articles.