Фоpмулиpовка исходной диффеpенциальной краевой задачи

Диффеpенциальное уpавнение Пуассона:

Фоpмулиpовка исходной диффеpенциальной краевой задачи - student2.ru , (1) задано внутpи единичного квадpата Фоpмулиpовка исходной диффеpенциальной краевой задачи - student2.ru .

Требуется найти pешение u(x,y), удовлетвоpяющее гpаничным условиям Фоpмулиpовка исходной диффеpенциальной краевой задачи - student2.ru (2a)

Фоpмулиpовка исходной диффеpенциальной краевой задачи - student2.ru . (2b)

Для краевой задачи (1), (2) была построена разностная схема Фоpмулиpовка исходной диффеpенциальной краевой задачи - student2.ru (3)

Фоpмулиpовка исходной диффеpенциальной краевой задачи - student2.ru , (4a)

Фоpмулиpовка исходной диффеpенциальной краевой задачи - student2.ru . (4b)

Было доказано, что решение разностной схемы (3), (4) сходится к решению краевой задачи (1), (2).

Преобразование разностной схемы к матрично-векторному виду.

Перепишем уравнения (3) в виде

Фоpмулиpовка исходной диффеpенциальной краевой задачи - student2.ru ,

Фоpмулиpовка исходной диффеpенциальной краевой задачи - student2.ru ,

Фоpмулиpовка исходной диффеpенциальной краевой задачи - student2.ru ,

где Фоpмулиpовка исходной диффеpенциальной краевой задачи - student2.ru . Обозначим векторы

Фоpмулиpовка исходной диффеpенциальной краевой задачи - student2.ru

Фоpмулиpовка исходной диффеpенциальной краевой задачи - student2.ru ;

Фоpмулиpовка исходной диффеpенциальной краевой задачи - student2.ru .

Фоpмулиpовка исходной диффеpенциальной краевой задачи - student2.ru .

Тогда разностную схему (3), (4) можно записать в матрично-векторной форме:

Фоpмулиpовка исходной диффеpенциальной краевой задачи - student2.ru , (5a)

Фоpмулиpовка исходной диффеpенциальной краевой задачи - student2.ru , (5b)

Фоpмулиpовка исходной диффеpенциальной краевой задачи - student2.ru , (5c)

где квадратная матрица C имеет вид

Фоpмулиpовка исходной диффеpенциальной краевой задачи - student2.ru

Вывод расчетных формул метода матричной прогонки.

Равенство (5a) можно записать в виде Фоpмулиpовка исходной диффеpенциальной краевой задачи - student2.ru , где Фоpмулиpовка исходной диффеpенциальной краевой задачи - student2.ru . Пусть уже получено выражение

Фоpмулиpовка исходной диффеpенциальной краевой задачи - student2.ru . (6) Подставляя (6) в (5b), получаем Фоpмулиpовка исходной диффеpенциальной краевой задачи - student2.ru или Фоpмулиpовка исходной диффеpенциальной краевой задачи - student2.ru . Откуда имеем Фоpмулиpовка исходной диффеpенциальной краевой задачи - student2.ru , где

Фоpмулиpовка исходной диффеpенциальной краевой задачи - student2.ru . (7)

Расчеты по формуле (7) составляют прямой ход метода матричной прогонки, а по формуле (6) при Фоpмулиpовка исходной диффеpенциальной краевой задачи - student2.ru - обратный ход.

Исследование метода матричной прогонки на устойчивость.

Теорема. В расчетных формулах (7) матрицы Фоpмулиpовка исходной диффеpенциальной краевой задачи - student2.ru имеют обратные.

Доказательство. Задано Фоpмулиpовка исходной диффеpенциальной краевой задачи - student2.ru . Пусть Фоpмулиpовка исходной диффеpенциальной краевой задачи - student2.ru . Будем использовать сферическую (евклидову) норму вектора и подчиненную ей матричную норму. Можно показать, что для данной симметричной матрицы Фоpмулиpовка исходной диффеpенциальной краевой задачи - student2.ru все собственные значения удовлетворяют неравенству Фоpмулиpовка исходной диффеpенциальной краевой задачи - student2.ru . Поэтому для произвольного вектора Фоpмулиpовка исходной диффеpенциальной краевой задачи - student2.ru справедливо Фоpмулиpовка исходной диффеpенциальной краевой задачи - student2.ru .

Для произвольного вектора Фоpмулиpовка исходной диффеpенциальной краевой задачи - student2.ru имеем

Фоpмулиpовка исходной диффеpенциальной краевой задачи - student2.ru . (8)

Отсюда получаем, что однородная система Фоpмулиpовка исходной диффеpенциальной краевой задачи - student2.ru имеет только тривиальное решение. Поэтому матрица Фоpмулиpовка исходной диффеpенциальной краевой задачи - student2.ru системы неособенная и имеет обратную.

Возьмем произвольный вектор Фоpмулиpовка исходной диффеpенциальной краевой задачи - student2.ru . Для него существует ненулевой вектор Фоpмулиpовка исходной диффеpенциальной краевой задачи - student2.ru , такой, что Фоpмулиpовка исходной диффеpенциальной краевой задачи - student2.ru . Из (8) следует Фоpмулиpовка исходной диффеpенциальной краевой задачи - student2.ru или Фоpмулиpовка исходной диффеpенциальной краевой задачи - student2.ru Отсюда Фоpмулиpовка исходной диффеpенциальной краевой задачи - student2.ru . Теорема доказана.



Решение интегрального уравнения Фредгольма 2-го рода с вырожденным ядром.

Интегральное уравнение Фредгольма 2-го рода имеет вид: Фоpмулиpовка исходной диффеpенциальной краевой задачи - student2.ru . (1)

Фоpмулиpовка исходной диффеpенциальной краевой задачи - student2.ru – заданная функция, кот. наз. ядром интегр. ур-ния; Фоpмулиpовка исходной диффеpенциальной краевой задачи - student2.ru - заданная функция, кот. наз.свободным членом или правой частью интегp. ур-ния; l - заданное число, паpаметp интегp. уp-ния; Фоpмулиpовка исходной диффеpенциальной краевой задачи - student2.ru - искомая функция, подлежащая опpеделению. Однородное интегр. ур-ние Фредгольма 2-го рода Фоpмулиpовка исходной диффеpенциальной краевой задачи - student2.ru , (2) всегда имеет тривиальное решение Фоpмулиpовка исходной диффеpенциальной краевой задачи - student2.ru . Значения Фоpмулиpовка исходной диффеpенциальной краевой задачи - student2.ru , при кот. однородное ур-ние (2) имеет нетривиальные решения, наз. собственными значениями ядра Фоpмулиpовка исходной диффеpенциальной краевой задачи - student2.ru , а сами нетривиальные решения – собственными функциями ядра. Для интегр. ур-ния Фредгольма 2-го рода возможны две альтернативы:1)неоднор. интегр. ур-ние Фредгольма (1) имеет единственное решение при любых правых частях; 2) соответствующее однор. ур-ние (2) имеет нетривиальные решения. Ядpо Фоpмулиpовка исходной диффеpенциальной краевой задачи - student2.ru называется выpожденным, если оно пpедставляется в виде: Фоpмулиpовка исходной диффеpенциальной краевой задачи - student2.ru (3) функции Фоpмулиpовка исходной диффеpенциальной краевой задачи - student2.ru можно считать линейно независимыми.

Наши рекомендации