Laboratorio 1105 3

M.Sc. Ing.
Freddy Zambrana Rodrguez

Fecha de publicacin: 13 de mayo de 2013
Fecha de Entrega: 31 de mayo de 2013
1.
MTODO DE PUNTO FIJO MULTIVARIABLE.
Al igual que en el mtodo de punto fijo de una ecuacin no lineal, el mtodo multivariable busca
transformar el sistema de ecuaciones no lineales de manera que la primera incgnita x1 se puede
despejar de la primera ecuacin, x2 de la segunda y as sucesivamente, llegndose a formar un
sistema de la forma:
x1 g ( x1 , x2 , x3 ,..., xn )
x2 g ( x1 , x2 , x3 ,..., xn )
xn g ( x1 , x2 , x3 ,..., xn )
Al igual que en el mtodo para resolver una ecuacin no lineal, puede predecirse el
comportamiento y las caractersticas de este mtodo. Como es de conocimiento, el
reordenamiento de las ecuaciones en el caso lineal afecta la convergencia, por lo que puede
esperarse que la convergencia del mtodo en estudio dependa de si se despeja x1 de f1, de f2 o de
otra ecuacin respectiva.
Para asegurar la convergencia del mtodo, tambin existe un criterio de convergencia equivalente
al existente para resolver una ecuacin no lineal por el mtodo de punto fijo, que puede aplicarse
antes de iniciar el proceso iterativo.
De cualquier forma, sea cual sea el sistema al que se haya llegado y se vaya a resolver con este
mtodo, puede aumentarse la velocidad de convergencia usando desplazamientos sucesivos en
lugar de los desplazamientos simultneos, es decir, se procedera al proceso iterativo mediante
x1k 1 g ( x1k , x2k ,..., xnk )
x2k 1 g ( x1k 1 , x2k ,..., xnk )
k 1
n
1.1
k 1
1
g(x
, x2k 1 ,..., xnk )
METODOLOGA DE TRABAJO.
La secuencia de pasos a realizar en la utilizacin del presente mtodo son:

1. Abra MATLAB haciendo clic sobre el icono
2. Escriba un archivo M de funcin cuyo formato se presenta a continuacin. Por ejemplo, sea
Criterio el nombre del archivo M que para un sistema formado por dos ecuaciones con dos
incgnitas debe presentar el siguiente formato:
function sol = Criterio(x,n)
if n == 1
g1 g1 g 2 g 2
,
;
,
;
x1 x 2 x1 x 2
sol sum( abs ( A));

else
sol g1 ; g 2 ;
end;
NOTA. Las ecuaciones que representan a cada derivada parcial, como a cada funcin,
deben ser escritas de la siguiente manera:
x1 x 22 10
se escribe como
x (1). * x( 2).^ 2 10
3. Una vez desarrollado el archivo M correspondiente, retorne a la pantalla principal de

MATLAB y proceda a verificar que el archivo no presente errores de sintaxis, ni lgicos.
4. En el indicador de MATLAB >>, llamar al programa PuntoFijoMultiple_L cuyos formatos
de utilizacin son:
PuntoFijoMultiple_L(MaxIte);
Para cualquiera de los formatos, es necesario sealar la informacin correspondiente,
donde:
MaxIte variable que seala el nmero mximo de iteraciones a realizar por el programa.
5. De ser necesario, introduzca la informacin solicitada por el respectivo programa con ayuda
del teclado.
6. Anote los resultados obtenidos.
1.2
PROBLEMA A RESOLVER.
Utilice el mtodo iterativo de Punto Fijo Multivariable para resolver el siguiente sistema de
ecuaciones no lineales:
x 2 3 y 2 15 0
x3 2 y 2 0
Determine la solucin absoluta con una aproximacin de 10-2
2
MTODO DE NEWTON RAPHSON MULTIVARIABLE.
Este mtodo es considerado, basndose en la solucin de ecuaciones no lineales, como un

mtodo de segundo orden. Consiste fundamentalmente en escribir el sistema de ecuaciones no
lineales problema, como un sistema de ecuaciones lineales que presenta la forma:
f1
f
f
h1 1 h2 ... 1 hn f1
x1
x2
xn
f 2
f
f
h1 2 h2 ... 2 hn f 2
x1
x2
xn
f n
f
f
h1 n h2 ... n hn f n
x1
x2
xn
Que en forma concisa se expresa como J h f . Donde J es la matriz Jacobiana, h es el nuevo
vector incgnita y -f el vector de las funciones con signo cambiado y evaluadas en X(k). La solucin
del sistema de ecuaciones lineales se realiza utilizando uno de los mtodos vistos en el laboratorio
anterior (sistemas de ecuaciones lineales). Con la solucin del sistema y la ecuacin:
xik 1 hi xik
se obtiene el siguiente vector de aproximacin X(k+1). Este procedimiento se repite hasta satisfacer
algn criterio de convergencia establecido.
2.1
La secuencia de pasos a seguir en la utilizacin del presente mtodo son:

