function [f,g]=devil(x)f=x(1)^2+x(2)^2;g(1)=x(1)-206.9014;g(2)=x(2)-0.7*x(1);g(3)=0.53*x(1)-x(2);g(4)=x(1)^2+x(1)*x(2)+x(2)^2-0.28*(x(1)+x(2)); x0=[0.1 0.06];opt(1)=1;>> x=constr('devil',x0,opt)不过你的约束条件有问题没有这样D,d同时满足你的三个条件吧