Приближенные методы решения дифференциальных уравнений
,
или, что тоже самое,
. (5.15)
Таким образом, мы показали, что всякое решение задачи Коши (5.2), (5.6) есть решение интегрального уравнения (5.15). С другой стороны, если y(x)- решение интегрального уравнения (5.15), то дифференцируя (5.15) по x, получаем, что y(x)- решение задачи Коши (5.2), (5.6).
Решение интегрального уравнения (5.15) будем искать с помощью метода последовательных приближений. Положим
y0(x)=y0, . (5.16)
— (5.17)
сжимающий [12], то последовательные приближения (5.16) сходятся к решению интегрального уравнения (5.15), а, следовательно и дифференциального уравнения y’ = f(x,y), удовлетворяющего условию y(x0) = y0. Желающие могут познакомиться с доказательством сжимаемости оператора (5.17) в [12].
Пример №1 . Найдём с помощью метода последовательных приближений решение уравнения y’ = y, удовлетворяющее условию y(0)=1. Подставляя y(0)=1 в (5.16), получаем
y0=1, …,
С другой стороны, решая исходную задачу Коши, имеем y = e x .
Таким образом, нами получено разложение функции e x в ряд Тейлора в нуле (ряд Маклорена).
Перейдём теперь к изложению численного метода Эйлера решения задачи Коши (5.2), (5.6). Разобьём отрезок [a,b], на котором мы ищем решение, на части точками x0 = a 2 )=yi+hf(xi,yi)+o(h 2 ).
Сравнивая с (5.17) видим, что погрешность формулы (5.17) равна o(h 2 ). К сожалению, метод Эйлера накапливает ошибку от шага к шагу. Поэтому на практике пользуются либо модификациями метода Эйлера, например методом прогноза и коррекции [14], либо другими методами, в частности методом Рунге-Кутта [14].
Метод последовательных приближений решения дифференциального уравнения
Пусть требуется найти решение дифференциального уравнения
Будем предполагать, что в некотором прямоугольнике для уравнения (1) выполнены условия а) и б) теоремы существования и единственности решения задачи (1)-(2).
Решение задачи (1)-(2) может быть найдено методом последовательных приближений , который состоит в следующем.
Строим последовательность функций, определяемых рекуррентными соотношениями
В качестве нулевого приближения можно взять любую функцию, непрерывную в окрестности точки , в частности — начальное значение Коши (2). Можно доказать, что при сделанных предположениях относительно уравнения (1) последовательные приближения сходятся к точному решению уравнения (1), удовлетворяющему условию (2), в некотором интервале , где
Оценка погрешности, получаемой при замене точного решения n-м приближением , даётся неравенством
где . Применяя метод последовательных приближений, следует остановиться на таком , для которого не превосходит допустимой погрешности.
Пример 1. Методом последовательных приближений найти решение уравнения , удовлетворяющее начальному условию .
Решение. Очевидно, что для данного уравнения на всей плоскости выполнены условия теоремы существования и единственности решения задачи Коши. Строим последовательность функций, определяемых соотношениями (3), приняв за нулевое приближение :
Ясно, что при . Непосредственной проверкой убеждаемся, что функция решает поставленную задачу Коши.
Пример 2. Методом последовательных приближений найти приближенное решение уравнения , удовлетворяющее начальному условию в прямоугольнике
Решение. Имеем , т. е. . За берем меньшее из чисел , т. е. . Последовательные приближения согласно (4) будут сходится в интервале . Составляем их
Абсолютная погрешность третьего приближения не превосходит величины
Замечание. Функция должна удовлетворять всем условиям теоремы существования и единственности решения задачи Коши.
Следующий пример показывает, что одной непрерывности функции недостаточно для сходимости последовательных приближений.
Пусть функция определена следующим образом:
На множестве , функция непрерывна и ограничена постоянной . Для начальной точки последовательные приближения при имеют вид:
Поэтому последовательность для каждого не имеет, предела, т. е. последовательные приближения не сходятся. Заметим также, что ни одна из сходящихся подпоследовательностей и не сходится к решению, поскольку
Если же последовательные приближения сходятся, то полученное решение может оказаться неединственным , как показывает следующий пример: .
Возьмем начальное условие ; тогда
Беря в качестве нулевого приближения функцию , будем иметь
так что все последовательные приближения равны нулю и поэтому они сходятся к функции, тождественно равной нулю. С другой стороны, функция представляет собой также решение этой задачи, существующее на полупрямой .
http://mathhelpplanet.com/static.php?p=metod-posledovatelnyh-priblizheniy