| 查看: 293 | 回复: 0 | ||
xiaohsu2006木虫 (著名写手)
|
[求助]
含泪跪求1stopt高版本跑个代码 已经写好了
|
|
代码1（550）： ConstStr r1=k1*CA, r2=k2*CB*CD, r3=k3*CB*CD, r4=k4*CB, r5=k5*CB, r6=k6*CB, r7=k7*CB, r8=k8*CC*CD, r9=k9*CC*CD, r10=k10*CD*CF, r10r=k10/3.64*CE*CH, r11=k11*CE*CF, r11r=k11/273.7*CD*CG; Variable t,CE,CF,CG,CH; InitialODEValue t=0, CA=0.730635793, CB=0, CC=1.498246463, CD=3.332429483, CE=0, CF=0, CG=0, CH=0; ODEFunction CA’=-r1; CB’=r1-r2-r3-r4-7*r5-4*r6-11*r7; CC’=-r8-r9; CD’=-r2-2*r3-r8-2*r9-r10+r10r+r11-r11r; CE’=r2+2*r3+6*r4+r8+2*r9+r10-r10r-3*r11+3*r11r; CF’=r2+3*r5+r8-r10+r10r-r11+r11r; CG’=3*r6+r11-r11r; CH’=r3+3*r7+r9+r10-r10r; Data; //t/s CE CF CG CH 0 0 0 0 0 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 代码2（700）： ConstStr r1=k1*CA, r2=k2*CB*CD, r3=k3*CB*CD, r4=k4*CB, r5=k5*CB, r6=k6*CB, r7=k7*CB, r8=k8*CC*CD, r9=k9*CC*CD, r10=k10*CD*CF, r10r=k10/1.62*CE*CH, r11=k11*CE*CF, r11r=k11/1.65*CD*CG; Variable t,CE,CF,CG,CH; InitialODEValue t=0, CA=0.730635793, CB=0, CC=1.498246463, CD=3.332429483, CE=0, CF=0, CG=0, CH=0; ODEFunction CA’=-r1; CB’=r1-r2-r3-r4-7*r5-4*r6-11*r7; CC’=-r8-r9; CD’=-r2-2*r3-r8-2*r9-r10+r10r+r11-r11r; CE’=r2+2*r3+6*r4+r8+2*r9+r10-r10r-3*r11+3*r11r; CF’=r2+3*r5+r8-r10+r10r-r11+r11r; CG’=3*r6+r11-r11r; CH’=r3+3*r7+r9+r10-r10r; Data; //t/s CE CF CG CH 0 0 0 0 0 30 0.004375416 0.003086408 0.003315778 0.015026091 60 0.019617178 0.003836818 0.007720506 0.019703998 180 0.038487661 0.003891664 0.011588798 0.031720487 300 0.048176609 0.003933202 0.015502672 0.036733943 600 0.062063433 0.004022325 0.016942726 0.044129082 1200 0.075194185 0.004214514 0.023919896 0.055351964 1800 0.076085076 0.004151429 0.020840457 0.053030061 待求参数范围均在[0,1]内，这个不知道怎么写到代码里，求大神帮忙，跪谢！用matlab和mandona搞了好久得不到好的结果，不知道1stopt咋样 |
回复此楼» 猜你喜欢
大豆异黄酮分离
已经有0人回复
-
湖南大学材料学院急招2026年博士生，临时增加一名博士联培指标
已经有10人回复
-
物理学I论文润色/翻译怎么收费?
已经有183人回复
-
天津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生1-2名
已经有1人回复
-
中国科学院物理研究所谌志国研究员团队招收2027年博士研究生
已经有5人回复
-
2026年中德博士后交流项目 - 新型量子和磁性材料：材料制备表征和中子散射研究
已经有12人回复
-
26申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有1人回复
-
如何从铁电相到顺电相。
已经有1人回复