-数值微分与积分.pptVIP

(c) 计算塔径（m） 例题 5-3 (本例选自《化工原理》) 计算方法 此处自变量为Y，故应将X用Y表示 (d) 计算传质单元数 NOG 数值积分法 F(Y) 塔的任一截面与塔顶间物料平衡 全塔物料平衡式 联立上两式求出 上式代入 (e) 计算传质单元高度和填料高度HOG （m） （m） 未知量X用Y表示 X （2）计算程序及结果 参考程序见教材p.59。 计算结果 出塔吸收液最大浓度 Xmax=0.0313 吸收剂最小用量 Lmin=0.0361 kmol/s 出塔吸收液实际浓度 x1=0.0188 吸收剂实际用量 L=0.0602 kmol/s 填料塔径 D=0.9404 填料层高度 Z=11.5691 m function y=f(x,Y1,Y2,M,x1,x2) y=1./(x-((x-Y2).*(x1-x2)./(Y1-Y2)+x2).*M); %被积函数F的表达式 必须是向量形式 clear KV=1.1111*10^(-4); Vt=0.5556;M=1.68;a0=190;u=0.8; Y1=5/95; %入塔气体丙酮比摩尔分率 Y2=(1-0.98)*Y1; %出塔气丙酮比摩尔分率 x2=0; %吸收剂入口浓度 V=273*(1-0.05)*Vt/293/22.4; %惰性气体流量 GA=V*(Y1-Y2); %吸收速率 X1m=Y1/M; %吸收液出口最大浓度 Lmin=GA/(X1m-x2); %吸收剂最小用量 x1=X1m*0.6; %吸收液出口实际浓度 L=GA/(x1-x2); %吸收剂实际用量 D=sqrt(4*Vt/(pi*u)); %塔径计算 Nog=quad(@f,Y2,Y1,[],[],Y1,Y2,M,x1,x2); %调用辛普生积分函数计算传质单元数 Hog=V/(KV*a0*(pi*D^2/4)); %计算传质单元高度 Z=Nog*Hog; %计算填料层高度 string=kmol/s; fprintf(\tXmax=%.4f\tLmin=%.4f\t%s\n\tx1=%.4f\t\tL=%.4f\t\t%s,X1m,Lmin,string,x1,L,string) fprintf(\n\tD=%.4f\t\tZ=%.4f\t\t%s,D,Z,m) 被积函数为列表函数时的积分问题如何计算？ x 1.0 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 f (x) 1.543 1.668 1.811 1.971 2.151 2.352 2.577 2.828 3.107 积分： 步长分别取(a)h=0.1, (b)h=0.2，(c)h=0.4 。 (a) h=0.1，用复化梯形公式 辛普生公式 积分： x 1.0 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 f (x) 1.543 1.668 1.811 1.971 2.151 2.352 2.577 2.828 3.107 (b) h=0.2，（5个点） 辛普生公式 (c) h=0.4, （3个点） 梯形公式 梯形公式 辛普生公式 计算结果比较 步长 h=0.1 h=0.2 h=0.4 梯形公式 1.7683 1.7728 1.7904 辛普生公式 1.7668 1.76693 1.7672 柯特斯公式 1.76679 步长相同时，辛普生公式 与用柯特斯公式的计算精度高于梯形公式 例题 5-4 用脉冲法测量某反应器的停留时间分布，由实验数据计算寿命分布密度与时间t的关系如下表 t 10 15 20 25 30 35 40 E(t)×103 0 1.92 9.0 20.95 31.10 32.20 30.60 求分布密度积分 解：根据所给数据，小区间分为6等分，步长h=5，节点个数为7，用辛普生复化求积公式计算。 本章要求 清楚数值积分的原理 掌握数值积分的梯形公式和辛普生公式、复化梯形公式和复化辛普生公式 了解龙贝格积分的主要思想 能用数值积分求解实际问题 分别用定步长和变步长辛普森积分计算 * 填料吸收塔 基于气相的总传质单元数（NTU）为 5.2 数值积分基础 定积分计算 并不是所有的积分都能求原函数和/或积分值，原因有： 被积函数形式复杂，找不到积分公式求原函数 原函数形式复杂，不便计算（如无穷级数） 被积函数以列表函数形式给出，无法用积分公式 问题：能否不通过求原函数来求积分值? 方法