Решить Ax = b для решения, которое минимизирует норму L2 x - y.
При заданном yy это означает минимизацию нормы L2 x - y, на которую наложено ограничение Ax = b.
По множителю Лагранжа я получаю следующее уравнение блочной матрицы.
Вектор u является половиной множителей Лагранжа.
/ A 0 \ / x \ / b \
| | | | = | |
\ I A^T / \ u / \ y /
Вероятно, не самый быстрый способ.