Étodos Uméricos: Integração Numérica

étodos
uméricos
INTEGRAÇÃO NUMÉRICA
Prof. Erivelton Geraldo Nepomuceno
PROGRAMA DE PÓS-GRADUAÇÃO EM ENGENHARIA ELÉTRICA
CENTRO FEDERAL DE EDUCAÇÃO

UNIVERSIDADE DE JOÃO DEL-REI TECNOLÓGICA DE MINAS GERAIS
PRÓ-REITORIA DE PESQUISA DIRETORIA DE PESQUISA E PÓS-GRADUAÇÃO
2016
Conteúdo
1. Fórmulas de Newton-Cotes.
2. Quadratura de Gauss-Legendre.
3. Comparação dos métodos de integ. simples.
4. Integ. dupla pelas fórmulas de Newton-Cotes.
5. Integ. dupla via fórmulas de Gauss-Legendre.
6. Comparação dos métodos para integ. Dupla.
Introdução
Seja uma função f(x) integrável no intervalo [a, b]:
 f ( x) dx  F (b)  F (a),
a
F ( x)  f ( x)
Quando a fórmula analítica for de difícil obtenção ou se forem

conhecidos somente valores discretos de f(x), se faz necessário o uso
de métodos numéricos para avaliar a integral de f(x).
Esses métodos consistem em aproximar f(x) por um polinômio

interpolador e determinar analiticamente a integral desse polinômio no
intervalo [a, b].
Formulas de Newton-Cotes
• Seja uma função f(x) aproximada por um PI, por exemplo, um
polinômio de Gregory-Newton:
Operador de Diferença Finita:
▪ Exemplo: Verificar a tabela de diferenças finitas:

• Para n = 1:
• Mudança de variável de x  ux e simplificando a notação ux  u:

• Usando a notação yi = f(xi):
• Integrando, analiticamente, o polinômio:
Que é conhecida como regra do trapézio.

• Aproximação de uma função f(x) por um PI P(x) de grau 1:
Exemplo: Calcular pela regra do trapézio a integral:
Polinômio de grau 1 passa pelos pontos a = x0 = 1 e b = x1 = 4

• Aproximando f(x) por polinômio P2(x) de grau 2:
• Mudança de variável:
• Equação de integração:
• Integrando analiticamente o polinômio:
Que é a regra de 1/3 de Simpson ou primeira regra de Simpson.

• Aproximação de uma função f(x) por um PI P(x) de grau 2:
Exemplo: Calcular usando a regra do 1/3 de Simpson a integral:
Para construir um polinômio de grau 2 são necessários 3 pontos. Assim

as abscissas por onde o polinômio irá passar são:
• Aproximando f(x) por polinômio P3(x) de grau 3:
• Mudança de variável:
• Equação de integração:
• Integrando analiticamente o polinômio:
Que é a regra dos 3/8 de Simpson ou segunda regra de Simpson.

• Aproximação de uma função f(x) por um PI P(x) de grau 3
Exemplo: Calcular usando a regra do 3/8 de Simpson a integral:
São necessários 4 pontos para construir um polinômio de grau 3. Assim

as abscissas são:
Sendo:
Resultado da integração melhora à medida que o grau do PI aumenta:
• Formulas de Newton-Cotes são da forma geral:
• Na prática é difícil o uso de um polinômio de grau superior a 3 para
integração numérica.
• O resultado é melhorado pela subdivisão do intervalo de integração
e aplicação de uma fórmula de Newton-Cotes em cada subintervalo.
1 – Regra do Trapézio (composta):
• Integração baseada em polinômio de grau 1:
• Subdividir [a, b] em m subintervalos iguais e aplicar a equação acima

a cada 2 pontos:
c0 = c m = 1
ci = 2, i = 1,2, ..., m-1
• Aplicável a qualquer número de subintervalos m.

• Integração da função f(x) utilizando 6 PI P(x) de grau 1:
Exemplo: Usando a regra do trapézio composta com m = 4
subintervalos, calcular:
Valor de h:
Dispositivo prático formado por quatro colunas, o com i = 0, 1, ..., m, xi =

a, a + h, a + 2h, ..., b, yi = f(xi) e ci sendo os coficientes de Cotes:
Exemplo: Usando a regra do trapézio composta com m = 5
subintervalos, calcular :
Valor de h:
2 – Regra de 1/3 de Simpson (composta):
• Subdividir [a, b] em m (múltiplo de 2) subintervalos iguais. Aplicar a

equação acima a cada 3 pontos:
Número de subintervalos m deve
ser múltiplo de 2, que é o grau do
PI usado.
c0 = c m = 1
ci = 4 se i impar
ci = 2 se i par
Exemplo: Usando a regra do 1/3 de Simpson composta com h = 0,25,
verificar que:
Valor de m:
Dispositivo prático:
Exemplo: Usando a regra do 1/3 de Simpson composta com m=6,
calcular:
Valor de h:
3 – Regra de 3/8 de Simpson (composta):
• Subdividir [a, b] em m (múltiplo de 3) subintervalos iguais. Aplicar a

