Constant
n=18,
S=[0.48,1.13,2.72,0.51,1.28,2.49,0.59,0.94,0.36,0.51,1.26,0.44,0.51,0.93,1.23,1.82,0.43,1.09],
A=[10.14,9.64,10.41,8.8,9.55,9.59,9.21,10.07,10.05,9.69,8.95,9.85,10.76,9.49,10.58,11.58,9.08,9.97],
V=[61.6,64.1,60.9,62.1,64.2,57.7,60.6,63.3,66.1,62.8,64.6,69.3,59.4,60.5,64.6,70.9,68.9,60.6],
aa=[4.45,3.65,5.8,2.85,4.7,7.1,6.85,4.7,4.45,4.7,6.85,5.8,5.8,3.85,7,4.45,2.85,7],
b=[2.9,2.5,4.5,2.8,3.5,10.9,6,3.5,2.9,3.5,6,4.5,4.5,2.5,8,2.9,2.8,8],
j=[1100,1007,1005,1061,1025,805,1245,1125,1000,1044,1087,1100,1241,1060,1095,1015,1045,1087];
ConstStr sumx=Sum(i=1:n)(x[i]), sumxa=Sum(i=1:n)(aa[i]*x[i]), sumxb=Sum(i=1:n)(x[i]/b[i]),
         SI=Sumxa/Sumx, CBI=(1-Sumx)/Sumxb;
Parameter x(n)=[0,0.2];
MinFunction Sum(i=1:n)(j[i]*x[i]);
            2.603+1.024*Sum(i=1:n)(A[i]*x[i]) <12.8;
            0.883*Sum(i=1:n)(S[i]*x[i]) <0.95;
            if(CBI>=1,3.601+0.236*CBI^2, if(CBI<=0.8, 5.58-1.65*CBI, 4.049+0.115*CBI^2))<=SI;
            sumx=1;

2楼: Originally posted by dingd at 2014-01-19 21:28:32
个人觉得你的模型有些问题：
1：约束条件：SI≥5.58-1.65CBI              ( CBI≤1)
                    SI≥3.601+0.236CBI2                 (CBI≥1)
                    SI≥4.049+0.115CBI2             ...

2楼: Originally posted by dingd at 2014-01-19 21:28:32
个人觉得你的模型有些问题：
1：约束条件：SI≥5.58-1.65CBI              ( CBI≤1)
                    SI≥3.601+0.236CBI2                 (CBI≥1)
                    SI≥4.049+0.115CBI2             ...