Число обусловленности матрицы
англ. Matrix condition number
Коэффициент усиления ошибки матрицей: κ(A) = ||A||·||A⁻¹|| показывает, во сколько раз относительная погрешность данных может вырасти в решении системы Ax = b.
Микрофон, который усиливает всё подряд, превращает шёпот докладчика в гул: сигнал слышен, но вместе с ним — стулья, вентиляция, шаги в коридоре. Матрица с большим числом обусловленности работает так же. Для системы полагают : во сколько раз относительная ошибка правой части может вырасти в решении. У симметричной матрицы в спектральной норме это отношение крайних собственных значений . Всегда ; идеал — единичная матрица, она ничего не усиливает. Общая постановка без матриц — в термине обусловленность задачи.
Живой пример: матрица со строками и имеет около . Сдвиньте второй элемент правой части с до — относительное возмущение в пять стотысячных — и решение переползёт с на : реакция в тысячи раз сильнее причины. Гильбертова матрица даёт — потеря около верных знаков из шестнадцати машинных. Правило счёта: съедает знаков. Крайние собственные значения для спектральной формы ищут численно — в уроке собственные значения: численные методы, сами системы — в материале численные методы СЛАУ.
Частые вопросы
Почему нельзя просто пересчитать точнее?
Потому что знаки теряются при самом умножении на плохо обусловленную матрицу, а не из-за ленивого процессора. При машинной точности около и числе в решении остаётся верных знаков: для от двойной точности остаются единицы, и никакой пересчёт их не вернёт — информация о решении физически не помещается в разрядную сетку. Лечат не арифметикой, а постановкой: масштабирование, регуляризация, переформулировка задачи.
Как оценить κ без обращения матрицы?
Для симметричной матрицы хватит крайних собственных значений: отношение — это и есть в спектральной норме, а крайние значения находят степенным методом и обратной итерацией. Для общей матрицы — через сингулярные числа. Косвенный сигнал — маленький ведущий элемент при разложении Гаусса: он почти всегда сопровождает плохую обусловленность. Практическая проверка: решить систему дважды с чуть возмущённой правой частью и посмотреть на разницу решений.
Решатель обратной матрицы
Найти обратную матрицу для κ(A)
Число обусловленности равно ‖A‖·‖A⁻¹‖: решатель посчитает обратную матрицу через алгебраические дополнения.
Открыть решатель