数据处理与数值分析论文管理.docxVIP

数据处理与数值计算结课论文学 院 班 级 年 级 姓 名 指导教师 2014年12月29日摘 要化学反应工程主要研究反应动力学和反应器，反应动力学是反应器设计与分析的基础。工业反应器中除了化学反应外，还伴随有反应组分的传递、反应热的吸收或排放、反应物料的流动与混合等物理过程，极其错综复杂，其研究方法以数学模拟为主。现在在解决化学反应工程实际问题的过程中常常遇到对数学模型求解和学学模型参数估值问题。常规方法难以或不能对其进行处理，通过数据处理与数值计算这门课程的学习，我基本掌握了插值，数值微分，数值积分，最小二乘法，（非）线性方程（组）解法，常微分方程（组）的初值和边值问题。并且对计算机软件在化工生产问题中的应用有所了解，现结合具体实例介绍如下。利用最小二乘法确定动力学方程环氧乙烷可杀灭细菌（及其内孢子）、霉菌及真菌，因此可用于消毒一些不能耐受高温消毒的物品。主要用于制造其他各种溶剂(如溶纤剂等),稀释剂,非离子型表面活性剂,合成洗涤剂、抗冻剂、消毒剂、增韧剂和增塑剂等。与纤维素发生羟乙基化可合成得水溶性树脂(其环氧乙烷含量约75%)。还可用作熏蒸剂、涂料增稠剂、乳化剂、胶黏剂和纸张上浆剂等。乙烯的深度氧化是乙烯环氧化生产环氧乙烷的工业生产中重要的副反应。对反应选择性有较大影响。在铂催化剂上，乙烯深度氧化的动力学方程可表示为r=kpApB/(1+KBpB)2式中pA、pB分别表示乙烯及氧的分压。在473K等温下的实验数据如表1:表1 实验数据 　PA*103MPaPB*103MPar*104mol/g*min18.990 3.230 0.672 214.220 3.000 1.072 38.860 4.080 0.598 48.320 2.030 0.713 54.370 0.890 0.610 67.750 1.740 0.834 77.750 1.820 0.828 86.170 1.730 0.656 96.130 1.730 0.694 106.980 1.560 0.791 112.870 1.060 0.418 试求该温度下的反应速率常数k和吸附平衡常数KB。【1】本题计划用最小二乘原理解决，最小二乘法的原理暨通过使数据的残差平方和最小寻找数据的最佳函数匹配。但是因为此模型公式非线性，不能直接套用线性最小二乘法，需要对公式进行线性化处理：Kr=kPAPB(1+KBPB)2 =1/+KBPB/直线方程形式为y=b+a*x。即有：y=x=PBa=KB/b=利用excel计算结果如下表2：表2最小二乘法数据表　xi*10^3yixi*yi*10^3xi^2*10^51 3.230 0.657 2.123 1.0433 2 3.000 0.631 1.892 0.9000 3 4.080 0.777 3.172 1.6646 4 2.030 0.487 0.988 0.4121 5 0.890 0.253 0.225 0.0792 6 1.740 0.402 0.700 0.3028 7 1.820 0.413 0.751 0.3312 8 1.730 0.403 0.698 0.2993 9 1.730 0.391 0.676 0.2993 10 1.560 0.371 0.579 0.2434 11 1.060 0.270 0.286 0.1124 SUM22.870 5.055 12.090 5.6875 b={- [（）*（）]}/[-*]a=-b由表中数据有：k=1/b2KB=a*有：乙烯深度氧化动力学方程可以用公式:r=表示对回归方程进行分析：如下表3表3 原始数据与回归数据误差　PA*103MPaPB*103MPar*104mol/g*minr^*104mol/g*minδ*104mol/g*min18.993.230.6720.67761506900.005615068969 214.2231.0721.12558237610.053582376119 38.864.080.5980.5666578595-0.031342140461 48.322.030.7130.82959876410.116598764077 54.370.890.610.5844587582-0.025541241806 67.751.740.8340.83424588030.000245880301 77.751.820.8280.8165868004-0.011413199621 86.171.730.6560.66595405450.009954054523 96.131.730.6940.6616366863-0.032363313740 106.981.560.7910.7888186693-0.0021813306