La matriz H es la matriz hessiana de segundas derivadas parciales. Encontramos a partir de la "diagonalización de la congruencia" que los valores propios de H son ++- \; \; . Por supuesto, las entradas diagonales de D no son los propios valores propios. Se trata de Ley de inercia de Sylvester .
La ley de inercia de Sylvester es un teorema del álgebra matricial sobre ciertas propiedades de la matriz de coeficientes de una forma cuadrática real que permanecen invariantes bajo un cambio de base. En concreto, si A es la matriz simétrica que define la forma cuadrática, y S es cualquier matriz invertible tal que D = SAS^T es diagonal, entonces el número de elementos negativos en la diagonal de D es siempre la misma, ya que todos esos S \; ; \; y lo mismo ocurre con el número de elementos positivos elementos positivos.
Si tuviera la energía para encontrar los vectores propios, podría escribir O^T H O = D_0 con la ortogonal O y la diagonal D_0, en cuyo momento las entradas diagonales de D_0 serían los valores propios. Sylvester dice que mi P^T H P = D da el mismo número de entradas diagonales positivas y el mismo número de negativas, también el mismo número de entradas nulas. Para este ejemplo, dos valores propios positivos y uno negativo.
P^T H P = D \left( \begin{array}{rrr} 1 & 0 & 0 \\ - \frac{ 4 }{ 3 } & 1 & 0 \\ - \frac{ 27 }{ 10 } & \frac{ 7 }{ 5 } & 1 \\ \end{array} \right) \left( \begin{array}{rrr} 6 & 8 & 5 \\ 8 & 14 & 2 \\ 5 & 2 & 4 \\ \end{array} \right) \left( \begin{array}{rrr} 1 & - \frac{ 4 }{ 3 } & - \frac{ 27 }{ 10 } \\ 0 & 1 & \frac{ 7 }{ 5 } \\ 0 & 0 & 1 \\ \end{array} \right) = \left( \begin{array}{rrr} 6 & 0 & 0 \\ 0 & \frac{ 10 }{ 3 } & 0 \\ 0 & 0 & - \frac{ 67 }{ 10 } \\ \end{array} \right)
E_j^T D_{j-1} E_j = D_j P_{j-1} E_j = P_j E_j^{-1} Q_{j-1} = Q_j P_j Q_j = I P_j^T H P_j = D_j Q_j^T D_j Q_j = H
H = \left( \begin{array}{rrr} 6 & 8 & 5 \\ 8 & 14 & 2 \\ 5 & 2 & 4 \\ \end{array} \right)
\==============================================
E_{1} = \left( \begin{array}{rrr} 1 & - \frac{ 4 }{ 3 } & 0 \\ 0 & 1 & 0 \\ 0 & 0 & 1 \\ \end{array} \right) P_{1} = \left( \begin{array}{rrr} 1 & - \frac{ 4 }{ 3 } & 0 \\ 0 & 1 & 0 \\ 0 & 0 & 1 \\ \end{array} \right) , \; \; \; Q_{1} = \left( \begin{array}{rrr} 1 & \frac{ 4 }{ 3 } & 0 \\ 0 & 1 & 0 \\ 0 & 0 & 1 \\ \end{array} \right) , \; \; \; D_{1} = \left( \begin{array}{rrr} 6 & 0 & 5 \\ 0 & \frac{ 10 }{ 3 } & - \frac{ 14 }{ 3 } \\ 5 & - \frac{ 14 }{ 3 } & 4 \\ \end{array} \right)
\==============================================
E_{2} = \left( \begin{array}{rrr} 1 & 0 & - \frac{ 5 }{ 6 } \\ 0 & 1 & 0 \\ 0 & 0 & 1 \\ \end{array} \right) P_{2} = \left( \begin{array}{rrr} 1 & - \frac{ 4 }{ 3 } & - \frac{ 5 }{ 6 } \\ 0 & 1 & 0 \\ 0 & 0 & 1 \\ \end{array} \right) , \; \; \; Q_{2} = \left( \begin{array}{rrr} 1 & \frac{ 4 }{ 3 } & \frac{ 5 }{ 6 } \\ 0 & 1 & 0 \\ 0 & 0 & 1 \\ \end{array} \right) , \; \; \; D_{2} = \left( \begin{array}{rrr} 6 & 0 & 0 \\ 0 & \frac{ 10 }{ 3 } & - \frac{ 14 }{ 3 } \\ 0 & - \frac{ 14 }{ 3 } & - \frac{ 1 }{ 6 } \\ \end{array} \right)
\==============================================
E_{3} = \left( \begin{array}{rrr} 1 & 0 & 0 \\ 0 & 1 & \frac{ 7 }{ 5 } \\ 0 & 0 & 1 \\ \end{array} \right) P_{3} = \left( \begin{array}{rrr} 1 & - \frac{ 4 }{ 3 } & - \frac{ 27 }{ 10 } \\ 0 & 1 & \frac{ 7 }{ 5 } \\ 0 & 0 & 1 \\ \end{array} \right) , \; \; \; Q_{3} = \left( \begin{array}{rrr} 1 & \frac{ 4 }{ 3 } & \frac{ 5 }{ 6 } \\ 0 & 1 & - \frac{ 7 }{ 5 } \\ 0 & 0 & 1 \\ \end{array} \right) , \; \; \; D_{3} = \left( \begin{array}{rrr} 6 & 0 & 0 \\ 0 & \frac{ 10 }{ 3 } & 0 \\ 0 & 0 & - \frac{ 67 }{ 10 } \\ \end{array} \right)
\==============================================
P^T H P = D \left( \begin{array}{rrr} 1 & 0 & 0 \\ - \frac{ 4 }{ 3 } & 1 & 0 \\ - \frac{ 27 }{ 10 } & \frac{ 7 }{ 5 } & 1 \\ \end{array} \right) \left( \begin{array}{rrr} 6 & 8 & 5 \\ 8 & 14 & 2 \\ 5 & 2 & 4 \\ \end{array} \right) \left( \begin{array}{rrr} 1 & - \frac{ 4 }{ 3 } & - \frac{ 27 }{ 10 } \\ 0 & 1 & \frac{ 7 }{ 5 } \\ 0 & 0 & 1 \\ \end{array} \right) = \left( \begin{array}{rrr} 6 & 0 & 0 \\ 0 & \frac{ 10 }{ 3 } & 0 \\ 0 & 0 & - \frac{ 67 }{ 10 } \\ \end{array} \right) Q^T D Q = H \left( \begin{array}{rrr} 1 & 0 & 0 \\ \frac{ 4 }{ 3 } & 1 & 0 \\ \frac{ 5 }{ 6 } & - \frac{ 7 }{ 5 } & 1 \\ \end{array} \right) \left( \begin{array}{rrr} 6 & 0 & 0 \\ 0 & \frac{ 10 }{ 3 } & 0 \\ 0 & 0 & - \frac{ 67 }{ 10 } \\ \end{array} \right) \left( \begin{array}{rrr} 1 & \frac{ 4 }{ 3 } & \frac{ 5 }{ 6 } \\ 0 & 1 & - \frac{ 7 }{ 5 } \\ 0 & 0 & 1 \\ \end{array} \right) = \left( \begin{array}{rrr} 6 & 8 & 5 \\ 8 & 14 & 2 \\ 5 & 2 & 4 \\ \end{array} \right)
1 votos
Si se trata de una tarea, ¿no se han tratado ya estos métodos en clase? (Diagonalización o completar los cuadrados, por ejemplo).