7.2. Разностный метод решения краевой задачи

Рассмотрим краевую задачу для дифференциального уравнения второго порядка следующего вида:

(7.25)

(7.26)

зададим шаг , n – целое. Точки , примем за узлы сетки, – неизвестные значения искомого решения в узлах. Выразим производную в узлах сетки по формуле численного дифференцирования

.

Пусть

.

Вместо дифференциальной краевой задачи (7.25)-(7.26) будем иметь разностную краевую задачу

(7.27)

(7.28)

где y_j – приближенное значение точного решения y(x_j) в узлах .

Перепишем систему линейных алгебраических уравнений (7.27)-(7.28) в виде

.

Эта система с трех диагональной матрицей при на [0, H] имеет решение, причем единственное, которое может быть получено методом прогонки, при этом условие гарантирует устойчивость прогонки. Дадим оценку этому решению.

Лемма 1. Пусть и числа таковы, что . Тогда для всех j. Пусть . Предположим, что . Следовательно, . Пусть q – наименьшее целое, для которого . Из определения d и q имеем: .

Тогда

- противоречие с .

Лемма 1 доказана.

Лемма 2. Если , то для любой системы чисел z_j выполняется неравенство

,

где .

Введем в рассмотрение функцию

Через обозначим

.

Очевидно, что .

Это многочлен второй степени. Для него конечная разность второго порядка , следовательно,

.

Отсюда следует, что

Очевидно, что

.

Числа удовлетворяют условиям леммы 1. Поэтому . Отсюда следует оценка

.

Имеем неравенство

.

Кроме того,

.

Поэтому имеем

.

Лемма 2 доказана.

Рассмотрим случай, когда функции P(x) и f(x) дважды непрерывно дифференцируемы. В курсе дифференциальных уравнений доказывается, что когда краевая задача (7.25)-(7.26) имеет единственное решение y(x), которое четырежды непрерывно дифференцируемо. Наша задача – оценить разность для . – это краевые условия.

Рассмотрим

Согласно дифференциальному уравнению (7.25) для любого j

Следовательно,

Левая часть этого равенства есть разность между приближенным значением второй производной в точке x_j, полученным по формуле численного дифференцирования, и точным значением этой производной y(x_j) и равна остаточному члену этой формулы

(7.29)

Согласно (7.27) имеем

. (7.30)

Вычтем (7.30) из (7.29)

т.е.

.

Воспользуемся леммой 2 для чисел . Имеем

.

Таким образом, при , т.е. неограниченном сгущении сетки, решение разностной задачи приближается к решению дифференциальной.

Разностный метод решения краевой задачи (7.25)-(7.26) используется также и при , хотя успешный результат заранее предвидеть трудно. Для оценки получаемого решения в этом случае нужно провести расчеты для различных значений шага h (не менее трех) и убедиться в том, что полученные значения функции в одних и тех же узлах близки между собой и их разность уменьшается, что говорит о стремлении решения к некоторому пределу при .

К оглавлению

Назад к разделу "7.1. Метод Рунге-Кутта"

Вперед к разделу "Литература"