Сформируйтесь наименьшие квадраты увеличили систему
S = spaugment(A,c)
S = spaugment(A)
S = spaugment(A,c)
создает разреженный, квадратный, симметричный неопределенный матричный S = [c*I A; A' 0]
. Матричный S
связан с задачей наименьших квадратов
r = b - A*x S * [r/c; x] = [b; 0]
Оптимальное значение остаточного масштабного коэффициента c
, включает min(svd(A))
и norm(r)
, которые являются обычно слишком дорогими, чтобы вычислить.
S = spaugment(A)
без заданного значения c
, использование max(max(abs(A)))/1000
.
Примечание
В предыдущих версиях MATLAB®, расширенная матрица использовалась разреженными решателями линейного уравнения, \и/, для неквадратных проблем. Теперь MATLAB вычисляет решение методом наименьших квадратов с помощью qr
факторизация A
.