Минимизация функции с использованием множителей Лагранжа в Matlab

Question

У меня проблема с моим кодом MATLAB, который я пишу, чтобы минимизировать эту функцию с двумя ограничениями (одно из них - неравенство, а другое - равенство) с множителями Лагранжа, которые используют условия KKT. Эта функция:

       MIN F = 2*x1^2 + 2*x2^2-6*x1*x2
        ineq=x1+2*x2<=3
        eq=3*x1+x2=9

Я лично подозреваю, что моя if Функция (оператор) имеет проблему, но не знаю, как это исправить .....

    clc
    warning off
    syms x1 x2 u s lambda
    f=2*x1^2+2*x2^2-6*x1*x2;
    g1=sqrt(x1^2*x2)+s^2-0.25;
    H1=x1*x2^2-0.1;
    Lagrange=f+u*g1+lambda*H1;
    Grad=gradient(Lagrange);
    S=solve(Grad);
    S=double([S.x1 S.x2 S.s S.u S.lambda]);M=size(S,1);
      for i=M:-1:1
    if imag(S(i,3))~=0 || S(i,1)<=0 || S(i,2)<=0
      end
    end
    x1=S(:,1);x2=S(:,2);s=S(:,3);u=S(:,4);lambda=S(:,5);
    out=table(x1,x2,s,u,lambda)
    x0=[1 0.5 -0.75];

Кажется, у меня проблема, потому что она дает мне как письменные, так и неправильные ответы. Можете ли вы помочь мне с этим?

Answer 1

Просто удалите соответствующую строку, когда if условие выполнено, используя S(i,:)=[];

---

Согласно вашей формулировке

Условие if должно выглядеть следующим образом

if imag(S(i,3))~=0 || S(i,1)<=0 || S(i,2)<=0
    S(i,:)=[];
end

Решение

out =

  1×5 table

      x1        x2       s      u       lambda 
    ______    _______    _    ______    _______

    0.3393    0.54288    0    3.6349    -2.6402

---

С моей точки зрения, ваши ограничения по равенству и неравенству формулировка неверна

ineq = x1 + 2*x2 <= 3 -------> g1 = x1 + 2*x2 - 3 + s^2

eq = 3*x1 + x2 = 9    -------> H1 = 3*x1 + x2 - 9

Проверка состояния if

• если u ( изменение g1) меньше нуля или нет

или

• , если s - реальное число или нет

for i=M:-1:1
      if imag(S(i,3))~=0 || S(i,4)<0 
      S(i,:)=[];
      end
  end

Весь код выглядит следующим образом

syms x1 x2 u s lambda
    f=2*x1^2+2*x2^2-6*x1*x2;
    g1 = x1+2*x2-3 +s^2 ;
    H1 = 3*x1+x2-9;

    Lagrange=f+u*g1+lambda*H1;
    Grad=gradient(Lagrange);

    S=solve(Grad)
 S=double([S.x1 S.x2 S.s S.u S.lambda]);M=size(S,1)

      for i=M:-1:1
          if imag(S(i,3))~=0 || S(i,4)<0 

             S(i,:)=[];

          end
      end

    x1=S(:,1);x2=S(:,2);s=S(:,3);u=S(:,4);lambda=S(:,5);
    out=table(x1,x2,s,u,lambda)

Решение

out =

  1×5 table

    x1    x2    s     u      lambda
    __    __    _    ____    ______

    3     0     0    13.2     -8.4 

---

Проверьте свой ответ, используя fmincon, как следует

% Check your answer using fmincon 
% Inequality constraint
%ineq = x1 + 2*x2 <= 3
A = [1, 2];
b = 3;
% Equality constraint
%eq = 3*x1 + x2 = 9
Aeq = [3, 1];
beq = 9;
%F = 2*x1^2 + 2*x2^2-6*x1*x2
F = @(x)2*x(1).^2 + 2*x(2).^2-6*x(1).*x(2);

x0 = [0,0];
[x_minimum, Feval] = fmincon(F, x0, A, b, Aeq, beq);

Решение

x_minimum =

    3.0000   -0.0000


Feval =

   18.0002