Resumen Disparo No Lineal

METODO DEL DISPARO NO LINEAL
*INTRODUCCION:

Las ecuaciones diferenciales ordinarias se dividen en dos grupos que
son las lineales y no lineales.

Las ecuaciones lineales son accesibles para encontrar una solucin
analtica mediante mtodos matemticos, en cambio las ecuaciones
no lineales con excepcin de algunas ecuaciones de primer orden,
son generalmente imposibles de resolver en trminos de las
funciones elementales comunes.

-EDO LINEAL:

Las EDO lineales tienen dos caractersticas importantes:

a. La variable dependiente y junto con todas sus derivadas son
de grado 1.
b. Cada coeficiente de y solo depende de x.

-EDO NO LINEAL:

Por ejemplo:

0 2 ' ' '
2
= + y x yy

En este ejemplo vemos que el coeficiente de la tercera derivada es
una funcin que no depende de x; y adems que el exponente de
y es dos.

( ) ( ) ( ) ( ) ) ( ...
0 1 1
1
1
x g y x
y
x
y
x
y
x
a
dx
d
a
dx
d
a
dx
d
a n
n
n n
n
n
= + + + +

0 ' arctan ' ' = +

e
x y y
y

En el ejemplo anterior vemos que es imposible encontrar una
solucin mediante funciones elementales de matemticas, entonces
surge la pregunta; como encontrar una solucin a este problema o
por lo menos una aproximacin a este caso.

Las ecuaciones diferenciales se pueden escribir de esta forma:

La cual es una EDO de orden n.

*ALGORITMO MATEMATICO:

Con el mtodo del disparo no lineal, nos centraremos en encontrar la
solucin o aproximaciones a problemas de contorno o problemas de
valor de frontera de ecuaciones diferenciales de orden 2.

Tenemos:

Con

Este mtodo es parecido al mtodo del disparo lineal, con la
diferencia de que la solucin a un problema lineal se puede escribir
por la combinacin lineal de las soluciones de dos problemas de
valor inicial, entonces se tiene que hacer una serie de soluciones
para encontrar el mejor parmetro de la pendiente y as
reemplazarlo en la condicin inicial, tomamos a ese parmetro como
t.

A continuacin el problema de condicin inicial:

0 ) ,..., ' ' , ' , , (
) (
= y y y y x F
n
) ' , , ( ' ' y y x f y = b x a s s
o = ) (a y | = ) (b y

Entonces se elige un primer parmetro t para iniciar el proceso.

Como el problema anterior es de valor inicial, la solucin de este
problema es

La solucin al problema anterior se puede hallar mediante cualquier
mtodo para EDO, por ejemplo el mtodo de Euler o el mtodo de
Runge-Kutta.

La finalidad es encontrar un parmetro t que al reemplazarlo en el
problema de valor inicial nos de una buena aproximacin al y(b)=.

Se empieza con un parmetro to, la cual determina la elevacin del
disparo, en otras palabras es la pendiente o la condicin inicial para
el y(a).

Si y(b,to) no esta lo suficientemente cerca de , corregimos nuestra
aproximacin tomando nuevas elevaciones o pendientes t1,t2,y as
) ' , , ( ' ' y y x f y = b x a s s
o = ) (a y t a y = ) ( '
) , ( t b y ) , ( t x y
sucesivamente hasta que y(b,tk), acierte en con la precisin
deseada o bajo un nivel de tolerancia.

Entonces el verdadero problema es determinar el parmetro t del
problema de valor inicial, con lo que queda:

Si hacemos:

Ya que b y son conocidos, queda una funcin que depende solo
del parmetro t.

Esto representa encontrar el t que satisfaga esta ecuacin:

Para esto recurrimos a dos mtodos bien conocidos:

>Mtodo de la Secante:

El mtodo de la secante es una forma recursiva para hallar nuevos
puntos, una forma recursiva que depende de dos puntos iniciales, la
forma general para f(x)=0:

Entonces elegimos dos parmetros, to y t1 para iniciar el proceso y
generar as los nuevos puntos para el G(t)=0.
0 ) , ( = | t b y
| = ) , ( ) ( t b y t G
0 ) ( = t G
) ( ) (
) )( (
2 1
2 1 1
1
x
f
x
f
x x x
f
x
x
k k
k k k
k k

=
Esto implica resolver dos problemas de valor inicial, porque
queremos el valor de G(to) y el valor de G(t1), luego reemplazarlo en
la siguiente formula:

Pero como sabemos:

Entonces la formula para hallar el tk es:

Luego se repite el proceso para hallar los siguientes tk, para obtener
la mejor aproximacin t.

>Mtodo de Newton:

El mtodo de Newton es una formula de recurrencia que depende de
un solo punto inicial, en casos generales tomamos la funcin f(x)=0,
entonces tenemos:

Este mtodo es mejor para hallar el parmetro t, ya que es una
forma de recurrencia que depende de un solo punto inicial to, por
) ( ) (
) )( (
2 1
2 1 1
1
t
G
t
G
t t t
G
t
t
k k
k k k
k k

=
| = ) , ( ) (
t
b y
t
G
k k
) , ( ) , (
) )( ) , ( (
2 1
2 1 1
1
t
b y
t
b y
t t t
b y
t
t
k k
k k k
k k

=
|
) ( '
) (
1
1
1
x
f
x
f
x
x
k
k
k k
=
otra parte presenta tambin un problema, la cual es no conocer la
derivada en to de la funcin G(t).

La formula queda as:

Se puede representar as:

El problema ahora es encontrar

Esto representa una dificultad ya que no sabemos una forma
explicita de y(b,t) para encontrar la derivada, en otras palabras solo
contamos con valores que provienen de puntos.

Para afrontar esta dificultad, hacemos que el problema de valor
inicial dependa de x y t, quedando:

Con valores iniciales:

La nueva funcin y(x,t), lo derivamos parcialmente con respecto a t,
por medio de la regla de la cadena.

Y al final nos queda:
) ( '
) (
1
1
1
t
G
t
G
t
t
k
k
k k
=
) , (
) , (
1
1
1
t
b
dt
dy
t
b y
t
t
k
k
k
k

=
|
) , ( t b
dt
dy
)) , ( ' ), , ( , ( ) , ( ' ' t x y t x y x f t x y =
o = ) , ( t a y t t a y = ) , ( '

Derivando las condiciones iniciales con respecto a t:

Las dos expresiones anteriores se asemejan a un problema de valor
inicial y para observarlo mejor hacemos un cambio de notacin:

Ahora la ecuacin diferencial queda:

Con condiciones iniciales:

Por consiguiente para usar el mtodo de Newton es necesario
resolver dos problemas de valor inicial en cada iteracin.

) , (
'
'
) , ( ) , (
' '
t x
t
y
y
f
t x
t
y
y
f
t x
t
y
c
c
c
c
+
c
c
c
c
=
c
c
0 ) , ( = t a
dt
dy
1 ) , (
'
= t a
dt
dy
) , ( ) , ( t x z t x
dt
dy
) , ( '
'
) , ( ) , ( ' ' t x z
y
f
t x z
y
f
t x z
c
c
+
c
c
=
0 ) , ( = t a z 1 ) , ( ' = t a z
Entonces la resolucin de la EDO NO LINEAL por el mtodo de
Newton es resolver:

Y tambin resolver:

Despus hallar el tk:

Y as podemos obtener los siguientes parmetros t.

*EJEMPLO:

Tenemos la EDO NO LINEAL

) ' , , ( ' ' y y x f y = b x a s s o = ) (a y
t
a y
k 1
) ( '

=
' ) ' , , ( ) ' , , ( ' '
'
z y y x f z y y x f z
y y
+ = b x a s s
0 ) , ( = t a z 1 ) , ( ' = t a z
) , (
) , (
1
1
1
t
b z
t
b y
t
t
k
k
k k

=
|
x
y y y
3
3
2 6 2 ' ' =
2 1 s sx
2 ) 1 ( = y 5 . 2 ) 2 ( = y

Lo primero es dar un valor inicial t y resolver el problema de valor
inicial

Y a la vez resolver:

Con condicin inicial:

Y as encontrar el mejor parmetro t con cierta tolerancia y luego
reemplazarlo en el problema de valor inicial, entonces se encuentra
las aproximaciones y(x), para x [1,2].

*ALGORITMO COMPUTACIONAL:

Usa el mtodo de R-K y el mtodo de Newton.

-Datos de entrada:

y; z: funciones
n: numero de particiones
tol: tolerancia
[a,b]: intervalo
y : condiciones de frontera
max: numero mximo de iteraciones

x
y y y
3
3
2 6 2 ' ' =
2 1 s sx 2 ) 1 ( = y t y = ) 1 ( '
' ) 2 6 2 (
'
) 2 6 2 ( ' '
3
3
3
3
z
x
y y
y
z
x
y y
y
z
c
c
+
c
c
=
z y z = ) 6 6 ( ' '
2
0 ) 1 ( = z 1 ) 1 ( ' = z

-Datos de salida:

y(x): aproximaciones de solucin

-Proceso:

h=(b-a)/n
t= (- )/(b-a)
ban=0

k=1;
while (k<=max)
y(a)=
y(a)=t
z(a,t)=0
z(a,t)=1

por medio de R-K resolver el problema de valor inicial de
y;
tambin el problema de valor inicial de z

if (y(n)- )<=tol
ban=1
k=max+1
else
t=t-(y(b)- )/z(b,t)
k=k+1
end
end

if ban==1
Imprimir y(x)
else
no convergi
end

*CONCLUSIONES:

Las ecuaciones diferenciales no lineales se presentan en la mayora
de casos experimentales en la fsica, por eso su importancia ya que
antes no se tena un computador para hacer clculos
extremadamente imposibles con lpiz y papel.

La solucin de estas ecuaciones solo se pueden aproximar, este
mtodo es sensible a errores de redondeo si es que la funcin y(x) es
muy creciente en el intervalo de [a,b]

*BIBLIOGRAFIA:

-Faires Burden, Richard Douglas-Mtodos Numricos.

-Zill, Dennos G.-Ecuaciones Diferenciales.

Resumen Disparo No Lineal

Cargado por

Copyright:

Formatos disponibles

Resumen Disparo No Lineal

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

Resumen Disparo No Lineal

Cargado por

Copyright:

Formatos disponibles

METODO DEL DISPARO NO LINEAL

0 ' arctan ' ' = +

También podría gustarte