МатВектор

Command Palette

Search for a command to run...

🔢 Численные методы

Число обусловленности матрицы

англ. Matrix condition number

Коэффициент усиления ошибки матрицей: κ(A) = ||A||·||A⁻¹|| показывает, во сколько раз относительная погрешность данных может вырасти в решении системы Ax = b.

Микрофон, который усиливает всё подряд, превращает шёпот докладчика в гул: сигнал слышен, но вместе с ним — стулья, вентиляция, шаги в коридоре. Матрица с большим числом обусловленности работает так же. Для системы полагают : во сколько раз относительная ошибка правой части может вырасти в решении. У симметричной матрицы в спектральной норме это отношение крайних собственных значений . Всегда ; идеал — единичная матрица, она ничего не усиливает. Общая постановка без матриц — в термине обусловленность задачи.

Живой пример: матрица со строками и имеет около . Сдвиньте второй элемент правой части с до — относительное возмущение в пять стотысячных — и решение переползёт с на : реакция в тысячи раз сильнее причины. Гильбертова матрица даёт — потеря около верных знаков из шестнадцати машинных. Правило счёта: съедает знаков. Крайние собственные значения для спектральной формы ищут численно — в уроке собственные значения: численные методы, сами системы — в материале численные методы СЛАУ.

граница усиления ошибки; для матрицы со строками (1, 1) и (1, 1.0001) число : сдвиг правой части на меняет решение с (2, 0) на (1, 1)

Частые вопросы

Почему нельзя просто пересчитать точнее?

Потому что знаки теряются при самом умножении на плохо обусловленную матрицу, а не из-за ленивого процессора. При машинной точности около и числе в решении остаётся верных знаков: для от двойной точности остаются единицы, и никакой пересчёт их не вернёт — информация о решении физически не помещается в разрядную сетку. Лечат не арифметикой, а постановкой: масштабирование, регуляризация, переформулировка задачи.

Как оценить κ без обращения матрицы?

Для симметричной матрицы хватит крайних собственных значений: отношение — это и есть в спектральной норме, а крайние значения находят степенным методом и обратной итерацией. Для общей матрицы — через сингулярные числа. Косвенный сигнал — маленький ведущий элемент при разложении Гаусса: он почти всегда сопровождает плохую обусловленность. Практическая проверка: решить систему дважды с чуть возмущённой правой частью и посмотреть на разницу решений.

пошагово

Решатель обратной матрицы

Найти обратную матрицу для κ(A)

Число обусловленности равно ‖A‖·‖A⁻¹‖: решатель посчитает обратную матрицу через алгебраические дополнения.

Открыть решатель