Teorema sobre la velocidad de convergencia en un método iterativo

(1/2) > >>

pabloN:
Supongamos que se quiere resolver el sistema (se asume no singular) mediante cierto método iterativo dado por:



en donde la matriz es invertible.

Si suponemos que la sucesión es convergente, se tendría:





Restando a se tiene que . Si se llama al error en el paso -ésimo queda:

. Llamando se tiene que el método es convergente si y sólo si (radio espectral menor a uno).

Hasta aquí un breve contexto en el que está inmerso el ejercicio. Dice así:

Sean los valores propios de la matriz asociada a cierto método iterativo, y vectores propios asociados a los .

Pruebe que si se cumple (o sea, el radio espectral es un indicador de la velocidad de convergencia).

Para la demostración lo que hice fue lo siguiente:

Sabemos que . Ahora al error en el paso 0 lo puedo expresar como combinación lineal de los vectores propios: . Entonces . Entonces llegamos a que:



Si se tiene la convergencia en la norma-1, se tiene la convergencia en cualquier otra norma ya que todas las normas en son equivalentes. Así que pasemos a estudiar:



Llamemos a la entrada -ésima del vector . Es decir,

. La entrada -ésima en virtud de lo dicho anteriormente será

Entonces (en este paso lo que hice fue dividir numerador y denominador entre )

Ésto fue para mostrar que para pero para que eso sea cierto debe cumplirse (desigualdad estricta). Eso podría concluirlo si la letra del problema fuera pero es en todos los casos... :-\

Mi pregunta es: ¿está bien lo que hice? ¿hay alguna manera de arreglarlo, o la demostración va por otro lado?

Desde ya, muchas gracias.

PD. Pido disculpas a los moderadores, en especial a administrador por darle problemas nuevamente, ya que creo que este hilo lo más adecuado sería que estuviera en el subforo de métodos numéricos, si bien es cierto que está directamente relacionado con álgebra lineal.

Editado

HernanV:
Hay algo que yo no entiendo. En este paso , ¿no estas intentando dividir vectores?.

pabloN:
Cita de: HernanV en 04/02/2012, 11:17:08 pm

Hay algo que yo no entiendo. En este paso , ¿no estas intentando dividir vectores?.


Hola HernanV, antes que nada muchísimas gracias por contestar. Está perfecto lo que dices; como estaba antes, estaba mal expresado. Fijate ahora a ver que te parece. Igual mi idea sigue siendo la misma. El problema es que no puedo concluir. ¿Cómo hago?

pabloN:
Cita de: pabloN en 04/02/2012, 08:10:37 pm

Si se tiene la convergencia en la norma-1, se tiene la convergencia en cualquier otra norma ya que todas las normas en   son equivalentes.


Eso es verdad. Si converge para alguna norma, entonces converge para cualquier otra norma (por equivalencia de normas). ¿Pero convergen a lo mismo? Creo que no tiene por qué. En dicho caso, tendría que haberse especificado la norma. Capaz que era la norma-2...

HernanV:
Te comento lo que veo, aunque no estoy seguro de que "el problema" este allí.

Sabes que el método es convergente si y sólo si el mayor autovalor de Q es en módulo menor que 1. O sea, eso lo sabes, lo podes dar por sentado que es así (¿verdad?).

Respecto de la norma, yo consideraría la norma infinito.

Navegación

[0] Índice de Mensajes

[#] Página Siguiente