之前你帮忙做出的结果已经很好了,但有个别数据太小,所以我做了替代数据,请帮忙再跑一下,万分感谢。。
下面是程序:
Title "test 1";
Parameters a1[0,],a2[0,],a3[0,],b1[0,],b2[0,],b3[0,],c1[0,],c2[0,],c3[0,],M[1.5,],C[y,],E[z,],N[1.5,];
Sharedmodel;
Variable t, x,y,z;
Function x'=a1*x*(y/M-1)-a2*y+a3*z;
y'=-b1*x+b2*y*(1-y/C)+b3*z*(1-z/E);
z'=c1*x*(x/N-1)-c2*y-c3*z;
RowData:
x=[1.0000,1.0777,1.1443,1.3497,1.5881,1.9814,2.1905,2.5799,2.9763,3.2445,3.5765,3.9315];
y=[1.0000,1.0597,1.0982,1.2035,1.3566,1.4610,1.4771,1.4862,1.6052,1.6094,1.8892,2.0334];
z=[1.0000,1.0920,1.2033,1.4031,1.6920,1.9382,2.2159,2.6186,2.9489,3.1536,3.5979,4.0233]; |