Jacobiano el nombre del archivo M para un sistema formado por dos ecuaciones con dos
function [matriz, func] = Jacobiano(x)
f 1 f 1 f 2 f 2
,
;
,
;
x1 x 2 x1 x 2
matriz
func [ f 1 , f 2 ];
NOTA . Las ecuaciones que representan a cada derivada parcial, como a cada funcin,
deben ser escritas de la siguiente manera:
2 x1 x 2 10
se escribe como
2 * x (1). * x( 2) 10

4. En el indicador de MATLAB >>, llamar al programa NewtonMultiple_L, cuyos formatos de
utilizacin son:
NewtonMultiple_L(MaxIte);
Para cualquiera de los formatos, es necesario enviar la informacin correspondiente,
donde:
MaxIte nmero mximo de iteraciones a realizar antes de abandonar el proceso.
del teclado.
2.2
Aplique el mtodo de Newton Raphson Multivariable para determinar una solucin absoluta del
sistema:
xe xy1 e y 2.5 0
2
x 2 y 2 0.5e xy 0
Con una aproximacin de 10-3
3
MTODO DE NEWTON RAPHSON MODIFICADO.
Este mtodo consiste en aplicar el mtodo de Newton Raphson modificado univariable

acondicionndolo para el caso de un sistema de n ecuaciones no lineales con n incgnitas, es
decir, aplicar la ecuacin una para cada variable de acuerdo a la siguiente expresin general:
x k 1 x k
f1 ( x k , y k , z k ,...)
f1
x xk , y k , z k ,...
Este mtodo converge a menudo si x0, y0, z0, estn muy cerca de x , y , z ,... y requiere la
evaluacin de solo 2n funciones por paso. Este mtodo puede ser diseado de dos maneras,
utilizando desplazamientos sucesivos y/o desplazamientos simultneos.
Este mtodo en muchos casos, no siempre puede producir convergencia ya que en algunos de los
arreglos a utilizar se producir divergencia. Es en este sentido que cuando n>3 las posibilidades
de despeje son varias n!, motivo por el cual, es imposible conocer cual de los diferentes arreglos
tiene viabilidad de convergencia y la eleccin de despeje se convierte en un proceso aleatorio.
Esta aleatoriedad es la mayor desventaja de este mtodo.
3.1
La secuencia de pasos a realizar en la utilizacin del mtodo son:

Modificado el nombre del archivo M para un sistema formado por dos ecuaciones con dos
function [matriz, func] = Modificado(x)
f 1 f 2
matriz
;
;
x y
func [ f 1 , f 2 ];
NOTA. Las ecuaciones que representan a cada derivada parcial, como a cada funcin,
deben ser escritas de acuerdo al siguiente formato:
12 x12 x 2 10 x1
se escribe como
12 x(1) 2 x ( 2) 10 x (1)
3. Una vez desarrollado el archivo M de funcin correspondiente, retorne a la pantalla

principal de MATLAB y proceda a verificar que el archivo no presente errores de sintaxis, ni
lgicos.
4. En el indicador de MATLAB >>, llamar al programa NewtonModificado_L, cuyos formatos

