Лабораторная работа № 1 решение систем линейных алгебраических уравнений методом гаусса с выбором главного элемента

ЦЕЛЬ РАБОТЫ

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

ТЕОРЕТИЧЕСКИЕ СВЕДЕНИЯ

Метод Гаусса для произвольной системы линейных алгебраических уравнений

основан на приведении матрицы системы к треугольной. Вычтем из второго уравнения системы (1) первое, умноженное на такое число, чтобы уничтожился коэффициент при . Затем таким же образом вычтем первое уравнение из третьего, четвертого и т.д. Тогда исключатся все коэффициенты первого столбца, лежащие ниже главной диагонали. Затем при помощи второго уравнения исключим из третьего, четвертого и т.д. уравнений коэффициенты второго столбца. Последовательно продолжая этот процесс, исключим из матрицы все коэффициенты, лежащие ниже главной диагонали.

Запишем общие формулы процесса. Пусть проведено исключение коэффициентов из столбца. Тогда остались такие уравнения с ненулевыми коэффициентами ниже главной диагонали:

Умножим -ю строку на число

и вычтем из -той строки. Первый элемент этой строки обратится в нуль, а остальные изменятся по формулам

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

После выполнения прямого хода получим треугольную систему

с матрицей

Треугольная система (5) легко решается обратным ходом по формулам

Замечания:

1. Исключение по формулам (3)-(4) нельзя проводить, если в ходе расчета на главной диагонали оказался нулевой элемент . Тогда в промежуточной системе (2) перестановкой строк необходимо переместить ненулевой элемент на главную диагональ и продолжить расчет.

2. Если элемент на главной диагонали мал, то эта строка

умножается на большие числа , что приводит к значительным ошибкам при вычитаниях. Чтобы избежать этого, каждый цикл всегда начинают с перестановки строк. Среди элементов столбца,, находят главный, т.е. наибольший по модулю в-том столбце, и перестановкой строк переводят его на главную диагональ, после чего выполняют исключения. В методе Гаусса с выбором главного элемента погрешность округления обычно невелика. Только для плохо обусловленных систем () устойчивость этого метода оказывается недостаточной.

3. Для контроля расчета полезно найти невязки:

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