matlab既約梯度法習題

既約梯度法的約束條件爲線性，目標函數可爲非線性。

驗證：

oldFmt = format("rat");
f=@(x) x(1)^2+x(2)^2-x(1)*x(2)-2*x(1)-3*x(2);
A = [1,1
    1,5];
b = [2,5];
Aeq=[];
beq=[];
lb=[0,0];
ub=[];
[x,fval] = fmincon(f,[0,0],A,b,Aeq,beq,lb,ub)
format(oldFmt);

x =

35/31 24/31

fval =

-111/31