Exact solution of linear least square problem

  • Thread starter Thread starter Demon117
  • Start date Start date
  • Tags Tags
    Linear Square
Click For Summary
SUMMARY

The exact solution to the linear least squares problem presented involves the matrix equation defined by the system of equations with the matrix A given as A = \left(\begin{array}{ccc}1 & 1 & 1 \\\epsilon & 0 & 0 \\0 & \epsilon & 0 \\0 & 0 & \epsilon \\ \end{array}\right) and the result vector \left(\begin{array}{c}1\\ 0\\ 0\\ 0\end{array}\right). The Cholesky factorization of A^{T}A yields A^{T}A = \left(\begin{array}{ccc} 1+\epsilon^{2} & 1 & 1 \\1 & 1+\epsilon^{2} & 1 \\1 & 1 & 1+\epsilon^{2}\end{array}\right). The final solution derived from the attempts is \vec{x} = \left(\begin{array}{c} 1/(3+\epsilon^{2}) \\1/(3+\epsilon^{2}) \\1/(3+\epsilon^{2}) \end{array}\right), indicating that the solution converges as epsilon approaches zero.

PREREQUISITES
  • Understanding of linear algebra concepts, specifically least squares problems.
  • Familiarity with matrix operations, including transposition and multiplication.
  • Knowledge of Cholesky factorization and its application in solving linear systems.
  • Basic calculus to comprehend the behavior of functions as epsilon approaches zero.
NEXT STEPS
  • Study the derivation of the least squares solution using normal equations.
  • Learn about the implications of matrix conditioning and how it affects solutions in linear algebra.
  • Explore numerical methods for solving linear systems, focusing on stability and accuracy.
  • Investigate the role of regularization techniques in least squares problems, particularly in ill-posed cases.
USEFUL FOR

Students and professionals in mathematics, engineering, and data science who are dealing with linear regression, optimization problems, or numerical analysis will benefit from this discussion.

Demon117
Messages
162
Reaction score
1

Homework Statement


What is the exact solution to the linear least square problem

[tex]\left(\begin{array}{ccc}1 & 1 & 1 \\\epsilon & 0 & 0 \\0 & \epsilon & 0 \\0 & 0 & \epsilon \\ \end{array}\right)[/tex][tex]\left(\begin{array}{c}x_{1}\\x_{2}\\x_{3}\end{array}\right)[/tex]=[tex]\left(\begin{array}{c}1\\ 0\\ 0\\ 0\end{array}\right)[/tex]

as a function of epsilon.


The Attempt at a Solution


I've tried numerous methods, including Cholesky factorization of [tex]A^{T}A[/tex], which is of course

[tex]A^{T}A[/tex] = [tex]\left(\begin{array}{ccc} 1+\epsilon^{2} & 1 & 1 \\1 & 1+\epsilon^{2} & 1 \\1 & 1 & 1+\epsilon^{2}\end{array}\right)[/tex]

The Cholesky factorization is a difficult one but in the end I get a vector that cannot possibly be a solution to this problem, what I get is

[tex]\vec{x}[/tex] = [tex]\left(\begin{array}{c} 1+\epsilon^{2} \\1 \\1 \end{array}\right)[/tex]

Would anyone be willing to give me a few pointers?
 
Physics news on Phys.org
matumich26 said:

Homework Statement


What is the exact solution to the linear least square problem

[tex]\left(\begin{array}{ccc}1 & 1 & 1 \\\epsilon & 0 & 0 \\0 & \epsilon & 0 \\0 & 0 & \epsilon \\ \end{array}\right)[/tex][tex]\left(\begin{array}{c}x_{1}\\x_{2}\\x_{3}\end{array}\right)[/tex]=[tex]\left(\begin{array}{c}1\\ 0\\ 0\\ 0\end{array}\right)[/tex]

as a function of epsilon.


The Attempt at a Solution


I've tried numerous methods, including Cholesky factorization of [tex]A^{T}A[/tex], which is of course

[tex]A^{T}A[/tex] = [tex]\left(\begin{array}{ccc} 1+\epsilon^{2} & 1 & 1 \\1 & 1+\epsilon^{2} & 1 \\1 & 1 & 1+\epsilon^{2}\end{array}\right)[/tex]

The Cholesky factorization is a difficult one but in the end I get a vector that cannot possibly be a solution to this problem, what I get is

[tex]\vec{x}[/tex] = [tex]\left(\begin{array}{c} 1+\epsilon^{2} \\1 \\1 \end{array}\right)[/tex]

Would anyone be willing to give me a few pointers?

Well, I've tried and tried again . . . but I keep ending up with something like this:

[tex]\vec{x}[/tex] = [tex]\left(\begin{array}{c} 1/(3+\epsilon^{2}) \\1/(3+\epsilon^{2}) \\1/(3+\epsilon^{2}) \end{array}\right)[/tex]
 

Similar threads

  • · Replies 2 ·
Replies
2
Views
2K
Replies
6
Views
3K
Replies
6
Views
1K
  • · Replies 5 ·
Replies
5
Views
2K
Replies
4
Views
2K
  • · Replies 9 ·
Replies
9
Views
2K
  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 4 ·
Replies
4
Views
1K
  • · Replies 4 ·
Replies
4
Views
2K
Replies
1
Views
2K