Найдем неотрицательное решение линейной задачи методом наименьших квадратов и сравним результат с решением задачи без ограничений.

Подготовьте C матрица и d вектор для проблемы $$\min ||Cx-d||$$.

C = [0.0372    0.2869
     0.6861    0.7071
     0.6233    0.6245
     0.6344    0.6170];
 
d = [0.8587
     0.1781
     0.0747
     0.8405];

Вычислите ограниченные и неограниченные решения.

x = lsqnonneg(C,d)

x = 2×1

         0
    0.6929

xunc = C\d

xunc = 2×1

   -2.5627
    3.1108

Все записи в x являются неотрицательными, но некоторые записи в xunc отрицательны.

Вычислите нормы остаточных значений для этих двух решений.

constrained_norm = norm(C*x - d)

constrained_norm = 0.9118

unconstrained_norm = norm(C*xunc - d)

unconstrained_norm = 0.6674

Неограниченное решение имеет меньшую норму невязки, потому что ограничения могут только увеличить норму невязки.

Установите Display опция к 'final' видеть выход когда lsqnonneg концы.

Создайте опции.

options = optimset('Display','final');

Подготовьте C матрица и d вектор для проблемы $$\min ||Cx-d||$$.

C = [0.0372    0.2869
     0.6861    0.7071
     0.6233    0.6245
     0.6344    0.6170];

d = [0.8587
     0.1781
     0.0747
     0.8405];

Вызовите lsqnonneg со структурой опций.

x = lsqnonneg(C,d,options);

Optimization terminated.

Вызовите lsqnonneg с выходными параметрами, чтобы получить решение, норму невязки и вектор невязок.

Подготовьте C матрица и d вектор для проблемы $$\min ||Cx-d||$$.

C = [0.0372    0.2869
     0.6861    0.7071
     0.6233    0.6245
     0.6344    0.6170];

d = [0.8587
     0.1781
     0.0747
     0.8405];

Получите решение и остаточную информацию.

 [x,resnorm,residual] = lsqnonneg(C,d)

x = 2×1

         0
    0.6929

resnorm = 0.8315

residual = 4×1

    0.6599
   -0.3119
   -0.3580
    0.4130

Проверьте, что возвращенная норма невязки является квадратом нормы возвращенного вектора невязок.

 norm(residual)^2

ans = 0.8315

Запросите все выходные аргументы исследовать процесс решения и решения после lsqnonneg концы.

Подготовьте C матрица и d вектор для проблемы $$\min ||Cx-d||$$.

C = [0.0372    0.2869
     0.6861    0.7071
     0.6233    0.6245
     0.6344    0.6170];

d = [0.8587
     0.1781
     0.0747
     0.8405];

Решите задачу, запросив все выходные аргументы.

[x,resnorm,residual,exitflag,output,lambda] = lsqnonneg(C,d)

x = 2×1

         0
    0.6929

resnorm = 0.8315

residual = 4×1

    0.6599
   -0.3119
   -0.3580
    0.4130

exitflag = 1

output = struct with fields:
    iterations: 1
     algorithm: 'active-set'
       message: 'Optimization terminated.'

lambda = 2×1

   -0.1506
   -0.0000

exitflag 1, указание на правильное решение.

x(1) = 0, и соответствующий lambda(1) $$\ne$$0 , правильно указывание на двойственность. Точно так же x(2) > 0, и соответствующий lambda(2) = 0.