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

计算机化工应用

课后习题答案

特别鸣谢 章奇明 材料支持

Work01

1. 已知某化工管道的真实长度为1000m，某次测量结果为1001m，其测量的绝

对误差，相对误差为多少？

答：绝对误差为1m，相对误差为0.1% 。

2. 已知e 的真实值为2.718281828459046···，取2.71828 作为e 的替代值时，具

有几位有效数字。取2.71827 作为e 的替代值时，有几位有效数字。

答：2.71828 具有 6 位有效数字。e − 2.71828 < 0.000001828 < 0.5×10−5 ，m-p=-5，

m=1 ，所以p=6 。2.71827具有5位有效数字，e − 2.71827 < 0.000011828 < 0.5×10−4 ，

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*106+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

Work02

Work03

Work04

Work05

Work06

Work07

Work08

Work09

Work10

Work11

1. 求解以下高阶常微分方程

1) 输出 y, 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. 求解下面的微分方程并作解函数图和相平面图

解：

function Excercise6_2

tspan=[0 30];

y0=[1 0.5];

ode45(@odefun,tspan,y0)

title('解函数图')

legend('x','y')

figure

opts=odeset('Outputfcn','odephas2');

sol=ode45(@odefun,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 表示甲苯。反应器中各