Алгоритм Бімана і Шофілда

2. Алгоритм Бімана і Шофілда

Інший корисний алгоритм, в якому немає нагромаження погрішностей округлення, як в алгоритмі Верле, належить Біміану і Шофілду. Запишемо алгоритм Бімана в наступному виді:

(25а)

і

(25б)

Підмітимо, що точність розразунку траєктторії по схемі (25) не вища, ніж в алгоритмі Верле. Її перевага заключається в тому, що просто вона краще зберігає енергію. Однак алгоритм Бімана не самостартуючий. Алгорим Бімана і алгоритм Верле в швидкісній вормулі викоритані в програмі BEEMAN

Завершимо наше обговорення оротким викладом двох методів, які зазвичай приводяться в підручниках по чисельному аналізу. Один приклад метода предиктор:

 

(26а)

Передбучуване значення коодринати дозволить оприділити прискорення

Тоді, використовуючи , отримаємо скоректироване значення v_n+1 і x_n+1

коректор:

(26б)

Скоректироване значення x_n+1 використовується для визначення нового передбачуваного значення а_n+1 і, значить, нових передбачених значень v_n+1 i x_n+r Ця процедура повторяється до тих пір, доки передбачення і скоррективонане значення x_n+1 відрізняються менше ніж на задану величину! Даний метод можна розробити на схемі більш високого порядку, які зв’язуються між собою не тільки  x_{n+1 ,}  x_n і v_n , але і так само значеннями v_n-1 і v_n-2. Замітимо, що метод предиктора-коректора не являється самостартуючим.

3. Метод Рунге-Кутта

Для пояснення методу Рунге-Кутта подивимось спочатку розв’язок диференціального рівняння першого порядку

(27)

Метод Рунге-Кутти другого порядку для розв’язку рівняння (27) модна, використовуючи стандартні значення, записати наступним чином:

(28)

Сенс формул (28) полягає у наступному: В методі Ейлера допускається, що для екстаполяції в наступну точку модна використовувати нахил кривої f(x_n,y_n)в точці (x_n,y_n) так чи однакше y_n+1=y_n+f(x_n,y_n)*∆x. Однак можна повисити почність оцінки нахилу, якщо методом Ейлера повести екстраполяцію в середню точку відрізку, а потім використати центральну похідну на всьому відрізку. Звідси оцінка нахилу в методі Рунге-Котти рівна

 де

Застосування методу Рунге-Кутти до рівнянь руху Ньютона дає

(29)

Оскільки методи Руиге-Кутти є такими, що самостартуючими, то їх часто використовують для вираховання декількох перших кроків для несамостартуючих алгоритмів.

4. Метод Рунге — Кутта 4-го порядку

Цей метод настільки широко розповсюджений, що його часто називають просто методом Рунге — Кутта.

Розглянемо задачу Коші для системи диференціальних рівнянь довільного порядку, що записується у векторній формі як:

Тоді значення невідомої функції в точці x_n+1 обчислюється відносно значення в попередній точці x_n по такій формулі:

де h— крок інтегрування, а коефіцієнти k _n розраховуються наступним чином:

Це метод 4-го порядку, тобто похибка на кожному кроці становить O(h⁵), а сумарна похибка на кінцевому інтервалі інтегрування є величиною O(h⁴) .

Прямі методи Рунге — Кутта

Група прямих методів Рунге — Кутта є узагальненням методу Рунге — Кутти 4-го порядку. Воно задається формулами

де

Конкретний метод визначається числом s і коефіцієнтами b_i,a_ij i c_i . Ці коефіцієнти часто впорядковують в таблицю

0

c₂ a₂₁

c₃ a₃₁ a₃₂

∙ ∙ ∙ ∙

∙ ∙ ∙ ∙

∙ ∙ ∙ ∙

c_s a_s1 a_s2 ∙ ∙ ∙ a_{s,s − 1}

b₁ b₂ b_{s − 1} b_s

Для коефіцієнтів методу Рунге — Кутта мають справджуватись умови

 

Якщо ми хочемо, щоб метод мав порядок p, то варто так само забезпечити умову  — наближення, отримане по методу Рунге — Кутти. Після багаторазового диференціювання ця умова перетвориться в систему поліноміальних рівнянь на коефіцієнти методу.