Modelos de Variable No Lineales (Varias Variable)

Descargar como pptx, pdf o txt
Descargar como pptx, pdf o txt
Está en la página 1de 17

Modelos de Variables No

Lineales

Varias Variables
Método de Iteración
Los métodos que se van a tratar para la solución de ecuaciones en
varias variables son: Método de Iteración y el Método de Newton.

Sea un sistema de ecuaciones no lineales de un tipo especial:


x1 = g1(x1, x2, . . . , xn),
x2 = g2(x1, x2, . . . , xn),
............
xn = gn(x1, x2, . . . , xn),
donde las funciones g1, g2, . . . , gn son reales, definidas y continuas
en un conjunto D. De forma más compacta el sistema se puede
escribir como:
x = G(x).
Método de Iteración

Definición. Se dice que una función G: D ⊂ Rn → Rn tiene un punto


fijo en p ∈ D si G(p) = p.

Para hallar la raíz vectorial p = (p1, p2, . . . , pn) de la ecuación,


frecuentemente resulta conveniente utilizar el método de iteración
x(k+1) = G(x(k)),
k = 0, 1, 2, . . ., donde la aproximación inicial x(0) ≈ p.

El método de iteración puede aplicarse también al sistema general


f1(x1, x2, . . . , xn) = 0,
f2(x1, x2, . . . , xn) = 0,
.........
fn(x1, x2, . . . , xn) = 0,
Método de Iteración

Usando notación vectorial para representar las variables x i, el


sistema asume la forma:
F(x) = 0.
donde F(x) es una función vectorial, definida y continua en la
vecindad D de una raíz separada p. Por ejemplo, escribamos
nuevamente este sistema de la forma
x = x + Λ F(x) = G(x),
donde Λ es una matriz no singular. El método ordinario de iteración
es directamente aplicable a esta ecuación. Si la función F(x) tiene
una derivada continua F’(x) en D, se sigue que

G’(x) = I + Λ F’(x) = I + ΛJ(x),


Método de Iteración

donde la matriz J(x) = F’(x) es la matriz Jacobiana, definida como


J(x) = Jij = (). Entonces elegiremos la matriz Λ de forma que

G’(x(0)) = I + Λ J(x(0)) = 0,

de donde, si la matriz J(x(0)) = F’(x(0)) es no singular, tendremos

Λ = −[J(x(0))]−1.

Si det[J(x(0))] = 0, debe elegirse una aproximación inicial x(0) diferente.


Es decir, es esencialmente el proceso de Newton modificado.
Método de Iteración
Ejemplo:

Utilícese el método de iteración para dar una solución aproximada del


sistema
x12 + x22 = 1,
x13 − x2 = 0.
A partir de una construcción gráfica, puede verse que el sistema
dado tiene dos soluciones que difieren únicamente en el signo. Nos
limitaremos a la solución positiva. Se puede tomar x (0) = (0.9, 0.5)t
para la aproximación inicial de la solución positiva. Estableciendo
Método de Iteración
Ejemplo:

tenemos

de donde

Como la matriz J(x(0)) es no singular, existe una matriz inversa, y de


este modo
Método de Iteración
Ejemplo:

Ahora hacemos

Usando esta última fórmula hallamos las aproximaciones sucesivas


siguientes para la solución del sistema dado (con aritmética de cuatro
dígitos significativos):
Método de Iteración
Ejemplo:
Método de Iteración
Ejemplo:

y así sucesivamente. Finalizando con la cuarta aproximación,


tenemos las raíces

x1 = 0.8261, x2 = 0.5636

y F(x) = (0.0000, 0.0001)t , con una exactitud de 10−4 .


Método de Newton Raphson
La forma general de un sistema de ecuaciones no lineales es:
f1(x1, x2, . . . , xn) = 0,
f2(x1, x2, . . . , xn) = 0,
.........
fn(x1, x2, . . . , xn) = 0,
Usando notación vectorial, el sistema asume la forma:
F(x) = 0.
Una de las ventajas del método de Newton–Raphson además de su
velocidad de convergencia, es que se puede generalizar fácilmente a
sistemas de ecuaciones no lineales. Esta generalización se conoce
como método de Newton.
Método de Newton Raphson

El método de Newton consiste en, dada la aproximación de la


solución x(k), tomar como nueva aproximación x(k+1) el valor de la
expresión anterior:

x(k+1) = x(k) − Df(x(k))−1f(x(k)), k = 0, 1, 2, . . .

donde x(0) es la aproximación inicial y:


Método de Newton Raphson

En la práctica no es necesario (ni conveniente) invertir la matriz Df(x


(k)
), sino que se utiliza un método menos costoso: resolver en cada
iteración el sistema de ecuaciones lineales siguiente:

Df(x(k))(x(k+1) − x(k)) = −f(x(k)).

Así, llamando δx(k) = x(k+1) − x(k); las iteraciones se basarán en:

x(k+1) = x(k) + δx(k)


Método de Newton Raphson
Ejemplo:

Resolver el siguiente sistema de ecuaciones con error menor que


tol = 10−5:

Las funciones a utilizar son:


Método de Newton Raphson
Ejemplo:
Método de Newton Raphson
Ejemplo:
Método de Newton Raphson
Ejemplo:

Resultados obtenidos:

También podría gustarte