中间变量有两种取值方法，怎么办？
Parameteruh=0.3,ul=0.2,ur=0.5,rH=2.2,rL=1.8,pai=5;LoopConstanttauH1=3.6:0.05:4.8;
PlotLoopDatatauH1,ObjFunction;
Constantk=1.4,etaC=0.9,etaT=0.9,tauL1=1,D=0.96,Cwf=1,UT=20;
ConstStrUUH=UT*uh,UUL=UT*ul,UUR=UT*ur;
ConstStrm=(k-1)/k;
ConstStrCH=rH*Cwf,CL=rL*Cwf;
ConstStrNR=UUR/Cwf,ER=NR/(NR+1);
ConstStrCHmin=min(Cwf,CH),CLmin=min(Cwf,CL);
ConstStrCHmax=max(Cwf,CH),CLmax=max(Cwf,CL);
ConstStrNH=UUH/CHmin,NL=UUL/CLmin;
//*******************************************************************
//*****************若CH不等于Cwf，则EH按下式计算***
ConstStrEH=(1-exp(-NH*(1-CHmin/CHmax)))/(1-(CHmin/CHmax)*exp(-NH*(1-CHmin/CHmax)));
//*****************若CH等于Cwf，则EH按下式计算***********************
ConstStrEH=1-exp(-NH);
//*****************若CL不等于Cwf，则EL按下式计算*********************
ConstStrEL=(1-exp(-NL*(1-CLmin/CLmax)))/(1-(CLmin/CLmax)*exp(-NL*(1-CLmin/CLmax)));
//*****************若CL等于Cwf，则EL按下式计算***********************
ConstStrEL=1-exp(-NL);
//*******************************************************************
//***********************以下部分不变********************************
ConstStrx=pai^m,y=D^m;
ConstStra1=x*y+(1-x*y)*etaT;
ConstStra2=-1+x+etaC;
ConstStra3=-2+2*x+etaC;
ConstStra4=-x*y+2*(-1+x*y)*etaT;
ConstStra5=(-1+x)*x*y+(1-x*y)*a3*etaT;
ConstStrk1=Cwf^2;
ConstStrb1=a2*CLmin*Cwf*EL+CHmin*EH*(a2*CLmin*EL*(-1+2*ER)+Cwf*(a2-a3*ER));
ConstStrb2=(-1+x*y)*a2*etaT;
ConstStrb3=(-1+x)*x*y;
ConstStrG=a1*b1+a2*a4*CLmin*Cwf*EL*ER+k1*(b2+b3*(-1+ER)+(1-x*y)*a3*ER*etaT);
ConstStrw=1/G*(x*y*etaC*(CHmin*EH*(Cwf+CLmin*EL*(-1+ER))*tauH1
+CLmin*Cwf*EL*tauL1)-a2*ER*(x*y*CHmin*EH*(Cwf-CLmin*EL)*tauH1
+a4*CLmin*Cwf*EL*tauL1)-a1*(-ER*(a3*CHmin*Cwf*EH*tauH1
+CLmin*(-Cwf+CHmin*EH)*EL*etaC*tauL1)+a2*(CLmin*Cwf*EL*tauL1
+CHmin*EH*(Cwf*tauH1+CLmin*EL*((-1+2*ER)*tauH1+ER*tauL1)))));
MaxFunctionw;
uh+ul+ur=1;
0<eta<0.5;如以上代码所示，中间变量EH和EL因CH和CL是否等于Cwf而有两种取值方法，请问这样的优化问题该怎么做？是不是要用到编辑模式？
请大家不吝赐教。 可用if语句：
//*****************若CH不等于Cwf，则EH按下式计算***
ConstStrEH=(1-exp(-NH*(1-CHmin/CHmax)))/(1-(CHmin/CHmax)*exp(-NH*(1-CHmin/CHmax)));
//*****************若CH等于Cwf，则EH按下式计算***********************
ConstStrEH=1-exp(-NH);
改为：
ConstStrEH=if(CH=Cwf, 1-exp(-NH), (1-exp(-NH*(1-CHmin/CHmax)))/(1-(CHmin/CHmax)*exp(-NH*(1-CHmin/CHmax))));
//*****************若CL不等于Cwf，则EL按下式计算*********************
ConstStrEL=(1-exp(-NL*(1-CLmin/CLmax)))/(1-(CLmin/CLmax)*exp(-NL*(1-CLmin/CLmax)));
//*****************若CL等于Cwf，则EL按下式计算***********************
ConstStrEL=1-exp(-NL);
改为：
ConstStrEH=if(CL=Cwf, 1-exp(-NL), (1-exp(-NL*(1-CLmin/CLmax)))/(1-(CLmin/CLmax)*exp(-NL*(1-CLmin/CLmax))));
shamohu 发表于 2012-5-12 15:09 static/image/common/back.gif
可用if语句：
//*****************若CH不等于Cwf，则EH按下式计算***
非常感谢！这下可以解决很多历史遗留问题了。
页:
[1]