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

Метод прогонки является модификацией метода Гаусса для частного случая систем уравнений с трехдиагональной матрицей вида:

Векторы:

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

Из 1-го уравнения системы (3.4) выражаем x₁:

где

Из 2-го уравнения системы (3.4) выражаем x₂:

Продолжая этот процесс при условии, что U₀=0, V₀=0, получим рекуррентную формулу для x_i:

где

После прямого хода система (3.4) примет вид:

Обратная прогонка состоит в последовательном вычислении неизвестных x_i. Сначала вычисляется неизвестное x_n из двух последних уравнений преобразованной системы (3.9):

Затем, используя выражение (3.7), последовательно вычисляются все остальные неизвестные

Алгоритм метода прогонки

1. Вычисляются прогоночные коэффициенты U_i, V_i по формулам (3.8) для i=1,2,…n-1 при условии, что U₀=0, V₀=0;

2. Определяются x_n из выражения (3.10);

3. По формуле (3.7) последовательно вычисляются все остальные неизвестные

3.3. Итерационные методы решения систем линейных алгебраических уравнений

В этой группе методов мы познакомимся с двумя старыми и простыми методами: методом Якоби и методом Гаусса – Зейделя [8, 14].

3.3.1. Метод Якоби (простых итераций)

Задана система линейных алгебраических уравнений

или в матричной форме

Для сходимости итерационного процесса, необходимо выполнения условия «преобладания диагональных элементов», т.е. диагональные элементы матрицы А должны удовлетворять условию:

Преобразуем систему (3.10) к эквивалентной, выражая неизвестное из каждого i-го уравнения:

Система (3.13) называется системой, приведенной к нормальному виду.

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

систему (3.13) можно записать в матричной форме

где

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

Таким образом, получили последовательность приближений (итераций):

Если эта последовательность имеет предел

то он является точным решением системы (3.11).

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

Критерий близости двух приближений может быть определен следующим образом:

Если условие (3.18) выполнено, то итерационный процесс прекращается. За приближенное решение системы (3.11) с заданной точностью  принимается (k)-е приближение, т.е.

.

Если условие (3.18) не выполняется, то итерационный процесс (3.17) необходимо продолжить до тех пор, пока условие не выполнится.

• За м е ч а н и я:

1. Начальный вектор может быть взят произвольным, так как сходимость итерационного процесса зависит только от свойств матрицы , и если процесс сходится при каком-нибудь начальном приближении, то он будет сходиться к тому же предельному вектору и при любом другом выборе этого начального приближения.

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

• Пример 3.2. Методом Якоби решить систему линейных алгебраических уравнений:

Решение: Условие преобладания диагональных коэффициентов матрицы системы выполнено. Приведем эту систему к нормальному виду:

В матричной форме систему (3.20) можно записать так:

За начальное (нулевое) приближение решения системы примем нулевой вектор, т.е.

Подставляя эти значения в правые части уравнения (3.20), получим первое приближение решения системы (первую итерацию):

Далее, подставляя это найденное приближение в систему (3.20), получим 2-ое приближение решения системы:

После новой подстановки будем иметь 3-е приближение:

Аналогично получим 4-ую итерацию:

Проверим выполнение условия «близости» двух итераций, т.е. условие (3.18):

Таким образом, за приближенное решение системы (3.19) с точностью ε=0,1 принимаем 4-ю итерацию

Чтобы получить решение СЛАУ (3.19) с точностью ε=0,001 , потребуется 8 итераций. Точное решение: х₁=1; х₂=2; х₃=1.

Решение данного примера с использованием электронных таблиц Excel приведено в разделе 3.6.3.