integração numérica com possível divisão por 'zero'

9

Estou tentando integrar

01t2n+2exp(αr0t)dt

que é uma simples transformação de

1x2nexp(αr0x)dx

usando porque é difícil aproximar numericamente integrais impróprias. No entanto, isso leva ao problema de avaliar o novo integrando próximo de zero. Será muito fácil obter o número adequado de nós de quadratura, visto que o intervalo tem apenas o comprimento 1 (para que o comparável possa ser muito pequeno), mas que tipo de considerações devo fazer ao integrar perto de zero?t=1xdt

Em algum nível, acho que simplesmente tomar é uma boa ideia onde é um número pequeno . No entanto, qual número devo escolher? Deveria ser máquina epsilon? A divisão por máquina epsilon é um número bem quantificado? Além disso, se a divisão epsilon da minha máquina (ou próximo a ela) fornecer um número incrivelmente grande, a obtenção de se tornará ainda maior.ϵ1t2n+2exp(αr0t)dtϵexp(1ϵ)

Como devo explicar isso? Existe uma maneira de ter uma integral numérica bem definida dessa função? Caso contrário, qual é a melhor maneira de integrar a função?

drjrm3
fonte
11
Você já pensou em usar Monte Carlo?
Faheem Mitha
Eu sinto que isso não resolveria o problema. A integração de Monte Carlo costuma ser reservada para integrais de alta dimensão. Eu teria exatamente os mesmos problemas com Monte Carlo, simplesmente teria menos controle sobre onde minha função está sendo avaliada.
drjrm3
Você pode estar certo.
Faheem Mitha
Eu acho que ainda seria bom ter uma resposta (talvez para uma pergunta separada e mais geral) explicando como se faz a integração numérica quando a função é divergente em um limite, para o caso geral em que não é possível fazer a integral analiticamente. Então, novamente, que poderia muito bem ser encontradas no Numerical Recipes ...
David Z
@ Faaem: "Monte Carlo é um método extremamente ruim; deve ser usado apenas quando todos os métodos alternativos forem piores." - Alan Sokal
JM

Respostas:

10

1xeax=1axeax11a1eax=eaa+eaa2=a+1a2ea
1xkeax=1axkeax1ka1xk1eax=eaa+ka1xk1eax
I(k)=eaa+kaI(k1)
I(0)=eaa
Matt Knepley
fonte
absolutamente nenhuma idéia de como eu ignorei isso. obrigado.
drjrm3
11
Substituições inteligentes e integração por partes devem sempre ser uma das primeiras coisas que você faz com integrais indisciplinadas.
JM
1x2nexp(αx)dx assuming n::nonnegint,α>0Γ(2n+1,α)α2n1
Erik P.
Na verdade, o Mathematica escolhe representar a resposta como ExpIntegralE [-2 n, ar]. Se você executar o FunctionExpand, ele fornecerá a mesma resposta que o Maple.
Searke
1

Dê uma olhada no QUADPACK . Possui rotinas para integração em domínios (semi-) infinitos.

GertVdE
fonte