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可以让我学习一下（我没搜到啊

）。谢谢

回复此楼

» 猜你喜欢

求取一些关于纳米材料和纳米技术相关的英文PPT。已经有0人回复
【复旦大学】二维材料方向招收2026年博士研究生1名已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有136人回复
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组2026级硕/博/博后招生已经有10人回复
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生已经有22人回复
反铁磁体中的磁性切换：两种不同的机制已成功可视化已经有0人回复
26申博推荐：南京航空航天大学国际前沿科学研究院光学方向招收博士生！已经有0人回复
求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复
新西兰Robinson研究所招收全奖PhD 已经有0人回复

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

动力学反应级数已经有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, 给个红包，谢谢回帖

引用回帖:

18楼: Originally posted by whqs8426212 at 2014-08-27 11:33:16
function k1k2k3k4
format long
clear all
clc
tspan = ;
C0 = ; %初值只有一个C0
k0 = ; %k0 E M n 参数的初值
lb = ; %上届
%ub = ; %下界

dataTP1=...
; %实验数据1
dataTP2=.. ...

能否帮我指导下，乳酸与正丁醇酯化反应的反应级数和反应速率常数
数据如代码

function KineticsEst1_int11
% 动力学参数辨识: 用积分法进行反应速率分析得到速率常数k和反应级数n
% Analysis of kinetic rate data by using the integral method
%
% Author: HUANG Huajiang
% Copyright 2003 UNILAB Research Center,
% East China University of Science and Technology, Shanghai, PRC
% $Revision: 1.0 $  $Date: 2003/07/27 $
%
% Reaction of the type -- rate = kCA^order
% order - reaction order
% rate -- reaction rate vector
% CA -- concentration vector for reactant A
% T -- vector of reaction time
% N -- number of data points
% k- reacion rate constant

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);

这个代码我照着书上弄的有错误