Численное решение обыкновенных дифференциальных уравнений

Численное решение обыкновенных дифференциальных уравнений

Все дифференциальные уравнения делятся на 2 большие категории: 1. 1. В ДУ входит функция только одной независимой переменной. 2. 2. В ДУ входит функция нескольких независимых переменных.

Обыкновенные дифференциальные уравнения : Общий вид ОДУ: F(x, y, y’’, …, y (n)(n) )=0)=0 Наивысший порядок производной — порядок уравнения. ОДУ часто представляют в виде: yy (n)(n) =f(x, y, y’, y’’, …, y (n(n -1 -1 )) )) уравнение, разрешенное относительно старшей производной.

Общее решение ОДУ — бесконечное множество функций в аналитическом виде. Точное решение ОДУ — только 1 функция (обычно в аналитическом виде). Для нахождения точного решения необходимо задать дополнительные условия (например, начальную точку). Получаем задачу Коши или краевую задачу. При использовании численных методов решения ОДУ решение представляет собой набор точек, которые лежат на кривой истинного решения или расположены вблизи нее (ввиду погрешности метода).

Численные методы решения ОДУ классифицируются следующим образом: 1. ОДНОШАГОВЫЕ и МНОГОШАГОВЫЕ. 2. ЯВНЫЕ и НЕЯВНЫЕ.

Также численные методы решения ОДУ характеризуются следующими показателями: 1. Точность 2. Устойчивость Точность – погрешность, с которой получается решение. Все методы характеризуется определенным порядком точности. Чем выше порядок – тем выше точность. Устойчивость метода характеризует возможность вообще получить достоверный результат.

Решение ОДУ 1 -го порядка Общий вид ОДУ 1 -го порядка: F(x, y, y’)=0 Уравнение разрешаем относительно старшей производной: y’=f(x, y) Начальная точка: xx 00 , y, y 00. . Количество точек на интегральной кривой решения nn Шаг hh , с которым будет изменяться значение аргумента.

