计算机应用基础-4-积分方程和应用.pptVIP

要使 则 M=5,取m=5, h=0.2 程序ME_3_2.m, S=3.8875, err=2.7286 4.3 数值积分问题 二 由给定数据进行梯形求积 4.3 数值积分问题 4.3 数值积分问题 利用trapz( )函数进行积分求解, 格式如下: S=trapz(x, y) 例 4-7 用定步长方法求解积分 画图： x=[0:0.01:3*pi/2, 3*pi/2]; y=cos(15*x); plot(x,y) 计算程序：ME_4_9 求理论值： 不同步距： syms x, A=int(cos(15*x), 0, 3*pi/2) h0=[0.1, 0.01, 0.001, 0.0001, 0.00001, 0.000001]; v=[ ]; For h=h0 x=[0:h:3*pi/2, 3*pi/2]; y=cos(15*x); I=trapz(x,y); v=[v; h,I, 1/15-I]; end 4.3 数值积分问题 syms x, A=int(cos(15*x), 0, 3*pi/2) 二 双重积分问题的数值解 4.3 数值积分问题 求解格式: (1) y=dblquad(Fun, xm, xM, ym, yM) % 矩形区域的双重积分 (2) y=dblquad(Fun, xm, xM, ym, yM, ?) % 限定精度的双重积分 【例4-8】求解 4.3 数值积分问题 f=inline (‘exp(-x.^2/2).*sin(x.^2+y)’, ‘x’, ‘y’); y=dblquad(f,-2,2,-1,1) 4.4 微积分的化工应用示例 一 微分方程应用示例—化学反应动力学速率的测定 动力学参数（级数和速率常数）常采用初始速率方法，对于每一个初始浓度，有初始速率： 例【4-7】 反应物A在一等温间歇反应器中进行如下反应：A B，实测A的浓度CA 随时间变化如下表： 时间 t/s 0 20 40 60 120 180 300 CA mol/L 10 8 6 5 3 2 1 反应动力学模型： 计算程序“ME_4_12 .m” 结果：k = 0.0053， n =1.3854 4.4 微积分的化工应用示例 二 积分方程应用示例—精馏塔理论板的计算 例【4-9】 氯仿-苯双组分精馏系统的起液平衡数据如下表所示，已知xf=0.4, xd=0.9, xw=0.15,精馏段与提留段回流比R=5，R‘=4，计算理论塔板数。 x 0.178 0.275 0.372 0.456 0.650 0.844 y 0.243 0.382 0.518 0.616 0.795 0.931 4.4 微积分的化工应用示例 计算程序”ME_4_11.m“ 精馏段理论板数：5 （块） 提馏段理论板数：6 （块） 共需理论板数为：11 （块） 4.3 微积分的化工应用示例 第四章 插值/拟合函数与积分求解 4.1 插值函数 （1）插值基础 （2）拉格朗日插值 （3）Matlab插值函数 4.2 拟合函数 （1）曲线拟合 （2） 最小二乘拟合 4.3 积分方程数值求解 一、插值基础 设函数y=f(x)在区间[a,b]上有定义，且已知在点 a ≤x0x1…xn ≤b上的值y0,y1,…,yn,若存在一简单函数P(x),使 P（xi）=yi (i=0,1,2,…,n) （4-1）成立，就称P(x)为f(x)的插值函数，点x0，x1，…，xn称为插值节点，包含插值节点的区间[a,b]称为插值区间，求插值函数P(x)的方法称为插值法。 若P(x)为插值多项式，即 P(x)=a0+a1x+…+anxn （4-2） 其中ai为实数，就称P(x)为插值多项式，相应的插值法称为多项式插值。若P(x)为分段多项式，就称分段插值。 4.1 插值函数  从几何图形上看，插值法就是求曲线y=P(x),使其通过给定的n+1个点（xi,yi),i=0,1,…,n,并用它近似已知曲线y=f(x)。 插值多项式存在的唯一性 设P(x)是形如（4-2）的插值多项式，则满足条件（4-1）的插值多项式（4-2）存在唯一。 4.1 插值函数 一 线性插值与抛物插值 通过平面的两点可以确定一条直线经过两点，这就是拉格朗日线