Diferencia finita

De Wikipedia, la enciclopedia libre Saltar a: navegacin, bsqueda Una diferencia finita es una expresin matemtica de la forma f(x + b) f(x +a). Si una diferencia finita se divide por b a se obtiene una expresin similar al cociente diferencial, que difiere en que se emplean cantidades finitas en lugar de infinitesimales. La aproximacin de las derivadas por diferencias finitas desempea un papel central en los mtodos de diferencias finitas del anlisis numrico para la resolucin de ecuaciones diferenciales.

Contenido
[ocultar]

1 Diferencias anterior, posterior y central 2 Relacin con las derivadas 3 Clculo de diferencias finitas 4 Derivadas de rdenes mayores 5 Mtodos de diferencias finitas 6 Vase tambin 7 Referencias 8 Bibliografa complementaria

[editar] Diferencias anterior, posterior y central

Diferencias finitas. Slo se consideran normalmente tres formas: la anterior, la posterior y la central. Una diferencia progresiva, adelantada o posterior es una expresin de la forma

Dependiendo de la aplicacin, el espaciado h se mantiene constante o se toma el limite h 0. Una diferencia regresiva, atrasada o anterior es de la forma

Finalmente, la diferencia central es la media de las diferencias anteriores y posteriores. Viene dada por

[editar] Relacin con las derivadas

La derivacin de la funcin f en un punto x est definida por el lmite

Si h tiene un valor fijado no nulo, en lugar de aproximarse a cero, el trmino de la derecha se convierte en

Por lo tanto, la diferencia anterior dividida por h aproxima a la derivada cuando h es pequeo. El error de esta aproximacin puede derivarse del teorema de Taylor. Asumiendo que f es continuamente diferenciable, el error es

La misma frmula es vlida en la diferencia posterior:

Sin embargo, la diferencia central lleva a una aproximacin ms ajustada. Su error es proporcional al cuadrado del espaciado (si f es dos veces continuamente diferenciable).

[editar] Clculo de diferencias finitas

La diferencia anterior puede considerarse un operador diferencial que hace corresponder la funcin f con f. El teorema de Taylor puede expresarse por la frmula

Donde D denota el operador derivada, que hace corresponder con su derivada decir, Formalmente, invirtiendo la exponencial,

, es

Esta frmula sigue siendo vlida en el sentido de que ambos operadores dan el mismo resultado cuando se aplican a un polinomio. Incluso para funciones analticas, las series de la derecha no convergen con seguridad, sino que puede tratarse de una serie asinttica. Sin embargo, pueden emplearse para obtener aproximaciones ms precisas de la derivada. Por ejemplo, Los dos primeros trminos de la serie llevan a:

El error de la aproximacin es del orden de h2. Las frmulas anlogas para los operadores posterior y central son

[editar] Derivadas de rdenes mayores

De forma anloga se pueden obtener aproximaciones en diferencias finitas para derivadas de orden mayor y operadores diferenciales. Por ejemplo usando la frmula de la diferencia central mostrada anteriormente con un espaciado de para y y aplicando la frmula de diferencia central a la derivada de en x, obtenemos la aproximacin de la diferencia central de la segunda derivada de f:

[editar] Mtodos de diferencias finitas

Otro aspecto importante es que las diferencias finitas aproximan cocientes diferenciales a medida que h se acerca a cero. As que se pueden usar diferencias finitas para aproximar derivadas. Esta tcnica se emplea a menudo en anlisis numrico, especialmente en ecuaciones diferenciales numricas ordinarias, ecuaciones en diferencias y ecuacin en derivadas parciales. Los mtodos resultantes reciben el nombre de mtodos de diferencias finitas.

Las aplicaciones habituales de los mtodos de diferencias finitas son en los campos de la computacin y reas de la ingeniera como ingeniera trmica o mecnica de fluidos.

Mtodo de las diferencias finitas

De Wikipedia, la enciclopedia libre Saltar a: navegacin, bsqueda Plantilla:Mate-esbozo En anlisis numrico, el mtodo de las diferencias finitas es un mtodo utilizado para calcular de manera aproximada las soluciones a las ecuaciones diferenciales usando ecuaciones diferenciales finitas para aproximar derivadas.

[editar] Ejemplo bsico de ecuacin de diferencias finitas en economa

