next up previous
Next: Orthogonal Projection of Up: No Title Previous: Example 1

Minimizing Ax-b

Let A be an tex2html_wrap_inline199 matrix. Suppose that Ax=b is an inconsistent system, we are interested in finding an x such that Ax is as close as possible to b.

Lets first look at the exercise 1. There are several ways to make your line ``close'' to given points, depending how we define ``closeness". Usual way is to add the square of d1, d2 d3 tex2html_wrap_i

nline209 dn, then minimize the sum of squares. See figure bellow.

tex2html_wrap_inline213 This method is called `` least square Approximation ".

We may also think that Ax in a tex2html_wrap_inline213 vector and b is another tex2html_wrap_inline213 vector. We want to minimize Ax-b. Assume that our vector space is an Inner Product Space with the usual Euclidean inner product, minimizing (ax-b) translates to minimizing the distance || Ax -b|| between the two vectors Ax and b . Note that || ax-b|| is the length of the vector Ax-b.

One thing which helps understanding the procedure of minimizing the distance between Ax and b is the fact that Ax is a vector in column space of A. (Why?)

Since Ax=b is not consistent, b is not in the column space of A. So we are looking for a vector, Ax, in the column space of A which is closest to the vector b. It can be proved ( see your text book) that such a vector is the orthogonal projection of b onto the column space of A. Now if Ax is the orthogonal projection of b onto col(A) then Ax-b is orthogonal to clo(A). ( why?) That is Ax-b is in tex2html_wrap_inline269 . So

tex2html_wrap_inline271

or

tex2html_wrap_inline273.

This system of linear equations is consistent ( Why?) and it is called Normal System. Moreover if the column of A are linearly independent, then tex2html_wrap_inline277 is invertible and tex2html_wrap_inline279

is the unique solution of the system tex2html_wrap_inline273.

If columns of A are orthonormal then tex2html_wrap_inline285 and we can easily compute tex2html_wrap_inline287, and the close st vector p in the column space of A to b is tex2html_wrap_inline295.




next up previous
Next: Orthogonal Projection of Up: No Title Previous: Example 1