| 查看: 1016 | 回复: 4 | ||
[求助]
这样的动力学方程能够有1stopt拟合出参数来么 已有1人参与
|
|
拟合出参数 k0 k1 k2 k3, 动力学方程如下 dz/dt=-2.938471412e-3*(1-y-z)^2*z*(1/2+3*log(1-x))/sqrt(1-3*log(1-x))*k0+1.10249053e-2*(1-y-z)^2*k1-2.938471412e-3*(1-y-z)*y*z*(1/2+3*log(1-x))/sqrt(1-3*log(1-x))*k2-2.938471412e-3*(1-y-z)*y*z*(1/2+3*log(1-x))/sqrt(1-3*log(1-x))*k3-2.606063248e-4*z-7.984680048e-4*y*z*(1/2+3*log(1-x))/sqrt(1-3*log(1-x))-6.945906685e3*z*z*(1/2+3*log(1-x))/sqrt(1-3*log(1-x)); dx/dt=2.938471412e-3*(1-x)*sqrt(1-3*log(1-x))*(1-y-z)^2*k0+2.938471412e-3*(1-x)*sqrt(1-3*log(1-x))*(1-y-z)*y+7.984680048e-4*(1-x)*sqrt(1-3*log(1-x))*y+6.945906685e3*(1-x)*sqrt(1-3*log(1-x))*z; dy/dt=(9.180837032e-3*(1-y-z)^2-2.938471412e-3*(1/2+3*log(1-x))/sqrt(1-3*log(1-x))*y*(1-y-z)^2)*k0+2.938471412e-3*(1/2+3*log(1-x))/sqrt(1-3*log(1-x))*y*(1-y-z)*y*k2+(27.21200097*y*(1-y-z)-2.938471412e-3*(1/2+3*log(1-x))/sqrt(1-3*log(1-x))*y*(1-y-z)*y)*k3-2.494699999e-3*y-7.984680048e-4*(1/2+3*log(1-x))/sqrt(1-3*log(1-x))*y*y-6.945906685e3*(1/2+3*log(1-x))/sqrt(1-3*log(1-x))*y*z; 0=<x<=1 0=<y<=1 0=<z<=1 t=0 时 x=0 y=0 z=0 实验给出的数据为（第二个值为整个表达式的值）： x [k0*1.189664539e-5*(1-y-z)^2+k2*1.189664539e-5*(1-y-z)*y+3.232663987e-6*y+28.1210797*z] 0.1 3.e-6 0.2 3.7e-6 0.3 4.1e-6 0.4 4.e-6 0.5 3.4e-6 0.6 2.3e-6 0.7 0.9e-6 0.8 0.7e-6 0.9 0.6e-6 请教各位大侠能够拟合出参数 k0 k1 k2 k3 (均需要大于0），不知道1stopt可以求解出来么 |
回复此楼» 猜你喜欢
【复旦大学】二维材料方向招收2026年博士研究生1名
已经有0人回复
-
北京纳米能源与系统研究所 王中林院士/曹南颖研究员课题组2026级硕/博/博后招生
已经有10人回复
-
物理学I论文润色/翻译怎么收费?
已经有200人回复
-
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生
已经有23人回复
-
反铁磁体中的磁性切换：两种不同的机制已成功可视化
已经有0人回复
-
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
-
新西兰Robinson研究所招收全奖PhD
已经有0人回复
-
石墨烯转移--二氧化硅衬底石墨烯
已经有0人回复
-
笼目材料中量子自旋液体基态的证据
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助--用1stopt拟合蓝宝石衬底的GaN薄膜红外反射光谱
已经有6人回复
- 求助1stopt拟合动力学参数
已经有4人回复
- 用1stOpt进行复数曲线拟合时，结果出错。急求教！！！
已经有5人回复
- 求助用matlab拟合动力学方程
已经有9人回复
- 如何用已知函数去拟合实验数据，得出参数
已经有10人回复
- 怎样拟合数据求出x=f(A,B,C),y=g(A,B,C,D)中的A,B,C,D参数
已经有18人回复
- origin拟合Birch-Murnaghan状态方程求B0，B0'和v0
已经有38人回复
- 用1stopt拟合非线性方程结果与其他软件拟合结果差异大
已经有5人回复
- 动力学方程参数估计方法
已经有14人回复
- 大家在使用1stOpt进行数据拟合的时候有没有遇到过这种情况？
已经有8人回复
- 1stopt或matlab如何求解以下的非线性方程并拟合出相应曲线？
已经有13人回复
- matlab 拟合反应动力学参数结果很差。大家帮忙看一下
已经有14人回复
- 如何用VTF方程拟合电导率与温度曲线
已经有8人回复
- 用origin或1stopt拟合一条隐函数方程曲线
已经有4人回复
- 怎么由测试数据得到拟合方程
已经有11人回复
- 【求助】拟合动力学方程求助
已经有13人回复
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.5小时
- 虫号: 291104
- 注册: 2006-10-28
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
cpl6618: 金币+50, ★★★很有帮助 2014-11-14 16:55:47
感谢参与，应助指数 +1
cpl6618: 金币+50, ★★★很有帮助 2014-11-14 16:55:47
|
参考下： 参数 最佳估算 -------------------- ------------- k0 20.1647513744796 k1 0.00302377880204425 k2 0.669576609638345 k3 0.238331929611515 ====== 结果输出 ====== 文件: 数据文件-1 No t 目标 x 计算 x 目标 k0*1.189664539e-5*(1-y-z)^2+k2*1.189664539e-5*(1-y-z)*y+3.232663987e-6*y+28.1210797*z 计算 k0*1.189664539e-5*(1-y-z)^2+k2*1.189664539e-5*(1-y-z)*y+3.232663987e-6*y+28.1210797*z 1 1 0.1 0.0998298792714663 3E-6 0.000393430027969079 2 2 0.2 0.199996948986132 3.7E-6 0.000392151136999175 3 3 0.3 0.3000918871667 4.1E-6 0.000402298980955307 4 4 0.4 0.40012978802504 4E-6 0.00042411850382553 5 5 0.5 0.500106755208013 3.4E-6 0.000461288107613004 6 6 0.6 0.600030007554687 2.3E-6 0.00052234991547497 7 7 0.7 0.699925254905772 9E-7 0.000627789264184172 8 8 0.8 0.799857809623034 7E-7 0.000838217623415368 9 9 0.9 0.900011169138298 6E-7 0.00144623022373588 |
2楼2014-11-14 15:58:22
3楼2014-11-14 16:23:55
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.5小时
- 虫号: 291104
- 注册: 2006-10-28
4楼2014-11-14 16:32:31
5楼2014-11-14 16:54:34