24小时热门版块排行榜    

查看: 1003  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 时间戳变了,能看出什么问题? +8 基诺咪客 2026-08-17 11/550 2026-08-18 17:13 by 悠然熊猫
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +50 医学老男孩 2026-08-13 109/5450 2026-08-18 16:56 by 医学老男孩
[基金申请] 重要消息,中午系统在维护 +9 yuleib84 2026-08-18 10/500 2026-08-18 16:52 by archvillain
[教师之家] 两姐妹,一个辛辛苦苦30岁工作当大学教师,另一个23岁就工作当开车教练 +3 瞬息宇宙 2026-08-11 7/350 2026-08-18 16:42 by LNP@mRNA
[基金申请] 明天放榜? +4 Shxjjxjkx 2026-08-18 4/200 2026-08-18 15:59 by abcabc1133
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 时间戳又变了8-15 +14 archvillain 2026-08-15 26/1300 2026-08-18 13:37 by phantomgost
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +13 Ldrop2023 2026-08-13 16/800 2026-08-18 12:25 by 淀粉搬运工
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 8/400 2026-08-18 11:47 by gltch
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] filecode=后面第一个是大写字母 +8 wangze12014 2026-08-14 10/500 2026-08-17 17:05 by xter9665
[基金申请] 今天系统多次维护,明天很可能放榜! +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
信息提示
请填处理意见