31/08/2011

Regra de 1/3 de Simpson Repetida

Seja o Polinômio de grau 2 que interpola f (x) nos pontos:

clip_image002

clip_image004

clip_image006

Das relações acima, podemos obter outras que serão úteis no desenvolvimento deste estudo:

clip_image008

clip_image010

clip_image012

clip_image014

clip_image016

clip_image018

Podemos utilizar o polinômio de Lagrange para estabelecer a fórmula de integração resultante da aproximação de f (x) por um polinômio P2(x):

clip_image020

Relembrando que o polinômio de Lagrange de grau 2 é dado por:

clip_image022

clip_image024

Assim, temos:

clip_image026

clip_image028

Substituindo as relações de (4) a (9) na integral acima, obtemos:

clip_image030

clip_image032

clip_image034

clip_image036

Para resolvermos as integrais, fazemos uma mudança de variáveis conveniente, de modo que:

clip_image038

clip_image040

clip_image042

Assim, para:

clip_image044

clip_image046

clip_image048

Desta forma, as integrais acima ficam:

clip_image050

clip_image052

clip_image054

clip_image056

clip_image058

clip_image060

clip_image062

clip_image064

clip_image066

A relação (19) é a fórmula de 1/3 de Simpson. Mas para determinadas funções, cometemos um erro ao aproximar f (x) por um polinômio P2(x):

clip_image068

A fim de contornar este problema, podemos aplicar a regra de Simpson repetida vezes no intervalo [a , b] = [x0 , xm]. Suponha que x0, x1, ... , xm são pontos igualmente espaçados, onde:

clip_image070

Uma condição necessária é que m seja para, pois cada parábola utiliza três pontos consecutivos e para cada par de intervalo, temos:

clip_image072

Então:

clip_image074

clip_image076

clip_image078

Assim, a fórmula de 1/3 de Simpson Repetida fica:

clip_image080

clip_image082

Exemplo: Calcular uma aproximação para a função abaixo utilizando a regra de 1/3 de Simpson repetida. Considere m = 10:

clip_image084

Desta forma, teremos:

clip_image086

Assim:

clip_image088

clip_image090

clip_image092

Se a integral definida na função, obteremos:

clip_image094

clip_image096

Comparando os resultados de (20) e (21), vemos que o erro aparece na sexta casa decimal. Se quisermos ter uma aproximação ainda melhor, basta aumentarmos o valor de m o quanto desejarmos.


Referências

[1] Cálculo Numérico – Aspectos Teóricos e Computacionais - Márcia Ruggiero
[2] Notas de Aula


Veja mais:

O Cálculo Integral
Regra dos Trapézios Repetida
Zeros Reais de Funções Reais - O Método de Newton – Raphson

COMO REFERENCIAR ESSE ARTIGO: Título: Regra de 1/3 de Simpson Repetida. Publicado por Kleber Kilhian em 31/08/2011. URL: . Leia os Termos de uso.


Siga também o blog pelo canal no Telegram.
Achou algum link quebrado? Por favor, entre em contato para reportar o erro.
Para escrever em $\LaTeX$ nos comentários, saiba mais em latex.obaricentrodamente.com.

7 comentários:

  1. Cara, acho que há um erro de digitação aí...

    e^1 é aprox. 2,7... não 1,7...

    ResponderExcluir
    Respostas
    1. Anônimo9/7/13 09:09

      cara, so que tem que subtrair por e^0 = 1

      Excluir
  2. Você diz no final do item (21)? Se for, neste caso não é $e^1$, mas sim:
    $e^1−e^0=2,7\cdots −1$
    e assim dá aproximadamente 1,7...

    ResponderExcluir
  3. Na parte do

    "Assim para:
    (16)
    (17)
    (18)"

    de onde você tira as conclusões?

    ResponderExcluir
  4. Veja que $x_0$, $x_1$ e $x_2$ são pontos definidos. Já $x$, é um valor qualquer. Sendo assim, tome as relações 13, 14, 15. Se $x = x_0$, pela relação 13, $z = 0$, já que $h$ não pode ser, pois $h = x1 - x0 = x2 - x2 = ...$
    Use um raciocínio análogo para 17 e 18.
    Abraços.

    ResponderExcluir
  5. Nessa parte==> Assim, a fórmula de 1/3 de Simpson Repetida fica:
    Aqui a fórmula, no caso faltou multiplicar o f(x1) por 4, ficaria assim f(x0) + 4*f(x1)....

    ResponderExcluir
    Respostas
    1. É verdade. Me parece que foi um erro na digitação. Agradeço pelo aviso. Em breve farei a correção. um abraço.

      Excluir

Whatsapp Button works on Mobile Device only

Pesquise no blog