| 查看: 2217 | 回复: 4 | ||
[求助]
求助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，结果点运行没反应，不知道是不是版本问题还是代码写的有问题。。 求各路牛人帮忙用高版本运行下，非常非常感谢~~ |
回复此楼» 本帖@通知
» 猜你喜欢
瑞典林雪平大学博士后招聘｜PEC Water Splitting 方向
已经有130人回复
-
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
物理学I论文润色/翻译怎么收费?
已经有127人回复
-
PVA溶解
已经有2人回复
-
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有11人回复
-
散金币祈福
已经有94人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求大神帮忙拟合一个非线性方程，求出模型参数
已经有15人回复
- 求1stopt拟合动力学参数
已经有7人回复
- 【求专家】根据一组数据，求其方程，寻最优拟合方程。要求截距为0
已经有3人回复
- 求助用matlab拟合动力学方程
已经有9人回复
- matlab 拟合方法
已经有4人回复
- 求大神帮忙用1stOpt 拟合复数数据
已经有13人回复
- 通过实验数据拟合，求解公式中的参数
已经有14人回复
- 微分方程组参数拟合的问题（多参数动力系统参数估计）
已经有11人回复
- 重金请教-如何用MAtlab对力学模型进行参数识别
已经有12人回复
- MATLAB微分方程参数拟合问题，求大神
已经有7人回复
- 大家在使用1stOpt进行数据拟合的时候有没有遇到过这种情况？
已经有8人回复
- 用MATLAB怎么打开点云并把它拟合到一个平面上？
已经有8人回复
- matlab数据拟合问题求助
已经有7人回复
- 动力学参数拟合
已经有26人回复
- matlab 拟合反应动力学参数结果很差。大家帮忙看一下
已经有14人回复
- matlab拟合拟均相动力学参数
已经有6人回复
- 如何用VTF方程拟合电导率与温度曲线
已经有8人回复
- matlab非线性参数拟合问题
已经有7人回复
- 【求助】用matlab最优化方法进行参数拟合
已经有17人回复
- 【求助】拟合动力学方程求助
已经有13人回复
- 【求助】使用MATLAB怎么实现拟合力场参数的程序？【已解决】
已经有15人回复
- 【求助】递推公式的参数拟合
已经有5人回复
2楼2014-04-01 09:45:56
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
【答案】应助回帖
感谢参与，应助指数 +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
|
真是太感谢了。不过这个拟合效果确实挺差的，不知道是不是模型不太合适。。 你能帮我再试试这组数据么？谢谢~ 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
5楼2015-12-07 10:55:40
20