|
|
@dingd，不好意思，楼主给出的方程中有错误，10楼的程序就不用看了，麻烦跑一下下面的程序吧：
Parameters K1,K2,K3,K4,K5,K6,K7,K8；
Variable t,T,D,M,G,W,F;
ODEFunction
T’= -K7*T*M+K8*D^2-K1*T*W+K2*D*F；
D’= 2*K7*T*M-2*K8*D^2-K1*T*W+K2*D*F-K3*D*W+K4*M*F；
M’=K7*T*M+K8*D^2+K3*D*W-K4*M*F-K5*M*W+K6*G*F；
W’=-K1*T*W+K2*D*F-K3*D*W+K4*M*F-K5*M*W+K6*G*F；
F’=K1*T*W-K2*D*F+K3*D*W+K4*M*F+K5*M*W+K6*G*F；
G’= K5*M*W-K6*G*F；
data；
0 99.04 0.81 0 0 0.15 0
1 68.34 18.69 0.96 3.9633 12.01 11.86
2 53.58 23.71 2.05 6.8178 20.66 20.51
3 46.46 26.96 3.02 7.7748 23.56 23.41
4 40.76 28.83 4.23 8.6394 26.18 26.03
5 38.11 29.21 4.52 9.2928 28.16 28.01
6 35.7 29.37 4.75 9.9594 30.18 30.03
8 32.51 30.29 4.80 10.692 32.4 32.25
10 31.04 30.78 4.99 10.9527 33.19 33.04
12 30.81 31.24 5.11 10.8372 32.84 32.69
24 28.84 29.08 6.09 11.8767 35.99 35.84 |
|
回复此楼
qinzhong6138