Решение задачи Дирихле методом конечных разностей

Рассмотрим задачу Дирихле в прямоугольной области (рис. 3.3.2). Тогда ее математическую формулировку можно представить в виде:

(3.3.8)

Здесь все величины заданы, кроме функции (искомая функция).

Заметим, что при решении задачи на ЭВМ естественным является задание направления оси вниз, поскольку это соответствует естественному порядку вывода (печати) на экран или принтер (вывод информации последовательно сверху вниз). Это особенно важно при «промежуточной» печати во время счета.

Рис. 3.3.2. Прямоугольная область, на которой рассматривается задача Дирихле.

Рис. 3.3.3. Дискретизация рассматриваемой области.

Разобьем исходную прямоугольную область на мелкие прямоугольники с помощью сетки с постоянным шагом по первому направлению (по координате ) и с постоянным шагом – по второму направлению (по координате ):

; , (3.3.9)

где и – число узлов сетки по 1-му и 2-му направлению, соответственно (рис. 3.3.3).

Будем рассматривать задачу (3.3.8) в узлах сетки – в точках пересечения линий разбиения области. Каждый узел сетки имеет номер, определяемый двумя величинами , где – номер узла по направлению оси , – номер узла по направлению оси .

Введем обозначения:

; ; . (3.3.10)

Вторые производные для внутренних узлов сетки по каждому направлению заменим вторыми разностями, т.е.

; . (3.3.11)

Тогда краевая задача в конечных разностях примет вид:

– для внутренних узлов:

, , ; (3.3.12)

– для граничных узлов:

(3.3.13)

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

С алгоритмической точки зрения наиболее простыми способами решения такой системы являются итерационные, в частности, метод Зейделя и метод простой итерации. Для этого, выделив из уравнений (3.3.12) диагональный член , получим:

. (3.3.14)

Далее, задав начальное приближение для внутренних точек и присвоив заданные значения для граничных точек, реализуем алгоритм метода Зейделя в виде

(3.3.15)

, , .

Здесь и в дальнейшем – номер итерации.

Счет ведется до тех пор, пока не будет выполнено условие

, (3.3.16)

где – заданная точность; – оценка погрешности на -м шаге итерации.

Заметим, что сходимость метода Зейделя для задач такого типа доказывается в более подробных курсах вычислительной математики. Она следует из, так называемой, положительной определенности оператора (матрицы) краевой задачи. Отметим, что в данном случае диагонального преобладания в матрице нет (и не требуется).

Дата добавления: 2014-01-20