1stopt吧 关注:351贴子:2,280
  • 1回复贴,共1

求助,1stopt代码检查修改

只看楼主收藏回复

Parameters X0, Xm, um,a,b,P0,Yx,Yp,m,S0;
Variables t, X, P,S;
SharedModel;
Function X= X0*Xm*exp(um*t)/(Xm-X0+X0*exp(um*t))
P=P0+a*(X0*Xm*exp(um*t)/(Xm-X0+X0*exp(um*t))-X0)+(b*Xm/um)*ln((Xm-X0+X0*exp(um*t))/Xm)
S=S0-(1/Yx+a/Yp)*(X-X0)-Xm/um*(b/Y+m)*ln((Xm-X0+X0*exp(um*t))/Xm);
Data;
t X P S
0 0.05 0 79
4 0.16 0 78
8 0.6515 0 76
12 1.5426 1.5 74
16 2.4047 4.9 69
20 3.1066 8.7 60
24 3.6693 12.4 50
28 3.8165 16.1 43
32 3.8874 17.2 40
36 3.9618 17.9 38
上面三个函数表达式单独运算正确,但是三个一起运算就有误,求大神指点,在线等


1楼2014-12-17 13:09回复

    Parameters X0,Xm,um,a,b,P0,Yx,Yp,m,S0;
    Variables t, X, P,S;
    SharedModel;
    Function X= X0*Xm*exp(um*t)/(Xm-X0+X0*exp(um*t));
    P=P0+a*(X0*Xm*exp(um*t)/(Xm-X0+X0*exp(um*t))-X0)+(b*Xm/um)*ln((Xm-X0+X0*exp(um*t))/Xm);
    S=S0-(1/Yx+a/Yp)*(X-X0)-Xm/um*(b/Y+m)*ln((Xm-X0+X0*exp(um*t))/Xm);


    2楼2015-04-03 17:33
    回复