24小时热门版块排行榜    

查看: 794  |  回复: 4
【悬赏金币】回答本帖问题,作者臣本布衣将赠送您 300 个金币

臣本布衣

新虫 (初入文坛)

[求助] 有偿求助,用Matlab的Levenberg−Marquardt计算反应速率常数 已有2人参与

有偿求助,用Matlab的Levenberg−Marquardt计算反应速率常数。
回复此楼
爱拼才会赢
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hzlhm

至尊木虫 (著名写手)

有问题可以先贴出来看一下,或发我私信

发自小木虫Android客户端
QQ:2120156492
2楼2022-05-25 19:37:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

臣本布衣

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by hzlhm at 2022-05-25 19:37:17
有问题可以先贴出来看一下,或发我私信

动力学方程为(L-H方程):
dCa/dt=-(k1+k3)KaCa/(1+KaCa+KbCb+KcCc)+k2KbCb/(1+KaCa+KbCb+KcCc)
dCb/dt=-k1KaCa/(1+KaCa+KbCb+KcCc)-(k2+k4)KbCb/(1+KaCa+KbCb+KcCc)
dCc/dt=-k3KaCa/(1+KaCa+KbCb+KcCc)+k4KbCb/(1+KaCa+KbCb+KcCc)

求反应速率常数k1,k2,k3,k4和吸附平衡常数Ka,Kb,Kc

实验数据:
t                        Ca                     Cb                  Cc
0                   0.084691                 0                   0
0.008333            0                0.065164         0.011266
0.01                    0                0.065409         0.012565
0.0125                0                0.060309         0.015934
0.016667            0                0.008297         0.067477
0.025                  0                0.005952         0.063515
0.05                    0                         0             0.0567

文献(J. Phys. Chem. C 2014, 118, 27961−27972)采用LM算法进行优化。

MATLAB小白,求大神赐教!
爱拼才会赢
3楼2022-05-30 15:09:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

引用回帖:
3楼: Originally posted by 臣本布衣 at 2022-05-30 15:09:06
动力学方程为(L-H方程):
dCa/dt=-(k1+k3)KaCa/(1+KaCa+KbCb+KcCc)+k2KbCb/(1+KaCa+KbCb+KcCc)
dCb/dt=-k1KaCa/(1+KaCa+KbCb+KcCc)-(k2+k4)KbCb/(1+KaCa+KbCb+KcCc)
dCc/dt=-k3KaCa/(1+KaCa+KbCb+KcCc)+k4KbCb/ ...

这种微分方程拟合问题还没有比1stOpt更方便强大的了吧,建议试试。
另外Ca数据没问题吗?起始值为0.084691,其余全为0,与模型也不符合吧
4楼2022-06-01 07:56:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

du660283

禁虫 (小有名气)

本帖内容被屏蔽

5楼2022-06-04 06:13:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 臣本布衣 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 本人考085602 化学工程 专硕 +7 不知道叫什么! 2026-03-15 8/400 2026-03-15 20:11 by 棒棒球手
[考研] 290求调剂 +3 孔志浩 2026-03-12 8/400 2026-03-15 15:30 by 孔志浩
[考研] 294求调剂 +3 Zys010410@ 2026-03-13 4/200 2026-03-15 10:59 by zhq0425
[考研] 311求调剂 +3 26研0 2026-03-15 3/150 2026-03-15 09:12 by JourneyLucky
[考研] 288求调剂 +4 奇点0314 2026-03-14 4/200 2026-03-14 23:04 by JourneyLucky
[考研] 本科南京大学一志愿川大药学327 +3 麦田耕者 2026-03-14 3/150 2026-03-14 20:04 by 外星文明
[考研] 289求调剂 +4 这么名字咋样 2026-03-14 6/300 2026-03-14 18:58 by userper
[考研] 中科大材料专硕319求调剂 +3 孟鑫材料 2026-03-13 3/150 2026-03-14 18:10 by houyaoxu
[考研] 297一志愿上交085600求调剂 +5 指尖八千里 2026-03-14 5/250 2026-03-14 17:26 by a不易
[考研] 复试调剂 +3 呼呼?~+123456 2026-03-14 3/150 2026-03-14 16:53 by WTUChen
[考研] 310 070300化学求调剂 +5 扑风铃的猫 2026-03-08 6/300 2026-03-14 03:22 by JourneyLucky
[考研] 0703化学求调剂 +5 很老实人 2026-03-09 5/250 2026-03-14 02:57 by JourneyLucky
[考研] 一志愿天大化工(085600)调剂总分338 +6 蔡大美女 2026-03-09 6/300 2026-03-14 02:46 by JourneyLucky
[考研] 307求调剂 +7 超级伊昂大王 2026-03-10 7/350 2026-03-14 00:49 by JourneyLucky
[考研] 0703化学调剂 +4 快乐的香蕉 2026-03-11 4/200 2026-03-13 22:41 by JourneyLucky
[考研] 泣血叩求调剂恩,愿以丹心报师恩 +6 Iuruoh 2026-03-11 6/300 2026-03-13 22:06 by JourneyLucky
[考研] 0703化学一志愿211 总分320求调剂 +5 玛卡巴卡啊哈 2026-03-11 5/250 2026-03-13 21:40 by JourneyLucky
[考研] 工科调剂 +4 Jiang191123! 2026-03-11 4/200 2026-03-13 15:15 by Miko19
[考研] 308求调剂 +3 是Lupa啊 2026-03-12 3/150 2026-03-13 14:30 by 求调剂zz
[考研] 0817化学工程与技术考研312分调剂 +3 T123 tt 2026-03-12 3/150 2026-03-13 10:49 by houyaoxu
信息提示
请填处理意见