| 查看: 285 | 回复: 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咋样 |
回复此楼» 猜你喜欢
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
-
新西兰Robinson研究所招收全奖PhD
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有149人回复
-
石墨烯转移--二氧化硅衬底石墨烯
已经有0人回复
-
笼目材料中量子自旋液体基态的证据
已经有0人回复
-
数学教学论硕士可以读数学物理博士吗？
已经有0人回复
-
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生
已经有4人回复
-
澳门大学 应用物理及材料工程研究院 潘晖教授课题组诚招博士后
已经有11人回复
-
求助NH4V4O10晶体的CIF文件
已经有0人回复
-
英国全奖博士招聘-深度学习与量子物理
已经有0人回复
-
间接带隙半导体有效质量求助
已经有0人回复