工程数学作业第十一次方健(1211)详解.docxVIP

姓名：方健 学号：652081701073 专业：化学工程 班级：化工1503 问题：利用MATLAB对与化学工程学科的实际问题进行解决，比如间歇反应器中的一级连串反应，在间歇反应器中进行液相反应制备产物E，反应如图1所示。反应可在180~260℃的温度范围内进行，各反应均为一级动力学关系：r＝－kC，式中已知参数：k01＝5.78052×1010，k02＝3.92317×1012，k03＝1.64254×104，k04＝6.264×108，Ea1=124670，Ea2=150386，Ea3=77954，Ea4=111528。初始浓度：CA=1kmol/m3，其余物质浓度为0。已知是产物B收率最大的最优反应温度为224.6℃。 试计算： 1)在最优反应温度下各组分浓度随时间的动态变化； 2)最优反应时间； 3)输出产物D对反应物浓度A的关系图。 A B C D E k1 k2 k3 k4 图1 间歇反应序列示意图 分析：根据以上问题，由于反应是一个一级反应，而且各个反应的速率常数可以根据方程计算得出，根据反应动力学可以得到以下数学模型： 根据上面的模型，编写微分方程组的m文件D_Equation.m如下： function dCdt =D_Equations(t,C,k0,Ea,R,T) T = 224.6 + 273.15; R = 8.31434; k0 = [5.