| 查看: 2278 | 回复: 11 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
matlab拟合反应动力学参数 已有2人参与
|
||
|
各位大佬，可以帮我拟合一下动力学参数吗，可以有偿。先谢过各位大佬。 t B C A 0 0.0024 0.0000 0.1894 6 0.0422 0.0000 0.1496 12 0.0838 0.0014 0.1066 18 0.1261 0.0061 0.0596 24 0.1470 0.0151 0.0297 30 0.1447 0.0361 0.0111 36 0.1279 0.0591 0.0047 42 0.1127 0.0767 0.0025 48 0.1009 0.0898 0.0011 54 0.0848 0.1070 0.0000 60 0.0721 0.1197 0.0000 ρ=0.0077 [H2]=4 r1 = k1*KA*[A]*KH2*[H2]/(1+KA*[A]+KB*[B]+KC*[C]+sqrt(KH2*[H2]))^3 r1 = k2*KB*[B]*KH2*[H2]/(1+KA*[A]+KB*[B]+KC*[C]+sqrt(KH2*[H2]))^3 1/ρ*d[A]/dt = -r1 1/ρ*d[B]/dt = r1-r2 1/ρ*d[C]/dt = r2 |
回复此楼» 猜你喜欢
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有3人回复
-
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有2人回复
-
物理学I论文润色/翻译怎么收费?
已经有115人回复
-
B区学生调剂-兰州交通大学材料科学与工程学院
已经有12人回复
-
山西大同大学物理学专业还有调剂名额，欢迎调剂！
已经有13人回复
-
桂林理工大学物理学专业招收调剂
已经有13人回复
-
VASP 的一组 GPU / CPU 基准测试记录
已经有0人回复
-
美国顶级物理期刊《应用物理快报》（APL）的编辑欺骗和歧视作者及AIP的官僚主义傲慢
已经有7人回复
-
津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生一名
已经有0人回复
-
【原创讨论】从电子约束到物质编辑：一套可迭代的环形磁场科技树
已经有0人回复
送红花一朵 
9楼2019-03-25 11:04:13
独孤神宇
版主 (知名作家)
- 应助: 490 (硕士)
- 贵宾: 0.008
- 金币: 31016.3
- 散金: 802
- 红花: 122
- 沙发: 1
- 帖子: 5600
- 在线: 856.5小时
- 虫号: 3522474
- 注册: 2014-11-06
- 性别: GG
- 专业: 机械动力学
- 管辖: 计算模拟
【答案】应助回帖
感谢参与，应助指数 +1
|
均方差(RMSE): 0.00213787148383087 残差平方和(SSR): 0.000137114834441316 相关系数(R): 0.998778073301579 相关系数之平方(R^2): 0.997557639708014 修正R平方(Adj. R^2): 0.995605193318058 确定系数(DC): 0.996987006259065 F统计(F-Statistic): 161.269285477097 参数 最佳估算 -------------------- ------------- k1 13.5564432263352 k2 9.47763860792964 ka 28.7183814838917 kb 4.13770213410864 kc 3.5066823001231 kh2 2.55638641506165 ====== 结果输出 ====== 文件: 数据文件-1 No t 目标 b 计算 b 目标 c 计算 c 目标 a 计算 a 1 6 0.0422 0.0424108511588022 0 0.000565879792438171 0.1496 0.14882326904876 2 12 0.0838 0.0837866483079356 0.0014 0.00281912430154442 0.1066 0.10519422739052 3 18 0.1261 0.121738165846101 0.0061 0.00847370300303641 0.0596 0.0615881311508623 4 24 0.147 0.145287462165307 0.0151 0.0200631340157014 0.0297 0.0264494038189916 5 30 0.1447 0.145967706327343 0.0361 0.0378459854381184 0.0111 0.00798630823453884 6 36 0.1279 0.132131188851013 0.0591 0.0577272070532937 0.0047 0.00194160409569314 7 42 0.1127 0.115034522234822 0.0767 0.0763300773851729 0.0025 0.0004354003800049 8 48 0.1009 0.0989887391692515 0.0898 0.0927163359628322 0.0011 9.49248679162901E-5 9 54 0.0848 0.0848715515617263 0.107 0.106908019877197 0 2.04285610764949E-5 10 60 0.0721 0.0726606420421092 0.1197 0.119134999676661 0 4.35828122956682E-6 |
2楼2019-03-23 20:10:57
3楼2019-03-24 10:18:32
独孤神宇
版主 (知名作家)
- 应助: 490 (硕士)
- 贵宾: 0.008
- 金币: 31016.3
- 散金: 802
- 红花: 122
- 沙发: 1
- 帖子: 5600
- 在线: 856.5小时
- 虫号: 3522474
- 注册: 2014-11-06
- 性别: GG
- 专业: 机械动力学
- 管辖: 计算模拟
【答案】应助回帖
|
这是用1stopt拟合的，采用通用全局优化算法，至于每次结果有些许差别那是正常的，只要总体趋势满足，拟合的统计参数符合要求就可以了。拟合只是寻找理想的最优解，结果往往只有更优，没有最优。 发自小木虫Android客户端 |
» 本帖已获得的红花（最新10朵）
4楼2019-03-24 13:03:01
10