opt28_2c.m

来自「《MATLAB实用指南》系列丛书DE源代码」· M 代码 · 共 21 行

function [c, ceq, K1,s] = mycon(x, s)
% 初始化取样区间
if isnan(s(1, 1)),
   s = [2 2];
end
% 样本集
w1x = 1:s(1,1):100;
w1y = 1: s(1, 2):100;
[wx,wy]=meshgrid(w1x,w1y);
%
% 半无限约束
K1=sin(wx*x(1)).*cos(wy*x(2))-1/1000*(wx-50).^2 -...
    sin(wx*x(3))-x(3)+sin(wy*x(2)).*cos(wx*x(1))-...
    1/1000*(wy-50).^2-sin(wy*x(3))-x(3)-1.5;
% 无有限非线性约束
c=[]; ceq=[];
% 网格图
mesh(K1)
title('Semi-infinite constraint')
drawnow