应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 程序代码 » matlab拟合反应动力学参数
121/2返回列表12下一页
查看: 2184  |  回复: 11
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

苏拉

铁虫 (小有名气)

[求助] 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
回复此楼

» 猜你喜欢
1楼 2019-03-23 17:23:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

独孤神宇

版主 (知名作家)

【答案】应助回帖

感谢参与，应助指数 +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
赞一下(1人) 回复此楼
数值计算
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的回帖

瑞aaa

新虫 (正式写手)
在拟合过程需要r1和r2吗

发自小木虫Android客户端
赞一下 回复此楼
5楼2019-03-24 13:19:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

苏拉

铁虫 (小有名气)
引用回帖:
5楼: Originally posted by 瑞aaa at 2019-03-24 13:19:11
在拟合过程需要r1和r2吗

r1和r2用浓度的微分dC/dt表示

发自小木虫Android客户端
赞一下 回复此楼
6楼2019-03-24 16:17:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

FMStation

至尊木虫 (知名作家)

【答案】应助回帖

感谢参与，应助指数 +1
https://ww2.mathworks.cn/help/matlab/ref/polyfit.html

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

https://zhuanlan.zhihu.com/p/21754490
赞一下 回复此楼
7楼2019-03-24 18:22:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

苏拉

铁虫 (小有名气)
引用回帖:
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时间比较急，也知道当伸手党不好，所以提出了可以有偿求助，不知我是否会错了意。
赞一下 回复此楼
8楼2019-03-25 09:51:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

苏拉

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

...

你好，可以私信发一份1stopt代码吗，我还有几组数据想自己算一下，麻烦你了。
赞一下 回复此楼
9楼2019-03-25 11:04:13
已阅   回复此楼   关注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的回帖
相关版块跳转 我要订阅楼主 苏拉 的主题更新
121/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 293求调剂 +3 是乐渝哇 2026-03-04 3/150 2026-03-04 23:12 by wutongshun
[考研] 274求调剂 +8 一个学习者 2026-03-04 8/400 2026-03-04 20:35 by kakakapanpan
[考研] 302材料工程求调剂 +6 Doleres 2026-03-01 7/350 2026-03-04 20:09 by caszguilin
[考研] 学硕材料275调剂 +9 路三三 2026-03-03 9/450 2026-03-04 17:02 by 梦天888
[考研] 0703化学 一志愿东北大学 279分求调剂 +4 聂乐平 2026-03-03 4/200 2026-03-04 16:31 by zhukairuo
[考研] 281求调剂 +3 Y?l?h 2026-03-04 3/150 2026-03-04 10:51 by 汪！？！
[考研] 一志愿西工大材料学，英一数二总分321分，求调剂。 +3 zz05zz 2026-03-02 5/250 2026-03-04 09:16 by lbsjt
[考研] 0854总分272 +5 打小就是老实人 2026-03-02 6/300 2026-03-04 01:41 by ouhaiyu
[考研] 0703化学学硕 288分 +4 旅行中的紫葡萄 2026-03-03 4/200 2026-03-03 23:41 by ch1979
[考研] 环境工程学硕288求助调剂 +7 多吃亿口芝士 2026-03-02 7/350 2026-03-03 23:00 by zhukairuo
[考研] 306求调剂 +7 张张张张oo 2026-03-03 7/350 2026-03-03 21:35 by L135790
[考研] 0856材料工程，初试313调剂 +7 卖个关子吧 2026-03-03 7/350 2026-03-03 21:33 by L135790
[考研] 化工专硕调剂 +4 利好利好. 2026-03-03 7/350 2026-03-03 21:30 by L135790
[考研] 276求调剂 +8 路lyh123 2026-02-28 10/500 2026-03-03 18:25 by xin吖
[考研] 求调剂 +4 Guo_yuxuan 2026-03-02 5/250 2026-03-03 14:39 by xiaomc_gzh
[考研] 284求调剂 +10 天下熯 2026-02-28 11/550 2026-03-02 11:03 by 无际的草原
[考研] 275求调剂 +3 L-xin? 2026-03-01 6/300 2026-03-02 10:22 by 热情沙漠
[考研] 279求调剂 +3 dua1 2026-03-01 4/200 2026-03-02 00:23 by 大脸蛋子
[考研] 311求调剂 +6 亭亭亭01 2026-03-01 6/300 2026-03-01 15:41 by 324616
[考研] 317一志愿华南理工电气工程求调剂 +6 Soliloquy_Q 2026-02-28 11/550 2026-03-01 11:14 by 歌liekkas
信息提示
关闭
请填处理意见
关闭