Una ecuacin sencilla en diferencias finitas

La solucin se ensaya por tanteo o aproximacin

Sustituyendo en la ecuacin inicial

La solucin ser

Resolvemos

Comprobamos si la solucin es correcta

Escribimos la solucin general

expresa una combinacin lineal de la solucin Si analizamos el Wronskiano de soluciones particulares obtendremos para t=0 y t=1

Si el Wronskiano es cero, no podemos determinar una solucin correcta. El mtodo para resolver

es idntico pero la solucin general se escribe en funcin del nmero e.

[editar] Bibliografa
METODO DE DIFERENCIAS FINITAS INTRODUCCION El mtodo de diferencias finitas es un clsica aproximacin para encontrar la solucin numrica de las ecuaciones que gobiernan el modelo matemtico de un sistema continuo. Es valioso familiarizarse con sta aproximacin porque tal conocimiento reforzar la comprensin de los procedimientos de elementos finitos. Bsicamente, en una solucin por diferencias finitas, las derivadas son reemplazadas por aproximaciones en diferencias finitas, convirtiendo entonces un problema de ecuaciones diferenciales en un problema algebraico fcilmente resoluble por medios comunes (especialmente matriciales). METODO DE EXPANSION DE TAYLOR El mtodo de expansin de Taylor es una forma alternativa de obtener aproximaciones de diferencia. Este mtodo no solo deduce las frmulas de diferencia sistemticamente, sino que tambin deduce los trminos de error. Para una derivada de p-simo orden, el nmero mnimo de puntos de datos requeridos para deducir una aproximacin de diferencia es , as por ejemplo una aproximacin de diferencia para la primera derivada de una funcin necesita por lo menos de dos puntos de datos. Consideremos la deduccin de la aproximacin de diferencia para en trminos de . La expansin de Taylor de alrededor de es (1). Resolviendo la ecuacin anterior para la primera derivada, tenemos

(2). Si ignoramos todos los trminos con excepcin del primero del miembro derecho de la ecuacin (2), obtendremos la aproximacin por diferencia hacia adelante. Los trminos que se ignoran constituyen el error de truncado, representado por el trmino inicial, . Los dems trminos desaparecen ms rpidamente que el inicial cuando disminuye. La aproximacin de diferencia hacia adelante, con el error de truncado incluido, se expresa como (3), dnde . El trmino indica que el error es aproximadamente proporcional al intervalo de la retcula . El error tambin es proporcional a la segunda derivada . De la misma manera podemos expandir alrededor de en la forma (4), y resolviendo nuevamente para la primera derivada, tenemos

y aqu de la misma manera (5), dnde . Esta aproximacin se denomina de diferencia hacia atrs. Tomemos ahora ambas aproximaciones y restemos (4) de (1):

(6), expresin de la cual se ha eliminado el trmino . Resolviendo para , obtenemos (7). Con el trmino de error incluido, la aproximacin de diferencia central se expresa como

(8), dnde .

Resulta interesante observar que gracias a la cancelacin del trmino , el error de la aproximacin es proporcional al cuadrado de y no a . Entonces, reduciendo reducimos el error con mayor rapidez que con las otras aproximaciones. Como ya se expuso, una aproximacin de diferencia de requiere al menos puntos de datos. Aumentando el nmero de puntos de datos puede obtenerse una aproximacin de diferencia mas exacta. Como ilustracin de lo anterior, deduciremos una aproximacin de diferencia para la primera derivada utilizando tres puntos de datos , de modo que tenemos un punto mas del mnimo requerido. Las expansiones para se escriben:

(9).

(10). Con stas dos ecuaciones es posible cancelar los trminos de la segunda derivada, de modo que el trmino inicial de los errores de truncado es el trmino de la derivada de tercer orden. Por otro lado, si se eliminaran los trminos de la tercera derivada de las ecuaciones (9) y (10) en lugar de los de la segunda derivada, la aproximacin de diferencia obtenida sera menos exacta porque el trmino del error inicial sera de segundo orden en lugar de ser de tercer orden. Multiplicado la (9) por 4 y restndole la (10), obtenemos

(11). Resolviendo para :

(12), dnde el trmino de error est dado por . La (12) es la aproximacin de diferencia hacia adelante de tres puntos. Su error es del mismo orden que el de la aproximacin por diferencia central de dos puntos. Anlogamente, la aproximacin de diferencia hacia atrs de tres puntos puede deducirse utilizando

