Saya memiliki fungsi sehingga
adalah terbatas, dan saya ingin mendekati integral ini.
Saya akrab dengan aturan kuadratur dan pendekatan integral monte carlo, tapi saya melihat beberapa kesulitan menerapkannya dalam domain tak terbatas. Dalam kasus monte carlo, bagaimana orang melakukan pengambilan sampel wilayah yang tidak terbatas (terutama jika wilayah yang berkontribusi lebih signifikan terhadap integral tidak diketahui)? Dalam kasus quadrature, bagaimana cara menemukan titik optimal? Haruskah saya memperbaiki wilayah sewenang-wenang besar yang berpusat di sekitar titik asal dan menerapkan aturan kuadratur yang jarang? Bagaimana saya bisa mendekati perkiraan integral ini?
sumber
Cara standar untuk melakukannya adalah mengekstraksi dari ekspresi untuk sebuah prefactor eksponensial, mentransformasikannya menjadi e - x 2 , dan kemudian menggunakan aturan kuadratur Gaussian (atau Gauss Kronrod) dengan ini sebagai bobot. Jika f mulus, ini biasanya memberikan hasil yang sangat baik.f(x) e−x2 f
Di , yang sama bekerja dengan berat e - | x | 2 , dan formula cubature yang sesuai dapat ditemukan, misalnya, dalam buku karya Engels, quadrature numerik dan cubature.R3 e−|x|2
Formula online ada di http://nines.cs.kuleuven.be/ecf/
sumber
Untuk quadrature satu dimensi, Anda dapat memeriksa buku tentang Quadpack (oldie emas tetapi masih sangat relevan dalam quadrature satu dimensi) dan teknik yang digunakan dalam algoritma QAGI, integrator otomatis untuk rentang tak hingga.
Teknik lain adalah rumus quadrature eksponensial ganda, diterapkan dengan baik untuk interval tak terbatas oleh Ooura .
Untuk cubature, Anda dapat berkonsultasi dengan Encyclopedia of formula cubature oleh Ronald Cools.
sumber
Jika Anda ingat bagaimana quadrature bekerja, maka Anda juga akan tahu satu cara bagaimana memperkirakan integral tak terbatas. Yaitu: untuk quadrature, Anda memperkirakan fungsinyaf( x ) Anda ingin mengintegrasikan dengan sesuatu yang serupa, katakanlah jumlahnya banyak f~( x ) (atau polinomial piecewise) yang Anda dapat menuliskan integralnya secara analitis. Anda mendapatkanf~ dari f dengan interpolasi pada titik interpolasi - yang kemudian akan menjadi titik quadrature Anda.
Untuk integral tak terbatas, satu pendekatan menggunakan pemikiran yang sama persis. Misalnya, Anda dapat mencoba perkiraanf( x ) pada seluruh baris menggunakan fungsi yang lebih sederhana, mis f~( x ) = e- x2p ( x ) dengan polinomial p ( x ) interpolasi itu f( x ) ex2 di sejumlah titik. Kemudian ada rumus sederhana untuk menghitung integral∫∞- ∞f~( x ) dx . Pilihan titik interpolasi mengikuti logika yang sama seperti yang dilakukan untuk derivasi quadrature biasa.
sumber
Jika Anda ingin menggunakan integrasi Monte Carlo, Anda bisa mulai dengan menggunakan sampling penting dengan sampler yang secara kasar mendekati integrand Anda. Semakin baik sampler Anda cocok dengan integrand Anda, semakin sedikit varians dalam perkiraan integral Anda. Tidak masalah jika domain Anda tidak terbatas selama sampler Anda memiliki domain yang sama.
sumber