Устойчивость решения задачи Коши

Наряду с уравнением (2.1) рассмотрим задачу

, (2.2)

где - возмущения”, вносимые в правую часть уравнения (2.1) и начальные условия.

Задача Коши (2.1) абсолютно устойчива [4], если такое, что из условия следует .

Иначе говоря, решение устойчиво, если “возмущенное” решение “не слишком сильно” отличается от исходного.

Рассмотрим условия устойчивости решения по начальному условию, то есть положим .

Оценим погрешность решения

. (2.3)

Для этого, вычитая уравнение (2.1) из выражения (2.2), получаем:

.

Воспользовавшись теоремой Лагранжа о среднем, для правой части последнего выражения получим:

.

Решаем полученное дифференциальное уравнение

разделяя переменные:

,

,

,

.

Очевидно, что погрешность z(t) не будет возрастать, если . Но это возможно в том случае, когда . Поскольку это требование должно выполняться для любого значения x, необходимо потребовать, чтобы выполнялось условие .

Пример 2.1. Расмотрим задачу: .

Точное решение этого уравнения имеет вид . В случае k < 0 условие устойчивости выполняется:

,

и различие между решениями уменьшается с увеличением значения аргумента (рис. 2.1а).

При k > 0 условие устойчивости не выполняется:

,

и различие между решениями увеличивается с ростом аргумента x (рис. 2.1б).

а б

Рис. 2.1. Примеры сходимости (а) и расходимости (б) решений действительного (нижние кривые) и возмущенного при начальных условиях (верхние кривые).

Метод Пикара [4]

Проинтегрируем уравнение (2.1) с помощью разделения переменных:

,

,

. (2.4)

Используя последнее соотношение, построим итерационный процесс:

, (2.5)

при этом в качестве “нулевого” приближения выбирается .

Интеграл в правой части процедуры (2.4) вычисляется точно или численно.

Пример 2.2. Найдем методом Пикара решение уравнения .

Начальное приближение, в соответствии с идеей метода, .

Первое приближение:

.

Аналогично определяются второе

,

третье

,

четвертое

приближения, и так далее. Становится очевидным, что получаемая последовательность решений представляет собой частичные суммы разложения в ряд Тэйлора экспоненты , являющейся точным решением поставленной задачи.

Рассмотрим условия сходимости метода Пикара.

Пусть определена ограниченная область , в которой правая часть уравнения (2.1) удовлетворяет условию Липшица

.

Вследствие ограниченности области выполнены соотношения:

.

Вычитая из соотношения (2.5) уравнение (2.4), определим погрешность приближенного решения:

.

Оценим модуль погрешности:

.

Поскольку (решение рассматривается внутри G), получаем:

,

,

,...,

.

В силу того, что

,

приходим к заключению, что погрешность решения уменьшается с увеличением числа итераций:

,

то есть метод Пикара сходится к точному единственному решению исходной задачи.


Понравилась статья? Добавь ее в закладку (CTRL+D) и не забудь поделиться с друзьями:  



double arrow
Сейчас читают про: