Курсовая работа: Решение параболических уравнений

.

Эти формулы имеет погрешность . В результате уравнение (1.10) заменяется разностным:

(1.13)

Перепишем (1.13) в виде:

. (1.14)

Данная вычислительная схема имеет следующую конфигурацию:

(1.15)

(1.16)

Система (1.14) – (1.16) представляет собой разностную задачу, соответствующую краевой задаче (1.10) – (1.12).

За величину мы положили .

(1.14) – (1.16) есть система линейных алгебраических уравнений с 3-диагональной матрицей, поэтому ее резонно решать методом прогонки, так как он в несколько раз превосходит по скорости метод Гаусса.

. (1.17)

Здесь , – некоторые коэффициенты, подлежащие определению. Заменив в (1.17) на будем иметь:

. (1.18)

Подставив уравнение (1.18) в (1.14) получим:

. (1.19)

Сравнив (1.17) и (1.19) найдем, что:

(1.20)

Положим в (1.14) и найдем из него :

,

.

(1.21)

Заметим, что во второй формуле (1.21) величина подлежит замене на согласно первому условию (1.15).

С помощью формул (1.21) и (1.20) проводим прогонку в прямом направлении. В результате находим величины

Затем осуществляем обратный ход. При этом воспользуемся второй из формул (1.15) и формулой (1.17). Получим следующую цепочку формул:

(1.22)

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