24小时热门版块排行榜    

查看: 608  |  回复: 0

xiaohsu2006

木虫 (著名写手)

[求助] 求高版本1stopt拟合参数,代码已写好

(k值均在[0,1]区间内)
代码如下:
ConstStr r1=k1*CA*CB,
r2=k2*CA*CB,
r3=k3*CA*CC,
r4=k4*CA*CC,
r5=k5*CA*CE,
r5r=k5/3.64*CD*CG,
r6=k6*CD*CE,
r6r=k6/273.7*CA*CF;

Variable t,CD,CE,CF,CG;
InitialODEValue t=0, CA=3.332429483, CB=0.012090459, CC=0.012090459, CD=0, CE=0, CF=0, CG=0;
ODEFunction CA’=-4*r1-4*r2-8*r3-16*r4-r5+r5r+r6-r6r;
CB’=-r1-r2;
CC’=-r3-r4;
CD’=4*r2+8*r3+16*r4+r5-r5r-3*r6+3*r6r;
CE’=4*r2+8*r3-r5+r5r-r6+r6r;
CF’=2*r1+r6-r6r;
CG’=2*r1+8*r4+r5-r5r;
Data;
//t/s                                CD                          CE                          CF                         CG
30                0.000149796        0.000760865        0.000217173        0.005495421
60                0.002059783        0.002324358        0.002416821        0.011683547
180                0.006129219        0.002600548        0.004684001        0.014443429
300                0.007712601        0.002375765        0.005132747        0.0146369
600                0.009297344        0.002233589        0.005363171        0.018179198
1200                      0.012671286        0.002038349        0.006919546        0.020292707
1800                      0.014069376        0.002154183        0.007150669        0.021269902
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

相关版块跳转 我要订阅楼主 xiaohsu2006 的主题更新
信息提示
请填处理意见