| 查看: 636 | 回复: 0 | ||
pangteng1990金虫 (正式写手)
|
[求助]
浓度一时间比法确定反应级数和计算速率常数 MATLAB程序出错
|
|
我按照文献中给出的程序运行的，可是总是出错，就是原始数据的输入，不知道怎么弄，老是出错，求大神帮忙 % 数据输入 edit t.m; %输入时间数据 edit c.m %输入浓度数据 %数据处理程序 fyjs.m load t.m; %调入时间数据 load c.m %调入浓度数据 N0=max(size(t)); C0=c(1,1); x=[]; tb=[]; x1=[]; for i=2:1:N0 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); %计算测定数据中毎两组数据的时间比 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:N0-1 if n==1 p=-log(1-x(1,k))./t(1,k); %计算反应级数n=1时的速率常数 else p=C0.^(1-n).*((1-x(1,k)).^(1-n)-1)./(n-1)./t(1,k); %计算反应级数n≠1时的速率常数 kc=[kc,p]; end end kc=sum(kc)./(N-1); %计算速率常数的平均值 T/min 0 120 190 250 312 370 500 c/g/L 16.9250 16.5639 16.4806 16.2028 16.0917 15.7306 15.4250 |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : 浓度_时间比法确定反应级数和计算速率常数_粟智.caj
- 附件 2 : ctb.m
2015-06-15 10:48:49, 143.42 K
2015-06-15 10:55:32, 1.32 K
» 猜你喜欢
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有4人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有5人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有7人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有7人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有7人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有6人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有8人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有10人回复
找到一些相关的精华帖子，希望有用哦~
matlab实现牛顿迭代法求非线性方程组程序运行出现问题
已经有8人回复
- MATLAB应用四阶龙格库塔法求解反应动力学问题?
已经有21人回复
- 如何求这个matlab程序画出的legend曲线的法矢量？
已经有6人回复
- 微分方程ode45求解，最小二乘法优化微分方程参数，程序运行求助
已经有5人回复
- MATLAB最大值输法
已经有5人回复
- Matlab同时拟合活化能，反应级数，速率常数
已经有26人回复
- 如何用最小二乘法及MATLAB解这个方程组？
已经有6人回复
- MATLAB 奇异值分解法svd
已经有3人回复
- 电磁场数值计算法与MATLAB实现
已经有307人回复
- 求用牛顿迭代法解非线性方程组的matlab程序
已经有4人回复
- 【课件】最优化方法及其Matlab程序设计
已经有390人回复
- MATLAB无法计算出这个含【级数和超几何函数】的式子
已经有6人回复
- 求助利用补偿最小二乘法辨识半参数模型的matlab程序
已经有3人回复
- 求助用matlab三次样条法算一个反应的反应速率
已经有13人回复
- 结构分析的有限元法与MATLAB 程序设计(徐荣桥).pdf
已经有321人回复
- 求助注水算法功率分配matlab程序和最小二乘法matlab程序
已经有4人回复
- 求助一个最小二乘法matlab拟合程序。
已经有7人回复
- 使用matlab和有限元差分法计算无限深方势?S的第一波函数
已经有8人回复
- matlab插值法的一个问题
已经有6人回复
- 求大神帮我改一下下面的程序！matlab迭代法的 不知道错在哪里！
已经有7人回复
- 【求助】催化反应动力学matlab计算各基元反应的速率常数时，该如何避免较小量被忽略？
已经有3人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
30