8.2. Задача Коши для оду

В зависимости от вида ДУ (1) задача Коши формируется следующим образом.

1. Если n = 1, то требуется найтиY =Y(x), удовлетворяющую уравнению:

(4)

и принимающую при x =x₀заданное значениеY₀:

Y(x₀) = Y₀ . (5)

Для определенности будем считать, что решение нужно получить для значений x >x₀. В качестве начального значения может быть произвольноеx, но чаще всего принимаютx₀= 0, что не влияет на разработку численного метода для (4). Заметим, что все численные методы разработаны для решения ОДУ именно первого порядка.

2. Задача Коши для ОДУ n-го порядка

; (6)

найти Y =Y(x), удовлетворяющую (6) и начальным условиям

,, …,; (7)

где – есть заданные числа.

3. Задача Коши для системы ДУ:

(8)

Задача Коши для системы (8) заключается в отыскании Y_i(x) (), удовлетворяющих (8) и начальным условиям:

;; … ;. (9)

Численные методы для решения ОДУ (4) и (5) применяются и для решения (8) и (9).

Дифференциальное уравнение n-го порядка (6) может быть приведено к системе (8) путем введения новых неизвестных функцийY_i(x),:

,, …,. (10)

Тогда (6) запишется следующим образом

Если удается найти общее решение для (4), (6), или системы (8), то задача Коши сводится к отысканию значений произвольных постоянных. Как правило, она решается приближенно.

8.3. Численные методы решения задачи Коши

Для решения задачи Коши (4) и (5) по технологии разностных методов введем последовательность точек x₀,x₁, ...,x_n и шагиh_i=x_i+1 –x_i(i = 0,1,...,n–1). В каждом узлеx_iвместо значений функцииY(x_i) вводятся числаy_i, как результат аппроксимации точного решенияY(x) на данном множестве точек. Функциюy, заданную в виде таблицы {x_i,y_i} называют сеточной функцией. Заменяя значение производной в уравнении (4) отношением конечных разностей осуществляем переход от дифференциальной задачи (4), (5) относительно функцииY(x) к разностной задаче относительно сеточной функции

; (11)

y₀ =Y₀ . (12)

Это разностное уравнение в общем виде, а конкретное выражение правой части для (11) зависит от способа аппроксимации производной. Для каждого численного метода получается свой вид уравнения (11).

Если в правой части уравнения (11) отсутствует y_i+1, т.е. значениеy_i+1вычисляется поkпредыдущим значениямy_i,y_i–1, ...,y_i–k+1, , то разностная схема называется явной. При этом имеет местоk-шаговый метод:k = 1 – одношаговый,k=2 – двухшаговый и т.д., т.е. в одношаговых методах для вычисленияy_i+1используется лишь одно найденное значение на предыдущем шагеy_i, в многошаговом – многие из них.

Если y_i+1входит в правую часть (11), то это будут неявные методы, реализация которых носит только итерационный характер.

8.3.1. Одношаговые методы решения задачи Коши

Простейшими численными методами для решения задачи Коши для ОДУ являются следующие методы.

1. Метод Эйлера

Этот метод основан на разложении искомой функции Y(x) в ряд Тейлора в окрестностях узлов системыx =x_i(i = 0, 1, 2, …,n), в котором отбрасываются все члены, содержащие производные второго и более высоких порядков. Как правило, используется равномерная сеткаx =x_i+1 –x_i =h=const(i =). Разложение запишем в виде

. (13)

Заменяя значение функции Y(x) в узлах сеткиx_i значениями сеточной функции и используя уравнение (4), получим

.

Тогда из (13) получим

. (14)

При i= 0, для узлаx =x₁:.

Далее по алгоритму (14)

;

. . .

.

Геометрическая интерпретация имеет вид:

EMBED Word.Picture.8

На рисунке линия «0» – точное решение, линии «1» и «2» – приближенные решения.

Искомая интегральная кривая y(x), проходящая через точку (x₀,y₀), заменяется ломаной с вершинами в точках (x_i,y_i). Каждое звено ломаной имеет направление, совпадающее с направлением интегральной кривой (4), которая проходит через точку (x_i,y_i).

Блок-схема алгоритма будет иметь следующий вид:

EMBED Word.Picture.8

Вывод полученных результатов выполняется на каждом шаге, но если необходимо сохранить результаты, то следует ввести массив значений y₀,y₁, ...,y_n.

Локальная погрешность метода Эйлера, как видно из (13), оценивается, как О(h²). Весь интервал [a,b] разбивается наnчастей, тогда общая погрешность

nО(h²) =О(h²) =О(h) – 1-й порядок.

Для оценки погрешности при машинном расчете пользуются двойным просчетом, т.е. на отрезке [x_i,x_i+1] расчет повторяют с шагомh/2 и погрешность более точного решенияу^*_i+1(при шагеh_i/2) оценивается как разность.