24小时热门版块排行榜    

Znn3bq.jpeg
查看: 2326  |  回复: 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
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

苏拉

铁虫 (小有名气)

引用回帖:
4楼: Originally posted by 独孤神宇 at 2019-03-24 13:03:01
这是用1stopt拟合的,采用通用全局优化算法,至于每次结果有些许差别那是正常的,只要总体趋势满足,拟合的统计参数符合要求就可以了。拟合只是寻找理想的最优解,结果往往只有更优,没有最优。

...

大哥,数据和模型发了,你抽空帮我用程序跑一下吧,非常感激。

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : Desktop.rar
  • 2019-03-25 15:59:38, 24.88 K
10楼2019-03-25 16:03:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 12 个回答

独孤神宇

版主 (知名作家)

【答案】应助回帖

感谢参与,应助指数 +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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

苏拉

铁虫 (小有名气)

引用回帖:
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计算的吗,可以的话还想求一份计算的代码,麻烦你了,我可以向你付费咨询的。
3楼2019-03-24 10:18:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (知名作家)

【答案】应助回帖

引用回帖:
3楼: originally posted by 苏拉 at 2019-03-24 10:18:32
你好,这里求出来的是全局最优解还是局部最优解,因为我用不同的初值拟合也能收敛,只是每次求出来的参数可能都不一样如(k1 = 17.932955894914, k2 = 801.597983738672, ka = 8.49463586176274, kb = 0.032424827 ...

这是用1stopt拟合的,采用通用全局优化算法,至于每次结果有些许差别那是正常的,只要总体趋势满足,拟合的统计参数符合要求就可以了。拟合只是寻找理想的最优解,结果往往只有更优,没有最优。

发自小木虫Android客户端

» 本帖已获得的红花(最新10朵)

数值计算
4楼2019-03-24 13:03:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 考博 +4 Eternitf 2026-04-29 4/200 2026-05-02 16:10 by yejenny
[找工作] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 qs8kylt1wr 2026-05-01 3/150 2026-05-02 16:06 by zahslv263y
[公派出国] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 qs8kylt1wr 2026-05-01 3/150 2026-05-02 15:51 by zahslv263y
[考博] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 qs8kylt1wr 2026-05-01 3/150 2026-05-02 15:51 by zahslv263y
[教师之家] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 qs8kylt1wr 2026-05-01 3/150 2026-05-02 15:36 by zahslv263y
[论文投稿] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 qs8kylt1wr 2026-05-01 3/150 2026-05-02 15:21 by zahslv263y
[博后之家] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 qs8kylt1wr 2026-05-01 3/150 2026-05-02 15:06 by zahslv263y
[公派出国] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 d8tbypao1q 2026-05-01 3/150 2026-05-02 14:21 by zahslv263y
[教师之家] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 d8tbypao1q 2026-05-01 3/150 2026-05-02 14:06 by zahslv263y
[公派出国] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 k17zww5nna 2026-05-01 4/200 2026-05-02 10:21 by ovidyohpyz
[找工作] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 k17zww5nna 2026-05-01 4/200 2026-05-02 10:06 by ovidyohpyz
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 k17zww5nna 2026-05-01 5/250 2026-05-02 09:21 by dz7y1p1i98
[考博] 药化及相关博士的申请 5+3 中药学渣 2026-04-29 3/150 2026-05-01 18:14 by babero
[论文投稿] 一篇MDPI论文改变了学习工作和生活 +3 helsinglee 2026-04-29 4/200 2026-05-01 11:45 by bobvan
[有机交流] 一个化合物的合成路线:CAS:367929-02-0 名称:8β-乙烯基雌二醇 50+4 zhquanbing 2026-04-28 4/200 2026-04-30 22:45 by zyqchem
[文学芳草园] 太白金星有点烦 +3 yexuqing 2026-04-27 3/150 2026-04-30 21:56 by myrtle
[考博] 中国地质大学(北京)博士招生补录,数理学院材料科学与工程专业和材料与化工专业 +3 liuxh329 2026-04-29 3/150 2026-04-30 15:42 by 顺利毕业zs
[基金申请] 收到国自然专家邀请后几年才会有本子送过来评 +3 barley5 2026-04-30 3/150 2026-04-30 14:58 by aspect3000
[论文投稿] 有没有快的中文核心比较快录用的,纳米材料光催化 5+3 2914439618 2026-04-27 4/200 2026-04-29 08:59 by 北京莱茵润色
[论文投稿] 有人投过CCC中国控制会议吗? 19+3 初九凡 2026-04-26 3/150 2026-04-27 08:58 by 北京莱茵润色
信息提示
请填处理意见