化学反应速率常数和反应级数的计算机求解.pdf

v01．25 第25卷第3期 烟台大学学报(自然科学与工程版) No．3 ofYantai Science{Illd JIll．2012 2012年7月 J0umal univer8ity(NaturalEngineeringEdi“on) 文章编号：l004．882012012)03_0232．03 化学反应速率常数和反应级数的计算机求解 黄雪征，姬鄂豫 (1．南阳理工学院土木工程学院，河南南阳473004；2．南阳理工学院生物与化学工程学院，河南南阳，473004) Basic编程工吴设计了一个计算程序。用于解决复杂的化学反应动力学 摘要：利用Vi蛐al 的反应级数和反应速率常数的计算问题．程序根据线性相关系数R的数值来判断计算反 应速率时的最优的时间变化量出。从而保证了计算结果的准确性，同时，程序具有强大的 绘图功能，所绘图形美观协调，可形象直观地表达出反应物质的浓度变化趋势以及反应速 率的对数值和物质浓度的对数值的线性拟舍程度．利用一个反应实例对程序进行了验证， 计算结果和文献参考值吻合良好． 关键词：化学反应速率常数；反应级数；光滑不等距插值；差商；一元线性回归 中图分类号：0643．11；71日12文献标志码：A 在反应动力学实验数据处理中，常用的方法有 浓度保持不变，作为常数并入速率常数中，则式(1) 微分法、积分法、半衰期法等．文献[1—3]分别采用可简化为 微分法和积分法进行了编程求解，但大部分程序没 移=一孥=kc：． t上’(2) 有绘图功能，部分程序在求解反应速率时计算步长 移2一百2‰“· 固定，对于某些时间跨度的实验数据计算误差往往 若有厅组实验数据，时间点分别为f．，如，…，￡。，采用 较大．本文针对上述情况，利用VB编程工具，设计 化学分析法，测定对应时刻A物质的浓度分别为 了一个交互性强、具有作图功能、且反应速率的计算 c，，乞，…，c。，可由光滑不等距插值公式H。计算A物 质的浓度． 步长可根据相关系数R来自行调整以寻求最优计 算结果的软件． 根据反应速率的定义和数值分析的知识D】，可 由下式求出各时间点A物质的反应速率． 1 化学反应动力学实验数据的计算原 功2一——————药r一——’ 理 ‘ 铲一趔亲趔，(i：2，3，…，n一1) 1．1化学分析法实验数据的计算 设有反应口A+6B+…呻P+mM+…。则有如q。一荪■—一，L‘2二’’，…’n．1， 霄n=一 下反应速率方程 口=_|}c：吐…． (1) (3) 改变其中一种组分如A物质的浓度，其他组分 式(3)中c。(屯+山)代表￡；+dt时刻A物质的浓度， 收稿日期：20ll嘲捌 基金项目：河南省科技厅项目(0524”∞15)， 作者简介：黄雪征(1978一)，河南南阳人，讲师，硕士，研究方向：化学过程的计算及计算机仿真． 万方数据 第3期 黄雪征，等：化学反应速率常数和反应级数的计算机求解 233 用物理化