Pendekatan naif untuk memecahkan persamaan diferensial stokastik (SDE) adalah:
- mengambil metode Runge – Kutta multi-langkah biasa,
- menggunakan diskritisasi yang cukup baik dari proses Wiener yang mendasarinya,
- buat setiap langkah metode Runge – Kutta dianalogikan dengan Euler-Maruyama.
Sekarang, ini gagal pada beberapa level dan saya mengerti mengapa. Namun, sekarang saya ditugaskan untuk meyakinkan orang-orang tentang fakta ini yang memiliki sedikit pengetahuan tentang metode Runge-Kutta dan persamaan diferensial stokastik untuk memulai. Semua argumen yang saya ketahui bukanlah apa-apa yang dapat saya komunikasikan dengan baik dalam konteks yang diberikan. Oleh karena itu, saya mencari argumen yang mudah dimengerti bahwa pendekatan di atas akan gagal.
runge-kutta
education
stochastic-ode
Wrzlprmft
sumber
sumber
Jawaban:
Mari kita ambil persamaan diferensial stokastik:
Berikut adalah beberapa argumen berbeda yang mengarah pada pemahaman intuitif mengapa matematika di balik metode tingkat tinggi diperlukan. Saya akan membahas dalam hal keteraturan yang kuat, yang sama dengan mengatakan "untuk gerak Brown yang diberikan , seberapa baik integral numerik memecahkan lintasan itu?"W( t )
Keteraturan Persamaan
Pertama-tama, metode yang diusulkan gagal untuk memperhitungkan fakta bahwa tidak terus menerus terdiferensiasi. Sebenarnya, Anda dapat menggunakan hasil Rossler untuk menunjukkan bahwa memperluas metode RK normal seperti yang Anda sarankan akan menghasilkan metode konvergen, tetapi mereka hanya akan memiliki urutan kuat 0,5. Alasannya adalah karena mereka diturunkan menggunakan kalkulus dengan X t menjadi terdiferensiasi dan memiliki serangkaian Taylor. Gerakan Brown tidak dapat dibedakan, dan sebaliknya memiliki kontinuitas Pemegang α < 0,5 asXt Xt α < 0,5
Namun, seperti dalam teori perturbasi, proses yang tidak cukup teratur tidak dapat diperluas dalam hal seri Taylor, tetapi dengan keteraturan Pemegang mereka dapat diperluas dalam bentuk seri Puiseux dengan syarat α , yaitu untuk gerak Brown ada ekstensi ke gagasan seri Taylor yang diperluas dalam hal seperti 1α α turunan. Seperti di kalkulus biasa, istilah pertama adalah "istilah linear", yaitu perubahandtuntukΔtdandWtuntukN(0,dt)dan Anda mendapatkan sesuatu yang tepat. Inilah sebabnya mengapa metode, termasuk hal-hal seperti Euler-Maruyama, bertemu dengan urutan yang kuat 0,5: mereka mendapatkan istilah pertama dalam seri Taylor yang benar. Namun, syarat orde yang lebih tinggi perlu memiliki koreksi untuk fakta bahwaXttidak dapat dibedakan secara terus menerus, itulah sebabnya metode normal gagal melakukannya.12 dt Δ t dWt N( 0 , dt ) Xt
Korelasi Seketika dan Integral Iterasi
Itu penjelasan heuristik yang cepat, tetapi ada sedikit lebih dari itu. Mari kita lihat beberapa detail lainnya. Serangkaian Taylor tidak hanya ekspansi dalam hal turunan, tetapi juga dapat dianggap sebagai jumlah yang lebih tinggi untuk diintegrasikan. adalah mengintegrasikan sekali. Tetapi jika Anda menambahkan d t 2 istilah, untuk mendapatkan unit yang tepat yang perlu Anda lakukan integral ganda. d t 2 mudah diintegrasikan dua kali, tetapi apa yang d W i t d W j tXt= X0+ Δ t f( t , Xt) dt2 dt2 dWsayatdWjt ? Ini adalah korelasi instan antara gerakan Brown. Anda perlu tahu ini untuk menghitung integral ganda. Jika Anda hanya melihat rata-rata, Anda dapat melakukannya. Tetapi dalam lintasan apa pun ada korelasi antara gerakan Brown yang berbeda dari sistem persamaan diferensial. Dengan asumsi tidak ada korelasi antara gerakan Brown adalah cara lain untuk mengkarakterisasi perluasan metode deterministik Maruyama, tetapi untuk mendapatkan istilah berikutnya dalam seri (istilah 1.0) Anda harus mendapatkan ini dengan benar. Koreksi Milstein justru menambahkan istilah korelasi ini. Ketika noise diagonal, ini setara dengan memahami bahwa tidak ada korelasi kecuali dengan dirinya sendiri, tetapi korelasi dengan diri sendiri hanyalah varian yang dt , dan harus ada koreksi vs d t , yaitu d W 2 - d t . Ketika ada noise non-diagonal, integral ganda ini harus didekati sehingga korelasi sesaat dari gerakan Brownian diperhitungkan, dan perkiraan umum di sini adalah perkiraan Wiktorsson yang kemudian membuat simulasi noise non-diagonal jadi rumit (Karena tidak ada solusi analitis bahkan untuk integral ganda).dW2t dt dW2- dt
Efek Rata-Rata Difusi
Kesimpulan
Tentu saja, dalam beberapa keadaan ada cara untuk menemukan generalisasi yang tepat yang memang memberikan metode urutan yang lebih tinggi, tetapi saya akan meninggalkan ini sebagai utas yang menggantung karena itulah satu poin dari makalah yang akan saya kirimkan segera. Semoga ini membantu.
sumber