GADS命令方式求解非线性约束问题

求解问题

$$\begin{gathered} max2.2x_1+3x_2+4.5x_3+3x_1^2+x_3^2; \\ 4x_1^2+2x_2^2+x_3-46\le 0; \\ {x}_1+x_1^2+x_2^2-5x_3-34\le 0; \\ 3x_1+2x_1^2+x_2+3x_3-67\le 0; \\ {x}_1+x_3^2=13; \\ {x}_1\le 0;x_2\le 10; \end{gathered}$$

主程序

clc; clear;
myFitnessFun=@(x) -(2.2*x(1)+3*x(2)+4.5*x(3)+3*x(1)^2+x(3)^2);
A=[-2 -4 -6;-1 0 0;0 1 0];
b=[-10;0;10];
Aeq=[ ];
beq=[ ];
lb=[];
ub=[];
init_pop=repmat([3 1 3],100,1);
options=gaoptimset('InitialPopulation',init_pop,'CreationFcn',@gacreationlinearfeasible, ...
                   'PopulationType','doubleVector' ,'PopulationSize',100,'EliteCount',1,...
                   'CrossoverFraction',0.75,'Generations',100,'StallGenLimit',500,...
                   'TolFun',1e-100,'PlotFcns',{@gaplotbestf,@gaplotbestindiv}); %参数设置t;
[x_best,fval,reason,output,final_pop]=ga(myFitnessFun,3,A,b,Aeq,beq,ub,lb,@myNolinCon,options);
display(x_best);
fprintf('最优适应度值:%.3f\n',-fval);

myNolinCon函数

%ex4_2非线性约束
function [nlneq, nleq]=myNolinCon(x)
nlneq=[4*x(1)^2+2*x(2)^2+x(3)-46; 
       x(1)+x(1)^2+x(2)^2-5*x(3)-34; 
       3*x(1)+2*x(1)^2+x(2)+3*x(3)-67];
nleq=x(1)+x(3)^2-13;
end

结果

x_best =

3.1546    1.2363    3.1377

最优适应度值:64.469