Y=@(X) -22.92*X(1)+52.32*X(2)+10.90*X(3)+2.62*X(1)*X(2)-11.41*X(1)*X(3)+17.93*X(2)*X(3)+1.58*X(1)^2-28.57*X(2)^2+18.11*X(3)^2+220.65;
x0=[1 1 1];lb=-1.682*[1 1 1];vb=1.682*[1 1 1];options = optimset('algorithm','active-set');
[X Y]=fmincon(Y,x0,[],[],[],[],lb,vb)
按上面语句试试