24小时热门版块排行榜    

CyRhmU.jpeg
查看: 691  |  回复: 4
【悬赏金币】回答本帖问题,作者臣本布衣将赠送您 300 个金币
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

臣本布衣

新虫 (初入文坛)

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

有偿求助,用Matlab的Levenberg−Marquardt计算反应速率常数。
回复此楼
爱拼才会赢
已阅   回复此楼   关注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的回帖
查看全部 5 个回答

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的回帖

du660283

禁虫 (小有名气)

本帖内容被屏蔽

5楼2022-06-04 06:13:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
信息提示
请填处理意见