linear least squares
Let be an matrix with and an matrix. We want to consider the problem
We want to find the vector which is closest to in the column space of .
Among the different methods to solve this problem, we mention Normal Equations, sometimes ill-conditioned, QR Decomposition, and, most generally, Singular Value Decomposition. For further reading, [Golub89], [Branham90], [Wong92], [Press95].
Example: Let us consider the problem of finding the closest point (vertex) to measurements on straight lines (e.g. trajectories emanating from a particle collision). This problem can be described by with
This is clearly an inconsistent system of linear equations, with more equations than unknowns, a frequently occurring problem in experimental data analysis. The system is, however, not very inconsistent and there is a point that lies “nearly” on all straight lines. The solution can be found with the linear least squares method, e.g. by QR decomposition for solving :
Originally from The Data Analysis Briefbook (http://rkb.home.cern.ch/rkb/titleA.htmlhttp://rkb.home.cern.ch/rkb/titleA.html)
S.S.M. Wong, Computational Methods in Physics and Engineering, Prentice Hall, 1992.
Gene H. Golub and Charles F. van Loan: Matrix Computations, 2nd edn., The John Hopkins University Press, 1989.
R.L. Branham, Scientific Data Analysis, An Introduction to Overdetermined Systems, Springer, Berlin, Heidelberg, 1990.
W.H. Press, S.A. Teukolsky, W.T. Vetterling, and B.P. Flannery, Numerical Recipes in C, Second edition, Cambridge University Press, 1995. (The same book exists for the Fortran language). There is also an Internet version which you can work from.
|Title||linear least squares|
|Date of creation||2013-03-22 12:05:50|
|Last modified on||2013-03-22 12:05:50|
|Last modified by||akrowne (2)|