24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

xiaohsu2006

木虫 (著名写手)

应助: 23 (小学生)
金币: 2520.4
散金: 207
红花: 4
帖子: 1362
在线: 1137.1小时
虫号: 1012350
注册: 2010-05-06
性别: GG
专业: 多相流热物理学

[求助] 求高版本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

回复此楼

» 猜你喜欢

比利时鲁文大学固体物理系招收CSC博士已经有0人回复
🌟 比利时新鲁汶大学（UCLouvain）诚邀CSC博士加入Pascal Gehring教授团队已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有79人回复
第一性原理计算方向2026级博士申请 PRB*1，四级484 已经有1人回复
求助VISSIM破解版软件已经有0人回复
求2026年在台湾举行的物理和材料领域国际学术会议信息已经有0人回复
求国际会议网站已经有1人回复
求取一些关于纳米材料和纳米技术相关的英文PPT。已经有0人回复
【复旦大学】二维材料方向招收2026年博士研究生1名已经有0人回复
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组2026级硕/博/博后招生已经有10人回复

» 本主题相关价值贴推荐，对您同样有帮助:

求大神指导1stopt 非线性最优化问题已经有3人回复
求高手给做一道题，用MATLAB或1stOpt，请把结果直接给出来，我的软件就不行，谢谢已经有5人回复
求高版本1stopt跑一个程序已经有7人回复
1stopt求最小值已经有37人回复
求大神帮忙用1stOpt 拟合复数数据已经有13人回复
1stopt做最小二乘法回归参数，请求大家帮忙指导下已经有3人回复
求前辈帮忙用1stOpt你个一下数据，写论文急用，先谢过已经有14人回复
求高版本1stopt拟合，已经有11人回复
1stopt中使用sharedmodel拟合曲线出现问题，求指示！！已经有6人回复
求大哥大姐帮跑一下1stopt~新手啊，版本低啊已经有9人回复
求大侠帮助跑一下1stopt 已经有12人回复
求高版本1stopt，解决一个问题，微分方程组参数拟合已经有3人回复
运算过程有复数出现，求高版本1stOpt 已经有6人回复
有1stOpt4.0版本的朋友帮忙拟合一下复数吧感激不尽呐已经有4人回复
学1stopt要多久？已经有12人回复
用1stopt进行数据拟合时，无法运行sharedmodel ,求助高手已经有7人回复

1楼2014-05-20 09:48:06

已阅关注TA 给TA发消息送TA红花 TA的回帖

xiaohsu2006

木虫 (著名写手)

应助: 23 (小学生)
金币: 2520.4
散金: 207
红花: 4
帖子: 1362
在线: 1137.1小时
虫号: 1012350
注册: 2010-05-06
性别: GG
专业: 多相流热物理学

送红花一朵

引用回帖:

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的回帖

查看全部 3 个回答