Page 31 - 6760
P. 31
ПРАКТИЧНЕ ЗАНЯТТЯ № 7
Тема: Метод Рунге-Кутта для систем диференціальних
рівнянь рівнянь.
Технічне забезпечення: ПЕОМ середовище програмування
Короткі теоретичні відомості
Для чисельного інтегрування систем диференціальних
рівнянь та диференціальних рівнянь вищих порядків
застосовуються ті ж методи і схеми, які наведені в роботі №7,
єдиною відмінністю є те, що вони записуються у векторній
формі; наприклад, для системи:
x f ;,tx x 0 (7.1)
x
0
схему четвертого порядку точності записується у вигляді:
k f x , ,
t
1 n n
k f x 5 . 0 ; tk 5 . 0 ,
2 n 1 n
k f x 5 . 0 ; tk 5 . 0 , (7.2)
3 n 2 n
k f x 5 . 0 ; tk 5 . 0 ;
4 n 3 n
x x k 2k 2k k . (7.3)
n 1 n 1 2 3 4
6
Кожен з виразів k буде складатись з такої кількості
i
розрахункових формул, яка дорівнює кількості невідомих в
системі рівнянь. Таким чином, для інтегрування систем из
кроком по часу одержується явна обчичлювальна схема.
Оскільки будь-яке рівняння вищого порядку може бути
зведене до системи диференціальних рівнянь з n невідомими (n –
порядок рівняння), то підхід (7.2 – 7.3) може бути застосований
до інтегрування таких рівняннь.
Приклад 1. Записати рівняння x 2 x x 2 x t f , 0 xx ,
0
0
x x у вигляді системи диференційних рівнянь. Вводимо нову
1
змінну x . В такому випадку рівняння записується у вигляді
y
системи:
yx
2
y xt f 2xy
31