应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 求助1stopt拟合动力学参数
51/1返回列表
查看: 1736  |  回复: 4
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

chen1215

新虫 (初入文坛)

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

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2014-03-31 17:13:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与，应助指数 +1
如果都要求大于0：

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

参数                  最佳估算
--------------------        -------------
a        2379421492.52982
b        5477871570.00215
c        4827348354.29146
d        0.000157698484986776
e        28.2125368960529
f        9.53786176032263E-8
g        333524832.481022

====== 结果输出 ======

文件: 数据文件-1
No        t        目标 x        计算 x        目标 y        计算 y        目标 z        计算 z        目标 l        计算 l        目标 m        计算 m        目标 n        计算 n
1        5        0.598029237        0.461228912051106        0.192347906        0.212216074803243        0.044174154        0.154698496352178        0.096131452        0.102539265793472        0.569090571        0.18430565779089        4.022277657        4.40706257020911
2        10        0.541233252        0.352387819112766        0.219628094        0.194209366200932        0.058131509        0.186948110939243        0.111689894        0.19713745274706        0.670960794        0.259898922792036        3.920407434        4.33146930520796
3        20        0.502127725        0.241106037906148        0.23819892        0.153262848683966        0.066640197        0.191130610932665        0.123715906        0.34518325147722        0.742627034        0.355979862816147        3.848741194        4.23538836518385
4        30        0.46364134        0.179231653280274        0.251976191        0.121526115166401        0.073187605        0.174250824742148        0.141877613        0.455674155811178        0.823984239        0.419251043228081        3.767383988        4.17211718477191
5        40        0.442517119        0.138636909389627        0.252181432        0.0975440616385669        0.084552378        0.153342396668236        0.15143182        0.541159381303569        0.875581648        0.465201935622196        3.71578658        4.1261662923778
6        60        0.404660654        0.0883552573120677        0.264979117        0.0646893622260779        0.093852475        0.114305811509068        0.167190503        0.663332317952786        0.954255576        0.527617939775541        3.637112651        4.06375028822446
赞一下 回复此楼
3楼2014-04-03 14:01:19
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答

终极爱好者

铜虫 (初入文坛)
想学 matlab 不知道 改怎么获取资料
赞一下 回复此楼
2楼2014-04-01 09:45:56
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chen1215

新虫 (初入文坛)
引用回帖:
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的回帖

蜗牛虫

银虫 (正式写手)
引用回帖:
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- ...

您好，您是怎么做到的？

肯求授以渔：）


非常感谢
赞一下 回复此楼
多思考勤记录嗯加油
5楼2015-12-07 10:55:40
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