LDL^T factorization of a tridiagonal matrix giving wrong result

  • Thread starter Thread starter stunner5000pt
  • Start date Start date
  • Tags Tags
    Factorization
Join the discussion
Registration is free. Ask a follow-up in this thread, or start your own.
3 replies · 7K views
stunner5000pt
Messages
1,447
Reaction score
5
Find the [itex]LDL^T[/itex] factorization of this matrix

[tex]\left(\begin{array}{ccc}{2&-1&0\\-1&2&-1\\0&-1&2\end{array}\right)[/tex]

now i can find the L matrix by gaussian elimination
that yields
[tex]L = \left(\begin{array}{ccc}{1&0&0\\\frac{-2}{3}&1&0\\0&\frac{-1}{2}&1\end{array}\right)[/tex]
[tex]D = \left(\begin{array}{ccc}{\frac{1}{4}&0&0\\0&\frac{1}{3}&0\\0&0&\frac{1}{2}\end{array}\right)[/tex]

i am pretty sure about the ansswer since i checked my working many times.
However this is not the answer at the back of the book! In fact i am not even close!
What am i doing wrong?? Can anyone please help me iwth this?
Thank you for your help!
 
Physics news on Phys.org
It's hard to tell what you did wrong when you did wrong when you don't tell us what you did! I did a quick "column reduction" to get L and didn't get any like you got.
 
well i got those answers by Gaussian Elimination
this is what i did

[tex]\left(\begin{array}{ccc}{2&-1&0\\-1&2&-1\\0&-1&2\end{array}\right)[/tex]

R3 + 2R2
[tex]\left(\begin{array}{ccc}{2&-1&0\\-2&3&0\\0&-1&2\end{array}\right)[/tex]

R2+3R1
[tex]\left(\begin{array}{ccc}{4&0&0\\-2&3&0\\0&-1&2\end{array}\right)[/tex]

and my textbook says that that the D matrix is formed by dividing the square terms of the lower matrix formed and multiply that by the elementary matrix yielding
[tex]D = \left(\begin{array}{ccc}{\frac{1}{4}&0&0\\0&\frac{ 1}{3}&0\\0&0&\frac{1}{2}\end{array}\right)[/tex]
 
can anyone tell me what i have done wrong? my answer is not even close to the tedxxt book's answer. However all my steps with the row reductions are correct, as you can see.

I was told that i was not supposed to use row reduction to get the lower matrix? SO what do i do then?