Студопедия

Главная страница Случайная страница

Разделы сайта

АвтомобилиАстрономияБиологияГеографияДом и садДругие языкиДругоеИнформатикаИсторияКультураЛитератураЛогикаМатематикаМедицинаМеталлургияМеханикаОбразованиеОхрана трудаПедагогикаПолитикаПравоПсихологияРелигияРиторикаСоциологияСпортСтроительствоТехнологияТуризмФизикаФилософияФинансыХимияЧерчениеЭкологияЭкономикаЭлектроника






случай. Метод Рунге - Кутта четвёртого порядка.






Этот метод также называют классическим методом Рунге - Кутта.

Пусть p =4, c 1=0, c 2= c 3= c 4=1, d 1= d 4= , d 2= d 3= . Из рекуррентных формул (1) получим алгоритм решения задачи Коши классическим методом Рунге - Кутта:

(3)

Графиком приближённого решения является ломаная, последовательно соединяющая точки Pi (xi, yi) (i =0, 1, …, m). С увеличением порядка численного метода, звенья ломаной приближаются к ломаной, образованной хордами интегральной кривой y=j (x), последовательно соединяющими точки (xi, j (xi)) на интегральной кривой.

Правило Рунге практической оценки погрешности решения для численного метода четвёртого порядка имеет вид:

Пример. Решить задачу Коши y'=x+y, классическим методом Рунге - Кутта на отрезке [0; 0, 4]. Найти решение на равномерной сетке с шагом 0, 1 в четырёх узловых точках.

Решение: Так как f (x, y)= x+y, то согласно формулам (3) получаем

;

;

Полагая x 0=0, y 0=1, последовательно находим при i =1:

0, 1(0+0, 05+1+0, 055)=0, 1105

0, 1(0+0, 1+1+0, 1105)=0, 121050

x 1=0+0, 1=0, 1; y 1=

при i =2:

=0, 1(0, 1+0, 05+1, 110342+0, 0605171)=0, 1320859

=0, 1(0, 1+0, 05+1, 110342+0, 06604295)=0, 1326385

=0, 1(0, 1+0, 1+1, 110342+0, 1326385)=0, 1442980

x 2=0, 1+0, 1=0, 2; y 2= .

Далее получаем при i =3, x 3=0, 3, y 3=1, 399717; при i =4, x 4=0, 4; y 4=1, 583648.

Погрешность полученного решения не превышает величины

 






© 2023 :: MyLektsii.ru :: Мои Лекции
Все материалы представленные на сайте исключительно с целью ознакомления читателями и не преследуют коммерческих целей или нарушение авторских прав.
Копирование текстов разрешено только с указанием индексируемой ссылки на источник.