next up previous


Least-Squares Method


Let us consider the solution of a system of linear equations
 equation2547
by the least-squares method in the case the condition of the Kronecker-Capelli theorem is not satisfied, i.e., the system has no solution in ordinary sense.

Example 4.1.1. Let the system be
displaymath8093
where tex2html_wrap_inline8109 and rank(A)=2. Let tex2html_wrap_inline8113 be the orthogonal projection of the vector tex2html_wrap_inline8115 onto the space tex2html_wrap_inline8117 Since the vector tex2html_wrap_inline8119 and rank(A)=2, the system tex2html_wrap_inline8123 has a unique solution. Taking into consideration that tex2html_wrap_inline8125 we get tex2html_wrap_inline8127 and tex2html_wrap_inline8129 or
 equation2590
The matrix ATA of the system (2) is regular since rank(A)=2. Therefore the system (2) is uniquely solvable on the given conditions and
 equation2596
By minimizing the square of the norm of discrepancy tex2html_wrap_inline8135
displaymath8094
( tex2html_wrap_inline8137), we obtain the same system (2), and hence the same solution tex2html_wrap_inline8139 determined by the formula (3), the least-square solution of the equation (1).

The line of reasoning given in example 4.1.1 can be realized also in a more general case.

Definition 4.1.1. If tex2html_wrap_inline6877 then system (2) is called the system of normal equations of system (1).

Proposition 4.1.1. If tex2html_wrap_inline6877 tex2html_wrap_inline8145 and suppose rank(A)=n, then the system of normal equations (2) of system (1) is uniquely solvable and the least-squares solution tex2html_wrap_inline8139 of the system (1) is given by (3).

Example 4.1.2.* Let us solve by the least-squares method the system of equations
displaymath8095

We form the system of normal equations tex2html_wrap_inline8153
displaymath8096
Thus,

displaymath8097

If tex2html_wrap_inline7037, tex2html_wrap_inline8145 and rank(A)<n, then the system of normal equations (2) has infinite number of solutions, which can be all expressed as
displaymath8098
where tex2html_wrap_inline8161 and tex2html_wrap_inline8163 From among the solutions tex2html_wrap_inline8139 we will find the one having the least norm, the so-called optimum solution tex2html_wrap_inline8167 From the orthogonality of the vectors tex2html_wrap_inline8169 and tex2html_wrap_inline8171 it follows that
displaymath8099
Since from tex2html_wrap_inline8171 tex2html_wrap_inline8175 it follows tex2html_wrap_inline8177 then
displaymath8100
and tex2html_wrap_inline8169 tex2html_wrap_inline8181 is the optimum solution tex2html_wrap_inline8183 of the equation tex2html_wrap_inline8185. Thus, tex2html_wrap_inline8187 tex2html_wrap_inline8189


next up previous