24小时热门版块排行榜    

查看: 312  |  回复: 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咋样
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xiaohsu2006 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 gy1nBQXYQJqL 2026-08-29 3/150 2026-08-29 21:32 by 4qO8zhfjr2AX
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 gy1nBQXYQJqL 2026-08-29 5/250 2026-08-29 21:22 by 4qO8zhfjr2AX
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 ASdOkHsho7FD 2026-08-28 11/550 2026-08-29 20:18 by 63GmZBFV29mG
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 有没有仍没收到信息的 +5 德尚中行 2026-08-27 6/300 2026-08-29 18:10 by manplx
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +6 ASdOkHsho7FD 2026-08-28 7/350 2026-08-29 17:31 by 4FFAWE8HcgUD
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 ASdOkHsho7FD 2026-08-28 6/300 2026-08-29 17:30 by 4FFAWE8HcgUD
[基金申请] 为什么到现在没收到通知? +4 tannykie 2026-08-29 4/200 2026-08-29 17:14 by lmz0216
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +5 ASdOkHsho7FD 2026-08-28 9/450 2026-08-29 14:03 by jCd0dEvKHShX
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 ASdOkHsho7FD 2026-08-28 6/300 2026-08-29 11:51 by jCd0dEvKHShX
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ASdOkHsho7FD 2026-08-28 6/300 2026-08-29 11:40 by jCd0dEvKHShX
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 ASdOkHsho7FD 2026-08-28 9/450 2026-08-29 11:18 by jCd0dEvKHShX
[基金申请] 中青基了要发朋友圈吗? +3 349506619 2026-08-28 3/150 2026-08-29 10:49 by 木风9012
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 ASdOkHsho7FD 2026-08-28 7/350 2026-08-29 09:26 by G6APbkg8SA6w
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
信息提示
请填处理意见