Estrategias para el método de Newton cuando el jacobiano en la solución es singular

12

Estoy tratando de resolver el siguiente sistema de ecuaciones para las variables y (todo lo demás son constantes): $P,x_1$ $x_2$

\frac{A (1 - P)}{2} - k_{1} x_{1} = 0 \frac{A P}{2} - k_{2} x_{2} = 0 \frac{(1 - P) (r_{1} + x_{1})^{4}}{L_{1}} - \frac{P (r_{1} + x_{2})^{4}}{L_{2}} = 0

$\frac{A(1-P)}{2}-k_1x_1=0 \\ \frac{AP}{2}-k_2x_2=0 \\ \frac{(1-P)(r_1+x_1)^4}{L_1}-\frac{P(r_1+x_2)^4}{L_2}=0$

Puedo ver que puedo convertir este sistema de ecuaciones en una sola ecuación de una sola variable resolviendo las ecuaciones 1 y 2 para y respectivamente y sustituyéndolas en la ecuación 3. Al hacerlo, puedo use el comando de matlab para encontrar la solución. Usando los parámetros , y , encontré que la verdadera solución es $(P)$ $x_1$ $x_2$ fzero $k_1=k_2=1$ $r_1=r_2=0.2$ $A=2$ . $P=x_1=x_2=0.5$

Sin embargo, cuando uso el método de newton aplicado al sistema de ecuaciones original de 3 variables - 3, las iteraciones nunca convergen a la solución, no importa cuán cerca comience a la solución verdadera . $x^*=(P^*,x_1^*,x_2^*)=(0.5,0.5,0.5)$

Al principio, sospeché que era un error en mi implementación del método de newton. Después de revisar varias veces, no encontré ningún error. Luego intenté usar una conjetura inicial , y he aquí: el jacobiano es singular. Sé que un jacobiano singular puede reducir el orden de convergencia, pero no creo que necesariamente evite la convergencia a la solución verdadera. $x_0=x^*$

Entonces, mi pregunta es, dado que el jacobiano del sistema en la verdadera solución es singular:

¿Qué otras condiciones son necesarias para demostrar que el método de newton no convergerá a la raíz?
¿Una estrategia de globalización (por ejemplo, búsqueda de línea) garantizaría la convergencia a pesar del singular jacobiano?

optimization convergence newton-method

— Paul
fuente

4

(1): Esto depende del comportamiento de los derivados del jacobiano (¡sic!) En el espacio nulo del jacobiano en la solución. En la práctica, nadie calcula estos derivados, y ni siquiera me molesté en recordar las condiciones precisas.

(2) funciona, aunque la convergencia es solo lineal.

Para obtener una convergencia superlineal (al menos en la mayoría de los casos), uno puede usar métodos tensoriales. Ver, por ejemplo,
https://cfwebprod.sandia.gov/cfdocs/CCIM/docs/SAND2004-1944.pdf
http://www.jstor.org/stable/10.2307/2156931
http://www.springerlink.com/ index / X5G827367G548327.pdf

— Arnold Neumaier
fuente

3

En lugar de una búsqueda de línea, puede intentar modificar el método de Newton al algoritmo Levenberg-Marquardt . La implementación es lo suficientemente simple si está usando Matlab.

— Leo Uieda
fuente