Dalam kursus lanjutan ini pada aplikasi teori fungsi kompleks pada satu titik dalam latihan integral yang sangat berosilasi
harus didekati untuk nilai menggunakan metode sadel di bidang kompleks.
Karena sifatnya yang sangat berosilasi, integral ini sangat sulit untuk dievaluasi menggunakan sebagian besar metode lain. Ini adalah dua fragmen grafik integrand untuk pada skala yang berbeda:
Perkiraan asimtotik urutan terkemuka adalah
dan penyempurnaan lebih lanjut (jauh lebih kecil) menambahkan istilah
Grafik nilai yang diperkirakan sebagai fungsi terlihat sebagai berikut:
Sekarang tiba pertanyaan saya: untuk melihat secara visual seberapa bagus perkiraannya, saya ingin membandingkannya dengan "nilai sebenarnya" dari integral, atau lebih tepatnya ke perkiraan yang baik untuk integral yang sama menggunakan algoritma independen. Karena kecilnya koreksi subleading, saya berharap ini akan menjadi sangat dekat.
Akhirnya saya mencoba keberuntungan saya dengan integrator Monte-Carlo menggunakan sampel penting yang saya terapkan, tetapi saya tidak berhasil mendapatkan hasil yang stabil juga.
sumber
Jawaban:
Gunakan teorema Plancherel untuk mengevaluasi integral ini.
sumber
Asimptotik
Sebagai output Anda mendapatkan sinus yang cukup bagus yang bertepatan dengan yang Anda dapatkan di atas.
Jika Anda ingin menemukan koefisien berikut, sedikit lebih canggih jika diperlukan. Gagasan kode di bawah ini adalah untuk mengambil beberapa nilai batas atas yang tinggi dan untuk "rata-rata" hasilnya.
Penjelasan
Contoh sederhana
Masalahmu
Kembali ke integral dari program Konstantin dan Yaroslav, Anda dapat melihat bahwa ia berperilaku persis sama dengan sinus - integral sebagai fungsi batas atas. Itu berarti Anda hanya perlu menghitung nilai dengan . Di bawah ini adalah plot dari beberapa nilai seperti itu dengan .Ix0(λ)=2∫x00cos(λcos(x))sinc(x)dx x0=πN+π2 λ=12π
Di sini Anda dapat melihat hasil dari metode percepatan lain. Saya mengatur ulang jumlah parsial dengan cara berikut dan mendapatkan urutan baru yang konvergennya jauh lebih cepat. Trik itu juga berguna jika Anda ingin mengevaluasi integral dengan presisi tinggi.S′N=12(SN+SN+1) S′N
sumber
Metode Ooura untuk integral sinus Fourier bekerja di sini, lihat:
Ooura, Takuya, dan Masatake Mori, Formula eksponensial ganda yang kuat untuk integral tipe Fourier. Jurnal matematika komputasi dan terapan 112.1-2 (1999): 229-241.
Saya menulis sebuah implementasi dari algoritma ini tetapi tidak pernah bekerja untuk membuatnya cepat (oleh, katakanlah caching node / bobot), tetapi tetap saja, saya mendapatkan hasil yang konsisten pada segala sesuatu di luar presisi float:
Berikut kodenya:
Anda tidak dapat benar-benar melihat perbedaan antara quadrature dan asimptotik karena mereka berada tepat di atas satu sama lain kecuali sebagai :λ→0
sumber