example3_3.m

来自「Matlab 6.X辅助优化计算与设计----源代码」· M 代码 · 共 22 行

%编写目标函数
function y=objfun(x)
y=exp(x(1))*(3*x(1)^2+2*x(2)^2+3*x(1)*x(2)+2*x(2)+1);

%编写返回约束值的函数
function [c,ceq]=confun(x)
%非线性不等式约束
c=[1.5+x(1)*x(2)-x(1)-x(2);-x(1)*x(2)-10];
%线性等式约束
ceq=[];

%初始点
x0=[-1,1]; 
%设置下界
lb=[0,0]; 
%无上界
ub=[]; 
%采用标准算法
options=optimset ('largescale','off');
[x,fval]=fmincon('objfun',x0,[],[],[],[],lb,ub,'confun',options)
[c,ceq]=confun(x)