Трёхдиагональной матрицей называют матрицу следующего вида:

Системы линейных алгебраических уравнений с такими матрицами встречаются при решении многих задач математики и физики. Краевые условия x1 и xn, которые берутся из контекста задачи, задают первую и последнюю строки. Так краевое условие первого рода F(x = x1) = F1 определит перую строку в виде C1 = 1, B1 = 0, а условие второго рода dF / dx(x = x1) = F1 будет соответствовать значениям C1 = − 1, B1 = 1.
Метод прогонки
Для решения систем вида используется метод прогонки, основанный на предположении, что искомые неизвестные связаны рекуррентным соотношением:


Используя это соотношение, выразим xi-1 и xi через xi+1 и подставим в i-e уравнение:

где Fi - правая часть i-го уравнения. Это соотношение будет выполняться независимо от решения, если потребовать


Отсюда следует:


Из первого уравнения получим:


После нахождения прогоночных коэффициентов α и β, используя уравнение (1), получим решение системы. При этом,

Ссылки
Алгоритм метода прогонки
Костомаров Д.П., Фаворский А.П. "Вводные лекции по численным методам"
0 件のコメント:
コメントを投稿