Reduction of quadratic form (principal axis)

xman
Messages
92
Reaction score
0
i keep getting nonzero off diagonal elements when i try to reduce to simple sum of squares, of the equation
[tex]2 x_{1}^{2}+2x_{2}^{2}+x_{3}^{2}+2x_{1}x_{3}+2x_{2}x_{3}[/tex]
what i have is
[tex]\left(\begin{array}{ccc} x_{1} & x_{2} & x_{3} \end{array}\right)<br /> \left(\begin{array}{ccc}<br /> 2 & 1 & 0 \cr<br /> 1 & 2 & 1 \cr<br /> 0 & 1 & 1 <br /> \end{array} \right) <br /> \left(\begin{array}{c} x_{1} \cr x_{2} \cr x_{3} \end{array} \right)[/tex]
so my thought was to calculate the eigenvalues of the coefficient matrix above, which yield complex solutions from the characteristic equation
[tex]1-6 \lambda+5 \lambda^{2}-\lambda^{3}=0[/tex]
From the complex eigenvalues I obtain complex eigenvectors, which i'll post if necessary, but are rather lengthy. From the eigenvectors I choose to use Gram-Schmidt orthogonalization to form an orthonormal basis set. From which I construct a matrix with the corresponding basis set, and use diagonalize the system I have the diagonalization matrix
[tex]D = \left(\mid n \rangle \langle m \mid \right)^{T} A \left( \mid n \rangle \langle m \mid \right)[/tex]
where the matrix
[tex]\left(\mid n \rangle \langle m \mid \right)[/tex]
is the orthonormal eigenvector matrix. When I'm done with all of this I'm not getting a diagonalized matrix. I was wondering if I am making a mistake in my approach, or if anyone else does get a diagonalized matrix equation.
 
on Phys.org
The equation 1-6x+5x^2-x^3=0 has all real solutions. You could have known this because your matrix is symmetric and symmetric matrices in the reals are orthogonally diagonalizable in the reals.
 
0rthodontist said:
The equation 1-6x+5x^2-x^3=0 has all real solutions. You could have known this because your matrix is symmetric and symmetric matrices in the reals are orthogonally diagonalizable in the reals.
That's interesting, I was bad and usinig Maxima to calculate the roots of the equation and getting imaginary components, interesting, when I plot and find the roots, you're correct the roots are all real. Thanks, for the tip, I didn't know that about symmetric matrices, either. Thanks again.
 
xman said:
That's interesting, I was bad and usinig Maxima to calculate the roots of the equation
:biggrin: I was using my TI-89.
 
Perhaps the -89 is superior to my cas, even mathematica, maple given imaginary components, though on order of [tex]10^{-16}[/tex] or so. I wonder why that is. Not to mention I'm glad I'm not the only one whose cheats on algebra parts of problems...well that may get you into trouble as I found out today. Hey thanks again.
 

Similar threads

  • · Replies 3 ·
Replies
3
Views
3K
  • · Replies 12 ·
Replies
12
Views
3K
Replies
2
Views
3K
  • · Replies 7 ·
Replies
7
Views
2K
Replies
0
Views
1K
  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 7 ·
Replies
7
Views
4K
Replies
6
Views
2K
Replies
2
Views
2K
  • · Replies 10 ·
Replies
10
Views
2K