Newton Raphson Multivariable
El método iterativo para sistema de ecuaciones converge linealmente.
Como en el método de una incógnita, pero puede crearse un método de convergencia
cuadrática; es decir, el método de newton–raphson multivariable. A continuación se
obtendrá este procedimiento para dos variables; la extensión a tres o más variables es viable
generalizando resultados.
Supóngase que se esta resolviendo el siguiente sistema:
Donde ambas funciones son continuas y diferenciables, de modo que puedan expandirse en la
serie de Taylor,esto es:
Donde f(x,y) se ha expandido alrededor del punto (a,b) y todas las derivadas parciales están
evaluadas en (x^k,y^k)
Expandiendo f1 alrededor de (x^k,y^k).
Donde todas las derivadas parciales están evaluadas en (x^k,y^k). de la misma forma puede
expandirese f2.
Donde todas las derivadas parciales están evaluadas en (x^k,y^k).
Ahora supóngase que x^(k+1) y y^(k+1) están cerca de la raíz buscada (x,y),que los lados
izquierdos de las dos ultimas ecuaciones son casi cero; además asúmase que x^k,y^k están tan
próximos de x^(k+1) que pueden omitirse los rerminos a partir de los que se encuentran
agrupados en los paréntesis rectangulares. Con esto las ecuaciones se simplifican a:
Para simplificar aun más, se cambia la notación con:
Y asi queda la (k+1)-esima iteración en términos de la k-esima,como se ve a continuación:
La sustitución de la ecuación y el arreglo dan como resultado:
La cual es un sistema de ecuaciones lineales con incógnitas h y j.
Este sistema de ecuaciones resultantes tiene solución única, siempre que el determinante de la
matriz de los coeficientes o matriz jacoviana J no sea cero, es decir:
Precisando:
El método de newton rapsón multivariable, consiste fundamentalmente en formar y resolver
el sistema de ecuaciones lineales, con la solución se obtiene una siguiente aproximación y
asi,una serie de iteraciones hasta a satisfacer algún criterio de convergencia establecido,
cuando converge este método lo hace con orden 2 y requiere que el vector (x0,y0) estén muy
cerca de la raíz buscada (x,y).
Representación gráfica de un sistema de ecuaciones:
grafica de la superficie f1(x,y).
plano tangente a ala superdicie f1(x,y) en el punto (1,1,1).
Intercesión del plano tangente y el plano x-y.
Generalización:
Para un sistema de n ecuaciones no lineales con n incognitas y retomando la notación vectorial
y matricial, las ecuaciones quedan:
De donde:
Y la matriz de derivadas parciales
(matriz jacoviana),ampliada en el vector de funciones queda:
Ejemplo:
Use el método de newton-raphon multivariable para encontrar una solución aproximada del
sistema:
Solución:
Primero se forma la matriz de derivadas parciales:
Y aumentada en el vector de funciones resulta:
Al evaluar en (x^0,y^0)=(0,0).
Tenemos:
Al resolver la matriz por eliminación gaussiana,tenemos los valores de h y j:
h= 0.8 j=0.88
Al sustituir en la ecuación:
Calculo de la distancia entre x^0 y x^1:
Segunda iteración:
Evaluando la matriz en (x^1.y^1).
Al resolver la matriz por eliminación gaussiana, tenemos los nuevos valores de h y j:
h =0.19179 j=0.11171
de donde:
calculo de la distancia entre x^1 y x^2:
Con la continuación de este proceso iterativo se obtienen los resultados siguientes:
