24小时热门版块排行榜

返回列表

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

chen1215

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 1351
帖子: 2
在线: 2.2小时
虫号: 2331211
注册: 2013-03-09

[求助] 求助1stopt拟合动力学参数已有1人参与

最近做反应动力学模型，需要拟合微分方程组的几个参数，之前用matlab，ode45求解，得到的数值解套用非线性最小二乘拟合，可是总也拟合不好。不知道是初值不合适还是模型有问题，但模型之前有文献报道过，有问题的可能性较小。
想求助用1stopt拟合试试，看看能不能得到比较好的效果~
我对1stopt 不太熟，试着写了下程序，不知道有没有问题，具体如下：
Parameter a,b,c,d,e,f,g;
Variable t,x,y,z,l,m,n;
Function x’=-a*x*n*(0.000384369/(1+d*x+e*y+f*z+g*m))；
y’=(a*x*n-b*y*n)*(0.000384369/(1+d*x+e*y+f*z+g*m))；
z’=(b*y*n-c*z*n)*(0.000384369/(1+d*x+e*y+f*z+g*m))；
l’=c*z*n*(0.000384369/(1+d*x+e*y+f*z+g*m))；
m’=(a*x+b*y+c*z)*(0.000384369/(1+d*x+e*y+f*z+g*m))；
n’= -m’；
RowData；
t=0 5 10 20 30 40 60;
x=0.930682749 0.598029237 0.541233252 0.502127725 0.46364134 0.442517119 0.404660654;
y=0    0.192347906 0.219628094 0.23819892 0.251976191 0.252181432 0.264979117;
z=0    0.044174154 0.058131509 0.066640197 0.073187605 0.084552378 0.093852475;
l=0    0.096131452 0.111689894 0.123715906 0.141877613 0.15143182 0.167190503;
m=0    0.569090571 0.670960794 0.742627034 0.823984239 0.875581648 0.954255576;
n=4.591368228 4.022277657 3.920407434 3.848741194 3.767383988 3.71578658 3.637112651;
就是已知x,y,z,l,m,n 随时间变化的值，想拟合出模型中的参数 a,b,c,d,e,f,g(这些参数的值都是大于等于0）
之前有下过低版本的1stopt，结果点运行没反应，不知道是不是版本问题还是代码写的有问题。。

求各路牛人帮忙用高版本运行下，非常非常感谢~~

回复此楼

» 本帖@通知

@dingd

» 猜你喜欢

求取一些关于纳米材料和纳米技术相关的英文PPT。已经有0人回复
【复旦大学】二维材料方向招收2026年博士研究生1名已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有66人回复
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组2026级硕/博/博后招生已经有10人回复
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生已经有22人回复
反铁磁体中的磁性切换：两种不同的机制已成功可视化已经有0人回复
26申博推荐：南京航空航天大学国际前沿科学研究院光学方向招收博士生！已经有0人回复
求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复
新西兰Robinson研究所招收全奖PhD 已经有0人回复

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

1楼2014-03-31 17:13:16

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

chen1215

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 1351
帖子: 2
在线: 2.2小时
虫号: 2331211
注册: 2013-03-09

引用回帖:

3楼: Originally posted by dingd at 2014-04-03 14:01:19
如果都要求大于0：

均方差(RMSE):0.291509496670823
残差平方和(SSE):3.05920031937395
相关系数(R): 0.41204372086458
相关系数之平方(R^2): 0.169780027903928
确定系数(DC): -32.2049580782807
F统计(F- ...

真是太感谢了。不过这个拟合效果确实挺差的，不知道是不是模型不太合适。。
你能帮我再试试这组数据么？谢谢~
t=0 2 5 10 15 20 30
x=0.97929461 0.848102002 0.791547773 0.742120992 0.665706681 0.627442016 0.567283612
y=0 0.111151416 0.154799598 0.178808192 0.203933465 0.230873341 0.255998614
z=0 0.029197641 0.031321393 0.058558032 0.078165633 0.1017788 0.122183142
l=0 0.000965082 0.003370178 0.008463798 0.012727551 0.019200453 0.030677511
m=0 0.142077349 0.222319921 0.295346436 0.454731222 0.492032299 0.601852621
n=1.929382737 1.787305387 1.707062816 1.634036301 1.474651515 1.437350438 1.327530115

赞一下

回复此楼

4楼2014-04-03 15:01:04

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

相关版块跳转我要订阅楼主 chen1215 的主题更新

返回列表