(13), dnde . Las aproximaciones de diferencia para la segunda derivada se deducen aplicando el mismo principio, el cual consiste en eliminar la primera derivada y el mayor nmero posible de derivadas de orden dos superior. Como ilustracin deduciremos la aproximacin de diferencia para en trminos de . Las expansiones de Taylor de y estn dadas por las ecuaciones (4) y (1) respectivamente. Sumando ambas obtenemos:

de forma equivalente . Entonces si truncamos despus del trmino y reacomodamos los trminos tendremos (14). La ecuacin anterior es la aproximacin de diferencia central para , dnde el error est representado por . Podemos deducir otra aproximacin de diferencia para en trminos de (el nmero mnimo de puntos de datos para es 3). Si multiplicamos por 2 la expansin de Taylor de y la restamos de , el resultado ser: . Resolviendo la anterior para la segunda derivada:

(15), en la cual el error est dado por . La ecuacin (15) es la aproximacin de diferencia hacia atrs para . El orden de su error de truncado es menor que el de la aproximacin de diferencia central,

dada por (14). De ste modo la mayor exactitud pertenece a la aproximacin de diferencia central. De forma similar podemos obtener aproximaciones de diferencia para derivadas superiores, pero la deduccin se hace cada vez mas laboriosa al aumentar tanto el nmero de trminos como el orden de la derivada. Sera til por lo tanto el desarrollo de algoritmos computacionales que permitan hallar automticamente la aproximacin de diferencia para un conjunto dado de datos. No obstante, seguidamente damos las expresiones de diferencias, cuyo uso es frecuente. Primera derivada Aproximaciones de diferencia hacia adelante

Aproximaciones de diferencia hacia atrs

Aproximaciones de diferencia centrales

Segunda derivada Aproximaciones de diferencias hacia adelante

Aproximaciones de diferencia hacia atrs

Aproximaciones de diferencia centrales

Tercera derivada Aproximaciones de diferencia hacia adelante

Aproximaciones de diferencia hacia atrs

Aproximaciones de diferencia centrales

APROXIMACION DE DIFERENCIA PARA DERIVADAS PARCIALES Las frmulas de aproximacin de diferencia para derivadas parciales de funciones multidimensionales son esencialmente iguales a las de diferenciacin de funciones unidimensionales. Consideremos una funcin bidimensional . La aproximacin de diferencia para la derivada parcial con respecto a , por ejemplo, puede deducirse fijando en un valor constante y considerando como una funcin unidimensional. Por tanto, las aproximaciones de diferencia hacia adelante, central y hacia atrs para stas derivadas parciales se pueden escribir, respectivamente:

(16). Las aproximaciones de diferencia central para las segundas derivadas de en estn dadas por:

(17). DIFERENCIAS FINITAS EN UNA DIMENSION Supongamos estar frente a un simple problema unidimensional de contorno, esto es, queremos determinar una funcin , la cual satisfaga una ecuacin diferencial dada en una regin , junto con condiciones de contorno apropiadas es y .

Como ejemplo, consideremos el anlisis de una barra uniforme (mdulo elstico longitudinal y rea de seccin transversal ) como la mostrada en la figura. La ecuacin diferencial que corresponde a la formulacin de ste problema es

(18), con las siguientes condiciones de contorno: Tomamos por simplicidad la funcin de carga longitudinal (variacin lineal).

Para resolver el problema va diferencias finitas, comenzamos por diferenciar la variable independiente , esto es, construimos un conjunto ( grilla malla) de puntos de grilla discretos, igualmente espaciados sobre el rango ( dominio) , dnde . El siguiente paso es reemplazar los trminos de la ecuacin diferencial que involucran diferenciacin por trminos que involucren solo operaciones algebraicas. Este proceso, necesariamente, involucra una aproximacin y puede lograrse haciendo uso de aproximaciones de diferencias finitas (deducidas anteriormente por medio de las expansiones de Taylor). Sustituyendo la aproximacin de diferencia central de la segunda derivada en un punto en (18), obtenemos:

(19), dnde es la carga en el punto de grilla y puede pensarse como la carga total aplicada sobre la estacin de diferencia finita. Tomando ahora las condiciones de contorno: (20)

