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

Материал из Алговики
Перейти к навигации Перейти к поиску
[непроверенная версия][непроверенная версия]
м
Строка 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 Литература

  1. Воеводин В.В. Вычислительные основы линейной алгебры. М.: Наука, 1977.
  2. Воеводин В.В., Кузнецов Ю.А. Матрицы и вычисления. М.: Наука, 1984.