equação acima a cada 4 pontos:
Número de subintervalos
m deve ser múltiplo de 3,
que é o grau do PI usado.
c0 = c m = 1
ci = 2 se i for múltiplo de 3
ci = 3 resto
Exemplo: Usando a regra dos 3/8 de Simpson com m = 6
Valor de h:
Exemplo: Usando usando a regra dos 3/8 de Simpson com passo de
integração h = 0,3, calcular:
Valor de m:
4 – Erro de Integração:
• Erro de truncamento de um polinômio de Gregory-Newton de grau n:
• Regra do trapézio (polinômio de grau n = 1)
• Erro de integração E1,1 cometido ao usar P1(x) no intervalo [x0, x1]

• Mudança de variável de x para
• Erro de integração global considerando os m subintervalos é:
• i é determinado em cada um dos m subintervalos.

• Se f”(x) for contínua no intervalo [a, b], então existe algum valor de
x =  [a, b] para o qual o somatório é igual a mf”().
• Considerando o passo de integração:
• Erro global de integração da regra do trapézio torna-se:
• Regra do 1/3 de Simpson:
• Regra do 3/8 de Simpson:
• Devido à dificuldade de determinar o valor de  ele é tomado como o

ponto no intervalo [a, b], no qual a derivada de f(ii,iv)(x) apresenta o
maior valor em módulo.
• As equações fornecem a cota máxima do erro de integração.
Exemplo: Utilizando a regra do 1/3 de Simpson com m = 2
Valor de h:
Erro de integração:
Resultado exato:
Exemplo: Usando as três primeiras fórmulas de Newton Cotes com m
= 6 subintervalos, calcular:
Regra do trapézio:
Regra do 1/3 de Simpson:
Regra dos 3/8 de Simpson :

Determinação de :
Erro de integração da regra do trapézio:
Erro de integração da regra dos 3/8 de Simpson :
Erro de integração da regra dos 1/3 de Simpson :

Erros de integração máximo e real cometidos:
Regra do 1/3 de Simpson produziu os menores erros máximo e erro

real.
Sinal negativo do erro En indica que a integração numérica foi por

excesso: In > Iexata.
Exemplo: Usando as três primeiras fórmulas de Newton-Cotes com E
< 10-2, calcular:
Determinação do valor m para cada fórmula:
Regra do trapézio:  = .
Regras de Simpson:  = /2.
Regra do trapézio:
Regra do 1/3 de Simpson:
Regra dos 3/8 de Simpson :

Passo de integração:
Pela regra do 1/3 de Simpson:
Verificação da exatidão:
Exemplo: Verificar o erro cometido no cálculo da integral a seguir
usando as seis primeiras fórmulas de Newton-Cotes, com m = 60:
A medida que o grau n do PI aumenta, o erro diminui.
Fórmulas utilizando grau par é melhor do que a de grau ímpar seguinte.

Quadratura de Gauss-Legendre
• Escolher pontos igualmente espaçados nas fórmulas de Newton-
Cotes simplifica os cálculos.
• Sem a imposição de espaçamento constante, podem ser obtidas

fórmulas que fornecem uma maior exatidão, usando o mesmo
número de pontos que Newton-Cotes.
1 – Fórmula para dois pontos:
• Com esse objetivo faz-se inicialmente a mudança de variável de x
para t, definida no intervalo [-1, 1]
• Derivando
• E definindo
• A integral
• Então deseja-se que,
• Expressão análoga à regra do trapézio
• Então deve-se encontrar valores de t1, t2, A1 e A2 que tornem a

exatidão a maior possível.
• Método construído de modo a ser exato para polinômios de grau até
3. Fazendo
• e impondo
• Para
• Para
• Para
• Para
• Sistema de equações não lineares de ordem 4
• Solução
Exemplo: Calcular:
Mudança de variável.
Dispositivo prático:
Resultado exato:
Exemplo: Calcular:
Valor exato:
Erro cometido:
Mais exato que pela regra do trapézio com m = 6 subintervalos,

equivalente a 7 pontos (30,8816).
2 – Fórmula geral:
• Determinar os valores dos pesos Ai, e das abscissas
• De modo que esta seja exata para polinômios de grau menor ou igual
a 2n - 1.
• Sabendo que:
• Impondo:
• Sistema de equações não lineares de ordem 2n

• Solução fornece os n pesos Ai e as n abscissas ti.
• A solução desse sistema linear pode ser evitada utilizando um
processo alternativo. Inicialmente sejam os polinômios de Legendre
definidos pela fórmula de recorrência:
• Com L0(x) = 1 e L1(x) = x.

• Por exemplo
• Para maior exatidão na fórmula de quadratura é suficiente que ti, i =
1, 2, ... . n sejam os zeros do polinômio de Legendre de grau n.
• Conhecidas as abscissas ti, o sistema não linear se reduz a um
sistema linear de ordem n
• Em vez de resolver este sistema via decomposição LU, os pesos Ai

podem ser obtidos por
• L’n(ti): derivada de Ln(x) na abscissa ti.
• L’n(ti): derivada de Ln(x) na abscissa ti.
Exemplo: Verificar com n=3 e n=4.que:
Para n = 3:
Para n = 4:
Exemplo: Com n=5 calcular:
Resultado exato:
Gauss-Legendre com n = 5 é mais exato que a regra do 1/3 de Simpson

