24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

heisekaishui

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 215
帖子: 52
在线: 59.2小时
虫号: 2445260

[交流] Matlab同时拟合活化能，反应级数，速率常数

各位，您好!
本人研究的方向是催化化学，最近在评估催化剂效果时想引入动力学的内容，对其中的一些地方有一些疑问。
受月之蓝版主的这个帖子http://muchong.com/bbs/viewthread.php?tid=6425538&target=self&page=1以及其它前辈的启发，对于反应级数以及反应速率常数的拟合有一定的了解。
但是如果我的动力学模型是这种形式：-dC/dt=k0*exp(-E/R*T)*C^m*P^n（只有一个反应物和氢气），其中P与T分别为压力与时间，就是说是一个加氢反应，试验中可以得到不同压力、不同温度下浓度随时间的变化值，应该如何拟合k0,E,m,n。
我的一个想法是先拟合不同压力下的k1（k1=k0*exp(-E/R*T)*P^n)（例如以上述帖子的代码方法），再根据k1=k2*P^n（k2=k0*exp(-E/R*T)）拟合出k2和n，此后拟合出不同温度下的k2，根据k2=k0*exp(-E/R*T)）拟合出E和k0
但是我觉得这种方法较为复杂，一个matlab里有好几个拟合程序，不知道是否有更简洁的CODE。

PS:不知道版里有没有相似的帖子有类似的CODE可以让我学习一下（我没搜到啊

）。谢谢

回复此楼

» 猜你喜欢

散金币祈福已经有187人回复
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布已经有1人回复
物理学I论文润色/翻译怎么收费? 已经有71人回复
澳门大学IAPME李海峰教授常年招收博士研究生已经有0人回复
澳门大学IAPME李海峰教授常年招收博士研究生已经有0人回复
澳门大学IAPME李海峰教授常年招收博士研究生已经有0人回复
深圳大学应用技术学院招聘凝聚态物理博士后已经有115人回复
深圳大学2027级光电信息工程/物理学研究生推免已经有70人回复
求助中医院临床中医中标的青年标书学习，谢谢已经有4人回复
浙大物院表面物理课题组长期招收普博及直博学生已经有257人回复
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry 已经有47人回复

» 本主题相关价值贴推荐，对您同样有帮助:

动力学反应级数已经有4人回复
求助用matlab三次样条法算一个反应的反应速率已经有13人回复
求助用matlab拟合动力学方程已经有9人回复
速率方程的问题求教！已经有4人回复
有关反应速率常数的计算已经有9人回复
SCR反应动力学方程求解的问题已经有4人回复
催化剂增大了速率常数但没有改变活化能-来自一篇经典文献，如何理解？已经有3人回复
动力学参数拟合已经有26人回复
matlab 拟合反应动力学参数结果很差。大家帮忙看一下已经有14人回复
根据活化能怎么判断一级反应还是零级反应已经有12人回复
焓变活化能键离解能能垒反应速率的巩固求教已经有14人回复
帮我看看聚合反应三个速率常数怎么算？已经有9人回复
【求助】拟合动力学方程求助已经有13人回复
【求助】催化反应动力学matlab计算各基元反应的速率常数时，该如何避免较小量被忽略？已经有3人回复
【求助】反应速率常数怎么算已经有5人回复
【求助】关于催化反应中的活化能已经有24人回复
【求助】使用Matlab拟合反应动力学方程问题已经有7人回复
【求助】使用Matlab预估动力学方程问题已经有13人回复

1楼 2014-08-25 17:30:13

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

guanwei

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 1075.1
帖子: 247
在线: 53.7小时
虫号: 1638561

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

引用回帖:

12楼: Originally posted by whqs8426212 at 2014-08-26 00:14:30
楼主方不方便给我一些数据，如：

某T，P（应该是多个吧？）
t（/s） C(mol)
0          XX
1          XX
2          XX

function KineticsEst1_int11

clear all
clc

t=[0 20 40 100 270 450 750 1410 1590 1810];%时间
rA=[0.0004765 0.0004507  0.0004250  0.0003477 0.001549 0.0001117 0.0000396 0.0000018  0.0000038 0.0000062];%反应速率
Pa=[0.042729739 0.041609737 0.040801805 0.038526838 0.033412812 0.030519156 0.028695168 0.02792077 0.027918742 0.0278887981]; %乳酸的摩尔分率
Pb=[0.123670897  0.122550895  0.121742963 0.119467996 0.11435397  0.111460314  0.109636326  0.108861928 0.1088599 0.108829956];%正T醇的摩尔分率
Pc=[0.833599364  0.834719366 0.835527299 0.837802265 0.842916291 0.845809947 0.847633933 0.848408333 0.848410361 0.848440305];%水的摩尔分率
Pd=[0.001120002 0.001927934 0.004202901 0.009316927 0.012210583 0.014034571 0.014808969 0.014810996  0.014840941  0.014850941];%乳酸正丁酯的摩尔分率
%线性拟合
P=2.337132745*Pa.*Pb-0.564655417*Pc.*Pd;y=rA';X=[ones(size(y))  P'];

b=X\y;k=b(2);
%非线性拟合
beta0=[k]
[beta,resnorm,residual,exitflag,output,lambda,jacobian] = ...
lsqnonlin(@ObjFunc,beta0,[],[],[],rA,Pa,Pb,Pc,Pd)

%拟合效果图(实验与拟合的比较)
figure(1);plot(t,rA,'.')
r_poly=beta(1)*(2.337132745*Pa.*Pb-0.564655417*Pc.*Pd);
hold on;plot(t, r_poly,'g')
figure(2);plot(Pb,rA,'.'); hold on;plot(Pb, r_poly,'g')

%-----------------------------------------------------------------------------一
Function f=ObjFunc(beta,rA,Pa,Pb,Pc,Pd)
f=rA-beta(1)*(2.337132745*Pa.*Pb-0.564655417*Pc.*Pd);

我用的这个程序拟合反应速率常数为什么有问题的呢