Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Metodo Taylor

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 4

Elaborado por J.C.

Alcerro
MÉTODO DE TAYLOR DE ORDEN SUPERIOR
1. Aplique el método de Taylor de orden 4 para aproximar las soluciones del siguiente problema
de valor inicial:
y  sen t  et , 0  t  1, y(0)  0, con h  0.5
Solución
Para obtener las soluciones se aplicará la ecuación (5.17) del libro de texto:
w0  
wi 1  wi  hT ( n ) (ti , wi ), i  0,1,..., N  1 (I)
donde
T ( n) (ti , wi )  f (ti , wi )  h2 f (ti , wi ) 
n1
 hn! f ( n1) (ti , wi ) (II)

con f (t , y)  sen t  et , h  0.5, n  4 y las n  1 derivadas implícitas


f (t , y )  cos t  e t
f (t , y )   sen t  e  t
f (t , y )   cos t  e t
Sustituyendo n  4 en (II) se obtiene
T (4) (ti , wi )  f (ti , wi )  h2 f (ti , wi )  h3! f (ti , wi )  h4! f (ti , wi )
2 3

y sustituyendo f y sus derivadas


T (4) (ti , wi )  sen ti  eti  0.5  ti  ti
2 [cos ti  e ]  3! [ sen ti  e ]  4! [ cos ti  e ]
0.52 0.53  ti

 0.24479167 cos ti  0.95833333sen ti  0.78645833eti

Luego, al sustituir T (4) (ti , wi ) y h en (I) se obtiene


w0  y (0)  0

wi 1  wi  0.5 0.24479167 cos ti  0.95833333sen ti  0.78645833eti  (III)

Con ti  a  hi  0  0.5i se obtiene que t0  0 , t1  0.5 y t2  1.

Seguidamente las aproximaciones se obtienen sustituyendo los valores ti y wi correspondientes


en (III):

Si t0  0 y w0  0  w1  w0  0.5 0.24479167 cos t0  0.95833333sen t0  0.78645833e t0 
 0  0.5  0.24479167 cos 0  0.95833333sen 0  0.78645833e  0

 0.515625
Si t1  0.5 y w1  0.515625


 w2  w1  0.5 0.24479167 cos t1  0.95833333sen t1  0.78645833e t1 

 0.515625  0.5 0.24479167 cos 0.5  0.95833333sen 0.5  0.78645833e0.5 
 1.09126773
Elaborado por J.C. Alcerro
En resumen se obtiene que las aproximaciones de las soluciones son

i ti wi  y(ti )
0 0 0
1 0.5 0.515625
2 1 1.09126773

2. Dado el problema de valor inicial


1 y
y  2   y 2 , 1  t  2, y (1)  1,
t t
1
con la solución exacta y (t )   .
t
a) Aplique el método de Taylor de orden dos con h  0.05 para aproximar la solución y
compárela con los valores reales de y.
b) Use las respuestas obtenidas en el inciso (a) y la interpolación lineal para aproximar
y(1.052) , y después compárelo con el valor real.
Solución
a) Para obtener las soluciones se aplicará la ecuación (5.17) del libro de texto:
w0  
wi 1  wi  hT ( n ) (ti , wi ), i  0,1,..., N  1 (I)
donde
T ( n) (ti , wi )  f (ti , wi )  h2 f (ti , wi ) 
n1
 hn! f ( n1) (ti , wi ) (II)

1 y
con f (t , y)  2
  y 2 , h  0.05, n  2 y las n  1 derivadas implícitas
t t
2 y y
f (t , y )   3  2   2 yy
t t t
2 y 1 
  3  2  y   2 y 
t t t 
 
2 y 1 y 1 
f (t , y )   3  2   2   y 2    2 y 
t t  t t  t



Sustituyendo n  2 en (II) se obtiene
T (2) (ti , wi )  f (ti , wi )  h2 f (ti , wi )
y sustituyendo f y sus derivadas
1 w  2 wi  1 wi 2  1 
T (2) (ti , wi )  2  i  wi 2  0.05 2   3
 2
  2
  wi   2 wi 
ti ti  ti ti  ti ti  ti 
Elaborado por J.C. Alcerro

