Я пытаюсь решить линейную систему, которая имеет несколько решений.Вот часть моего кода:
from scipy.sparse.linalg import lsqr
solution = lsqr(M, b)[0]
Теперь, если матрица M такая:
[1 1 1 1]
[0 0 0 0]
[0 0 0 0]
[0 0 0 0]
[0 0 0 0]
[0 0 0 0]
и b это:
[1 0 0 0 0 0]
Решениеэтот код задается следующим образом:
[1/4 1/4 1/4 1/4]
Однако для моих целей я хотел бы получить в качестве решения вектор с максимально возможным числом нулей, поэтому в этом случае это будет так:
[1 0 0 0]
Есть ли способ сделать это?Я в порядке, используя пакеты, кроме scipy.sparse.linalg тоже.Спасибо!