Метод декомпозиции
Сначала исходная матрица A раскладывается на две треугольные матрицы B и C таким образом, что A = BC. Формулы для получения элементов матриц B и C:
Диагональные элементы матрицы C равны 1, остальные элементы матриц B и C нулевые:
Важен порядок вычисления элементов матриц B и C. Сначала вычисляется первый столбец матрицы B, затем первая строка матрицы C, затем второй столбец B, затем вторая строка C и т.д. После этого сначала решается СЛАУ By = d, а затем – СЛАУ Cx = y. По аналогии с (2.2.8), для решения этих систем можно записать
Определитель исходной матрицы A можно вычислить по формуле
Метод декомпозиции обладает тем же недостатком, что и метод Гаусса. В формуле (2.2.11) происходит деление на диагональные элементы матрицы B. Если в процессе решения требуемый диагональный элемент получится равным нулю, то этот метод также даст сбой. Аналогично, тогда может помочь только перестановка строк исходной СЛАУ, но делать этого, в рамках данной практической работы, мы не будем.
|
(2.2.10)
(2.2.11)
(2.2.12)
(2.2.13)
(2.2.14)