de utilizacin son:
NewtonModificado_L(MaxIte);
para cualquiera de los formatos, es necesario enviar la informacin correspondiente,
donde:
MaxIte l nmero mximo de iteraciones a realizar antes de abandonar el proceso.
5. De ser necesario introduzca la informacin solicitada por el respectivo programa con ayuda
del teclado.
3.2
Encuentre una solucin absoluta ms prxima al origen del siguiente sistema de ecuaciones no
lineales:
x 2 x y 2 1
y sin( x 2 ) 0
Emplee el mtodo de Newton Raphson modificado con un error de 10-2
4
MTODO DE NEWTON RAPHSON CON OPTIMIZACION DE T.
La utilizacin de la ecuacin vectorial de aceleracin permite estudiar como mejorar los mtodos
existentes. Es en este sentido, que la utilizacin de esta ecuacin vectorial en el mtodo de
Newton Raphson Multivariable permite generar el mtodo de Newton Raphson con optimizacin de
t. Mtodo que se obtiene a travs de la transformacin de la ecuacin matricial del mtodo de
Newton Raphson multivariable a la forma vectorial, cuya forma para un sistema de 2x2 es:
x ( k 1)
x (k )
y ( k 1) y ( k )
f 1
x
1
f 2
x
f 1
y
f 2
y
f1 ( x (k ) , y ( k ) )
(k )
(k )
f 2 (x , y )
Este mtodo utiliza en forma explicita la informacin de las derivadas para generar algoritmos
eficientes que localicen el ptimo, ya que distintos valores de t llevara a distintos vectores de xk+1,
algunos ms cercanos a la raz x que los dems. Para seleccionar los valores de t es necesario
seleccionar en principio un intervalo de bsqueda [a, b], para lo cual se calculan valores de t dentro
de este intervalo utilizando las ecuaciones
t a
ba
F
t b
ba
F
donde F son los trminos de la serie de Fibonacci. El proceso iterativo se realiza hasta que el valor
de funcin suma de residuos al cuadrado
2
S f 1 ( x, y ) f 2 ( x, y )
sea cero o prximo a cero, ya que esta funcin es el indicativo de la cercana de x k con respecto a
la raz
4.1
La secuencia de pasos a realizar son:

Optimizado el nombre del archivo M para un sistema formado por dos ecuaciones con dos
incgnitas debe presentar la forma:
function [matriz, func] = Optimizado(x)
f 1 f 1 f 2 f 2
matriz
,
;
,
;
x y x y
func f 1 , f 2 ;
NOTA . Las ecuaciones que representan a cada derivada parcial, como a cada funcin,
deben ser escritas de acuerdo al siguiente formato:
12 x12 x 2 10 x1
se escribe como
12 x(1) 2 x ( 2) 10 x (1)

4. En el indicador de MATLAB >>, llamar al programa NewtonOpT_L cuyos formatos de
utilizacin son:
NewtonOpT_L (MaxIte);
Para cualquiera de los formatos, es necesario enviar la informacin correspondiente,
donde:
MaxIte l nmero mximo de iteraciones a realizar antes de abandonar el proceso.
del teclado.
4.2
Resuelva el siguiente sistema de ecuaciones no lineales para obtener la solucin cerca de (2.5,
0.2, 1.6). Con una aproximacin de 10 -3 por el mtodo de Newton con optimizacin de t utilizando
como intervalo de bsqueda los valores de a = -0.3 y b = -2.5
x 2 y 2 z 2 9
xyz 1
x y z 2 0
PROBLEMAS QUE NO PRESENTAN SOLUCIN.
1
NO PRESENTA SOLUCIN.
PROBLEMAS QUE NECESITAN DE UN ERROR Y VECTOR INICIAL
1
ERROR = 0.001
VECTOR INICIAL =[1,2,3]
SOLUCIN = 1.2345
= -0.9987
= 9.8765
El trabajo a presentar debe encontrarse en formato texto (Block de Notas de Windows). El trabajo en
cuestin, se debe encontrar almacenado en el disco de trabajo como
GUIA_3.TXT

Laboratorio 1105 3

Cargado por

Información del documento

Título original

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

Laboratorio 1105 3

Cargado por

Copyright:

Formatos disponibles

M.Sc. Ing.

Freddy Zambrana Rodrguez

MTODO DE PUNTO FIJO MULTIVARIABLE.

, x2k 1 ,..., xnk )

La secuencia de pasos a realizar en la utilizacin del presente mtodo son:

sol sum( abs ( A));

3. Una vez desarrollado el archivo M correspondiente, retorne a la pantalla principal de

MTODO DE NEWTON RAPHSON MULTIVARIABLE.

Este mtodo es considerado, basndose en la solucin de ecuaciones no lineales, como un

La secuencia de pasos a seguir en la utilizacin del presente mtodo son:

3. Una vez desarrollado el archivo M correspondiente, retorne a la pantalla principal de

MTODO DE NEWTON RAPHSON MODIFICADO.

Este mtodo consiste en aplicar el mtodo de Newton Raphson modificado univariable

La secuencia de pasos a realizar en la utilizacin del mtodo son:

3. Una vez desarrollado el archivo M de funcin correspondiente, retorne a la pantalla

4. En el indicador de MATLAB >>, llamar al programa NewtonModificado_L, cuyos formatos

MTODO DE NEWTON RAPHSON CON OPTIMIZACION DE T.

La secuencia de pasos a realizar son:

3. Una vez desarrollado el archivo M correspondiente, retorne a la pantalla principal de

PROBLEMAS QUE NO PRESENTAN SOLUCIN.

También podría gustarte