Разложение Холецкого

Материал из Википедии — свободной энциклопедии

(Перенаправлено с Теорема Холецкого)
Перейти к: навигация, поиск

Разложение Холецкого — представление симметричной положительно-определённой матрицы A в виде A = LLT (или A = UTU, U = LT), где Lнижняя треугольная матрица со строго положительными элементами на диагонали. Разложение Холецкого всегда существует и единственно для таких матриц A.

Этот вид разложения может применяться для численного решения системы линейных уравнений Ax = b, если матрица A симметрична и положительно-определена. Выполнив разложение A = LLT, решение x получается последовательным решением двух треугольных систем уравнений: Ly = b и LTx = y. Такой способ решения иногда называется методом квадратных корней. По сравнению с более общими методами, такими как метод Гаусса или LU-разложение, он устойчивее численно и требует примерно вдвое меньше арифметических операций.

Существует также обобщение этого разложения на случай комплекснозначных матриц. Если A — положительно-определённая эрмитова матрица, то существует разложение A = LL * , где L — нижняя треугольная матрица, а L * — её эрмитово дополнение.

[править] Формулы

Элементы матрицы L можно вычислить, начиная с верхнего левого угла матрицы, по формулам:

\textstyle L_{ii} = \sqrt{A_{ii} - \sum_{k=1}^{i-1} L_{ik}^2}

\textstyle L_{ij} = \frac{1}{L_{jj}} \left(A_{ij} - \sum_{k=1}^{j-1} L_{ik} L_{jk} \right), если j < i.

[править] Реализация в математических пакетах программ

В системах MATLAB и Octave разложение выполняется командой U = chol(A).

[править] Ссылки