Parameters x1[0.56,0.8],x2[0.001,0.031],x3[600,1100];
MaxFunction    -1793444.01*x1*x1*x1-0.9269*x3*x3-10308400*x2*x2+3746080.65*x1*x1+333.10356*x3+877480.696*x2-2158343.84*x1+7.1376*x1*x3+391.7883*x2*x3-982145*x1*x2+416476.113;

2楼: Originally posted by 月只蓝 at 2013-10-22 15:03:01
1stopt软件计算代码和结果：
Parameters x1,x2,x3;
MaxFunction    -1793444.01*x1*x1*x1-0.9269*x3*x3-10308400*x2*x2+3746080.65*x1*x1+333.10356*x3+877480.696*x2-2158343.84*x1+7.1376*x1*x3+391.7883*x2*x3 ...

2楼: Originally posted by 月只蓝 at 2013-10-22 15:03:01
1stopt软件计算代码和结果：
Parameters x1,x2,x3;
MaxFunction    -1793444.01*x1*x1*x1-0.9269*x3*x3-10308400*x2*x2+3746080.65*x1*x1+333.10356*x3+877480.696*x2-2158343.84*x1+7.1376*x1*x3+391.7883*x2*x3 ...

function max_Fmincon
format short

clear all
clc

x0 = [0.6; 0.01;800];

[x,fval] = fmincon(@ObjFunc,x0,[],[],[],[],[],[],@NlinCons)
% ------------------------------------------------------------------
ff=-ObjFunc(x)

function f = ObjFunc(x)
x1=x(1);
x2=x(2);
x3=x(3);
f = -(0.00038596*x3*x3*x3-1793444.01*x1*x1*x1-...
    0.9269*x3*x3-10308400*x2*x2+3746080.65*x1*x1+...
    333.10356*x3+877480.696*x2-2158343.84*x1+...
    7.1376*x1*x3+391.7883*x2*x3-982145*x1*x2+416476.113);
% ------------------------------------------------------------------

function [c,ceq] = NlinCons(x)
c(1) = x(1)-0.8;
c(2)= 0.56-x(1);
c(3) = x(2)-0.031;
c(4)= 0.001-x(2);
c(5) = x(3)-1100;
c(6)= 600-x(3);
ceq=[];