Курсовая работа: Численные методы
что эквивалентно умножению (13) на элементарную нижнюю треугольную матрицу
Таким образом, для рассмотренного примера процесс исключения Гаусса с выбором главного элемента по столбцу записывается в
виде
(15)
По построению матрица
(16)
является верхней треугольной матрицей с единичной главной диагональю.
Отличие от обычного метода Гаусса состоит в том, что в качестве сомножителей в (16) наряду с элементарными треугольными матрицами могут присутствовать элементарные матрицы перестановок .
Покажем еще, что из (16) следует разложение
PA=LU , (17)
где L - нижняя треугольная матрица, имеющая обратную, P - матрица перестановок.
Для этого найдем матрицу
(18)
По свойству 2) матрица получается из матрицы перестановкой второй и третьей строк,
Матрица согласно свойству 3) получается из перестановкой второго и третьего столбцов
т.е. -нижняя треугольная матрица, имеющая обратную.
Из (18), учитывая равенство , получим
(19)
Отсюда и из (16) видно, что
где обозначено . Поскольку Р -матрица перестановок и L -нижняя треугольная матрица, свойство (17) доказано. Оно означает, что метод Гаусса с выбором главного элемента по столбцу эквивалентен обычному методу Гаусса, примененному к матрице РА , т.е. к системе, полученной из исходной системы перестановкой некоторых уравнений.
3. Общий вывод. Результат, полученный ранее для очень частного примера, справедлив и в случае общей системы уравнений (1).
А именно, метод Гаусса с выбором главного элемента по столбцу можно записать в виде
(20)
где - элементарные матрицы перестановок такие, что
и -элементарные нижние треугольные матрицы.