(21), en la cual hemos tomado las estaciones extremas y la aproximacin de diferencia central para la primera derivada. El punto de grilla en slo se coloca con el fin de imponer la condicin de contorno. Para la solucin por diferencias finitas aplicamos (19) a todas las estaciones y utilizando las condiciones de contorno anteriores, obtenemos:

(22) (los elementos no mostrados de la matriz son nulos). Aqu y . De modo matricial podemos escribir (22) de la forma (23), dnde evidentemente:

. La ecuacin (22) es idntica a la que se hubiera derivado utilizando una serie de n elementos de resorte, cada uno de rigidez . Las cargas en los puntos de grilla correspondientes a se obtendran usando el valor de carga distribuida en el punto de grilla l y multiplicando ese valor por la longitud de contribucin (h para los puntos de grilla internos y para el punto de grilla final). Tal vez con ste ejemplo no se aprecie la utilidad de las diferencias finitas, pues la naturaleza de la formulacin diferencial hace que su resolucin analtica sea viable por mtodos de uso comn. No obstante, lo importante de recalcar y que es conclusin general es que hemos reemplazado un problema de determinacin de una funcin continua desconocida por un problema de resolucin de una ecuacin matricial para un conjunto de valores discretos . Esta es la esencia del mtodo. Debe recordarse que la solucin slo aproxima a la solucin exacta del problema porque hemos reemplazado derivadas por diferencias.

La solucin exacta corresponde a: . Se deja como ejercicio plantear el problema con nmeros crecientes de puntos de grilla y ver como evoluciona el error comparando los resultados obtenidos con los que se obtienen a partir de la solucin exacta. Es evidente que el error decrece a medida que se aumenta el nmero de puntos de grilla. Esto es conclusin inmediata de la formulacin de las aproximaciones de diferencia por medio de las expansiones de Taylor. DIFERENCIAS FINITAS EN MAS DE UNA DIMENSION El problema de aproximacin de ecuaciones diferenciales en dos ms variables independientes es obviamente un poco ms comprometido, aunque los principios utilizados son idnticos a los de una dimensin. Consideremos un problema de torsin elstica de una barra prismtica (regin rectangular) , regido por la ecuacin diferencial siguiente:

(24). Aqu es el mdulo elstico transversal , dnde es el mdulo elstico longitudinal y es la relacin de Poisson; es el ngulo de torsin de cada seccin y es la funcin de tensin que satisface la condicin en los contornos. El momento torsor est dado por y la tensin tangencial en una direccin cualquiera en la seccin se obtiene a partir de . Para aplicar el mtodo de diferencias finitas en sta situacin, procedemos exactamente de la misma manera que en el caso unidimensional. A tal fin, construimos un conjunto de puntos de grilla , igualmente espaciados en el rango con , y tambin un conjunto de puntos de grilla igualmente espaciados sobre el rango

, con . La regin en la cual se requiere la solucin est entonces cubierta por una grilla rectangular de diferencias finitas, a travs del trazado de lneas paralelas al eje a travs de cada punto ; y de la misma forma, trazando paralelas al eje a travs de cada punto . Un punto tpico de grilla est dado entonces por las coordenadas . El mtodo de diferencias finitas es ahora aplicable a la ecuacin (24), lo que significa que nuevamente reemplazaremos los trminos que involucran ahora derivadas parciales por sus correspondientes aproximaciones de diferencias finitas.

Aplicaremos a la resolucin de la siguiente barra prismtica, utilizando la grilla que vemos a continuacin:

Por condiciones de simetra, la solucin necesita ser obtenida slo para una cuarta parte de la seccin, como se muestra en la figura anterior. Utilizaremos una malla de tamao . Notamos que el valor de la funcin de tensin debe ser proporcional a la constante , y por simplicidad tomamos . Utilizando las aproximaciones por diferencia (17), tenemos que, para un punto como el : (25). El uso de las condiciones de simetra requiere que a lo largo del eje y que similarmente a lo largo del eje . Aplicando sta condicin por ejemplo en el punto como el , tenemos que la aproximacin de diferencia central de la primera derivada era

(despreciando el trmino de error). Entonces las condiciones son:

