24小时热门版块排行榜    

Znn3bq.jpeg
查看: 932  |  回复: 2

xiaohsu2006

木虫 (著名写手)

[求助] 求高版本1stopt跑个程序,代码写好了

哪位大哥帮我跑个代码,我用matlab拟合了好久,初值老取不好,谢谢!代码写好了,如下
ConstStr r1=k1*CA,
r2=k2*CA,
r3=k3*CC*CB,
r4=k4*CC*CB,
r5=k5*CC,
r6=k6*CC,
r7=k7*CC,
r8=k8*CC,
r9=k9*CC,
r10=k10*CC,
r11=k11*CD*CB,
r12=k12*CD*CB,
r13=k13*CD,
r14=k14*CD,
r15=k15*CD,
r16=k16*CD,
r17=k17*CD,
r18=k18*CD,
r19=k19*CH*CB,
r19r=k19/3.64*CJ*CG,
r20=k20*CH*CG,
r20r=k20/273.7*CI*CB;
Variable t,CG,CH,CI,CJ;
InitialODEValue t=0, CA=0.016789972, CB=3.332429483, CC=0, CD=0, CE=0, CF=0, CG=0, CH=0, CI=0, CJ=0;
ODEFunction CA’=-r1-r2;
CB’=-5*r3-11*r4-5*r11-11*r12-r19+r19r+r20-r20r;
CC’=r1-r3-r4-r5-r6-r7-r8-r9-r10;
CD’=r2-r11-r12-r13-r14-r15-r16-r17-r18;
CE’=r5+r13;
CF’=r6+r14;
CG’=6.5*r3+12.5*r4+r7+6.5*r11+12.5*r12+r15+r19-r19r-3*r20+3*r20r;
CH’=6*r3+r8+6*r11+r16-r19+r19r-r20+r20r;
CI’=r9+r17+r20-r20r;
CJ’=6*r4+r10+6*r12+r18+r19-r19r;
Data;
//t/s        CG        CH          CI          CJ
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的回帖

dingd

铁杆木虫 (职业作家)

内容已删除

» 本帖已获得的红花(最新10朵)

2楼2014-05-20 14:16:53
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

xiaohsu2006

木虫 (著名写手)

送红花一朵
引用回帖:
2楼: Originally posted by dingd at 2014-05-20 14:16:53
参数没范围限制吗?下面是一组,供参考,计算很费时间:
均方差(RMSE):0.000514570470417657
残差平方和(SSE):7.41391753272375E-6
相关系数(R): 0.971702664574647
相关系数之平方(R^2): 0.944206068341469
确 ...

谢谢你,很不好意思,我把参数限制忘了放上去了,因为我没在这个软件上跑过,参数限制为:
k1 [0,3]
k2 [0,3]
k3 [0,1]
k4 [0,1]
k5 [0,1]
k6 [0,1]
k7 [0,1]
k8 [0,1]
k9 [0,1]
k10 [0,1]
k11 [0,1]
k12 [0,1]
k13 [0,1]
k14 [0,1]
k15 [0,1]
k16 [0,1]
k17 [0,1]
k18 [0,1]
k19 [0,1]
k20 [0,1]
3楼2014-05-20 14:54:02
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xiaohsu2006 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 26/27申博自荐 10+4 ZXW0202 2026-05-22 9/450 2026-05-24 08:47 by bjvtcliu
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 hvkbtfonbv 2026-05-23 3/150 2026-05-24 08:01 by 9ps9vgkqva
[基金申请] 河北省自然科学基金 +5 Peterchao 2026-05-18 7/350 2026-05-23 16:57 by articlefan
[基金申请] 揭秘青基评审内幕:几个A才能顺利中标 +3 国自然国社科中 2026-05-23 4/200 2026-05-23 15:37 by 2000zf36392
[基金申请] 青B发送上会通知了吗 +5 chemBioBro 2026-05-22 7/350 2026-05-23 12:35 by zhuifengzhy
[考博] 博士申请 +3 焦晓明 2026-05-21 3/150 2026-05-23 11:26 by mlc840311
[教师之家] 某211大学教师把个人教师官方主页改成:我跑了我跑了我跑了!官宣跑路! +3 zju2000 2026-05-21 4/200 2026-05-23 08:33 by xli1984
[基金申请] 今年审到国自然15份,谈谈感受 +27 国自然国社科中 2026-05-17 28/1400 2026-05-22 16:25 by hunter无悔
[论文投稿] 投稿求助,期刊 +4 希冀,有书读 2026-05-20 8/400 2026-05-22 10:16 by 希冀,有书读
[文学芳草园] 献血感触 +7 呀呀好傻 2026-05-19 13/650 2026-05-21 20:15 by 呀呀好傻
[基金申请] 面上本子正文33页,违规吗?会被低分嘛? +14 1234567wang 2026-05-17 16/800 2026-05-21 17:58 by 脆脆的饼干
[基金申请] 国自然评分 +4 无名者登山 2026-05-20 5/250 2026-05-21 16:35 by swuq
[基金申请] 国自然上会要求 +7 无名者登山 2026-05-18 11/550 2026-05-21 15:50 by draco1987
[基金申请] 提交了我也来说说感想 +9 fummck 2026-05-20 10/500 2026-05-21 14:17 by draco1987
[基金申请] 评审有感 +15 popular289 2026-05-18 26/1300 2026-05-21 10:35 by 西葫芦炒鸡蛋
[考博] 云南大学材料与能源学院解琳课题组钙钛矿博士招生 +4 光伏爱好者 2026-05-17 6/300 2026-05-21 01:37 by 17520222914
[有机交流] 反应很差,大量原料没有反应 5+3 Mr.Zot 2026-05-19 8/400 2026-05-20 22:19 by Equinoxhua
[考博] 如果工作了想读博,可以边工作边读全日制嘛? 30+3 铁达火车 2026-05-18 5/250 2026-05-20 09:33 by tfang
[考博] 博士申请 +5 星…… 2026-05-18 6/300 2026-05-18 23:49 by 糊糊涂涂好
[硕博家园] 我在等一个没有答案的答案 +3 Love_MH 2026-05-17 3/150 2026-05-18 02:22 by 竹林孤影
信息提示
请填处理意见