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

姓名：方健 学号：652081701073 专业：化学工程 班级：化工1503 问题：利用MATLAB对与化学工程学科的实际问题进行解决，比如使用正交配置法求解有轴向扩散的固定床反应器中的催化剂的温度和浓度分布，柱形固定床反应器中的催化的温度和浓度方程为： -?T?r=1=BiwT1,z-Twz, -?c?r=1=0 Tr,0=T0, cr,0=c0 其中R(c,T)为催化反应的速率方程，其形式为： 分析：从上面的问题可以知道，以上是一个偏微分方程的数值求解问题，由于该方程含有初值条件，还有左右边界条件，所以应该使用pdepe函数进行求解，于是先学习了pdepe函数的具体用法如下： @pdefun(函数格式描述)：[c,f,s]=pdefun(x,t,u,du) @pdebc（边界条件描述）：[pa,qa,pb,qb]=pdebc(x,t,u,du) l:左边界，r:右边界 调用格式： sol=pdepe(m, pdepe, icfun, bcfun, xmesh, tspan) sol=pdepe(m, pdepe, icfun, bcfun, xmesh, tspan, options) 输入参数： 1. m表示定解问题的对称性，m=0, 1, 2分别代表平板、圆柱和球形。 2. pdefun是描述PDE问题的函数，输出变量为标准形式中的c