Analisis Numerico II Parcial
Analisis Numerico II Parcial
Analisis Numerico II Parcial
DIFERENCIACION NUMRICA o Primera derivada o Segunda derivada INTEGRACIN NUMERICA o Frmulas de Newton-Cotes: Mtodo del Trapecio Mtodo de Simpson o Cuadratura Gaussiana ECUACIONES DIFERENCIALES ORDINARIAS (EDOS) o EDOs de Primer orden o Sistemas de EDOs o EDOs de orden Superior ECUACIONES DIFERENCIALES PARCIALES (EDPS) o Elptica o Parablica o Hiperblica.
Pg. 1
FEPOL 2010-2012
DIFERENCIACION NUMRICA A pesar que derivar una funcin es bastante sencillo, se plantearon mtodos numricos para la diferenciacin con el objeto de resolver ecuaciones diferenciales con derivadas parciales. Los mtodos de diferenciacin de ecuaciones se basan en el desarrollo de Taylor de la funcin a diferenciar. Aproximacin Numrica de la Primera derivada Para el clculo de la primera derivada de una funcin las frmulas se derivan de la siguiente expresin:
( ) ( )[ ( )( ) ] ( )[ ( )( ) ] ( )[ ( )( ] )
Con un valor de
y h(
( )
Con un valor de
y h(
Con un valor de
y h(
( )
Aproximacin Numrica de la Segunda derivada Para el clculo de la segunda derivada de una funcin se emplea la siguiente expresin:
( ) ( ) ( ) ( )
La cual es una frmula de diferencias centradas que considera valores de h (incrementos para la variable independiente)
Considere la tabla de datos de la posicin de una partcula en funcin del tiempo t en segundos x x(t) 0.1 3.46 0.2 4.42 0.3 4.99 0.4 4.65 0.5 3.72
Pg. 2
FEPOL 2010-2012
Utilizando todos los datos y con el menor error posible determine la aceleracin de la partcula en cada uno de los instantes dados Solucin: La aceleracin es la segunda derivada de la funcin. As para los valores centrales los resultados sern los siguientes:
( ) ( ) ( ) ( ( ( ( ) ) ) ) ( )
( ( (
) ) )
Para el clculo de la aceleracin en t=0.1 y t=0.5 seg. es necesario encontrar los valores de la primera derivada para luego volverlos a derivar
( ) ( ) ( ( ) ( ( ) ) ) ( )
) ( ( ) ) ( ( ) ) ( ( ) ( ( ) )
) (
Para la cual aplicamos nuevamente la primera derivada y as obtener la aceleracin a t=0.1 y t=0.5 seg.
( ) ( ( ) ( ) ( ) ) ( ) ) ( )
( ( )
Pg. 3
FEPOL 2010-2012
Considere la tabla de datos de la posicin de una partcula en funcin del tiempo t en segundos x x(t) 1.2 3.46 1.29 4.42 1.30 4.99 1.31 4.65 1.40 3.72
Utilizando todos los datos y con el menor error posible determine la velocidad de la partcula en cada uno de los instantes dados Solucin: La velocidad es la segunda derivada de la funcin. Pero no se pueden usar las formulas simplificadas debido a que los incrementos para la variable independiente (valores de h) son distintos: Es necesario aplicar entonces la frmula general para as encontrar la primera derivada.
( ) ( )[ ( )( ) ] ( )[ ( )( ) ] ( )[ ( )( ] )
( )[
[ ( (
(
) )(
)(
)
)
]
( )[
[ (
(
( )
)(
)(
]
]
( )[
(
( ( )
)(
)(
]
]
( )
( )[
[ ( (
(
)
)(
)(
)
)
( )[
[ (
(
( )
)(
)(
]
]
( )[
(
( (
)(
) )(
]
]
( )
( (
) )(
] )
( (
) )(
] )
( (
) )(
] )
Pg. 4
FEPOL 2010-2012
( ( ) )( ( ) )( ( ( ) )(
] )
[ (
] )
] )
(
) )(
)(
)
)
]
( )[
(
( )
)(
)(
( )[
(
( ( )
)(
)(
[ (
( )
x v(t)
1.2 -31.03
1.29 52.38
1.3 11.5
1.31 -31.63
1.4 10.96
INTEGRACION NUMRICA Suele ser complicado encontrar el valor de una integral ( ) por medio de tcnicas analticas, de ah la importancia de los mtodos numricos para integrar. Frmulas de Newton-Cotes Estos mtodos reemplazan a la funcin ( ) por el polinomio de interpolacin de Lagrange. De acuerdo al desarrollo que se le d a la aproximacin resultante se derivan los mtodos de: Trapecios, Simpson, Regla de los 3/8, entre otros. En general no conviene generar mtodos de Newton-Cotes de grados ms altos debido a la caracterstica oscilante de polinomio de interpolacin de Lagrange. Los mtodos ms usados son el del Trapecio y el de Simpson, pero con una ligera variacin, que consiste en subdividir el intervalo sobre el cual se va a integrar en pequeos sub-intervalos para aplicar en stos el mtodo. As se tienen las llamadas frmulas generalizadas. La frmula generalizada de los trapecios es la siguiente: ( ) ( ( ) ( ) ( )+
Pg. 5
FEPOL 2010-2012
( )
( ( )
( )
) ,
Aproximar el permetro de la regin ubicada en el primer Cuadrante, acotada por los ( ) ejes coordenados y la curva C: ; t * + Utilice El mtodo de los ( ) trapecios con n=8 Solucin: La longitud de curva de una funcin paramtrica est dada por: (
( ))
( ))
( )
( )
( )
( )
( )
( ( )
( )
( )+
Para realizar la aproximacin se debe conocer el valor de h y se tabulan los datos a emplear en la aproximacin.
i 0 1 2 3 4 5 6 7 8 ( ) 1.732050 1.743003 1.776822 1.818971 1.870828 1.921286 1.963046 1.990462 2
Pg. 6
FEPOL 2010-2012
[ ( )
[
( )
(
( ( )
( )
( )
( )
( )
( )
( ))]
)]
)]
, la integral impropia:
Estimar tambin el error del resultado obtenido. Solucin: Para problemas de integrales indefinidas antes de proceder con la aplicacin del mtodo se deben eliminar las indeterminaciones (evaluando los lmites en la funcin a integrar), por lo que suele ser necesario aplicar un cambio de variable y particionar el intervalo sobre el cual se va a integrar. Si slo se realiza el cambio de variable la integral resultante es la siguiente:
( )
Evaluando los lmites en esta nueva integral se puede notar que la indeterminacin aun persiste, por ello se dividen el intervalo de integracin de la forma mostrada
Pg. 7
FEPOL 2010-2012
Para la primera integral no hay problemas de indeterminacin, en cambio en la segunda an est presente este problema, a esta integral se le aplica entonces el cambio de variable anteriormente propuesto, obtenindose ahora:
i 0 1 2 3 4
)]
Haciendo una analoga de la integral calculada recientemente con la primera integral, se puede notar que son idnticas, por esto su aproximacin ser la misma.
Pg. 8
FEPOL 2010-2012
El error es la diferencia entre el valor exacto y el valor aproximado. | El valor exacto se lo encuentra resolviendo la integral. | |
Cuadratura Gaussiana Las frmulas de Newton-Cotes requieren de puntos igualmente espaciados lo cual puede ser una desventaja debido a que no aprovecha las caractersticas de la funcin. Por lo cual se han desarrollado frmulas con el mismo nmero de puntos (pero no igualmente espaciados) para encontrar el valor numrico de la integral de una manera ms precisa, stas frmulas se denominan cuadratura Gaussiana. Para el empleo de ste mtodo se utiliza el siguiente cambio de variable:
( )
n 2 3 1,1
Pg. 9
FEPOL 2010-2012
Calcule el valor de la siguiente integral usando el mtodo de cuadratura gaussiana con n=3.
( (
))
( )
( )
( (
))
) [
( (
),
) ]
[ ( (
))
)]
( ( [
),
) ]
))
( (
))
Integracin Doble. Los mtodos numricos empleados para el clculo de integrales simples pueden generalizarse sin ningn inconveniente para integrales mltiples. Se pueden emplear los mtodos nombrados anteriormente de forma nica o combinados, ) para el clculo sobre regiones (
Pg. 10
FEPOL 2010-2012
Calcular la siguiente integral, con el algoritmo de la integral doble de Simpson: Donde R es la regin acotada por:
Usar n=m=4
Solucin: Para problemas de integrales dobles se debe conocer la regin sobre la cual se est integrando. En este caso esa regin es una circunferencia cuyo radio es 3 y que se encuentra centrada en el origen. Planteando nuevamente la integral
La variable y en cambio estar acotada por la parte superior e inferior de la circunferencia. As para cada valor de habr un incremento de y ( ( Evaluando para ( ) ; ( ) Evaluando para ( ) ( Evaluando para ( ) ; ( )
Pg. 11
) ( )
( )
; )
FEPOL 2010-2012
) ( )
( (
)+
( (
( (
))
( (
)))
( (
) [ [ [
( ( ( ( ( ( ( ( (
( ( ) ) ) ) ( ( ( (
) ) ) ) )
( ( ( ( ( ( ( ( (
)) ) ) ) )
( ( ( ( ( (
Pg. 12
FEPOL 2010-2012
y -3 -1.5 0 1.5 3 27 3.375 0 3.375 27 27 6.084367 0 6.084367 27 27 6.75 0 6.75 27 27 6.084367 0 6.084367 27 27 3.375 0 3.375 27
[ ((
( ( ( ( (( [
)) ( ) ( ))
( )) ( ) ( ( ))] ) ] ( ))
( (
( )
ECUACIONES DIFERENCIALES ORDINARIAS (EDOS) EDOs de Primer orden Se trata de resolver el problema de la forma: ( ( ) Se particiona [ ] con ; ; ) en [ ]
Los mtodos a utilizar para resolver este tipo de ejercicios son de dos tipos: Mtodos de Taylor de orden n Mtodos de Runge-Kutta de orden n
Pg. 13
FEPOL 2010-2012
Mtodos de Taylor de orden n: Desarrollando la funcin desconocida y en serie de Taylor alrededor de manipular la expresin y despreciando el error se llega a: ( ) ( ) ( ) , despus de
La cual se denomina frmula de Taylor de orden n. En la misma cada aproximacin se obtiene a partir de la anterior y se desprecia el residuo en cada paso, es previsible que el error vaya aumentando con cada aproximacin. Particularizando lo anteriormente planteado se tiene: Con n=1 (Mtodo de Euler). ( Con n=2 (Mtodo de Taylor de orden 2). ( ) ( ) )
Nota: En general no es recomendable desarrollar mtodos de Taylor de orden mayor a 3, puesto que la frmula respectiva involucrara el clculo de derivadas de rdenes mayores para la funcin f, lo cual puede tornarse complicado al tratar con dos variables.
a) Establecer el algoritmo correspondiente para la ecuacin dada b) Escribir la tabla de resultados para Solucin: Desarrollando el polinomio de Taylor para n=1 (Euler): ( Se establece la funcin ( ) despejando
( )
) , teniendo entonces: ( )
Pg. 14
FEPOL 2010-2012
( )
Siguiendo con el proceso iterativo y resumiendo estos valores en una tabla, se logra:
i 0 1 2 3 4 5 6 7 8 9 10 1 0.8 0.69184075 0.60974725 0.53078359 0.44716690 0.35773352 0.26459861 0.17151673 0.08289747 0.00312236
Pg. 15
FEPOL 2010-2012
No olvide establecer todos los pasos necesarios para establecer el algoritmo Solucin: Desarrollando el polinomio de Taylor para n=2: ( Siendo: ( ) ( ) ( ) ) ( )
Si se particulariza el algoritmo ( ( ) ( ) * (
( ( * ( * ) ( ( * ( * ) ( )
Pg. 16
FEPOL 2010-2012
( ( ( * ( * ) (
( )
( ) )
( )(
( )
( ) )
( ( ))
i 0 1 2 3 4 5
Mtodos de Runge-Kutta de orden n: El mtodo de Runge-Kutta de orden n permite resolver una ecuacin deferencial con una precisin equivalente al mtodo de Taylor del mismo orden, pero con la ventaja de que no involucra el clculo de derivadas de la funcin. Entre los ms utilizados s encuentran los mtodos de Runge-Kutta de orden dos y de orden cuatro. El algoritmo del mtodo de Runge-Kutta de cuarto orden es el siguiente: ( ( ( ( ( ) ) * * )
2y
4 x2 dx xydy 0 , y(1) 2 , 1 x 2
Usando el mtodo de Runge-Kutta de cuarto orden: a) Escriba el algoritmo para la funcin especfica f ( x, y ) . b) Escriba la tabla de resultados para h 0.2 .
Pg. 17
FEPOL 2010-2012
) )
( ( )
( ( ( ( ( ( ( ( (
) )( ) )( ) )( (
( ) ( )
) )
) )
) ) )
( ( ( ( ( ( ) )(
) )( (
( ) ) ) )
) )
( (
Pg. 18
FEPOL 2010-2012
( ( ( ) )( ( ( ) ( ) ) ( ( ( ) ( )( ( ( ) )) ( )) ) ) )
i 0 1 2 3 4 5
1 x y ' xy x 1 x , y(0) 2 , 0 x 1 2
2 2
Usando el mtodo de Runge-Kutta de cuarto orden: Escriba el algoritmo para la funcin especfica f ( x, y ) . Escriba la tabla de resultados para h 0.1 . Solucin: Convirtiendo el problema anterior a la forma conocida, se tiene: ( ) ( (
( )
) )
FEPOL 2010-2012
( (
( (
)( ( ) )
) )
( ( ((
( ( ( (
)( ( )( ( ) ) ) )
) )
i 0 1 2 3 4 5
( ) en [ ]
Pg. 20
FEPOL 2010-2012
Los algoritmos para la resolucin de este tipo de problema son los siguientes: Euler:
Taylor 2:
Runge-Kutta de orden 4:
( ( ( ( ( )
) ) )
FEPOL 2010-2012
( ( ( ( ( )
) ) )
Para j=0
( Para i=1
( ) ( )
Para i=2
( ) ( ( ))
Para i=1
( *( ( )* ( ( )* ) ( )+
Para i=2
( ) (( )( ( )*)
Para i=1
( *( ( )* ( ( ) )* ( )+
Para i=2
Pg. 22
FEPOL 2010-2012
((
)(
)*)
Para i=1
( [( ( )) ( ( ) )) ( )]
Para i=2
( ) (( )( ( )))
Valores de
Iteracin 1 2 3 4 5 j 1 -0.2 0 -0.28791093 -0.00970586 -0.39046479 -0.04645867 -0.50355755 -0.12325638 -0.61915251 -0.25419927 2 -0.24 -0.002 -0.33573144 -0.02319612 -0.44486540 -0.07759658 -0.56158767 -0.17904814 -0.67564784 -0.34169791 3 -0.2438 -0.0024 -0.33916446 -0.02463074 -0.44719168 -0.08031661 -0.56181150 -0.18311025 -0.67254750 -0.34678249 4 -0.28828 -0.009752 -0.39081767 -0.04654490 -0.50383981 -0.12335101 -0.61929780 -0.25423168 -0.72430551 -0.4522589
Pg. 23
FEPOL 2010-2012
( ) La estrategia es convertir a esta ecuacin diferencial de orden superior en un sistema de ecuaciones de primer orden. Se aplica el cambio de variable:
Finalmente el problema se torna en la resolucin de un sistema de ecuaciones diferenciales de primer orden, de la forma: ( ) ( )
( )
La ecuacin de un movimiento angular est dada por y+10sen(y)=0; 0t1; y(0)=0; y(0)=0.1. Empleando el mtodo de Runge-Kutta de 4 orden generalizado y un paso de 0.20, aproximar la solucin de la ecuacin en t =0.40 Solucin: Aplicando el cambio de variable se tiene:
Pg. 24
FEPOL 2010-2012
Con el problema convertido en un sistema de edos de primer orden se emplea lo anteriormente explicado:
( ( ( ( ( )
) ) )
Para j=0
( Para i=1
( ) ( )
Para i=2
( ) ( ( ))
Para i=1
( ) * ( )+
Para i=2
( ) ( ( ( ))*
Pg. 25
FEPOL 2010-2012
Para i=1
( ) * ( )+
Para i=2
( ) ( ( ( ))*
Para i=1
( ) [ ( )]
Para i=2
( ) ( ( ( )))
Valores de
Iteracin 1 2 j 1 0.02 0 0.01613344 -0.03733120 2 0.02 -0.01999966 0.01240031 -0.05346044 3 0.01800003 -0.01999966 0.01078739 -0.04972857 4 0.01600006 -0.03599812 0.00618772 -0.05889965
Nota: Las aproximaciones de la ecuacin diferencial de segundo orden son las respuestas de la aproximacin de es decir 0.01866668 y 0.03011612
Pg. 26
FEPOL 2010-2012
Resolucin numrica de problemas con valores en la frontera Se plantea u problema de la siguiente forma: ( ) ( ) ( ) Estos problemas estn asociados a mltiples aplicaciones fsicas, en particular la deflexin de vigas sujetas en los extremos sobre las cuales se ejerce una carga uniforme. Se particiona [ ] con ; ( ) ( ) ; [ ]
Usando las aproximaciones de la primera derivada y despreciando el residuo se puede ) ecuaciones con ( ) incgnitas, cuyo planteamiento llegar a un sistema de ( matricial es el siguiente:
Siendo la matriz A:
( ( )) ( ) ( ( ) ( )) ( ) ( ( ) ( )) ( [ ( ) ) ( ( ( ) )) ( ( ( ) ))]
( ) ( ) ( ) ( [ ( ) )]
( )
( )) ( ) ( )
( [ ( )
(
) ( )) ]
Pg. 27
FEPOL 2010-2012
Solucin: Se debe llevar el problema a la forma conocida para hacer una semejanza con la forma general e identificar ( ) ( ) ( ). Las condiciones de frontera son y , valores ) respectivamente. que para el presente ejercicio son -1 y (
( )
( )
( )
, ( ( ) )
) )
( ))
( ) (
( )
( ))
( ) (
( )
( ))
( ) (
( )
( ))]
( (
((
(
))
)) (
)) )
((
( (
))
)) (
)) )
((
( (
))
)) (
)) )
((
))]
Pg. 28
FEPOL 2010-2012
( )
( )) ( ) ( )
( ( (
( )
( )) ]
( ) ) (
))(
))(
)]
ECUACIONES DIFERENCIALES PARCIALES EDP Elptica Se trata de resolver el problema de distribucin de calor en una regin plana ( ( La funcin ( ) ( ) ) ( ( ) )
La principal aplicacin fsica es que esta ecuacin permite determinar la temperatura ( ) en cualquier punto ( ) interior de una placa, conociendo cual es la temperatura en los puntos frontera de la placa. Se particiona la regin (la cual no necesariamente ser rectangular) en intervalos iguales para la variable y para la variable Para los incrementos se usan las siguientes variables: . y
Pg. 29
FEPOL 2010-2012
Usando el mtodo de DIFERENCIAS CENTRADAS se denota el siguiente sistema de ( ) ( ) ( ) ecuaciones con ( ) incgnitas * ( * + ( ) ( * ( ) ( )
( (
) ) (
( ) )
( ) ( )
( (
) )
En el eje x tomar h=1/4, y en el eje y k=1/3. Solucin: Observando los datos se puede notar que la regin sobre la cual se trabaja es una ) regin cuadrada y la funcin (
Los puntos donde se desea conocer los resultados son obtenidos as: , ( * ( * ( *
Pg. 30
FEPOL 2010-2012
( *
) ec (1)
FEPOL 2010-2012
) ec (2)
) ec (3)
) ec (4)
) ec (5)
Pg. 32
FEPOL 2010-2012
) ec (6)
[ Finalmente se obtiene:
La placa plana mostrada en la figura est construida con cierto metal, y se ha determinado que la temperatura en los bordes de la placa es la que se indica en la figura. Adems se tiene que el trmino no homogneo asociado a la ecuacin elptica ( ), respectiva es f(x,y)=20. El problema consiste en determinar la
temperatura en los puntos del interior de la malla que se muestra en la figura. Determinar el algoritmo en diferencias finitas que resuelve el problema. Plantear el sistema de ecuaciones lineales que resuelve el problema. Utilice el mtodo de Gauss para resolver el sistema de ecuaciones generado.
Pg. 33
FEPOL 2010-2012
2.5
1.5
La temperatura se incrementa linealmente de 0C a 100C
Solucin: Ha pesar de no tratarse de una placa de forma rectangular el planteamiento para la solucin de este problema sigue siendo el mismo. Se desea saber la temperatura en los nodos interiores de la placa:
Pg. 34
FEPOL 2010-2012
Pg. 35
FEPOL 2010-2012
( [ ] [ ] [
La matriz A es diagonal estrictamente dominante, as que el sistema puede ser resuelto por un mtodo iterativo. Utilizando el mtodo de Gauss-Seidel y con se tiene:
k 0 1 2 3 4 5 6 7
1 65.666666 75.145833 81.173177 82.303304 82.515202 82.554933 82.562383 1 34 43.479166 49.506510 50.636637 50.848536 50.888267 50.895716 1 38.916666 63.026041 67.546549 68.394144 68.553068 68.582867 68.588454 1 78.479166 84.506510 85.636637 85.848536 85.888267 85.895716 85.897113
[ [ ]
Pg. 36
FEPOL 2010-2012
En las ecuaciones parablicas se tiene una variable espacial y otra temporal, en la denominada ecuacin de difusin.
( (
) )
( )
La funcin (
de la barra en el instante .
Esta ecuacin representa la difusin de calor en una barra de longitud, en el tiempo. En los extremos de la barra en cualquier instante la temperatura de la barra es la misma. En el instante inicial se conoce exactamente la temperatura de la barra. Se particiona con y el tiempo con
Donde:
, [ y,
( ) ( ) ( [ ( ) )]
con
Pg. 37
FEPOL 2010-2012
Dada la ecuacin:
( ( Siendo, ( ) Aproximar (
) )
( )
Solucin: Los puntos que se desea evaluar se los obtiene de la siguiente manera:
Donde ejercicio
, para este
FEPOL 2010-2012
( *
) ( ) ( ) + [ ]
][
], para t=0.1seg
][
], para t=0.2seg
][
], para t=0.3seg
EDP Hiperblica Al igual que en las ecuaciones parablicas se tiene una variable espacial y otra temporal, a este tipo de ecuacin se le conoce tambin como ecuacin de la onda.
( ( ( La funcin (
) ) )
( ) ( )
Los extremos de la cuerda estn sujetos en cualquier instante t. En el instante inicial se conoce exactamente la forma de la cuerda. Se particiona con y el tiempo con
Pg. 39
FEPOL 2010-2012
Donde: ( ) ( ) ( [ teniendo,
( ) ( ) ( [ ( ) )] y ( [ ( ( ) ( ) ) )]
) ( ) ( )]
con
( ( ( Aproximar (
) ) )
( (
) ) ( )
Solucin: La cuerda tiene longitud 1m y los puntos en los que se desea saber la posicin en cada instante se los obtiene de la siguiente manera: Al ser h el incremento para la variable espacial (x) y conocer las condiciones de frontera los valores restantes resultan de , teniendo asi: ( ( ) ) ( ( ) )
Pg. 40
FEPOL 2010-2012
Antes de avanzar con el procedimiento se deben reconocer cada una de las funciones, as: ( ) ( ) ( ) ( )
( ( ( (
( ( ( (
)) )) )) ))]
+ [
( ( ( (
*
( ( ( (
)) )) )) ))]
+
( ( [ ( (
) ) ] ) )
( ( [ ( (
) ) ] ) )
+ , para t=0.2seg
Con lo realizado se puede plantear el algoritmo para encontrar cada posicin a diferentes instantes de tiempo.
este ejercicio
FEPOL 2010-2012
]*
+, para t=0.4seg
]*
+, para t=0.6seg
]*
+, para t=0.8seg
Estos valores de representan las posiciones de cada uno de los puntos que interesan en el instante t = 0.8 segundos
Pg. 42