| 查看: 321 | 回复: 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咋样 |
回复此楼» 猜你喜欢
散金币祈福
已经有186人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
-
物理学I论文润色/翻译怎么收费?
已经有100人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
深圳大学应用技术学院招聘凝聚态物理博士后
已经有114人回复
-
深圳大学2027级光电信息工程/物理学研究生推免
已经有70人回复
-
求助中医院临床中医中标的青年标书学习，谢谢
已经有4人回复
-
浙大物院表面物理课题组长期招收普博及直博学生
已经有257人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有39人回复