Чи є ярлики для чисельно наближених систем звичайних диференціальних рівнянь, коли автономні?


10

Існуючі алгоритми вирішення ODE обробляють функції , де . Але у багатьох фізичних системах диференціальне рівняння є автономним, тому , , при цьому залишилося. За допомогою цього спрощеного припущення, які покращення можна побачити в існуючих числових методах? Наприклад, якщо , проблема перетворюється на і ми переходимо до зовсім іншого класу алгоритмів інтегрування одновимірних інтегралів. При максимально можливе поліпшення - це зменшення розмірностіyRndydydt=f(y,t)yRnyRntn=1t=dydydt=f(y)yRntn=1 n>1yt=dyf(y)n>1yна 1, оскільки залежний від часу випадок може бути змодельований додаванням до , змінивши домен від на .y y R n R n + 1tyyRnRn+1

Відповіді:


2

Я б сказав, що одне суттєве поліпшення полягає в тому, що в межах підходів, що крокують на час, де ви поширюєте використовуючи карту рішення , ви можете визначте пропагатор (або принаймні його частини) один раз, а потім повторно використовуйте його на кожному кроці.ynyn+1=U(yn)U

Наприклад, у лінійному випадку у вас буде , де - матриця. Оператор рішення складається в основному з експоненціалу матриці. Для автономних систем ця затратна матрична експоненціальна оцінка потрібна лише один раз для повного розповсюдження - на відміну від системи, залежної від часу, де ви повинні виконувати цю оцінку на кожному кроці часу.ty=AyAU(y)=exp(AΔt)y

Для нелінійних систем це не так просто, але залежно від алгоритму можуть бути використані певні дорогі оцінки.

Використовуючи наш веб-сайт, ви визнаєте, що прочитали та зрозуміли наші Політику щодо файлів cookie та Політику конфіденційності.
Licensed under cc by-sa 3.0 with attribution required.