化工计算机技术课后答案,2014计算机化工应用课后作业答案(全)-隋志军

《2014计算机化工应用课后作业答案(全)-隋志军》由会员分享，可在线阅读，更多相关《2014计算机化工应用课后作业答案(全)-隋志军(38页珍藏版)》请在技术文库上搜索。

1、dC=zeros(5,1)； dC(1)=-k1*C(1)； dC(2)= k1*C(1)- k2*C(2)； dC(3)= k2*C(2) - k2*C(3)； dC(4)= k3*C(3)- k4*C(4)； dC(5)= k4*C(4)； tx=313:10:363； ma=653 549.2 469.7 406.0 355.0 314.8； mb=393.0 365.4 339.6 315.3 292.4 270.9； m1=pchip(tx,ma,T)； m2=pchip(tx,mb,T)； Miu=(m1(1/3)*xA+m2(1/3)*xB)3； disp(The calcula。

2、ted viscosity is ,num2str(Miu) 七、 ( 20 分) function Test22014C07 T=0:5:60； C0=1 0 0 0 0； T,C=ode45(MemEq,T,C0)； plot(T,C(:,1),-o,T,C(:,2),-+,T,C(:,3),-d,T,C(:,4),-*,T,C(:,5), -) xlabel(Reaction Time(s) ylabel(Concentration (mol/L) legend(CA,CB,CC,CD,CE) function dC=MemEq(t,C) k1=0.04； k2=0.05； k3=0.1。

3、； k4=0.08； tao=17.5 40.1 69.2 103.7 144.2 186.3； y=tao./q； x=ones(6,1),q； beta=regress(y,x)； K=1/beta(2) qe=beta(1)*K/2 plot(q,y,bo) hold on tcal=0.01:0.001:0.06； ycal=beta(1)+beta(2)*tcal； plot(tcal,ycal,k-) xlabel(tao s) ylabel(tao/q) legend(Experiment,Simulation) 六、 ( 15 分) function Miu=CalM(T) i。

4、f T363 warning(The input temperature is not in the range, the results may not be right) return end xA=0.45；xB=0.55； 模拟练习答案 四、 ( 10分) function Test20142C04 P=(T) exp(23.8296-3480.3./T-5.0811*log10(T)； Tcal=100:500； Pcal=P(Tcal)； plot(Tcal,Pcal,k-) xlabel(Temperature) ylabel(Pressure) TP=(T) 101.325-e。

5、xp(23.8296-3480.3./T-5.0811*log10(T)； tp=fzero(TP,700 900)； fprintf(The boiling T is %.2fn,tp) 五、 function Test20142C05 q=0.01:0.01:0.06； 模拟题 t,y=ode45(Hydro,0 0.5,CH0,CM0,CX0)； plot(t,y(:,1),b+-,t,y(:,2),rp:,t,y(:,3),k-.d) legend(CH,CM,CX) xlabel(Residence Time s) ylabel(Concentration mol/m3) topt=。

6、t(y(:,3)=max(y(:,3)； fprintf(The optimum time for Xylene production is %.4fn,topt) function dC=Hydro(t,C) k1=55.20；k2=30.20； CH=C(1)；CM=C(2)；CX=C(3)； dCH=-k1*CH0.5*CM-k2*CX*CH0.5； dCM=-k1*CH0.5*CM； dCX=k1*CH0.5*CM-k2*CX*CH0.5； dC=dCH；dCM；dCX； 反应器进料中含 66.7的三甲基苯和和 33.3的氢气，总摩尔流率为 0.036。反 应速率常数 k1 55.20。

7、， k2 30.20，停留时间 00.5。试编写一个 MATLAB 函 数实现以下任务： 1) 求出氢气，三甲基苯和二甲苯浓度随反应时间的变化曲线，将计算结果输出 如下图所示； 2) 计算二甲苯产量最高时的停留时间，采用 fprintf 函数结果输出在屏幕上。 function Work11_3B CH0=0.036*0.667；CM0=0.036*0.333；CX0=0； y0=1 0.5； ode45(odefun,tspan,y0) title(解函数图 ) legend(x,y) figure opts=odeset(Outputfcn,odephas2)； sol=ode45(ode。

8、fun,tspan,y0,opts)； title(相平面图 )； t=1.5 2 2.5； function f=odefun(t,y) dy1=-y(1).3-y(2)； dy2=y(1)-y(2).3； f=dy1；dy2； 3管式反应器计算 在管式反应器中进行 1,3,5-三甲基苯加氢脱烷基生成间二甲苯反应。反应器内发 生的反应如下： 其中 M 表示三甲基苯， H 表示氢气， X 表示二甲苯， Me 表示甲苯。反应器中各 物质浓度( CH，氢气； CM，三甲基苯； CX，二甲苯)随停留时间 的变化可由 以下常微分方程组表示： Work11 1. 求解以下高阶常微分方程 1) 输出 y,。

9、 y和 y”与 x 的关系图； 2) 输出 x=1.2 时 y, y和 y”的值。 解： function CACE6_13 x=1:0.05:1.5； opt=odeset(outputfcn,odeplot) x,y=ode45(odefun,x,1 10 30,opt)； yout=y(find(x=1.2),:) function dy=odefun(x,y) dy(1)=y(2)； dy(2)=y(3)； dy(3)=(2*y(3)+3*y(2)+3*exp(2*x)./(x.3)； dy=dy(1)；dy(2)；dy(3)； 2. 求解下面的微分方程并作解函数图和相平面图 解： f。

10、unction Excercise6_2 tspan=0 30； Work10 Work09 Work08 Work07 Work06 Work05 Work04 Work03 Work02 Work01 1. 已知某化工管道的真实长度为 1000m，某次测量结果为 1001m，其测量的绝 对误差，相对误差为多少？ 答：绝对误差为 1m，相对误差为 0.1%。 2. 已知 e 的真实值为 2.718281828459046，取 2.71828 作为 e 的替代值时，具 有几位有效数字。取 2.71827 作为 e 的替代值时，有几位有效数字。 答： 2.71828 具有 6 位有效数字。 e 。

11、2.71828 0.000001828 0.5105， m-p=-5， m=1，所以 p=6。 2.71827具有 5位有效数字， e 2.71827 0.000011828 0.5104， m-p=-4， m=1，所以 p=5。 3. 假定使用一台十进制计算机，它表示的浮点数具有 7 位尾数和 1 位指数，超 过计算机存储位数的数字均被舍去，试分别计算以下表达式的值，写出运算过程。 1) 123456.7+101.3254 2) 3563.212/12.53221 答： 1) 123456.7+101.3254 0.1234567*106+0.1013254*103 0.1234567*10。

12、6+0.0001013*106 0.1235580*106 2) 3563.212/12.53221 0.3563212*104/0.1253221*102 0.2843243*103 4. 以下浮点数运算采用 IEEE 双精度格式， eps 表示从 1 到一个较大浮点数的距 离；试写出其计算结果并简要说明计算过程。 1) 2+eps-2 2) 2-2+eps 3) 2.0e308-1.1e308 4)1.5*(1.8e308+0.5e308) 5) 0/(1+1e-16)-1) 6) 1.8e308-1.8e308 7) (1 + 3e-16) 1 答： 1) 0； eps 被舍掉 2) eps， 3) Inf 无穷减有限实数，因此结果为无穷 4) Inf 有限实数乘以无穷也为无穷 5) NaN 0/0 结果为非数 6) NaN 无穷减无穷结果为非数 7) eps 或 2.22e-16 距离 1 最近的实数是 1+eps，下个不小于 1+2eps，因此 1+3e-16 实际 被认为是 1+eps 计算机化工 应用 课后习题答案 特别 鸣谢 章奇明 材料支持 。