Сформируйте дополненную систему методом наименьших квадратов
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.