example5_14.m

来自「MATLAB7辅助控制系统设计与仿真 第5章」· M 代码 · 共 16 行

function [sol,eval]=f552(sol,options)
x1=sol(1);
x2=sol(2);
r1=0.1;
r2=0.8;
%约束条件
g1=x1-2*x2+1;
g2=x1.^2/4-x2.^2+1;
%加惩罚项的适值
if (g1>=0)&(g2>=0)
  eval=(x1-2).^2+(x2-1).^2;
else
  eval=(x1-2).^2+(x2-1).^2+r1*g1+r2*g2;
  eval=-eval;
end