Convergencia no monotónica en problema de punto fijo


13

Antecedentes

Estoy resolviendo una variante de la ecuación de Ornstein-Zernike de la teoría líquida. En resumen, el problema puede representarse como la resolución del problema de punto fijo , donde es un operador integroalgebraico y es la función de solución (la función de correlación directa de OZ). Estoy resolviendo por iteración Picard, donde proporciono una solución de prueba inicial y genero nuevas soluciones de prueba mediante el esquema donde es un parámetro ajustable que controla la mezcla de yA Δ j + 1 d r | c j + 1 ( r ) - c j ( r ) | < ϵ . A λ A c = cUNC(r)=C(r)UNC(r)C0 0(r)

cj+1=α(Acj)+(1α)cj ,
αCACutilizado en la próxima solución de prueba. Para esta discusión, supongamos que el valor de no es importante. Repito hasta que la iteración converja dentro de la tolerancia deseada, : En mi variante del problema, depende de un parámetro , y mi pregunta es acerca de cómo la convergencia de depende de este parámetro.αϵ
Δj+1rerEl |Cj+1(r)-Cj(r)El |<ϵ .
UNλUNC=C

Para un amplio rango de valores para , el esquema de iteración anterior converge exponencialmente rápido. Sin embargo, a medida que disminuyo , eventualmente alcanzo un régimen en el que la convergencia no es monotónica, como se muestra a continuación. λλλaparición de convergencia no monótona

Preguntas clave

En soluciones iterativas a problemas de punto fijo, ¿la convergencia no monotónica tiene algún significado especial? ¿Indica que mi esquema iterativo está al borde de la inestabilidad? Lo más importante , ¿debería la convergencia no monótona hacerme sospechar que la solución "convergente" no es una buena solución para el problema de punto fijo?

Respuestas:


1

Supongamos que es la variable independiente desconocida en la solución de , entonces el método de punto fijo convergerá desde un punto siempre que el Jacobian , donde es una constante . En general, no es un único punto, sino el dominio atravesado por el esquema iterativo.x = f ( x ) x fXX=F(X)XFX(X)αα<1X

  1. Su solución es convergente, aunque no de forma monotónica. Verifique en su jacobiano varios valores de y la variable de solución para ver si pasa de satisfacer los criterios de convergencia a no satisfacerlos, lo que podría explicar lo que está viendo.λ

  2. Si su solución ha convergido dentro de una tolerancia relativa establecida adecuadamente que también representa números pequeños, entonces lo ha hecho.


¿Puedes aclarar tu segundo punto?
Endulum

La diferenciaentre dos iteraciones sucesivas podría compararse con donde es una tolerancia relativa. El |Xj+1-XjEl |El |XjEl |ϵϵ
NameRakes
Al usar nuestro sitio, usted reconoce que ha leído y comprende nuestra Política de Cookies y Política de Privacidad.
Licensed under cc by-sa 3.0 with attribution required.