Метод Эйлера y’=f(x, y) Представляем производную функции y’y’ в виде конечной разности. За величину приращения аргумента xx принимается величина шага hh : : Тогда: Аналогично: В общем виде: — — формула Эйлераh yy x y y x 01 0 lim’ ), ( 0001 yxfhyy yxf h yy hxx yxfhyy ii iiii 1 1 1112 ), (

Графическая интерпретация метода Эйлера Предположим, что мы имеем истинное решение ОДУ: y=y(x). . xy y=y(x) x 0 x 1 x 2 h hy 0 y 2 y 1 От начальной точки (x(x 00 , y, y 00 )) проводим касательную к графику до пересечения с линией x=xx=x 11. . Получаем новую точку (x(x 11 , y, y 11 )). . Из графика видно, что с увеличением количества шагов погрешность возрастает. Если величину шага hh уменьшить, то результат получится более точным.

Блок-схема метода Эйлера НАЧАЛО ввод f(x, y), x 0 , y 0 , n, h i=1, n x=x 0 +hy=y 0 +h*f(x 0 , y 0 ) КОНЕЦ вывод x, y x 0 =x y 0 =y

Характеристика метода Эйлера Метод Эйлера является одношаговым методом, то есть для расчета последующей точки необходимо знать только координаты предыдущей. Данный метод использует явную схему. В правой части формулы Эйлера стоят все известные величины. Метод Эйлера не является устойчивым методом, поэтому он применяется только для ОДУ, решением которых являются достаточно гладкие функции. Этот метод характеризуется первым порядком точности (точность низкая).

Метод Эйлера-Коши Данный метод является модификацией метода Эйлера. Он основан на том, что половину шага совершается с тангенсом угла наклона касательной в предыдущей точке, а вторую – с тангенсом угла наклона в последующей точке. xy y=y(x) x 0 x 1 h/2 y 0 y 1 Решение, найденное методом Эйлера-Коши Решение, найденное методом Эйлера

Таким образом, формула для решения ОДУ методом Эйлера-Коши будет следующей: — — формула Эйлера-Коши В левой и в правой части формулы Эйлера-Коши стоит неизвестная искомая величина yy i+1 i+1 — метод является неявным. . Для его реализации находится приближенное значение yy i+1 i+1 методом Эйлера, подставляется в правую часть формулы Эйлера-Коши и находится уточненное значение yy i+1 i+1. . Метод является одношаговым, устойчивым, 2 -го порядка точности (более точным, по сравнению с методом Эйлера)), ( 2 111 iiiiiiyxfyxf h yy )~ , (), (2/ ), (~ 1111 1 iiiiii ii iiii yxfyxfhyy hxx yxfhyy

Блок-схема метода Эйлера — Коши НАЧАЛО ввод f(x, y), x 0 , y 0 , n, h i=1, n x=x 0 +hy 1 =y 0 +h*f(x 0 , y 0 ) КОНЕЦ вывод x, yy=y 0 +h/2*(f(x 0 , y 0 )+ f(x, y 1 )) x 0 =x y 0 =y

Метод Рунге-Кутта 4 порядка Данный метод 4 -й порядок точности, является одношаговым, имеет явную схему, но не всегда устойчив. Для реализации этого метода используются достаточно громоздкие формулы Рунге-Кутта: Несмотря на сложные формулы, данный метод является самым распространенным. hxx kkkkyy kyhxfhk k yh xfhk yxfhk ii ii ii 1 43211 34 2 3 1 22 61 , 2, 2), (

Блок-схема метода Рунге-Кутта 4 порядка НАЧАЛО ввод f(x, y), x 0 , y 0 , n, h i=1, n x=x 0 +h КОНЕЦ вывод x, yk 1 =h*f(x 0 , y 0 ) k 2 =h*f(x 0 +h/2, y 0 +k 1 /2) k 3 =h*f(x 0 +h/2, y 0 +k 2 /2) k 4 =h*f(x 0 +h, y 0 +k 3 ) y=y 0 +1/6*(k 1 +2 k 2 +2 k 3 +k 4 ) x 0 =x y 0 =y

Метод Адамса Многошаговый метод (4 -шаговый) Для расчета последующей точки необходимо знать координаты четырех предыдущих точек. Как правило в задачах обычно известна только одна начальная точка. Поэтому три последующие точки вычисляются с использованием одношаговых методов, а затем используется 4 -шаговый метод Адамса. Данный метод имеет 4 -й порядок точности, явную схему, но не всегда устойчив. hxx yxfyxfh yy ii iiiii 1 3322111 ), (9), (37), (59), (

Решение систем обыкновенных дифференциальных уравнений Рассмотрим систему из двух дифференциальных уравнений 1 -го порядка: Оба уравнения необходимо разрешить относительно старшей производной: Пусть заданы начальные условия: xx 00 , y, y 00 , z, z 00. . 0)’, ‘, , , ( 2 1 zyzyx. F ), , (‘2 1 zyxfz zyxfy

Данную систему можно решить любым методом, применимым для решения единичных ОДУ. Метод Эйлера: Метод Эйлера-Коши: hxx zyxfhzz zyxfhyy ii iiiii 1 21 11 ), , ( )~ , ~ , (), , ( 2 ), , (~ 111221 1111111 21 11 iiiiiiii ii iiiii zyxfh zz zyxfh yy hxx zyxfhzz zyxfhyy

Метод Рунге-Кутта 4 -го порядка: )22( 61 ), , ( ) 2, 2, 2( ), , ( 43211 3324 3314 22 23 22 13 11 22 11 12 21 11 llllzz kkkkyy lzkyhxfhl lzkyhxfhk l zk yh xfhl l zk yh xfhk zyxfhl zyxfhk ii ii iii iii

Блок-схема метода Рунге-Кутта 4 порядка для системы двух ОДУ НАЧАЛО ввод f 1 (x, y, z), f 2 (x, y, z), x 0 , y 0 , z 0 , n, h i=1, n x=x 0 +h КОНЕЦ вывод x, y, zk 1 =h*f 1 (x 0 , y 0 , z 0 ) L 1 =h*f 2 (x 0 , y 0 , z 0 ) k 2 =h*f 1 (x 0 +h/2, y 0 +k 1 /2, z 0 +L 1 /2) L 2 =h*f 2 (x 0 +h/2, y 0 +k 1 /2, z 0 +L 1 /2) k 3 =h*f 1 (x 0 +h/2, y 0 +k 2 /2, z 0 +L 2 /2) L 3 =h*f 2 (x 0 +h/2, y 0 +k 2 /2, z 0 +L 2 /2) k 4 =h*f 1 (x 0 +h, y 0 +k 3 , z 0 +L 3 ) L 4 =h*f 2 (x 0 +h, y 0 +k 3 , z 0 +L 3 ) y=y 0 +1/6*(k 1 +2 k 2 +2 k 3 +k 4 ) z=z 0 +1/6*(L 1 +2 L 2 +2 L 3 +L 4 ) x 0 =x y 0 =y z 0 =z

Решение обыкновенных дифференциальных уравнений высших порядков Любое дифференциальное уравнение высшего порядка можно привести к системе дифференциальных уравнений 1 -го порядка путем замены переменных. Рассмотрим дифференциальное уравнение 2 -го порядка: Заданы начальные условия: xx 00 , y, y 00 , y’ 00 Разрешим уравнение относительно старшей производной: Заменим первую производную y’y’ функцией zz. Тогда y’’=z’ , а y’y’ 00 = = zz 00 Получим систему: Решаем полученную систему известными методами. ), , (‘ ‘ zyxfz zy 0)», ‘, , (yyyx. F )’, , (»yyxfy

Пример решения ОДУ 2 -го порядка: 2. 3, 5. 4, 2 2 75, 1′ ‘ : ‘. 2 ‘ 75, 1» : 2. 3)2(‘, 5. 4)2(, 0’5, 3»2 000 43 43 4 zyx x y x z z zy систему. Получим zyзамену. Произведем x y y йпроизводностаршейноотносительуравнение. Разрешим yyyxyyx