Calcular el polinomio característico de la matriz dispersa real

9

Dada una matriz dispersa genérica con m << n (corrección: m n 2 ) elementos distintos de cero (típicamente m O ( n ) ). A es genérico en el sentido de que no tiene propiedades específicas (por ejemplo, definición positiva), y no se supone ninguna estructura (por ejemplo, bandas).ARn×nmn2mO(n)A

¿Cuáles son algunos de los buenos métodos numéricos para calcular el polinomio característico o el polinomio mínimo de ?A


fuente
3
Parece que quieres calcular todos los valores propios. ¿Por qué quieres el polinomio y cómo quieres que se exprese? La base monomial está extremadamente mal condicionada, por lo que es probable que los coeficientes no puedan calcularse de manera estable en aritmética de precisión finita.
Jed Brown
@JedBrown más de una contemplación. En mi respuesta a esta pregunta , di un método algebraico para invertir una matriz, que es bien conocida en álgebra computacional (por ejemplo, matrices sobre anillos y campos conmutativos). Quiero saber si podría usarlo para matrices numéricas. Tenga en cuenta que, a los fines de esta pregunta, estoy interesado en métodos numéricos para encontrar el polinomio característico / mínimo en lugar de inverso.

Respuestas:

1

Si la complejidad de no es un obstáculo, es posible que desee ver el método Danilevskii. Es bastante conocido en la literatura rusa sobre álgebra lineal numérica, pero no hay mucha información en inglés. Puedes comenzar desde este enlace .O(norte3)

La idea es bastante directa: la matriz se reduce gradualmente a la forma normal de Frobenius mediante transformaciones de similitud "de eliminación gaussiana". Si no encuentra la información, puedo hacer que el algoritmo sea más elaborado.

faleichik
fuente
1

Puede usar un método numérico como Factorización QR o Método de potencia y sus realives (potencia inversa, etc.) para calcular los valores propios de su matriz genérica. Después de eso, puede calcular su polinomio característico por factorización como: (λ-λ1) (λ-λ2) ... (λ-λn) = 0 donde λi son los valores propios calculados. Aquí hay una breve presentación sobre el poder y los métodos QR:

QR-Power

chemeng
fuente
0

metroO(norte2)metroO(norte)norteO(metro)

Wolfgang Bangerth
fuente
metronorte2,metroO(norte)