com m = 6 (30,4337).
3 – Erro de integração:
• Erro de integração da fórmula de Gauss-Legendre
• : abscissa, na qual a derivada f2n(x) apresenta o maior valor em

módulo no intervalo [a, b].
• Cota máxima do erro de integração da fórmula de Gauss-Legendre.

Exemplo: Com n = 2 e o respectivo erro de integração calcular:
Para n = 2:
Determinação de .
Cálculo do erro máximo:
Valor exato:
Erro real:
Comparação dos Métodos
As comparações são realizadas por meio dos exemplos a seguir:
Foram utilizadas regras simples de Newton Cotes, isto é o grau do

polinômio é igual ao número de subintervalos. O número de ponto de
Gauss Legendre é igual a m+1, sendo m o número de intervalos de
Newton-Cotes, de forma a ter o mesmo número de pontos avaliados.
Quadratura de Gauss-Legendre mais vantajosa !!!

Referencias Bibliográficas
1. Frederico Ferreira Campos Filho, Algoritmos Numéricos.

Étodos Uméricos: Integração Numérica

Enviado por

Direitos autorais:

Formatos disponíveis

Étodos Uméricos: Integração Numérica

Enviado por

Dados do documento

Título original

Direitos autorais

Formatos disponíveis

Compartilhar este documento

Compartilhar ou incorporar documento

Opções de compartilhamento

Você considera este documento útil?

Este conteúdo é inapropriado?

Direitos autorais:

Formatos disponíveis

Étodos Uméricos: Integração Numérica

Enviado por

Direitos autorais:

Formatos disponíveis

étodos

PROGRAMA DE PÓS-GRADUAÇÃO EM ENGENHARIA ELÉTRICA

CENTRO FEDERAL DE EDUCAÇÃO

Quando a fórmula analítica for de difícil obtenção ou se forem

Esses métodos consistem em aproximar f(x) por um polinômio

▪ Exemplo: Verificar a tabela de diferenças finitas:

• Mudança de variável de x  ux e simplificando a notação ux  u:

• Integrando, analiticamente, o polinômio:

Que é conhecida como regra do trapézio.

Polinômio de grau 1 passa pelos pontos a = x0 = 1 e b = x1 = 4

Que é a regra de 1/3 de Simpson ou primeira regra de Simpson.

Para construir um polinômio de grau 2 são necessários 3 pontos. Assim

Que é a regra dos 3/8 de Simpson ou segunda regra de Simpson.

São necessários 4 pontos para construir um polinômio de grau 3. Assim

• Subdividir [a, b] em m subintervalos iguais e aplicar a equação acima

• Aplicável a qualquer número de subintervalos m.

Dispositivo prático formado por quatro colunas, o com i = 0, 1, ..., m, xi =

• Subdividir [a, b] em m (múltiplo de 2) subintervalos iguais. Aplicar a

• Subdividir [a, b] em m (múltiplo de 3) subintervalos iguais. Aplicar a

• Regra do trapézio (polinômio de grau n = 1)

• Erro de integração E1,1 cometido ao usar P1(x) no intervalo [x0, x1]

• Erro de integração global considerando os m subintervalos é:

• i é determinado em cada um dos m subintervalos.

• Regra do 1/3 de Simpson:

• Regra do 3/8 de Simpson:

• Devido à dificuldade de determinar o valor de  ele é tomado como o

Regra do 1/3 de Simpson:

Regra dos 3/8 de Simpson :

Erro de integração da regra do trapézio:

Erro de integração da regra dos 3/8 de Simpson :

Erro de integração da regra dos 1/3 de Simpson :

Erros de integração máximo e real cometidos:

Regra do 1/3 de Simpson produziu os menores erros máximo e erro

Sinal negativo do erro En indica que a integração numérica foi por

Determinação do valor m para cada fórmula:

Regra do 1/3 de Simpson:

Regra dos 3/8 de Simpson :

A medida que o grau n do PI aumenta, o erro diminui.

Fórmulas utilizando grau par é melhor do que a de grau ímpar seguinte.

• Sem a imposição de espaçamento constante, podem ser obtidas

• Expressão análoga à regra do trapézio

• Então deve-se encontrar valores de t1, t2, A1 e A2 que tornem a

Mais exato que pela regra do trapézio com m = 6 subintervalos,

• Determinar os valores dos pesos Ai, e das abscissas

• Sistema de equações não lineares de ordem 2n

• Com L0(x) = 1 e L1(x) = x.

• Em vez de resolver este sistema via decomposição LU, os pesos Ai

Gauss-Legendre com n = 5 é mais exato que a regra do 1/3 de Simpson

• : abscissa, na qual a derivada f2n(x) apresenta o maior valor em

• Cota máxima do erro de integração da fórmula de Gauss-Legendre.

Cálculo do erro máximo:

Foram utilizadas regras simples de Newton Cotes, isto é o grau do

Quadratura de Gauss-Legendre mais vantajosa !!!

Você também pode gostar