15.05.2013 Views

Cálculo Numérico - Engenharia Civil UEM

Cálculo Numérico - Engenharia Civil UEM

Cálculo Numérico - Engenharia Civil UEM

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

CAPÍTULO 11. INTEGRAÇÃO NUMÉRICA 322<br />

c) a função é dada simplesmente através de uma tabela-conjunto de pares ordenados obtidos como<br />

resultados de experiências. Aí não se conhece a expressão analítica da função em termos do argumentos<br />

x.<br />

As fórmulas de integração são de manejo fácil e prático e nos permite, quando a função f(x) é conhecida,<br />

ter uma idéia do erro cometido na integração numérica, como veremos mais adiante.<br />

Consideraremos integrais da forma:<br />

onde ω(x) ≥ 0 e contínua em [a, b].<br />

b<br />

a<br />

ω(x)f(x)dx ,<br />

A função ω(x) é chamada função peso e é igual a zero somente num número finito de pontos.<br />

Usaremos Fórmulas de Quadratura para aproximar a integral.<br />

Fórmulas de quadratura são aquelas que aproximam a integral usando combinação linear dos valores<br />

da função, isto é:<br />

Seja R(f) o erro, isto é:<br />

b<br />

a<br />

R(f) =<br />

ω(x)f(x)dx <br />

b<br />

a<br />

n<br />

Akf(xk) .<br />

k=0<br />

ω(x)f(x)dx −<br />

n<br />

Akf (xk) .<br />

Definição 11.1 - O grau de precisão de uma fórmula de quadratura é por definição o maior inteiro<br />

m tal que R x k = 0, k = 0, 1, . . . , m e R x m+1 = 0 .<br />

Observe que isto é equivalente a dizer que a fórmula de quadratura tem grau de precisão m se é exata<br />

para todo polinômio de grau ≤ m e é não exata para polinômios de grau m + 1.<br />

11.2 Fórmulas de quadratura interpolatória<br />

Sejam x0, x1, . . . , xn, n + 1 pontos distintos em [a, b] e sejam f0, f1, . . . , fn, n + 1 valores de uma<br />

função y = f(x) sobre x0, x1, . . . , xn.<br />

Seja Pn(x) o polinômio de interpolação da função y = f(x) sobre os n + 1 pontos. Pela fórmula de<br />

Lagrange, (Capítulo 10), temos que:<br />

Pn(x) =<br />

k=0<br />

n<br />

fkℓk(x) .<br />

k=0

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!