integrasi numerik dengan kemungkinan pembagian dengan 'nol'

9

Saya mencoba mengintegrasikan

01t2n+2exp(αr0t)dt

yang merupakan transformasi sederhana dari

1x2nexp(αr0x)dx

menggunakan karena sulit untuk memperkirakan secara integral integral secara numerik. Ini, bagaimanapun, mengarah pada masalah mengevaluasi integrand baru mendekati nol. Akan sangat mudah untuk mendapatkan jumlah node quadrature yang tepat karena intervalnya hanya panjang 1 (sehingga sebanding dapat dibuat sangat kecil), tetapi pertimbangan seperti apa yang harus saya buat ketika mengintegrasikan mendekati nol? dtt=1xdt

Pada tingkat tertentu, saya berpikir bahwa dengan mengambil adalah ide yang baik di mana adalah sejumlah kecil . Namun, nomor apa yang harus saya pilih? Haruskah itu mesin epsilon? Apakah pembagian berdasarkan mesin epsilon merupakan angka yang terkuantifikasi dengan baik? Selain itu, jika pembagian epsilon mesin saya (atau dekat dengan itu) memberikan angka yang sangat besar, maka mengambil akan menjadi lebih besar.ϵexp(1ϵ1t2n+2exp(αr0t)dtϵexp(1ϵ)

Bagaimana saya harus menjelaskan ini? Apakah ada cara untuk memiliki integral numerik yang terdefinisi dengan baik dari fungsi ini? Jika tidak, apa cara terbaik untuk mengintegrasikan fungsi?

drjrm3
sumber
1
Sudahkah Anda mempertimbangkan untuk menggunakan Monte Carlo?
Faheem Mitha
Saya merasa itu tidak akan memperbaiki masalah. Integrasi Monte Carlo sering dicadangkan untuk integral dimensi tinggi. Saya akan mengalami masalah yang sama persis dengan Monte Carlo, saya hanya akan memiliki sedikit kontrol di mana fungsi saya sedang dievaluasi.
drjrm3
Kamu mungkin benar.
Faheem Mitha
Saya pikir itu masih akan baik untuk memiliki jawaban (mungkin untuk pertanyaan yang terpisah, lebih umum) menjelaskan bagaimana seseorang melakukan integrasi numerik ketika fungsinya berbeda pada satu batas, untuk kasus umum di mana tidak mungkin untuk melakukan integral analitik. Kemudian lagi, itu bisa juga ditemukan dalam Numerical Recipes ...
David Z
@Faheem: "Monte Carlo adalah metode yang sangat buruk; harus digunakan hanya ketika semua metode alternatif lebih buruk." - Alan Sokal
JM

Jawaban:

10

Ini dapat dilakukan dengan integrasi oleh bagian-bagian: dan melanjutkan dengan induksi sehingga dan .

1xe-Sebuahx=-1Sebuahxe-Sebuahx1--1Sebuah1e-Sebuahx=e-SebuahSebuah+e-SebuahSebuah2=Sebuah+1Sebuah2e-Sebuah
1xke-Sebuahx=-1Sebuahxke-Sebuahx1--kSebuah1xk-1e-Sebuahx=e-SebuahSebuah+kSebuah1xk-1e-Sebuahx
saya(k)=e-SebuahSebuah+kSebuahsaya(k-1)
saya(0)=e-SebuahSebuah
Matt Knepley
sumber
sama sekali tidak tahu bagaimana saya mengabaikan ini. Terima kasih.
drjrm3
1
Penggantian dan integrasi yang cerdas oleh bagian harus selalu menjadi salah satu hal pertama yang Anda lakukan dengan integral yang sulit dikendalikan.
JM
Sering kali ide yang baik untuk menanyakan sistem aljabar komputer setiap kali Anda memiliki integral seperti ini. Maple mengevaluasi " " segera ke ; Saya yakin Mathematica melakukan hal yang sama. (Masih ide yang bagus untuk memverifikasinya secara numerik, tentu saja, yang biasanya dapat dilakukan oleh orang-orang ini.)1x2nexp(-αx)dx asumsi n::nonnegint,α>0Γ(2n+1,α)α-2n-1
Erik P.
Sebenarnya Mathematica memilih untuk mewakili jawaban sebagai ExpIntegralE [-2 n, ar]. Jika Anda menjalankan FunctionExpand di atasnya, maka itu memberikan jawaban yang sama seperti Maple.
Searke
1

Lihatlah QUADPACK . Ini memiliki rutinitas untuk integrasi lebih dari (semi-) domain tak terbatas.

GertVdE
sumber