Parte 2 de BSOR
Cargado por
Joel Abel Pin EspinalParte 2 de BSOR
Cargado por
Joel Abel Pin EspinalA medida que el radio espectral se vuelve más pequeño, la tasa de convergencia
aumenta, De hecho, −log10 ρG es una indicación del número de dígitos decimales por
los que La solución mejora con cada iteración convergente. Los radios espectrales de
la iteración de Jacobi, Gauss Seidel y SOR las matrices están relacionadas entre sí.
La relación entre los radios espectrales de la matriz de iteración SOR y la iteración de
Jacobi matriz es
1 1 1 /2
ρ1SOR = ω ρJ + [ ω2 ρ2J −4 ( ω−1 ) ]
/2
7.133
2 2
Cuando w = 1, PSOR se reduce a la iteración de Gauss-Seidel y Eq.7.133 se reduce a
la relación entre los radios espectrales del Matrices de iteración de Gauss-Seidel y
Jacobi.
ρ1GS/2= ρJ 7.134
Debido a que el logaritmo del radio espectral da una indicación de la tasa de
convergencia, se puede deducir de la Ec. 7.134 de la iteración Gauss-Seidel converge
dos veces más rápido que la iteración de Jacobi. iteración simplemente porque
log 10 ρG=2 log 10 ρJ 7.135
Un examen detallado de la ecuación. 7.133 revela que ningún valor de ω hace
ρ SOR<1 if | ρJ|>1. Por lo tanto, la iteración SOR diverge si la iteración de Jacobi diverge.
Puede ser muestra que se obtiene un valor mínimo de ρ SOR cuando
ω 2 ρ2J −4 ( ω−1 )=0 7.136
Evalúa ω opt es sustituyendo ρ J
2
ω opt = 7.138
( 1+ √ 1−ρGS )
La matriz se puede determinar a partir de un único cálculo de prueba con ø = 1
mediante el uso de la disminución de la diferencia entre aproximaciones sucesivas. La
relación de las diferencias máximas (mejoras) de iteraciones sucesivas satisface
|d (k+1)|
(k )
≈ ρGS=ρ2J 7.139
|d |
para valores suficientemente grandes de k. Una vez que el cociente de la Ec. 7.139
varía a un valor constante, se utiliza como ρGS en la ecuación. 7.138 para calcular ω opt ,
y las iteraciones se pueden continuar con la aceleración.
Con fines académicos, es posible realizar una serie de ejercicios numéricos en los que
se resuelve un conjunto de ecuaciones lineales con valores de o entre 1 y 2. La figura
7.35 muestra la naturaleza cualitativa de la relación entre el número de iteraciones
necesarias para la convergencia y el valor de ω.
