%首先編寫目標(biāo)函數(shù)的.m文件
function  f = myfun (x, s)
f = sum( (x - 0.5).^2 );

%然后編寫描述非線性等式和不等式約束以及半無窮約束的函數(shù)
function [c, ceq, K1, K2, s] = mycon (X, s)
 if isnan (s(1,1))
   s = (0.2, 0 ;0.2, 0); %初始化取樣步長
 end

 %取樣值
 w1 = 1 : s(1, 1) : 100;
 w2 = 1 : s(2, 1) : 100;
 %半無窮約束
 K1=sin(w1*X(1)).*cos(w1*X(2))-1/1000*(w1-50).^2-sin(w1*X(3))-X(3)-1;
 K2=sin(w2*X(2)).*cos(w2*X(1))-1/1000*(w2-50).^2-sin(w2*X(3))-X(3)-1;
 c = [ ]; ceq = [ ];
 %畫出半無窮約束的圖形
 plot(w1, K1, '-', w2, K2, '+');
 title('semi-infinite constraint');

 %調(diào)用函數(shù)fseminf
 x0 = [0.5; 0.2; 0.3]; %初值
 [x, fval] = fseminf (@myfun, x0, 2, @mycon);