Материал: Моделирование технологических процессов. учеб. пособие. Арсентьев А.В., Строгонов А.В

Подставим числовые значения известных величин:

Трех диагональная матрица

Из первого уравнения системы выражаем y₁:

Подставим y₁ во второе уравнение и выразим y₂:

Подставим y₂ в третье уравнение и выразим y₃:

Подставим y₃ в четвертое уравнение и найдем y₄:

Далее находим значения y₃, y₂, y₁:

Сложив 2-ое уравнение с 3-м получим:

Подставляя во второе уравнение получим:

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

Значение производной в крайней левой точке:

Значение производной в крайней правой точке найдём из системы:

Сложив 2-ое уравнение с 3-м получим:

Подставляя во второе уравнение получим:

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

Значение производной в крайней правой точке:

Найдём значение второй производной из следующей системы:

Общий вид системы отличается от системы (1.5) только правой частью.

Находим из данной системы значения коэффициентов. Сложив 2 уравнение с 3-м получим:

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

Итоговое выражение:

В качестве примера рассмотрим систему:

Выведем формулы для пятиточечного шаблона. В качестве неосновных функций выберем полиномы 1-ой, 2-ой, 3-ей, 4-ой и 5-ой степени.

Применяя метод интерполяции функции полиномами, получаем следующую систему:

1.5. Метод прогонки

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

Рассмотрим линейное дифференциальное уравнение

с двухточечными линейными краевыми условиями

где непрерывны на отрезке [a, b].

От линейного дифференциального уравнения обычным приемом перейдем к конечно-разностному уравнению. Для этого разобьем отрезок [a, b] на n равных частей с шагом Полагая (i = 0, 1, 2, …, n) и вводя обозначения: , для внутренних точек i = 1, 2, …, n – 1) отрезка [a, b] вместо линейного дифференциального уравнения получаем систему конечноразностных уравнений

…………………………………………

………………………………………….

Линейная система состоит из (n + 1) уравнений первой степени относительно неизвестных. Эту систему можно решить обычным способом. Однако мы сейчас укажем другой, более короткий путь, получивший название метода прогонки.

Метод прогонки осуществляется в два этапа: прямая прогонка, при которой находятся вспомогательные коэффициенты и обратная прогонка при которой находятся уже значения самой функции.

Прямая прогонка. Из первого уравнения находим:

,

где

Подставим во второе уравнение значение :

Так как , то получаем следующее:

,

где

В общем случае:

,

где

Обратная прогонка. При , получаем, что

Так как , то

Отсюда находим значение

Алгоритм прогонки

Вычисление прогоночных коэффициентов

При получаем:

При

При

Таким образом для решения системы трехдиагональных уравнений методом прогонки вначале необходимо найти прогоночные коэффициенты – прямой ход прогонки. Для начала счета необходимо принять, что при i = 0 γ_i = 0 и β_i = 0.

Затем выполняется обратный ход прогонки, который заключается в расчете вектора неизвестных u_i с учетом только что рассчитанных прогоночных коэффициентов.

Рассмотрим в качестве примера следующую систему уравнений: