请问1stopt最优值的精度问题

stextinction

借问
我用1stopt计算以下的问题的最大值
因为是按照外国的论文算个结果啊
所以知道答案的
然后我下的是1.5pro
选择不同的算法计算
最大也只是算到了0.4左右
差了许多
这个问题是模拟移动床的问题
有四个核心的非线性方程
我把这些方程都当作约束了
现在算出来约束倒是满足的不错
但是最大值有点小了
不知道是不是因为没有用2.0专业版？
请问
2.0对1.5的优势大么？
我考虑如果能帮我算出最大值
我就去买个正版的用用
可是我就害怕买了也算不出来
请问有没有大虾能帮我运行下？
要是结果好
我就买个正版的去
穷学生
买个正版一个半月的生活费就没有啦
所以问问能不能谁帮我算算看看准不准
所以在此求助
谢谢！


问题如下


Title "smb" ;
Parameter Cf[5,40],dc[40,80], dp[5,100],Lc[10,100,0],Nc1[2,14,0],Nc2[2,14,0],Nc3[2,14,0],
Nctot=8[8,20,0],v[0.02,150],u01[0.02,200], u02[0.02,200],u03[0.02,200],u04[0.02,200], Cs01[0.02,200];

Constant Y1=0.99,Y2=0.99,Pur_reqE=0.99,Pur_reqR=0.99,epsb=0.24,ep=0.6,mu=0.318/60,rou=0.79,of_bed=600,
Dinf1=6.12*0.0001,Dinf2=6.12*0.0001,Dop1=2*0.00001,Dop2=2*0.00001,a1=1.12,b1=0.0168,a2=2.12,b2=0.0383,
DV_rate=5.73*0.01,ts_min=0.5,deltaPmax=2000,P=(1-epsb)/epsb;

Conststr
Cs02=((Cs01*a2*b1*(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)+Cs01*a1*b2-(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)*a1+(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)*a2)/((((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)*(a2*b1*(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)+a1*b2))),
gamaplusf=(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf),
Cp2=((a2-a1)/(a1*b2+a2*b1*(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf))),
gamaminus= (((-a1*(1+b2*((Cs01*a2*b1*(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)+Cs01*a1*b2-(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)*a1+(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)*a2)/((((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)*(a2*b1*(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)+a1*b2))))+a2*(1+b1*Cs01))-((-a1*(1+b2*((Cs01*a2*b1*(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)+Cs01*a1*b2-(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)*a1+(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)*a2)/((((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)*(a2*b1*(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)+a1*b2))))+a2*(1+b1*Cs01))^2+(4*a1*a2*b1*b2*Cs01*((Cs01*a2*b1*(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)+Cs01*a1*b2-(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)*a1+(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)*a2)/((((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)*(a2*b1*(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)+a1*b2)))))^0.5)/2/a2/b1/((Cs01*a2*b1*(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)+Cs01*a1*b2-(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)*a1+(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)*a2)/((((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)*(a2*b1*(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)+a1*b2)))),
Cp1=((a1-a2)/(a2*b1+a1*b2/(((-a1*(1+b2*((Cs01*a2*b1*(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)+Cs01*a1*b2-(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)*a1+(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)*a2)/((((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)*(a2*b1*(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)+a1*b2))))+a2*(1+b1*Cs01))-((-a1*(1+b2*((Cs01*a2*b1*(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)+Cs01*a1*b2-(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)*a1+(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)*a2)/((((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)*(a2*b1*(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)+a1*b2))))+a2*(1+b1*Cs01))^2+(4*a1*a2*b1*b2*Cs01*((Cs01*a2*b1*(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)+Cs01*a1*b2-(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)*a1+(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)*a2)/((((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)*(a2*b1*(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)+a1*b2)))))^0.5)/2/a2/b1/((Cs01*a2*b1*(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)+Cs01*a1*b2-(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)*a1+(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)*a2)/((((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)*(a2*b1*(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)+a1*b2)))))),
delta1=(ep+(1-ep)*a2+DV_rate/(1-epsb)),delta2=(ep+(1-ep)*(a1/(1+b2*((a2-a1)/(a1*b2+a2*b1*(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)))))+DV_rate/(1-epsb)),
delta3=(ep+(1-ep)*(a2/(1+b1*Cs01+b2*((Cs01*a2*b1*(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)+Cs01*a1*b2-(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)*a1+(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)*a2)/((((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)*(a2*b1*(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)+a1*b2)))))+DV_rate/(1-epsb)),delta4=(ep+(1-ep)*(a1/(1+b1*((a1-a2)/(a2*b1+a1*b2/(((-a1*(1+b2*((Cs01*a2*b1*(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)+Cs01*a1*b2-(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)*a1+(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)*a2)/((((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)*(a2*b1*(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)+a1*b2))))+a2*(1+b1*Cs01))-((-a1*(1+b2*((Cs01*a2*b1*(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)+Cs01*a1*b2-(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)*a1+(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)*a2)/((((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)*(a2*b1*(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)+a1*b2))))+a2*(1+b1*Cs01))^2+(4*a1*a2*b1*b2*Cs01*((Cs01*a2*b1*(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)+Cs01*a1*b2-(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)*a1+(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)*a2)/((((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)*(a2*b1*(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)+a1*b2)))))^0.5)/2/a2/b1/((Cs01*a2*b1*(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)+Cs01*a1*b2-(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)*a1+(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)*a2)/((((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)*(a2*b1*(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)+a1*b2))))))))+DV_rate/(1-epsb)),
Eb1=(u01*epsb*dp*0.0001/(0.2+0.011*(rou*u01*epsb*dp*0.0001/60/mu)^0.48)),Eb2=(u02*epsb*dp*0.0001/(0.2+0.011*(rou*u02*epsb*dp*0.0001/60/mu)^0.48)),Eb3=(u03*epsb*dp*0.0001/(0.2+0.011*(rou*u03*epsb*dp*0.0001/60/mu)^0.48)),Eb4=(u04*epsb*dp*0.0001/(0.2+0.011*(rou*u04*epsb*dp*0.0001/60/mu)^0.48)),
kl1=(1.09*(u01)^(1/3)*(epsb*dp*0.0001/Dinf2)^(-2/3)),kl2=(1.09*(u02)^(1/3)*(epsb*dp*0.0001/Dinf1)^(-2/3)),
kl3=(1.09*(u03)^(1/3)*(epsb*dp*0.0001/Dinf1)^(-2/3)), kl4=(1.09*(u04)^(1/3)*(epsb*dp*0.0001/Dinf1)^(-2/3)),
ku1=(1/((dp*0.0001)^2/60/ep/Dop2+(dp*0.0001/6/(1.09*(u01)^(1/3)*(epsb*dp*0.0001/Dinf2)^(-2/3))))),
ku2=(1/((dp*0.0001)^2/60/ep/Dop1+(dp*0.0001/6/(1.09*(u02)^(1/3)*(epsb*dp*0.0001/Dinf1)^(-2/3))))),
ku3=(1/((dp*0.0001)^2/60/ep/Dop2+(dp*0.0001/6/(1.09*(u03)^(1/3)*(epsb*dp*0.0001/Dinf1)^(-2/3))))),
ku4=(1/((dp*0.0001)^2/60/ep/Dop1+(dp*0.0001/6/(1.09*(u04)^(1/3)*(epsb*dp*0.0001/Dinf1)^(-2/3))))),
beta1=(ln(u01*(u03-u04)*Y2/u04/(u01-u02)/(1-Y2))), beta2=(ln((u01-u02)*(u03*Cs01-(u03-u02)*Cf)/u02/(u03-u02)/Cf/(1-Y1))),
beta3=(ln(((Cs01*a2*b1*(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)+Cs01*a1*b2-(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)*a1+(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)*a2)/((((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)*(a2*b1*(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)+a1*b2)))*(u03-u04)/(u03-u02)/Cf/(1-Y2))), beta4=(ln (Y1*u04*(u01-u02)/(u03-u04)/u01/(1-Y1)));



MaxFunction pi*dc^2/4*epsb*(u03-u02)*Cf*Y1/of_bed/((Nc1+Nc2+Nc3+(Nctot-Nc1-Nc2-Nc3))*Lc*pi*(dc^2)/4)*1440  ;
Cs01<Cf;
((Cs01*a2*b1*(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)+Cs01*a1*b2-(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)*a1+(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)*a2)/((((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)*(a2*b1*(((-a1*(1+b2*Cf)+a2*(1+b1*Cf))+((-a1*(1+b2*Cf)+a2*(1+b1*Cf))^2+(4*a1*a2*b1*b2*Cf*Cf))^0.5)/2/a2/b1/Cf)+a1*b2)))<Cf;
Cf=Cf;
u03<u01;
u02<u03;
u04<u01;
u02<u01;
u04<u03;
u03-u02<=(P*(delta3-delta2)^2/4/(beta3*delta3^2/ku3/Nc3/Lc+beta2*delta2^2/ku2/Nc2/Lc)-(beta3*Eb3/Nc3/Lc+beta2*Eb2/Nc2/Lc));
-u03*Cs01+(u03-u02)*Cf<0;
(u03-u02)/(u03-u04)*Cf*Y1=Cp1;
Cf*Y2/(Cf*Y2+Cf*(1-Y1))>=Pur_reqE;
Cf*Y1/(Cf*Y1+Cf*(1-Y2))>=Pur_reqR;
v-u01/(1+P*delta1)=-beta1/(1+P*delta1)/Nc1/Lc*(Eb1+P*v^2*delta1^2/ku1);
v-u02/(1+P*delta2)=-beta2/(1+P*delta2)/Nc2/Lc*(Eb2+P*v^2*delta2^2/ku2);
v-u03/(1+P*delta3)=beta3/(1+P*delta3)/Nc3/Lc*(Eb3+P*v^2*delta3^2/ku3);
v-u04/(1+P*delta4)=beta4/(1+P*delta4)/(Nctot-Nc1-Nc2-Nc3)/Lc*(Eb4+P*v^2*delta4^2/ku4);
Lc-0.5*v>=0;
(((150*mu*u01*Nc1*Lc/dp^2*((1-epsb)/epsb)^2*(1000000)/6+1.75*rou*(u01)^2*Nc1*Lc/dp*(1-epsb)/epsb/3.6)*14.7/1.013/100000)+((150*mu*u02*Nc2*Lc/dp^2*((1-epsb)/epsb)^2*(1000000)/6+1.75*rou*(u02)^2*Nc2*Lc/dp*(1-epsb)/epsb/3.6)*14.7/1.013/100000)+((150*mu*u03*Nc3*Lc/dp^2*((1-epsb)/epsb)^2*(1000000)/6+1.75*rou*(u03)^2*Nc3*Lc/dp*(1-epsb)/epsb/3.6)*14.7/1.013/100000)+((150*mu*u04*(Nctot-Nc1-Nc2-Nc3)*Lc/dp^2*((1-epsb)/epsb)^2*(1000000)/6+1.75*rou*(u04)^2*(Nctot-Nc1-Nc2-Nc3)*Lc/dp*(1-epsb)/epsb/3.6)*14.7/1.013/100000))<=deltaPmax;
(a2*b1*gamaplusf^2+a1*b2*gamaplusf)*Cs02-(a2*b1*gamaplusf+a1*b2)*Cs01+(a1-a2)*gamaplusf=0;

Nctot-Nc1-Nc2-Nc3>=2;
Nctot-Nc1-Nc2-Nc3<=10;

shamohu

问题长的吓人啊！
1：“因为是按照外国的论文算个结果啊”，模型建的是否合理正确？表达式都对吗？有个约束“cf=cf”，不奇怪吗？笔误？其它地方还有吗？问题背景不明，又这么长，只有你自己检查了！
2：“所以知道答案的”，将你知道的答案代回模型验算一下：各约束是否满足？目标函数多少、与已知答案相同吗？从侧面检验你建的模型是否正确。
3：希望有人化几个小时或更多时间在这个问题上...。

stextinction

啊，cf=cf确实是笔误
因为我开始是两个变量Cf1和Cf2
约束则是CF1=Cf2
化简约束的时候直接用一个变量CF就代表两个了
这个答案我带回去了，约束满足，目标函数也是一致的
至少证明数学公式都是正确的啊
模型肯定是没有错
里面牵扯的有几个非线性方程
一般方法没法子定初值的
所以老困惑的～

shamohu

答案都知道还没法定初值吗？
不知道本计算的目的，是想验证老外的结果，还是想求一个更好的结果？把已知的结果放上来看看！