1 wi  2 wi  1 2wi wi 2wi 2 wi 2 
 2
  wi
2
 0.025  3  2   3  2  2    2wi 3  
ti ti  t i t i  ti ti ti ti ti 
1 wi 0.075 wi 2
T (2) (ti , wi )    0.05 wi
3
 wi
2
  0.075
ti 2 ti ti 3 ti

Luego, al sustituir T (2) (ti , wi ) y h en (I) se obtiene


w0  y (1)  1
1 w 0.075 w2 
wi 1  wi  0.05  2  i  0.05wi 3  wi 2  3  0.075 i  (III)
 ti ti ti ti 

Con ti  a  hi  1  0.05i se obtiene que t0  1 , t1  1.05, t2  1.1, t3  1.15, t4  1.2, t5  1.25,


t6  1.3, t7  1.35, t8  1.4, t9  1.45, t10  1.5, t11  1.55, t12  1.6, t13  1.65, t14  1.7, t15  1.75, t16  1.8,
t17  1.85, t18  1.9, t19  1.95 y t20  2 .

Seguidamente las aproximaciones se obtienen sustituyendo los valores ti y wi correspondientes


en (III):
 1 w 0.075 w0 2 

Si t0  1 y w0  1  w1  w0  0.05 2  0
 0.05w0  w0  3  0.075
3 2

 t0 t0 t0 t0 
 
 1 1 0.075 ( 1) 2 
 1  0.05  2   0.05(1)  (1)  3  0.075 
3 2

1 1 1 1 

w1  0.9525
1 w 0.075 w2 
Si t1  1.05 y w1  0.9525  w2  w1  0.05  2  1  0.05w13  w12  3  0.075 1 
 t1 t1 t1 t1 

 1 0.9525
 0.9525  0.05    0.05(0.9525)3
 1.052 1.05

0.075 (0.9525) 2 
 (0.9525) 2   0.075  
1.053 1.05 

w2  0.9093138

De igual manera se obtienen w3 , w4 ,..., w20 . En resumen se obtiene que las aproximaciones de las
soluciones, las soluciones reales y los errores reales son
Elaborado por J.C. Alcerro
1
i ti wi  y(ti ) y (ti )   wi  y (ti )
t
0 1 1 −1 0
1 1.05 0.9525 −0.95238095 0.000119
2 1.1 0.9093138 −0.90909091 0.0002238
3 1.15 0.8698799 −0.86956522 0.0003147
4 1.2 0.8337302 −0.83333333 0.0003147
5 1.25 0.80047127 −0.8 0.0003969
6 1.3 0.76977013 −0.76923077 0.00047127
7 1.35 0.74134301 −0.74074074 0.00060227
8 1.4 0.7149466 −0.71428571 0.00066089
9 1.45 0.6903711 −0.68965517 0.00071593
10 1.5 0.66743464 −0.66666667 0.00076798
11 1.55 0.64597879 −0.64516129 0.0008175
12 1.6 0.6258649 −0.625 0.0008649
13 1.65 0.60697107 −0.60606061 0.00091046
14 1.7 0.58918976 −0.58823529 0.00095447
15 1.75 0.57242569 −0.57142857 0.00099712
16 1.8 0.55659418 −0.55555556 0.00103862
17 1.85 0.54161965 −0.54054054 0.00107911
18 1.9 0.52743452 −0.52631579 0.00111873
19 1.95 0.51397809 −0.51282051 0.00115758
20 2 0.50119576 −0.5 0.00119576

b) y(1.052)  ?

Nótese que t  1.052 está contenido en el intervalo [1.05, 1.1] . Entonces considerando que
t0  1.05 y t1  1.1 , cuyas imágenes aproximadas son f (1.05)  0.9525 y f (1.1)  0.9093138
(tomados de la tabla anterior), y construyendo un polinomio de Lagrange de grado uno, se
obtiene que
P1 (t )  L0 (t )  f (t0 )  L1 (t )  f (t1 )
 1.052  1.1   1.052  1.05 
y (1.052)  P1 (1.052)    (0.9525)   (0.9093138)
 1.05  1.1   1.1  1.05 
   
 0.9144  0.03637255
y (1.052)  0.9507725

También podría gustarte