用普朗尼法估算化学动力学参数及程序设计.pdfVIP

第8卷 第4期 《计算机 与应用化学 Vol8．Ⅳo4 991年 】0月 Computers and AppLied Chemistry 0ctober l99 用普朗尼法估算化学动力学 参数及程序设计 姚 小 民 (西北工业大学) 摘 要 本文提 出了用普 朗尼 (Prony)法 ，根据离散数据表估算线性反应系统动力学参数 的 方法 并介绍了该方法的程序设计，本文方法简单 明确．计算复杂度较低，程序通用性强。 关键词 ：昔 朗尼法 动力学参数 程序设计 许多有工业意义的反应都可近似地看作一级或拟一级反应 … 。因而，线性常微分方 程组的参数估值在化学反应工程等方面具有重要 的意义 目前 ，微分方程法的解法有 ：曲 线拟合法 ，特征向量法 ，优化方法 ，拉普拉斯方法 等 。 本文提 出利用普 朗尼 (Prony)方法估算动力学参数的方法，本文方法简单、明确 ， 较之其它方法，其计算复杂度较低 ，如：曲线拟合法在求解时，需要选取适 当的权 函数， 并要对实验数据进行数值积分处理：特征向量法，其公式繁琐，迭代搜索过程复杂；而拉 普拉斯方法，在处理 由于拖尾积分 引起 的误差时，需要写 出常微分方程组的解析表达式， 这对于阶数较高的系统 ，会带来一定量的数学推导 。而用普朗尼方法估算动力学参数，整 个求解过程 由一系列求线性方程组的解和多项式的根组成，便于在计算机上实现，且程序 通用性强 二、计算方法 一 级反应系统的常微分方程组 的表示形式为： dX (t) — - ： KX() (1) o f 其 中 f： 白变量，可表示反应时间、反应器体积或长度等； fJ_N维 向量 ，X (f) = (，()，x ，()，…，h (f)) ，可表示各组分 的浓 度、分压或转化率等； ：Ⅳ阶矩阵，表示参数矩 阵 根据文献 3【L 参数矩阵 可表示为 K =BAB一 (2) 荩 中 ，，屯… 为系统的特征值，八为 i．．一， 构成的对角矩阵， 为 ／＼阵的特 奉文 f1990年 10月23n收到 第 4期 姚小民：用普朗尼法估算化学动力学参数及程序设计 ’279 征向量矩阵 由 (2)式可知，已知 与八，则 可唯一确定 2．1利用普 朗尼方法计算 ／、阵 普 朗尼方法 ，是在瞬变 电磁场 的研究 中，为从 目标 的瞬态 响应 中提取极点和计算 留数而提供 的一种行之有效的方法 它主要是把具有衰减振荡特性 的时间响应数据 (计算 的或测量的)用复指数和 的函数形式逼近 对于化学反应过程 中的线性反应系统，因参与 反应的各组分 的浓度 ，分压或转化率等具有某种衰减特性 ，则可表示成它们与时间数据 的 实指数和的函数形式 。因此，通过普 朗尼法可求 出反应系统的特征值 设系统的特征值有 Ill^l… liI，由于反应系统中的置 (f)(l，2， “，M都具有相同的解结 构，即相 同的 ^】，^2，…， ，则我们可 以从一个 t (f)的数据表来确定 zI，i2，…， 取某 一个 (f) 为 I (t)，设 ，(f)=∑A e