الگوریتم حل AWE

الگوریتم حل AWE

الگوریتم حل AWE
فرض کنید می خواهید معادله ای از فرم را حل کنید
که در آن A یک ماتریس مربع، x یک کمیت مجهول، y یک کمیت شناخته شده، و k یک پارامتر است. به طور کلی، اجزای معادله به k بستگی دارد . ارزیابی شکل موج مجانبی (AWE) یک روش اساسی مبتنی بر بسط Padé یا بسط تیلور است که می تواند برای سرعت بخشیدن به حل چنین معادلاتی برای متغیر k به طور قابل توجهی استفاده شود.
الگوریتم مورد استفاده برای حل کننده AWE از شرح بخش 13.4 Ref پیروی می کند. 21 . آن ارائه، به نوبه خود، از مقالات اصلی پیروی می کند ( مرجع 22 و مرجع 23 ). شکل کلی مسائلی که حل کننده AWE برای آن در نظر گرفته شده است
که در آن وابستگی به k می تواند غیر خطی باشد. یک بسط تیلور کوتاه شده x ( k ) حول مقدار پارامتر 0 را می توان به صورت نوشتاری
(20-21)
بردارهای ضریب مجهول
معمولاً در ادبیات به لحظاتی اشاره می شود. با تمایز مکرر نسبت به k و ارزیابی در 0 ، گشتاورها را می توان بر حسب A ( 0 ) (و مشتقات آن) و y ( 0 ) (و مشتقات آن) به صورت بیان کرد.
و
جایی که ≥ 1 . ظاهراً سیستم خطی باید برای چند سمت راست حل شود. علاوه بر این، مشتقات A ( k ) و y ( k ) باید محاسبه شوند. ممان‌ها نیز به انتخاب 0 بستگی دارند و از آنجا که انبساط‌ها احتمالاً باید در اطراف چندین نقطه انجام شوند، ممکن است چندین مرحله حل مورد نیاز باشد. هنگامی که ممان ها در دسترس هستند، می توان از آنها برای نشان دادن هر جزء، l ( k ) از x ( k ) بر حسب تقریب Padé استفاده کرد.
(20-22)
که در آن L برابر با 1، 2 یا 3 خواهد بود. اگر = 3 و = 6، j می تواند با دستکاری معادله 20-21 و معادله 20-22 ، راه حلی برای
i سپس توسط داده می شود
به طور معمول، شما به پاسخ پارامتر مقداری علاقه مند هستید ، در یک بازه، . رویکرد تکراری زیر به دقت معقولی برای کل بازه دست می یابد:
1
فرض کنید min 1 و max 2 .
2
بسط های حول min و max را محاسبه کنید . نشان دادن تقریب داده شده توسط آن بسط با و ، به ترتیب.
3
یک یا چند نقطه، k را در بازه انتخاب کنید و محاسبه کنید و .
4
اگر تکرارها انجام شود. در غیر این صورت، فاصله اصلی را نصف کنید و این کار را برای هر دو بازه تکرار کنید.