Exact solution of linear least square problem

  • Thread starter Thread starter Demon117
  • Start date Start date
  • Tags Tags
    Linear Square
Join the discussion
Registration is free. Ask a follow-up in this thread, or start your own.
1 reply · 2K views
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]