To'g'ridan-to'g'ri ko'p tortishish usuli - Direct multiple shooting method

Hududida matematika sifatida tanilgan raqamli oddiy differentsial tenglamalar, to'g'ridan-to'g'ri ko'p tortishish usuli a raqamli usul ning echimi uchun chegara muammolari. Usul yechim izlanadigan intervalni bir necha kichikroq oraliqlarga ajratadi, har bir kichik intervallarda boshlang'ich qiymat masalasini echadi va butun intervalda yechim hosil qilish uchun qo'shimcha mos keladigan shartlarni qo'yadi. Usul chiziqsizlikni taqsimlashda sezilarli yaxshilanishni tashkil etadi raqamli barqarorlik bitta tortishish usullari.

Yagona tortishish usullari

Shunga o'xshash chegara muammolarini (BVP) hal qilish uchun tortishish usullaridan foydalanish mumkin

{ displaystyle y '' (t) = f (t, y (t), y '(t)), quad y (t_ {a}) = y_ {a}, quad y (t_ {b}) = y_ {b},}

unda vaqt belgilanadi t_a va t_b ma'lum va biz izlayapmiz

{ displaystyle y (t), quad t in (t_ {a}, t_ {b}).}

Yagona tortishish usullari quyidagicha davom etadi. Ruxsat bering y(t; t₀, y₀) boshlang'ich qiymat muammosining echimini (IVP) belgilang

{ displaystyle y '' (t) = f (t, y (t), y '(t)), quad y (t_ {0}) = y_ {0}, quad y' (t_ {0}) ) = p}

Funktsiyani aniqlang F(p) orasidagi farq sifatida y(t_b; p) va belgilangan chegara qiymati y_b: F(p) = y(t_b; p) − y_b. Keyin har bir echim uchun (y_a, y_b) bizda mavjud chegara muammosi y_a=y₀ esa y_b a ga to'g'ri keladi ildiz ning F. Ushbu ildizni har qanday kishi hal qilishi mumkin ildiz topish usuli ma'lum bir uslubga bog'liq bo'lgan old shartlar qondirilishini hisobga olgan holda. Buning uchun ko'pincha taxmin qilish kerak bo'ladi y_a va y_b. Odatda, analitik ildiz topish mumkin emas va shunga o'xshash takroriy usullar Nyuton usuli ushbu vazifa uchun ishlatiladi.

Chegaraviy masalalarni raqamli echimi uchun bitta tortishishning qo'llanilishi bir nechta kamchiliklarga duch keladi.

Berilgan dastlabki qiymat uchun y₀ IVP eritmasi aniq intervalda mavjud bo'lishi kerak [t_a,t_b] funktsiyani baholashimiz uchun F uning ildizi izlanmoqda.

Yuqori chiziqli yoki beqaror ODElar uchun bu dastlabki taxminni talab qiladi y₀ haqiqiy, ammo noma'lum echimga juda yaqin bo'lish y_a. Haqiqiy echimdan ozgina tanlangan dastlabki qiymatlar o'ziga xosliklarga yoki ODE hal qiluvchi usulining buzilishiga olib kelishi mumkin. Bunday echimlarni tanlash, ammo takroriy ildiz topish usulida muqarrar.

Sonli aniq raqamlar butun vaqt oralig'ida ODE ni echishga imkon beradigan boshlang'ich qiymatlarni topishni umuman imkonsiz qilishi mumkin.
ODE ning nochiziqliligi samarali ravishda nochiziqlikka aylanadi Fva chiziqli bo'lmagan tizimlarni echishga qodir bo'lgan ildizlarni aniqlash usulini talab qiladi. Bunday usullar odatda sekinroq birlashadi, chunki chiziqsizliklar yanada og'irlashadi. Chegaraviy muammolarni hal qilishning samaradorligi bundan aziyat chekadi.
Hatto barqaror va yaxshi konditsionerlangan ODElar ham beqaror va shartli bo'lmagan BVPlarni yaratishi mumkin. Dastlabki qiymat taxminining ozgina o'zgarishi y₀ ODE echimida juda katta qadam yaratishi mumkin y(t_b; t_a, y₀) va shuning uchun funktsiya qiymatlarida F uning ildizi izlanmoqda. Ildizni analitik bo'lmagan usullar kamdan-kam hollarda ushbu xatti-harakatga dosh bera oladi.

Ko'p marta otish

To'g'ridan-to'g'ri tortishish usuli intervalni ajratadi [t_a, t_b] qo'shimcha katakchalarni kiritish orqali

{ displaystyle t_ {a} = t_ {0}

.

Usul qandaydir qiymatlarni taxmin qilish bilan boshlanadi y barcha tarmoq nuqtalarida t_k 0 with bilan k ≤ N - 1. Ushbu taxminlarni quyidagicha belgilang y_k. Ruxsat bering y(t; t_k, y_k) dan chiqadigan eritmani belgilang kpanjara nuqtasi, ya'ni boshlang'ich qiymat masalasining echimi

{ displaystyle y '(t) = f (t, y (t)), quad y (t_ {k}) = y_ {k}.}

Ushbu echimlarning barchasi birlashtirilib, agar qiymatlar bo'lsa, doimiy traektoriyani hosil qiladi y tarmoq nuqtalarida o'yin. Shunday qilib, chegara masalasining echimlari quyidagi sistemaning echimlariga to'g'ri keladi N tenglamalar:

{ displaystyle { begin {aligned} & y (t_ {1}; t_ {0}, y_ {0}) = y_ {1} & qquad qquad vdots & y (t_ {N-1} ; t_ {N-2}, y_ {N-2}) = y_ {N-1} & y (t_ {N}; t_ {N-1}, y_ {N-1}) = y_ {b} . end {hizalangan}}}

Markaziy N−2 tenglamalar - mos kelish shartlari, birinchi va oxirgi tenglamalar - shartlar y(t_a) = y_a va y(t_b) = y_b chegara muammosidan. Ko'p tortishish usuli ushbu tenglamalar tizimini echish orqali chegara masalasini hal qiladi. Odatda, ning modifikatsiyasi Nyuton usuli oxirgi vazifa uchun ishlatiladi.

Bir nechta tortishish va o'z vaqtida parallel usullar

Qabul qilish uchun bir nechta tortishish qabul qilingan parallel uchun hal qiluvchilar dastlabki qiymat muammolari.^[1]Masalan, Parareal parallel ravishda o'z vaqtida integratsiya qilish usuli ga yaqinlashganda ko'p tortishish algoritmi sifatida olinishi mumkin Jacobian.^[2]

Adabiyotlar

^ Kiehl, Martin (1994). "Dastlabki qiymat masalalarini echish uchun parallel ko'p martalik tortishish". Parallel hisoblash. 20 (3): 275–295. doi:10.1016 / S0167-8191 (06) 80013-X.
^ Gander, Martin J.; Vandewalle, Stefan (2007). "Parareal vaqtni tahlil qilish ‐ parallel vaqt ‐ integratsiya usuli". Ilmiy hisoblash bo'yicha SIAM jurnali. 29 (2): 556–578. CiteSeerX 10.1.1.92.9922. doi:10.1137 / 05064607X.

Stoer, Yozef; Bulirsch, Roland (2002), Raqamli tahlilga kirish (3-nashr), Berlin, Nyu-York: Springer-Verlag, ISBN 978-0-387-95452-3. 7.3.5-bo'limlarga va undan keyin qarang.
Bok, Xans Georg; Plitt, Karl J. (1984), "Optimal boshqarish muammolarini to'g'ridan-to'g'ri hal qilish uchun tortishish algoritmi", 9-IFAC Jahon Kongressi materiallari (PDF), Budapesht
Morrison, Devid D.; Riley, Jeyms D .; Zancanaro, Jon F. (1962 yil dekabr), "Ikki nuqtali chegara masalalari uchun bir nechta tortishish usuli" (PDF), Kommunal. ACM, 5 (12): 613–614, doi:10.1145/355580.369128

[1] Kiehl, Martin (1994). "Dastlabki qiymat masalalarini echish uchun parallel ko'p martalik tortishish". Parallel hisoblash. 20 (3): 275–295. doi:10.1016 / S0167-8191 (06) 80013-X.

[gander2007-2] Gander, Martin J.; Vandewalle, Stefan (2007). "Parareal vaqtni tahlil qilish ‐ parallel vaqt ‐ integratsiya usuli". Ilmiy hisoblash bo'yicha SIAM jurnali. 29 (2): 556–578. CiteSeerX 10.1.1.92.9922. doi:10.1137 / 05064607X.

[1]

[2]