24小时热门版块排行榜    

查看: 858  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 301求调剂 +3 A_JiXing 2026-03-16 3/150 2026-03-17 08:21 by 无际的草原
[考研] 283求调剂 +3 听风就是雨; 2026-03-16 3/150 2026-03-17 07:41 by 热情沙漠
[考研] 341求调剂 +5 捣蛋猪猪 2026-03-11 5/250 2026-03-17 06:56 by Muziwl
[考研] 机械专硕325,寻找调剂院校 +3 y9999 2026-03-15 5/250 2026-03-16 19:58 by y9999
[考研] 材料专硕306英一数二 +4 z1z2z3879 2026-03-16 6/300 2026-03-16 19:38 by z1z2z3879
[考研] 化学调剂0703 +8 啊我我的 2026-03-11 8/400 2026-03-16 17:23 by 我的船我的海
[考研] 一志愿985,本科211,0817化学工程与技术319求调剂 +5 Liwangman 2026-03-15 5/250 2026-03-16 17:10 by 我的船我的海
[考研] 一志愿苏州大学材料工程(085601)专硕有科研经历三项国奖两个实用型专利一项省级立项 +3 大火山小火山 2026-03-16 5/250 2026-03-16 16:54 by barlinike
[考研] 304求调剂 +3 曼殊2266 2026-03-14 3/150 2026-03-16 16:39 by houyaoxu
[考研] 0856求调剂 +3 刘梦微 2026-03-15 3/150 2026-03-16 10:00 by houyaoxu
[考研] 0856专硕279求调剂 +5 加油加油!? 2026-03-15 5/250 2026-03-15 11:58 by 2020015
[考研] 本科南京大学一志愿川大药学327 +3 麦田耕者 2026-03-14 3/150 2026-03-14 20:04 by 外星文明
[考研] 328求调剂 +3 5201314Lsy! 2026-03-13 6/300 2026-03-14 15:31 by hyswxzs
[考研] 26考研调剂 +3 ying123. 2026-03-10 3/150 2026-03-14 00:18 by JourneyLucky
[考研] 0703化学调剂 +4 快乐的香蕉 2026-03-11 4/200 2026-03-13 22:41 by JourneyLucky
[考研] 材料与化工求调剂一志愿 985 总分 295 +8 dream…… 2026-03-12 8/400 2026-03-13 22:17 by 星空星月
[考研] 295求调剂 +3 小匕仔汁 2026-03-12 3/150 2026-03-13 15:17 by vgtyfty
[考研] 308求调剂 +3 是Lupa啊 2026-03-12 3/150 2026-03-13 14:30 by 求调剂zz
[考博] 26申博求助 +3 跳跃饼干 2026-03-10 4/200 2026-03-10 21:15 by Tntcnn
[考研] 一志愿:武汉理工,材料工程,英二数二 总分314 +3 2202020125 2026-03-10 4/200 2026-03-10 13:54 by xiongyaxuan
信息提示
请填处理意见