小木虫

当前位置：首页 > 计算模拟 >matlab拟合反应动力学参数

matlab拟合反应动力学参数

作者苏拉: 来源: 小木虫 550 11 举报帖子

各位大佬，可以帮我拟合一下动力学参数吗，可以有偿。先谢过各位大佬。

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 返回小木虫查看更多

今日热帖

精华评论

独孤神宇

均方差(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楼: Originally posted by 独孤神宇 at 2019-03-23 20:10:57
均方差(RMSE): 0.00213787148383087
残差平方和(SSR): 0.000137114834441316
相关系数(R): 0.998778073301579
相关系数之平方(R^2): 0.997557639708014
修正R平方(Adj. R^2): 0.995605193318058
确定系数(DC): ...

你好，这里求出来的是全局最优解还是局部最优解，因为我用不同的初值拟合也能收敛，只是每次求出来的参数可能都不一样如(k1 = 17.932955894914， k2 = 801.597983738672, kA = 8.49463586176274, KB = 0.0324248271839244, kc = 2.11449075468859, kH2 = 0.537868687305816）。另外，可以加上带5%误差分布的拟合参数吗，是用MATLAB计算的吗，可以的话还想求一份计算的代码，麻烦你了，我可以向你付费咨询的。
瑞aaa

在拟合过程需要r1和r2吗
FMStation

https://ww2.mathworks.cn/help/matlab/ref/polyfit.html

https://ww2.mathworks.cn/products/curvefitting/features.html

https://zhuanlan.zhihu.com/p/21754490
苏拉

引用回帖:
7楼: Originally posted by FMStation at 2019-03-24 18:22:24
https://ww2.mathworks.cn/help/matlab/ref/polyfit.html

https://ww2.mathworks.cn/products/curvefitting/features.html

https://zhuanlan.zhihu.com/p/21754490

不是这样的，目前需要尽快出拟合结果，我用其他方法也算过了，但怀疑自己算的精度不够，想再用MATLAB计算验证一下，从头学MATLAB时间比较急，也知道当伸手党不好，所以提出了可以有偿求助，不知我是否会错了意。