Компактная схема метода Гаусса для трёхдиагональной матрицы и её модификации: различия между версиями
[непроверенная версия] | [непроверенная версия] |
Frolov (обсуждение | вклад) м |
Frolov (обсуждение | вклад) |
||
Строка 12: | Строка 12: | ||
:<math> | :<math> | ||
− | |||
u_{11} = a_{11} \\ | u_{11} = a_{11} \\ | ||
u_{i i+1} = a_{i i+1}, \quad i \in [1, n-1] \\ | u_{i i+1} = a_{i i+1}, \quad i \in [1, n-1] \\ | ||
l_{i+1 i} = a_{i+1 i} / u_{i i} , \quad i \in [1, n-1] \\ | l_{i+1 i} = a_{i+1 i} / u_{i i} , \quad i \in [1, n-1] \\ | ||
u_{ii} = a_{ii} - l_{i i-1} u_{i-1 i}, \quad i \in [2, n] \\ | u_{ii} = a_{ii} - l_{i i-1} u_{i-1 i}, \quad i \in [2, n] \\ | ||
− | |||
</math> | </math> | ||
+ | Подробно её свойства описаны на [[Компактная схема метода Гаусса для трёхдиагональной матрицы, последовательный вариант|соответствующей странице]] | ||
=== <math>LDU</math>-разложение === | === <math>LDU</math>-разложение === |
Версия 14:46, 1 июля 2015
Основные авторы описания: А.В.Фролов
Содержание
1 Компактная схема метода Гаусса для трёхдиагональной матрицы и её модификации
Треугольное разложение трёхдиагональной матрицы базируется на формулах этого же разложения для матриц общего вида. Обычно встречающиеся в прикладных задачах трёхдиагональные матрицы имеют свойство диагонального преобладания. Поэтому нужды в перестановках не возникает. Использование свойств, базирующихся на формуле Бине-Коши для миноров произведения двух матриц, показывает, что для трёхдиагональной матрицы [math]LU[/math]-разложение будет содержать две двухдиагональные матрицы. Поэтому формулы компактной схемы метода Гаусса и её аналогов существенно упрощаются.
1.1 [math]LU[/math]-разложение
1.1.1 Компактная схема метода Гаусса для трёхдиагональной матрицы
Компактная схема метода Гаусса вычисляет такое [math]LU[/math]-разложение, в котором матрица [math]L[/math] имеет на диагонали только единицы. Формулы метода следующие:
- [math] u_{11} = a_{11} \\ u_{i i+1} = a_{i i+1}, \quad i \in [1, n-1] \\ l_{i+1 i} = a_{i+1 i} / u_{i i} , \quad i \in [1, n-1] \\ u_{ii} = a_{ii} - l_{i i-1} u_{i-1 i}, \quad i \in [2, n] \\ [/math]
Подробно её свойства описаны на соответствующей странице
1.2 [math]LDU[/math]-разложение
2 Компактная схема метода Гаусса для трёхдиагональной эрмитовой матрицы и её модификации
2.1 [math]LL^T[/math]-разложение
2.2 [math]LDL^T[/math]-разложение
3 Разложения для трёхдиагональной эрмитовой матрицы специального вида
4 Существующие реализации алгоритма
5 Литература
- Воеводин В.В. Вычислительные основы линейной алгебры. М.: Наука, 1977.
- Воеводин В.В., Кузнецов Ю.А. Матрицы и вычисления. М.: Наука, 1984.