| 查看: 1172 | 回复: 0 | ||
pangteng1990金虫 (正式写手)
|
[求助]
浓度一时间比法确定反应级数和计算速率常数 附带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
- 附件 2 : fanyingjishu.m
2015-06-16 07:50:38, 143.42 K
2015-06-16 07:52:29, 1.33 K
» 猜你喜欢
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组招收2026硕士研究生1名
已经有0人回复
-
山西大同大学物理学硕士研究生招收调剂生
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有201人回复
-
[调剂信息]211智能人工感知方向国家青年特聘专家课题组招收调剂研究生
已经有0人回复
-
［调剂信息］211智能人工感知方向国家青年特聘专家课题组招收调剂研究生
已经有0人回复
-
[调剂信息]211智能人工感知方向国家青年特聘专家课题组招收调剂研究生
已经有0人回复
-
一志愿物理所凝聚态物理352求调剂
已经有0人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有2人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有0人回复
-
070200求调剂，一志愿某211，288分
已经有17人回复
找到一些相关的精华帖子，希望有用哦~
MATLAB应用四阶龙格库塔法求解反应动力学问题?
已经有21人回复
- 如何求这个matlab程序画出的legend曲线的法矢量？
已经有6人回复
- 微分方程ode45求解，最小二乘法优化微分方程参数，程序运行求助
已经有5人回复
- Matlab同时拟合活化能，反应级数，速率常数
已经有26人回复
- 求用牛顿迭代法解非线性方程组的matlab程序
已经有4人回复
- MATLAB无法计算出这个含【级数和超几何函数】的式子
已经有6人回复
- 求助用matlab三次样条法算一个反应的反应速率
已经有13人回复
- 求助一个最小二乘法matlab拟合程序。
已经有7人回复
- 求助啊！求一程序，用matlab程序做，用蒙特卡洛方法模拟
已经有10人回复
- 【求助】催化反应动力学matlab计算各基元反应的速率常数时，该如何避免较小量被忽略？
已经有3人回复
- 【求助】MATLAB编写最大熵法计算概率密度的程序
已经有8人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
