24小时热门版块排行榜

返回列表

pangteng1990

金虫 (正式写手)

应助: 4 (幼儿园)
金币: 695.8
散金: 19
红花: 2
帖子: 312
在线: 72.3小时
虫号: 2328422
注册: 2013-03-08
性别: GG
专业: 无机化工

[求助] 浓度一时间比法确定反应级数和计算速率常数附带MATLAB程序出错，求大神教我

edit t.m %输入时间数据
edit c.m %输入浓度数据、
%数据处理程序 fyjs.m
load t.m %调入时间数据
load c.m %调入浓度数据
NO=max(size(t));
CO=c(1,1);
x=[];
tb=[];
for i=2:1:NO
a=(c(1,1)-c(1,i))./c(1,1); %将浓度数据转化为转化率
x=[x,a];
end
t(1:1)=[];
N=max(size(t));
sn=[];
for i=1:1:N-1
for j=i+1:1:N
      b=t(1,j)./t(1,i); %计算测定数据中毎两组数据的时间比
      tb=[tb,b];
end
end
st=sum(tb); %计算总时间比
tn=[];
for n=-3:0.001:3; %逐步逼近搜索，从-3级到3级，步长为0.001
xn=[];
for i=1:1:N-1
      for j=i+1:1:N
         if n==1
            w=log(1-x(1,j))./log(1-x(1,i)); %反应级数n=1时的浓度比
         else
            w=((1-x(1,j)).^(1-n)-1)./((1-x(1,i)).^(1-n)-1); %反应级数n≠1时的浓度比
         end
         xn=[xn,w];
      end
end
tn=[tn,n];
sn=[sn,sum(xn)]; %建立不同级数下总浓度比集合
end
d=[];
for k=1:1:max(size(sn))
d=abs([d,sn(1,k)-st]); %建立总浓度比和总时间比的差的集合
end
for k=1:1:max(size(sn))
if abs(sn(1,k)-st)==min(d); %搜索最小差值的位置
      n=tn(1,k); %记录最小差值的反应级数
end
end
kc=[];
for k=1:1:NO-1
if n==1
      p=-log(1-x(1,k))./t(1,k); %计算反应级数n=1时的速率常数
else
      p=CO.^(1-n).*((1-x(1,k)).^(1-n)-1)./(n-1)./t(1,k); %计算反应级数n≠1时的速率常数
      kc=[kc,p];
end
end
k=sum(kc)./(N-1); %计算速率常数的平均值
n
k

回复此楼

» 本帖附件资源列表

欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
附件 1 : 浓度_时间比法确定反应级数和计算速率常数_粟智.caj

2015-06-16 07:50:38, 143.42 K

附件 2 : fanyingjishu.m

2015-06-16 07:52:29, 1.33 K

» 猜你喜欢

澳门大学IAPME李海峰教授常年招收博士研究生已经有0人回复
澳门大学IAPME李海峰教授常年招收博士研究生已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有287人回复
深圳大学应用技术学院招聘凝聚态物理博士后已经有123人回复
求助中医院临床中医中标的青年标书学习，谢谢已经有5人回复
浙大物院表面物理课题组长期招收普博及直博学生已经有259人回复
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry 已经有47人回复
25届硕士申请，材料计算方向已经有2人回复
课题组招1-2名硕士已经有1人回复
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上已经有0人回复

一个人的武林

1楼 2015-06-16 07:51:47

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

智能机器人

Robot (super robot)

我们都爱小木虫

找到一些相关的精华帖子，希望有用哦~

MATLAB应用四阶龙格库塔法求解反应动力学问题? 已经有21人回复
如何求这个matlab程序画出的legend曲线的法矢量？已经有6人回复
微分方程ode45求解，最小二乘法优化微分方程参数，程序运行求助已经有5人回复
Matlab同时拟合活化能，反应级数，速率常数已经有26人回复
求用牛顿迭代法解非线性方程组的matlab程序已经有4人回复
MATLAB无法计算出这个含【级数和超几何函数】的式子已经有6人回复
求助用matlab三次样条法算一个反应的反应速率已经有13人回复
求助一个最小二乘法matlab拟合程序。已经有7人回复
求助啊！求一程序，用matlab程序做，用蒙特卡洛方法模拟已经有10人回复
【求助】催化反应动力学matlab计算各基元反应的速率常数时，该如何避免较小量被忽略？已经有3人回复
【求助】MATLAB编写最大熵法计算概率密度的程序已经有8人回复

点击这里搜索更多相关资源

科研从小木虫开始，人人为我，我为人人

相关版块跳转我要订阅楼主 pangteng1990 的主题更新