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

Deber 5

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

ESCUELA POLITECNICA NACIONAL

FACULTAD DE INGENIERIA MECANICA


ELEMENTOS FINITOS

Nombre: Lizeth Carolina Barragán Machado


Paralelo: GR3
Fecha: 22/06/2017

DEBER N*5

POLINOMIOS DE INTERPOLACIÓN: DIFERENCIAS DIVIDIDAS DE NEWTON Y


DE LAGRANGE

En análisis numérico, la interpolación polinómica es una técnica de interpolación


de un conjunto de datos o de una función por un polinomio. Es decir, dado cierto
número de puntos obtenidos por muestreo o a partir de un experimento se
pretende encontrar un polinomio que pase por todos los puntos.

Dada una función f de la cual se conocen sus valores en un número finito de


abscisas 𝑋0, 𝑋1, . . . , 𝑋𝑚, se llama interpolación polinómica al proceso de hallar un
polinomio Pm(x) de grado menor o igual a m, cumpliendo 𝑃𝑚(𝑥𝑘 ) = 𝑓(𝑥𝑘 ), ∀ 𝑘 =
0, 1, . . . , 𝑚.

A este polinomio se le llama Polinomio interpolador de grado m de la función 𝑓.

La interpolación polinómica es un método usado para conocer, de un modo


aproximado, los valores que toma cierta función de la cual sólo se conoce su
imagen en un número finito de abscisas. A menudo, ni siquiera se conocerá la
expresión de la función y sólo se dispondrá de los valores que toma para dichas
abscisas.

El objetivo será hallar un polinomio que cumpla lo antes mencionado y que


permita hallar aproximaciones de otros valores desconocidos para la función con
una precisión deseable fijada. Por ello, para cada polinomio interpolador se
dispondrá de una fórmula del error de interpolación que permitirá ajustar la
precisión del polinomio.

CALCULO DEL POLINOMIO INTERPOLADOR

Se dispone de varios métodos generales de interpolación polinómica que


permiten aproximar una función por un polinomio de grado m. El primero de
estos polinomios es el método de las diferencias divididas de Newton. Otro de
los métodos es la interpolación de Lagrange:

METODO DE LAS DIFERENCIAS DIVIDIDAS DE NEWTON

Sea 𝑓𝑛 una variable discreta de n, elementos y sea Xn otra variable discreta de


n, elementos los cuales corresponden, por parejas, a la imagen u ordenada y
abscisa de los datos que se quieran interpolar, respectivamente, tales que:
Este método es muy algorítmico y resulta sumamente cómodo en determinados
casos, sobre todo cuando se quiere calcular un polinomio interpolador de grado
elevado. El polinomio de grado n-1 resultante tendrá la forma:

definiendo como

y definiendo como:

Los coeficientes son las llamadas diferencias divididas.


Una vez se hayan realizado todos los cálculos, nótese que hay (muchas) más
diferencias divididas que coeficientes . El cálculo de todos los términos
intermedios debe realizarse simplemente porque son necesarios para poder
formar todos los términos finales. Sin embargo, los términos usados en la
construcción del polinomio interpolador son todos aquellos que involucren a .
Estos coeficientes se calculan mediante los datos que se conocen de la
función : [1]

INTERPOLACIÓN DE LAGRANGE

Sea la función a interpolar, sean las abscisas conocidas de y


sean los valores que toma la función en esas abscisas, el
polinomio interpolador de grado de Lagrange es un polinomio de la forma
donde son los llamados polinomios de Lagrange, que se calculan de este
modo:

Nótese que en estas condiciones, los coeficientes están bien definidos y


son siempre distintos de cero. [2]

CUADRATURA DE GAUSS- LEGENDRE

Las cuadraturas de Gauss - Legendre es el nombre de una clase de técnicas que


aplica tal estrategia para obtener una aproximación más precisa de la integral.

El objetivo de la cuadratura de Gauss - Legendre es determinar las abscisas


𝑋1 𝑦 𝑋2 y dos coeficientes w1 y w2 de manera que la fórmula:

sea exacta para polinomios cúbicos de la forma f(x) = a3x3 + a2x2 + a1x +a0.
Como hay que determinar cuatro números w1, w2, x1 y x2 en la expresión
anterior, se deben seleccionar cuatro condiciones que deben cumplirse. Usando
el hecho de que la integración es aditiva, será suficiente con exigir que la integral
anterior sea exacta para las cuatro funciones f(x) = 1, x, x2, x3. Por lo tanto, las
cuatro condiciones de integración son:

De esta manera, el sistema de ecuaciones no lineales que se debe resolver es:


La solución del sistema anterior está dada por:

Así, se ha encontrado los nodos y los coeficientes o pesos con los que se
construye la cuadratura de Gauss - Legendre. En consecuencia, si f es continua
en [-1;1], resulta:

La cuadratura de Gauss - Legendre con dos nodos G2(f) tiene grado de precisión
n=3 y si f Є C4[-1;1], entonces:

Cuadratura de Gauss - Legendre con más puntos

Además de la fórmula de dos puntos descripta en la sección anterior, es posible


desarrollar versiones con más puntos en la forma general:

donde n es la cantidad de nodos que se toma.

Los nodos xn,k y los pesos wn,k que hay que usar, en cada caso, se encuentran
tabulados. En la siguiente tabla, se muestran los valores correspondientes para
las cuadraturas de Gauss - Legendre con hasta cinco nodos, así como la forma
de los términos de error En(f) correspondientes a las aproximaciones:
Los nodos que se usan para construir una cuadratura de n puntos son las raíces
del polinomio de Legendre de grado n y los pesos, como ya se explicó
anteriormente, se obtienen resolviendo un sistema de ecuaciones no lineales.
Como se puede observar en la tabla, el término de error de una cuadratura de
tres puntos es proporcional a la derivada de orden seis de la función f(x). Esto
implica que ese tipo de cuadratura será exacta si el integrando es un polinomio
de grado cinco o menor. De mismo modo, una cuadratura de cuatro puntos será
exacta si la función a integrar es un polinomio de grado siete o menor. Esto se
debe a que su término de error involucra una derivada de orden ocho de la
función f(x). En general, una cuadratura de Gauss - Legendre de n puntos será
exacta para funciones polinomiales de grado menor o igual que 2n - 1.

Traslación de la Cuadratura de Gauss - Legendre

Para aplicar la cuadratura de Gauss - Legendre en un intervalo [a;b], se debe


efectuar el cambio de variable:
De esta manera,

Por lo tanto, la fórmula de cuadratura está dada por:

[3]

BIBLIOGRAFIA:

[1] www.upv.es/mattel/asig/labormat/ejercici/practi07.doc

[2] https://es.wikipedia.org/wiki/Interpolación_polinómica

[3] http://www.frsn.utn.edu.ar/GIE/AN/IN/Cuadraturas_Gauss_Legendre.html

También podría gustarte