Жорстке диференціальне рівняння

Жорстке диференціальне рівняння — це диференціальне рівняння, для якого складно отримати розв'язок за допомогою прямих чисельних методів типу Адамса чи Рунге-Кутти. Строгого загальноприйнятого математичного визначення жорстких диференціальних рівнянь та жорстких систем диференціальних рівнянь немає. Жорсткість — це більше концепція для чисельного розв'язку звичайних диференціальних рівнянь. До жорстких систем відносяться задачі хімічної кінетики, нестаціонарні процеси у складних електронних схемах, задачі теплопровідності та дифузії та багато інших.

Диференціальне рівняння порядку n має у загальному випадку n сталих часу, і, коли деякі з них значно відрізняються за величиною, або одна з них достатньо мала порівняно з інтервалом часу, на якому шукають розв'язок, то задача стає «жорсткою» і її практично неможливо розв'язати звичайними чисельними методами. В таких випадках крок повинен бути дуже малим, щоб можна було враховувати приріст тих складових розв'язку, які найшвидше змінюються, навіть після того, як їх внесок стане практично непомітним. Але зменшення кроку веде до збільшення часу обчислень і накопичення похибок. Класичні прямі (явні) методи типу Адамса чи Рунге-Кутти вимагають як правило для розв'язку неприйнятно дрібного кроку. Тому для цих рівнянь застосовують неявні методи.

Мотиваційні приклади

Перший приклад

При використанні для інтегрування жорсткого звичайного диференціального рівняння, явні чисельні методи проявляють нестійкість

Розглянемо задачу Коші

$\,y'(t)=-15y(t),\quad t\geq 0,y(0)=1.\qquad \qquad \qquad \quad \quad (1)$

Точний розв'язок (блакитним) такий

$y(t)=e^{-15t}\,$ with $y(t)\to 0$ as $t\to \infty .\qquad \qquad \qquad \qquad (2)$

Ми шукаємо чисельний розв'язок, який показує таку саму поведінку.

Зображення (праворуч) ілюструє вади різноманітних методів застосованих до цього рівняння.

1. Метод Ейлера з кроком h = 1/4 дико коливається і швидко виходить за межі графу (показано червоним).
2. Метод Ейлера із вдвічі меншим кроком, h = 1/8, призводить до розв'язку у межах графу, але коливається навколо нуля (показано зеленом).
3. Інтегрування за правилом трапецій (тобто, двокроковий метод Адамса-Мултона) заданий як

   $y_{n+1}=y_{n}+{\frac {1}{2}}h\left(f(t_{n},y_{n})+f(t_{n+1},y_{n+1})\right).\qquad (3)$

   Застосування цього метода замість метода Ейлера дає значно кращий результат (синій). Чисельний розв'язок монотонно спадає до нуля, так само як і точний розв'язок.

Другий приклад

Припустимо ми маємо частинку з координатами $(x(t),y(t)),$ і, припустимо що ми хочемо щоб $y-$координата завжди була нулем. Одним зі способів зробити це є додати складову $-ky(t)$ до $y'(t),$ де $k$ це велика додатна стала. Якщо $k$ достатньо велика, тоді частинка ніколи не відійде далеко від $y(t)=0.$ Однак, припустимо що також існує обмеження на $x-$ координату. Отже, нехай ми маємо таке диференціальне рівняння:

${\dot {\mathbf {Y} }}={\frac {d}{dt}}{\begin{pmatrix}-x(t)\\-ky(t)\end{pmatrix}}.$

(Ми також припустимо, що частинка стартує не з положення коли $y_{0}=0.$) Таке рівняння має місце коли частинка притягується до лінії $y=0$ і не так сильно до $x=0.$ Якщо ми розв'яжемо ЗДР достатньо далеко в часі, то ми очікуємо, що частинка прибуде до $(0,0).$

Розглянемо застосування метода Ейлера до цього рівняння. Якщо взяти крок $h,$ отримуємо

$\mathbf {Y} _{new}=\mathbf {Y} _{0}+h{\dot {\mathbf {Y} }}(t_{0})={\begin{pmatrix}x_{0}\\y_{0}\end{pmatrix}}+h{\begin{pmatrix}-x_{0}\\-ky_{0}\end{pmatrix}}.$

Що дає

$\mathbf {Y} _{new}={\begin{pmatrix}(1-h)x_{0}\\(1-hk)y_{0})\end{pmatrix}}.$

Якщо ми розглянемо $y$ складову цього рівняння, ми бачимо, що якщо $|1-hk|>1$ тоді отримане $y_{new}$ матиме абсолютне значення більше ніж $|y_{0}|.$ Інакше кажучи, якщо $|1-hk|>1,$ метод Ейлера не сходиться до відповіді. Отже, нам краще щоб $1-hk>-1$ або $hk<2.$ Найбільший крок який ми можемо взяти мусить бути меншим ніж $2/k.$

Тепер, якщо $k$ велика стала, тоді нам доведеться робити дуже маленькі кроки. Тобто частинка зісковзує у напрямку $(0,0)$ дуже повільно. Навіть якщо частинка зовсім близько до $y_{0}=0,$ нам потрібно робити настільки малі кроки, що її поступ уздовж осі $x$ майже непомітний. У цьому випадку жорсткість з'являється через дуже велике $k.$

Див. також

• Чисельні методи