Kuadratur Adaptif

Kuadratur adaptif secara otomatis membagi interval integrasi di mana integran sulit, menggunakan estimasi galat lokal untuk memenuhi akurasi yang diminta dengan evaluasi fungsi sesedikit mungkin.

Definition

Kuadratur adaptif adalah strategi integrasi numerik apa pun yang menggunakan estimasi galat aproksimasi lokal untuk memutuskan di mana dan seberapa halus domain integrasi harus dibagi sehingga toleransi galat keseluruhan yang ditentukan tercapai secara efisien.

Scope

Topik ini mencakup estimasi galat lokal dengan membandingkan aturan-aturan dengan orde atau tingkat penghalusan yang berbeda, pembagian dua interval secara rekursif (Simpson adaptif dan Gauss-Kronrod adaptif), anggaran galat global dan kriteria penghentian, penanganan singularitas dan fitur tajam, serta desain integrator otomatis produksi seperti yang ada di pustaka QUADPACK.

Core questions

• Bagaimana galat lokal dari estimasi kuadratur dihitung tanpa mengetahui integral yang tepat?
• Bagaimana pembagian rekursif memusatkan upaya di mana integran paling bervariasi?
• Kriteria penghentian apa yang secara andal mencapai toleransi yang diminta sambil menghindari pekerjaan yang sia-sia?
• Bagaimana singularitas dan diskontinuitas yang dapat diintegrasikan dideteksi dan ditangani secara kuat?

Key theories

Estimasi galat lokal dan pembagian

Membandingkan estimasi kasar dan yang lebih halus (atau berorde lebih tinggi) pada subinterval menghasilkan estimasi galat lokal; jika melebihi bagian toleransi yang dialokasikan untuk subinterval tersebut, subinterval dibagi dan prosedur berulang, jika tidak, kontribusinya diterima.

Strategi adaptif global

Alih-alih memperlakukan subinterval secara independen, integrator adaptif global menyimpan antrean subinterval yang diurutkan berdasarkan galat estimasi dan selalu menyempurnakan yang terburuk, yang menangani singularitas terlokalisasi secara efisien dan mendasari rutin QUADPACK.

Mechanisms

Pada setiap subinterval, integrator mengevaluasi sepasang aturan tertanam — misalnya pasangan Gauss-Kronrod atau dua estimasi Simpson pada penghalusan yang berbeda — yang perbedaannya mengestimasi galat lokal. Metode adaptif lokal berulang dengan membagi dua subinterval mana pun yang galat estimasinya terlalu besar. Metode adaptif global mempertahankan antrean prioritas subinterval yang dikunci berdasarkan galat estimasi dan berulang kali membagi subinterval terburuk saat ini hingga estimasi galat yang dijumlahkan memenuhi toleransi. Ekstrapolasi dan penanganan bobot khusus ditambahkan untuk mengatasi singularitas titik akhir dan integran osilasi.

Clinical relevance

Kuadratur adaptif adalah apa yang diandalkan oleh rutin integrasi tujuan umum dalam perangkat lunak ilmiah untuk memberikan hasil dengan akurasi yang ditentukan pengguna tanpa pengguna harus menganalisis integran; ini penting untuk integran dengan puncak, perilaku lapisan batas, atau singularitas terintegralkan yang akan mengalahkan aturan tetap, dan ini mendasari integrator otomatis dalam paket numerik dan statistik yang banyak digunakan.

History

Integrasi otomatis yang dikendalikan galat berkembang pesat pada tahun 1970-an dan awal 1980-an, yang berpuncak pada paket QUADPACK (1983), yang rutin Gauss-Kronrod adaptifnya dengan ekstrapolasi menjadi standar de facto dan kemudian diadopsi, diporting, atau diimplementasikan kembali dalam banyak sistem perangkat lunak numerik dan statistik.

Key figures

• Robert Piessens
• Philip J. Davis
• Philip Rabinowitz

Seminal works

• davis1984
• piessens1983

Frequently asked questions

Bagaimana integrator adaptif mengetahui galat jika tidak mengetahui jawabannya?

Ini mengestimasi galat lokal dengan membandingkan dua aproksimasi dengan akurasi yang berbeda pada subinterval yang sama — misalnya aturan orde lebih tinggi dan orde lebih rendah. Perbedaannya mengaproksimasi galat dan memandu di mana harus menyempurnakan, meskipun integral sebenarnya tidak diketahui.

Kapan kuadratur adaptif mengalami kesulitan?

Ini dapat disesatkan oleh integran yang mulus pada titik-titik sampel tetapi memiliki fitur tersembunyi di antaranya, oleh integran yang sangat berosilasi, atau oleh singularitas yang tidak dapat diintegrasikan. Aturan khusus, transformasi, atau metode integrasi osilasi kemudian diperlukan.