Solución iterativa a una ecuación no lineal.


8

Me disculpo de antemano si esta pregunta es tonta. Necesito calcular la raíz de

uf(u)=0

Donde es un vector real yf ( u ) es una función con valor de vector real. Comencé con el método de Newton (que funcionó), pero luego me di cuenta de que un método mucho más simple sería una solución iterativauf(u)

ui+1=f(ui)

Esto es mucho más rápido y aparentemente tan preciso / estable como el método de Newton.

Ahora las preguntas:

  • ¿Es este el enfoque correcto o debería usar un método diferente?
  • ¿Hay algo que se pueda decir sobre su tasa de convergencia, estabilidad, acc, etc.?
  • ¿Es globalmente convergente?

Gracias de antemano a todos por la atención.


3
Esto se conoce como iteración de punto fijo. No estoy bien versado en el tema, pero al menos debería darle algunas palabras nuevas para lanzar a Google. Si no recuerdo mal, aparecen puntos fijos para una amplia variedad de funciones con una amplia variedad de puntos de partida.
Godric Seer

1
f(u)

44
f()f(u)=u(g)1g(u)ui+1=f(ui) g(u)=0

Respuestas:



5

El fractal de Feigenbaum es un buen ejemplo de cuán extraña puede ser la iteración del punto de fijación:

http://en.wikipedia.org/wiki/Feigenbaum_fractal

http://en.wikipedia.org/wiki/File:Logistic_Bifurcation_map_High_Resolution.png

El segundo enlace traza el comportamiento de la iteración de punto fijo aplicada al mapa logístico ya que uno de los parámetros varía. Para ciertos valores converge, aunque solo linealmente. Para otros valores converge a un ciclo de longitud variable. Para otra clase de valores, se comporta de manera completamente caótica.

En otras palabras, el comportamiento de la iteración de punto fijo depende completamente de la función en cuestión. Incluso las funciones que parecen similares pueden exhibir un comportamiento radialmente diferente.

Nota : Como señala Jed, la iteración de Newton puede ser igualmente extraña .


Para ser justos, muchos fractales populares son los conjuntos de Julia de la iteración de Newton en una ecuación simple.
Jed Brown

4

El teorema de punto fijo de Banach describe la situación estándar cuando una iteración de punto fijo es globalmente convergente. Especialmente la parte de unicidad del teorema indica que solo puede esperar convergencia local si la solución no es única.

fn=fff


2

Puede considerar útil esta referencia: una homotopía para resolver problemas de puntos fijos grandes, dispersos y estructurados. R. Saigal. Matemáticas de Investigación de Operaciones, vol. 8, núm. 4 (noviembre de 1983), págs. 557-578.


1

Este método es correcto y se llama "sustitución sucesiva". Por favor, mire la página 189 de esta referencia para más detalles.

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.