Metodos Numericos

4- DIFERENCIACIÓN E INTEGRACIÓN NUMÉRICA
Una integral definida tiene la forma,
Donde a, b son el límite inferior y superior de integración, respectivamente, f(x) es la

función a integrar y dx es la diferencial de x.
La integral, geométricamente, representa el área delimitada por el lugar geométrico de la

función, f(x), el eje de la abscisas, y la dos rectas verticales x = a y x = b, ver la figura,
Cuando, la integral definida es muy difícil o de plano imposible de resolver, existen

estrategias que permiten tener una solución aproximada. La integral definida permite que,
geométricamente, se determine una aproximación a través de trapecios. Otra forma de
obtener una solución aproximada a la integral definida es tratar de ajustar un polinomio al
lugar geométrico de la función a integrar, y así en vez de integrar la función f(x), se integra
el polinomio Pn(x),
Al igual que la integral, la derivada tiene un significado geométrico que permite una
aproximación numérica. Así, la derivada representa, geométricamente, la pendiente de una
recta tangente,
La recta es tangente al lugar geométrico de la función f(x). Toman dos puntos alrededor del
punto de tangencia es posible representar la pendiente a través de,
La idea anterior puede retomarse como base para determinar una aproximación a la
derivada de una función, junto con estrategias que mejoren dicha aproximación.
4.1. Ecuaciones de diferencias divididas finitas para datos

uniformemente distribuidos.
Este método consiste en una aproximación de las derivadas parciales por expresiones
algebráicas con los valores de la variable dependiente en un limitado número de
puntos seleccionados. Como resultado de la aproximación, la ecuación diferencial
parcial que describe el problema es reemplazada por un número finito de ecuaciones
algebraicas, en términos de los valores de la variable dependiente en puntos
seleccionados.
El valor de los puntos seleccionados se convierten en las incógnitas. El sistema de

ecuaciones algebraicas debe ser resuleto y puede llevar un número largo de
operaciones aritméticas.
Método de Diferencias Divididas
El método de NEWTON de diferencias divididas es otra forma de obtener el polinomio

interpolador.
Dado cierto número de puntos obtenidos pormuestreo o a partir de un experimento se

prentede 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 x1, x1,
....., xm, se llama interpolación polinómica al proceso de hallar un polinomio pm(x) de grado
menor o igual a m.
Sea fn una variable discreta de n elementos y se xn otra variable discreta de n elementos los
cuales corresponden por parejas a la imagen u ordenada y abcisa de los datos que se
quieran interpolar, respectivament, tales que:
Este método es muy algorítmico, por lo que resulta muy cómodo, más cuando se quiere
calcular un polinomio interpolador de grado elevad. El polonimio resultante tendrá la
forma.
Los coeficientes aj son las llamadas diferencias dividias.

4.2. Ecuaciones para derivar datos irregularmente
espaciados.
4.3. Ecuación de integración de Newton-Cotes.
Las fórmulas de Newton - Cotes son los tipos de integración numérica más
comunes. Se basan en la estrategia de reemplazar una función complicada o
datos tabulados por un polinomio de aproximación que es fácil de integrar:
donde fn(x) es un polinomio de la forma
donde n es el grado del polinomio. Por ejemplo, en la Figura 1 se utiliza un

polinomio de primer grado como una aproximación, mientras que en la Figura
2, se emplea una parábola con el mismo propósito.
La integral también se puede aproximar usando un conjunto de polinomios

aplicados por pedazos a la función o datos, sobre segmentos de longitud
constante. Así, en la Figura 3, se usan tres segmentos de línea recta para
aproximar la integral.
Existen formas cerradas y abiertas de las fórmulas de Newton - Cotes. En
esta sección sólo se analizarán las formas cerradas. En ellas, se conocen los
datos al inicio y al final de los límites de integración.
4.4. Aplicaciones de la diferenciación e integración

numérica.
Grado de precisión de una fórmula de integración numérica
El grado de precisión de una fórmula de integración numérica es el

