Register to reply 
Diagonalizing a 3x3 second derivative matrix 
Share this thread: 
#1
Mar105, 10:25 PM

Emeritus
Sci Advisor
PF Gold
P: 2,977

I've been working on this problem lately where I've been looking at the second derivatives of 2D and 3D density fields. Now, the second derivatives of the field can be represented in a matrix, which can be thought of as an Ndimensional ellipse with the principal axes aligned along some angle in which the second derivative matrix is diagonal.
Anyway, I've solved the full problem in 2D (angle and all), but in 3D, I've only been able to diagonalize the matrix. I haven't yet figured out how to determine the angle that represents. Does an easy analytic solution exist for this or will I have to resort to numerical guesstimations? 


#2
Mar205, 06:15 AM

Sci Advisor
HW Helper
PF Gold
P: 4,119

Did you determine the eigenvectors of this matrix?
Knowing those, you can use the dotproduct to determine the various direction cosines. 


#3
Mar305, 07:20 AM

Emeritus
Sci Advisor
PF Gold
P: 2,977




#4
Mar305, 11:10 AM

Sci Advisor
HW Helper
P: 9,453

Diagonalizing a 3x3 second derivative matrix
the answer he gave you was for a general situation. the eigenvalues are roots of the characteristic polynomial, which has a general formula, the same for all matrices.



#5
Mar305, 05:51 PM

Emeritus
Sci Advisor
PF Gold
P: 2,977




#6
Mar305, 10:07 PM

Sci Advisor
HW Helper
P: 9,453

i still do not understand how a general method does not give you a general result. do you mean you wanted someone to write down the characteristic polynomial for you?
if so, it is just the determinant of the matrix [A  XId], where A is the matrix of second partials, or do you also want me to write down the determinant? 


#7
Mar405, 01:39 AM

Emeritus
Sci Advisor
PF Gold
P: 2,977

However, to find the eigenvectors, I then need to solve a linear system of equations which is usually singular. Since I can't do this with a simple inverse, writing a routine to do it is tricky. I can do it with gaussjordan elimination, but that's not so simple to program, so I needed the fortran library to do that. I was hoping there was a simple general solution for the rotation angle in this problem (assuming the eigenvectors make up the rotation matrix), but I have yet to find one in 3D. 


#8
Mar505, 03:56 PM

Sci Advisor
HW Helper
P: 9,453

gosh how could I have been so careless. Please forgive me, i just lost my head.
now that i know you want eigenvectors, perhaps (in case your matrix is invertible) the formula involving determinants could help you. I guess you knolw there is a formula for the, solution of an invertible linear system called cramers rule. but i wouldn't want to waste your time, so just feel free to ignore this if it is not helpful. In case it m,ay be useful, here is a other short characterizatioon of an eigenvector for a symmetric matrix A, which however requires you to solve for thew points on the unit sphere where a derivative is zero: let v be a unit vector such that the function Ax.x is minimized over the unit sphere at x = v. Then v is an elgenvector for A. 


Register to reply 
Related Discussions  
Help! Diagonal matrix similar to upper triangular matrix?  Linear & Abstract Algebra  1  
Diagonalizing a 3x3 matrix  Precalculus Mathematics Homework  4  
Replacing total derivative with partial derivative in Griffiths' book  Advanced Physics Homework  3  
The stiffness matrix in matrix analysis of strutures  Classical Physics  0  
Matrix, making R2 to R3, finding standard matrix A! where did i mess up?  Calculus & Beyond Homework  1 