De igual manera se aplican en todos los puntos situados en el contorno de la regin , llegndose a condiciones similares en todos los casos. Planteando ecuaciones de tipo (25) para todos los puntos interiores de la regin, tenemos:

Las anteriores, han sido planteadas sistemticamente, siguiendo solamente la regla de la aproximacin de diferencia de la segunda derivada. Aplicando los criterios de simetra, y las condiciones de contorno sobre los lmites , las anteriores se reducen al siguiente conjunto:

Disponiendo las anteriores de forma matricial, con , tenemos:

(26). Nuevamente tenemos un sistema del tipo , el que puede resolverse para el vector incgnita por cualquier mtodo adecuado. Resolviendo, entonces se tiene que como solucin del sistema de ecuaciones (26). Para evaluar el momento torsor, se utiliza la regla trapezoidal en un dominio bidimensional, obtenindose , valor que puede compararse con la solucin exacta . Similarmente, la mxima pendiente est en el punto y una posible aproximacin al valor absoluto de la mxima tensin de corte es

, utilizando la aproximacin de diferencia hacia atrs para la derivada . Nuevamente podemos comparar con el valor exacto dado por (error: ). La aproximacin obtenida por el uso de la frmula de diferencia hacia atrs es de menor orden de exactitud que la aproximacin utilizada para la formulacin principal del problema. Podemos mejorar nuestra aproximacin de utilizando tres valores de la funcin de tensin sobre la seccin central como sigue. Denotando el punto , y , podemos escribir, segn la expansin de Taylor:

, dnde D est en la lnea AB.

, dnde E est en AC. De ambas es posible eliminar el trmino de , obteniendo entonces:

. Utilizando el primer trmino del lado derecho como aproximacin, el error cometido es entonces del mismo orden que el cometido en la aproximacin de la ecuacin que gobierna el problema. Insertando adecuadamente los valores, tenemos (error: ). Obviamente este resultado presenta mayor exactitud que el obtenido con la aproximacin de diferencia hacia atrs. APROXIMACION Y CONVERGENCIA Las soluciones uni y bidimensionales para ecuaciones diferenciales parciales ordinarias derivadas anteriormente por procedimientos numricos de diferencias finitas, ilustran las posibilidades de la discretizacin. El aparentemente inabordable ( a lo sumo matemticamente dificultoso) problema de resolucin de ecuaciones diferenciales en derivadas parciales ha sido reemplazado por un problema puramente algebraico en el cual debe resolverse un cierto nmero de ecuaciones simultneas. Para problemas pequeos es viable una resolucin manual, pero es conveniente desarrollar algoritmos computacionales para automatizar las operaciones de clculo. Lo importante es comprender que existe una posibilidad de solucin, aunque sta involucre una aproximacin. Hemos mostrado ya que el error en las aproximaciones de diferencias finitas decrece incrementando la densidad del mallado. Para aplicar el proceso a una situacin en la cual no disponemos de la solucin exacta, es necesario estudiar la convergencia del mtodo de acuerdo al refinamiento de la malla, en un intento por estimar la magnitud de los errores ocurridos al producirse una aproximacin. Si, por ejemplo, el error de una aproximacin es del orden de , entonces los resultados de dos soluciones sobre grillas de espaciado pueden extrapolarse como se detalla a continuacin. Digamos que y corresponden a las soluciones para las grillas anteriores 1 y 2 respectivamente y que corresponde a la solucin exacta en el punto que estamos considerando. De sta forma, an cuando no conocemos la magnitud del error, podemos escribir:

, de la cual podemos extraer la solucin exacta.

Esta relacin se conoce como extrapolacin de Richardson, y proporciona un mtodo para mejorar la solucin a partir de los resultados obtenidos para dos grillas de distinto tamao de espaciado. Es aplicable tambin a casos bidimensionales y tridimensionales. Bibliografa O.C. Zienkiewicz; El Mtodo de Elementos Finitos; Revert; 1982. O.C. Zienkiewicz & K. Morgan; Finite Elements and Approximation; John Wiley & Sons; 1983. K-J. Bathe; Finite Element Procedures; Prentice Hall; 1996. J.N. Reddy & M.L. Rasmussen; Anlisis Matemtico Avanzado; Limusa; 1992. S. Nakamura; Anlisis Numrico y Visualizacin Grfica con MatLab; Prentice Hall; 1997.