número natural n que verifica que el error de truncamiento E[P i]=0 para
todos los polinomios Pi(x) de grado i ≤ n, y existe un polinomio Pn+1(x)
de grado n+1 tal que E[Pn+1]≠0.
A continuación, se explicarán las fórmulas que se obtienen cuando
el grado del polinomio de aproximación es:
 uno (Regla del trapecio)
 dos (Regla de Simpson)
 tres (Regla 3/8 de Simpson
Regla del trapecio
La regla del trapecio es la primera de las fórmulas cerradas de Newton - Cotes.

Corresponde al caso donde el polinomio de aproximación es de primer grado:
Teniendo en cuenta que la ecuación de la recta que pasa por los puntos (a;f(a)) y
(b;f(b)) es:
el área bajo esta línea recta en el intervalo [a;b] está dada por:
Esta integral constituye una aproximación de la integral de f(x) en dicho intervalo. El resultado
de la integral anterior es:
que se denomina regla del trapecio.

Geométricamente, la regla del trapecio consiste en aproximar el área debajo de la curva
definida por f(x), por el área bajo la recta que une los puntos (a, f(a)) y (b, f(b)). Recuerde que la
fórmula para calcular el área de un trapecio es la altura por el promedio de las bases.
Error de la regla del trapecio
Cuando se emplea la integral bajo un segmento de línea recta para aproximar la

integral bajo una curva, obviamente se tiene un error que puede ser importante. Una
estimación del error de truncamiento E de la regla del trapecio es:
donde ξ Є (a;b). La expresión anterior indica que si la función a integrar es lineal, la regla del
trapecio será exacta. Es decir, la regla del trapecio tiene grado de precisión n = 1. Además,
sólo es posible aplicar la regla del trapecio si f(x) es de clase C 2[a;b].
La regla compuesta del trapecio
Una forma de mejorar la precisión de la regla del trapecio consiste en dividir el

intervalo de integración [a;b] en varios segmentos, y aplicar el método a cada uno de
ellos. Las áreas de los segmentos se suman después para obtener una aproximación
de la integral en todo el intervalo. Las expresiones resultantes se llaman fórmulas de
integración, de aplicación múltiple o compuestas.
La siguiente figura muestra el formato general y la nomenclatura que se usará
para obtener integrales de aplicación múltiple.
Si hay n+1 puntos igualmente espaciados, existen n segmentos del mismo
ancho:
Si a y b se designan como x0 y xn, respectivamente, la integral completa se

representará como:
Sustituyendo la regla del trapecio en cada integral se obtiene:
o, agrupando términos:
El error que se comete al aplicar la regla compuesta del trapecio está dado por:
Esto significa que el error es de orden O(h2). Además, cuando las derivadas de f(x) se
conocen, es posible estimar el número de subintervalos necesarios para alcanzar la
precisión deseada.
Regla de Simpson
Además de aplicar la regla del trapecio con una

segmentación más fina, otra forma de obtener una estimación
más precisa de una integral consiste en usar polinomios de grado
superior para unir los puntos. Por ejemplo, si se toma el punto
medio del intervalo de integración [a;b], los tres puntos se
pueden unir con una parábola. La fórmula que resulta de tomar
la integral bajo ese polinomio se conoce como regla de Simpson.
Es decir, la regla de Simpson se obtiene cuando el polinomio de
aproximación es de segundo grado.
Si a = x0, b = x2, se llama x1 al punto medio del
intervalo [a;b] y f(x) se representa por un polinomio de
Lagrange de segundo grado, la integral se transforma en:
Después de integrar y trabajar algebraicamente se obtiene:
donde h=(b-a)/2. Esta expresión se conoce como regla de Simpson y es la segunda

fórmula de integración cerrada de Newton - Cotes.
Error de la regla de Simpson
Se puede demostrar que la aplicación de un solo segmento de la regla de

Simpson tiene un error de truncamiento dado por la fórmula:
donde ξ Є (a;b). Así, la regla se Simpson es más precisa que la regla del trapecio.
Además, en lugar de ser proporcional a la tercera derivada, el error es proporcional a
la cuarta derivada. Esto es porque el término del coeficiente de tercer grado se hace
cero durante la integración de la interpolación polinomial. En consecuencia, la regla de
Simpson alcanza una precisión n = 3 aún cuando se base en sólo tres puntos. En otras
palabras, la regla de Simpson da resultados exactos para polinomios cúbicos aún
cuando se obtenga de una parábola.
Cabe destacar que es posible aplicar esta regla si f(x) es de clase C4[a;b].
La regla compuesta de Simpson
Así como en la regla del trapecio, la regla de Simpson se mejora al dividir el

intervalo de integración en varios segmentos de un mismo tamaño:
La integral total se puede representar como:

Al sustituir la regla de Simpson en cada integral se obtiene:
Una forma más compacta de escribir la expresión anterior es:
El error que se comete al aplicar la regla compuesta de Simpson está dado por:
Esto significa que el error es de orden O(h 4). También, cuando se conocen las
derivadas de f(x), es posible estimar el número de subintervalos necesarios para
alcanzar la precisión deseada.
La regla 3/8 de Simpson
De manera similar a la obtención de la regla del trapecio y

Simpson, es posible ajustar un polinomio de Lagrange de tercer
grado a cuatro puntos e integrarlo:
para obtener:
donde h=(b-a)/3. Esta expresión se llama regla 3/8 de

Simpsondebido a que h se multiplica por 3/8. Ésta es la tercera
fórmula de integración cerrada de Newton - Cotes.
Error de la regla 3/8 de Simpson
Se puede demostrar que la aplicación de un solo segmento de la regla 3/8 de

Simpson tiene un error de truncamiento E:
donde ξ Є (a;b). Como el denominador de la expresión del error de truncamiento, en

este caso, es menor que el de la regla de Simpson, se puede concluir que la regla 3/8
de Simpson es menos precisa que la regla de Simpson.
En general, se prefiere la regla de Simpson, ya que alcanza una precisión n = 3
con tres puntos en lugar de los cuatro puntos requeridos en la versión 3/8. No
obstante, la regla 3/8 es útil cuando el número de segmentos es impar.
Es posible aplicar la regla 3/8 de Simpson si f(x) es de clase C 4[a;b].
Tecnológico nacional de mexico
Instituto tecnológico de cerro azu
}

