Sistema de Ecuaciones No Lineales

Sistema de ecuaciones no
lineales
En esta página, vamos a obtener las raíces de un sistema de ecuaciones transcendentes
aplicando el procedimiento de Newton-Raphson. Finalmente, utilizaremos la función
MATLAB fsolve para este propósito.
Sea el sistema de ecuaciones
{x2−y2+2y=02x+y2−6=0{x2−y2+2y=02x+y2−6=0
Representamos gráficamente las funciones implícitas mediante ezplot
>> hold on
>> ezplot('x^2-y^2+2*y',[-8,8,-8,8])
>> ezplot('2*x+y^2-6',[-8,8,-8,8])
>> hold off
>> title('Funciones implícitas')
>> grid on
Observamos que el sitema de ecuaciones tiene cuatro raíces, vamos a obtener estas
raíces aplicando el procedimiento de Newton-Raphson
Procedimiento de Newton-Raphson
Para aplicar este procedimiento para un sistema de n ecuaciones con n incógnitas,
representamos la variable x por x1 y la variable y por x2. El sistema de dos ecuaciones se
escribe de una forma más general
{f1(x1,x2)=0f2(x1,x2)=0{f1(x1,x2)=0f2(x1,x2)=0
Supongamos que en la etapa k de proceso de cálculo partimos de un punto (x1, x2)

cualesquiera y nos movemos a otro muy próximo (x1+Δx1, x2+Δx2). Los valores de las
funciones son f1 y f2 en dicho punto son aproximadamente
f1(x1+Δx1,x2+Δx2)≈f1(x1,x2)+∂f1∂x1Δx1+∂f1∂x2Δx2f2(x1+Δx1,x2+Δx2)≈f2(x1,x2)+∂f2∂
x1Δx1+∂f2∂x2Δx2f1(x1+Δx1,x2+Δx2)≈f1(x1,x2)+∂f1∂x1Δx1+∂f1∂x2Δx2f2(x1+Δx1,x2+Δx2)≈f2(x
1,x2)+∂f2∂x1Δx1+∂f2∂x2Δx2
Si el punto (x1+Δx1, x2+Δx2) es una solución del sistema de ecuaciones, entonces
f1(x1,x2)+∂f1∂x1Δx1+∂f1∂x2Δx2=0f2(x1,x2)+∂f2∂x1Δx1+∂f2∂x2Δx2=0f1(x1,x2)+∂f1∂x1Δx1+∂f
1∂x2Δx2=0f2(x1,x2)+∂f2∂x1Δx1+∂f2∂x2Δx2=0
Escribimos el sistema de ecuaciones en forma matricial para despejar Δx1 y Δx2
(f1f2)+⎛⎝⎜∂f1∂x1∂f2∂x1∂f1∂x2∂f2∂x2⎞⎠⎟(Δx1Δx2)=0(f1f2)+(∂f1∂x1∂f1∂x2∂f2∂x1∂f2∂x2)(Δx
1Δx2)=0
Denominamos vector x al vector (x1,x2), el vector función F está formado por dos
elementos que son las funciones (f1,f2) y la matriz cuadrada de dimensión dos es el
Jacobiano J. Despejamos Δx1 y Δx2 del sistema de ecuaciones o el vector Δx.
F(x)+JΔx=0Δx=−J−1FF(x)+JΔx=0Δx=−J−1F
J-1 es la matriz inversa de J y Δx es el vector diferencia entre el vector que nos da las
coordenadas del nuevo punto xk+1, conocidas las del punto previo xk
xk+1=xk−J−1Fxk+1=xk−J−1F
Hemos obtenido una expresión similar al procedimiento de Newton-Raphson que

utilizamos para calcular una raíz de la ecuación f(x)=0
Para un sistema de n ecuaciones
⎧⎩⎨⎪⎪⎪⎪⎪⎪⎪⎪f1(x1,x2...xn)=0f2(x1,x2...xn)=0...fn(x1,x2...xn)=0{f1(x1,x2...xn)=0f
2(x1,x2...xn)=0...fn(x1,x2...xn)=0
El procedimiento se escribe
⎛⎝⎜⎜⎜⎜x′1x′2...x′n⎞⎠⎟⎟⎟⎟=⎛⎝⎜⎜⎜x1x2...xn⎞⎠⎟⎟⎟−⎛⎝⎜⎜⎜⎜⎜⎜⎜∂f1
∂x1∂f2∂x1...∂fn∂x1∂f1∂x2∂f2∂x2...∂fn∂x2............∂f1∂xn∂f2∂xn...∂fn∂xn⎞⎠⎟⎟⎟⎟⎟⎟⎟−1⎛⎝⎜⎜⎜
⎜f1(x1,x2...xn)f2(x1,x2...xn)...fn(x1,x2...xn)⎞⎠⎟⎟⎟⎟(x1'x2'...xn')=(x1x2...xn)−(∂f1∂x1∂f
1∂x2...∂f1∂xn∂f2∂x1∂f2∂x2...∂f2∂xn............∂fn∂x1∂fn∂x2...∂fn∂xn)−1(f1(x1,x2...xn)f2(x1,x2...xn)...
fn(x1,x2...xn))
o bien, X=X-J\F, utilizando el operador MATLAB, división por la izquierda \. Se
obtiene el nuevo punto, el vector X de la izquierda, a partir del punto X previo, a la
derecha de la igualdad
Para el sistema de dos ecuaciones que hemos planteado al principio de esta sección
{x21−x22+2x2=02x1+x22−6=0{x12−x22+2x2=02x1+x22−6=0
El Jacobiano es
J=(2x12−2x2+22x2)J=(2x1−2x2+222x2)
El código MATLAB para calcular las raíces es similar al empleado para calcular una
raíz de la ecuación f(x)=0
F=@(x) [x(1)^2-x(2)^2+2*x(2);2*x(1)+x(2)^2-6];
J=@(x) [2*x(1),-2*x(2)+2;2,2*x(2)];
x=[-4.6;-3.8]; %punto inicial
i=0;
while i<100
y=-J(x)\F(x);
if sqrt(norm(y)/norm(x))<0.001
disp('Solución')
disp(x)
break;
end
x=x+y;
i=i+1;
end
if i>=100
disp('Se ha soprepasado el número de iteracciones');
end
En la primera línea se define el vector columna F de las funciones. En la segunda, la
matriz cuadrada J que representa el Jacobiano. En la tercera, las coordenadas del punto
inicial. En la representación gráfica de las dos funciones (figura más arriba), utilizamos
el icono denominado Data cursor del menú de la ventana gráfica, para conocer el valor
aproximado de una de las raíces (0,5, 2.2), que tomaremos como vector inicial de
partida para el procedimento de cálculo.
La función norm calcula el módulo un vector. El procedimiento concluye cuando los

puntos xmy xm+1, están muy próximos, es decir, cuando el error relativo |xm+1-xm|/|xm| sea
menor que una cantidad prefijada
Si el procedimiento no convergiera, el programa entraría en un bucle indefinido.

Limitando el número de iteracciones se asegura que siempre va a terminar.
0.6252
2.1794
>> F(x)
ans = 1.0e-09 *
0.3991
0.0536
Cuando introducimos el valor de la raíz buscada x en la función F, vemos que se
obtiene un valor próximo a cero
Con Data cursor buscamos un valor próximo a la segunda raíz (-4.6,-3.8).
Modificamos en la tercera línea de código el vector inicial que tomamos para calcular la
raíz y obtenemos al correr el script modificado
-4.8642
-3.9659
>> F(x)
ans = 1.0e-09 *
0.2845
0.1122
El problema más importante para aplicar este procedimiento es la elección del vector
inicial, cuando el número de ecuaciones es tres o más.
La función fsolve
MATLAB dispone de una función denominada fsolve que realiza un cálculo similar
>> F=@(x) [x(1)^2-x(2)^2+2*x(2);2*x(1)+x(2)^2-6];

>> fsolve(F,[0.5,2.2])
ans = 0.6252 2.1794
>> fsolve(F,[-4.6,-3.8])
ans = -4.8642 -3.9659
Ejemplo
Calcular las raíces del sistema de dos ecuaciones:
f1(x,y)=2x2-xy-5x-1=0
f2(x,y)=x+3log10x-y2=0
En la figura vemos el punto de intersección y mediante Tools/Data Cursor sus

coordenadas aproximadas..
x=linspace(2.5,5.5,50);
y1=(2*x.^2-5*x-1)./x;
y2=sqrt(x+3*log10(x));
plot(x,y1,'b',x,y2,'r')
Vamos a utilizar la función fsolve de MATLAB para obtener el punto de intersección.
Definimos la función F(x), donde x es un vector, x(1) representa a la variable x y x(2)

representa la variable y. La primera ecuación, es el primer elemento del vector F y la
segunda ecuación, es el segundo elemento de dicho vector
F=@(x) [2*x(1)^2-x(1)*x(2)-5*x(1)-1, x(1)*3*log10(x(1))-x(2)^2];

x0 =[3 2]; %valor inicial
[x,fval] = fsolve(F,x0);
fprintf('La solución es x=%1.3f, y=%1.3f\n',x(1),x(2));
fprintf('Valores de la función = %g\n',fval)
A fsolve tenemos que pasarle la función que hemos definido y la aproximaxión inicial
(x0, y0) que hemos encontrado anteriormente de forma gráfica. Esta función nos
devuelve las coordenadas (x,y) del punto de intersección buscado y los valores de las
funciones f1(x,y) y f2(x,y) en dicho punto.
La solución es x=3.958, y=2.664

Valores de la función = -1.33692e-010
Valores de la función = -3.07305e-009
Ejercicio
Resolver el sistema de ecuaciones no lineales
sin(xy)+exp(−xz)−0.9=0zx2+y2−−−−−−√−6.7=0tan(yx)+cosz+3.2=0sin(xy)+exp(−xz
)−0.9=0zx2+y2−6.7=0tan(yx)+cosz+3.2=0
Tomando x0=1, y0=2 y z0=2 como aproximación inicial

Grado en Ingeniería de Energías Renovables
Angel Franco García, Copyright © 2016
Código En Matlab.
Inicialmente, al igual que en los programas de este tipo que hemos
hecho hasta este momento, nuestra función en Matlab, solo recibirá
como argumento de entrada, una función que depende de una sola
variable (definida con anticipación). Llamaremos a nuestra función
‘pfalse’
A continuación se realiza una representación gráfica de la ecuación

para proporcionar un intervalo adecuado, en esta sección se piden un
valor inferior y uno superior desde donde gráficar la ecuación.
posteriormente se procede como tal con el desarrollo y ejecución del

me´todo de posición falsa, que es prácticamente idéntico al de
bisección, solo que cambiando el paso 2, donde aquí aplicamos la
formula de posición falsa en la linea 18 del código. (clic en la imagen
para verla mejor).
En este caso, también es de resaltar la suma importancia del calculo
del error, ya que es este el que nos va a proporcionar nuevamente la
condición de finalización para nuestro programa, igualmente sin el
correcto calculo del error relativo este método no tendría buen
desempeño a nivel computacional (prácticamente sería inservible).
El modo de uso de esta función es la siguiente, se define en la ventana

de comandos la variable simbólica (en este caso X) y se introduce una
función (para este caso Y) después de ubicar correctamente el current
folder de Matlab, se llama la función ‘pfalse‘ y se le pasa como
argumento la funcio ‘Y’, posteriormente se ingresan los limites donde
queremos ver la gráfica y finalmente se introduce el intervalo donde esta
la raíz.
EJEMPLO: halle la raíz positiva de la siguiente función.
procedemos a ejecutar el siguiente código en Matlab.
inicialmente vemos la siguiente representación, que es lo que

buscamos ver cada vez que realizamos el análisis gráfico.
Con esta podemos elegir los limites del intervalo coo Xl=0.5 y Xu=1.5
Como sabemos por simple inspección para este casó, la raíz de la

ecuación es 1 y vemos que el programa no ha dado una respuesta
aproximada, esta precisión la podemos mejorar variando el umbral de
error, dentro del While.
Este es un método bastante poderoso, lo es para ecuaciones no

lineales, que son las que regularmente encontramos en ingeniería y
para las cuales utilizamos con gran frecuencia los métodos numéricos,
es por esto que es de suma importancia tener un buen dominio y
comprensión tanto de la parte teórica de su funcionamiento, ya que
todos los métodos tienen sus falencias o por lo menos en algunos
casos unos se comportan mejor que otros, así, teniendo un buen
conocimiento de cada uno, sabremos elegir el más apropiado para
nuestras necesidades.
El código de este programa (archivo.m) lo puedes descargar de este

enlace, cualquier pregunta que tengas acerca de este archivo, puedes
comentar esta entrada.
Espero que este artículo haya sido de ayuda, y si lo fue puedes ver
demás métodos y códigos en este blog, saludos.

Sistema de Ecuaciones No Lineales

Cargado por

Copyright:

Formatos disponibles

Sistema de Ecuaciones No Lineales

Cargado por

Información del documento

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Sistema de Ecuaciones No Lineales

Cargado por

Copyright:

Formatos disponibles

Sistema de ecuaciones no

Sea el sistema de ecuaciones

Representamos gráficamente las funciones implícitas mediante ezplot

Supongamos que en la etapa k de proceso de cálculo partimos de un punto (x1, x2)

Si el punto (x1+Δx1, x2+Δx2) es una solución del sistema de ecuaciones, entonces

Escribimos el sistema de ecuaciones en forma matricial para despejar Δx1 y Δx2

Hemos obtenido una expresión similar al procedimiento de Newton-Raphson que

Para un sistema de n ecuaciones

La función norm calcula el módulo un vector. El procedimiento concluye cuando los

Si el procedimiento no convergiera, el programa entraría en un bucle indefinido.

>> F=@(x) [x(1)^2-x(2)^2+2*x(2);2*x(1)+x(2)^2-6];

Calcular las raíces del sistema de dos ecuaciones:

En la figura vemos el punto de intersección y mediante Tools/Data Cursor sus

Definimos la función F(x), donde x es un vector, x(1) representa a la variable x y x(2)

F=@(x) [2*x(1)^2-x(1)*x(2)-5*x(1)-1, x(1)*3*log10(x(1))-x(2)^2];

La solución es x=3.958, y=2.664

Resolver el sistema de ecuaciones no lineales

Tomando x0=1, y0=2 y z0=2 como aproximación inicial

A continuación se realiza una representación gráfica de la ecuación

posteriormente se procede como tal con el desarrollo y ejecución del

El modo de uso de esta función es la siguiente, se define en la ventana

procedemos a ejecutar el siguiente código en Matlab.

inicialmente vemos la siguiente representación, que es lo que

Como sabemos por simple inspección para este casó, la raíz de la

Este es un método bastante poderoso, lo es para ecuaciones no

El código de este programa (archivo.m) lo puedes descargar de este

También podría gustarte

Pfad - The Proxy pFad of © 2024 Garber Painting. All rights reserved.

>> F=@(x) [x(1)^2-x(2)^2+2x(2);2x(1)+x(2)^2-6];

F=@(x) [2x(1)^2-x(1)x(2)-5x(1)-1, x(1)3*log10(x(1))-x(2)^2];