Metodos Numericos

Cargado por

Copyright:

Formatos disponibles

Metodos Numericos

Cargado por

Información del documento

Descripción original:

Título original

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

Metodos Numericos

Cargado por

Copyright:

Formatos disponibles

4- DIFERENCIACIÓN E INTEGRACIÓN NUMÉRICA

Una integral definida tiene la forma,

Donde a, b son el límite inferior y superior de integración, respectivamente, f(x) es la

La integral, geométricamente, representa el área delimitada por el lugar geométrico de la

Cuando, la integral definida es muy difícil o de plano imposible de resolver, existen

4.1. Ecuaciones de diferencias divididas finitas para datos

El valor de los puntos seleccionados se convierten en las incógnitas. El sistema de

El método de NEWTON de diferencias divididas es otra forma de obtener el polinomio

Dado cierto número de puntos obtenidos pormuestreo o a partir de un experimento se

Los coeficientes aj son las llamadas diferencias dividias.

donde fn(x) es un polinomio de la forma

donde n es el grado del polinomio. Por ejemplo, en la Figura 1 se utiliza un

La integral también se puede aproximar usando un conjunto de polinomios

4.4. Aplicaciones de la diferenciación e integración

Grado de precisión de una fórmula de integración numérica

El grado de precisión de una fórmula de integración numérica es el

La regla del trapecio es la primera de las fórmulas cerradas de Newton - Cotes.

que se denomina regla del trapecio.

Cuando se emplea la integral bajo un segmento de línea recta para aproximar la

La regla compuesta del trapecio

Una forma de mejorar la precisión de la regla del trapecio consiste en dividir el

Si a y b se designan como x0 y xn, respectivamente, la integral completa se

Sustituyendo la regla del trapecio en cada integral se obtiene:

Además de aplicar la regla del trapecio con una

Después de integrar y trabajar algebraicamente se obtiene:

donde h=(b-a)/2. Esta expresión se conoce como regla de Simpson y es la segunda

Error de la regla de Simpson

Se puede demostrar que la aplicación de un solo segmento de la regla de

Así como en la regla del trapecio, la regla de Simpson se mejora al dividir el

La integral total se puede representar como:

Una forma más compacta de escribir la expresión anterior es:

La regla 3/8 de Simpson

De manera similar a la obtención de la regla del trapecio y

donde h=(b-a)/3. Esta expresión se llama regla 3/8 de

Se puede demostrar que la aplicación de un solo segmento de la regla 3/8 de

donde ξ Є (a;b). Como el denominador de la expresión del error de truncamiento, en

También podría gustarte