1、第1页end一、引例：数学家生日蛋糕一、引例：数学家生日蛋糕某数学家学生要送一个特大蛋糕来庆贺他90岁生日。为了纪念他提出口腔医学悬链线模型，学生们要求蛋糕店老板将蛋糕边缘半径作成以下悬链线函数 r=2-(exp(2h)+exp(-2h)/5,0h1(单位：米)。蛋糕成本取决于蛋糕重量和表面积（底面除外），问怎样计算重量和表面积？第2页end解 设高为H，半径 r,比重为k若蛋糕是单层圆盘，则蛋糕重量和表面积分别为：W=kHr2 S=2Hr+r2若蛋糕是双层，每层高H/2，下层半径r1,上层半径r2，则 W=kH(r12+r22)/2 S=H(r1+r2)+r22假如蛋糕是n层，每层高H/n，
2、半径分别r1,rn,则rHr1r2第3页end若蛋糕边缘是曲线r=r(h),0hH，各层半径近似为ri=r(i-1/2)H/n),i=1,n,那么当n，第4页end二、数学理论复习：积分二、数学理论复习：积分函数f(x)在区间a,b上积分定义为其中 a=x0 x1xn=b,xi=xi-xi-1,i(xi-1,xi),i=1,2,n若在a,b上,F(x)=f(x),则二重积分定义为第5页end三、数值积分：梯形法和重积分三、数值积分：梯形法和重积分 1、梯形法 在xi-1,xi上f(x)近似为一直线，用弦线代替，则 设f(x)在a,b上大于0，a=x0 x1xn=b,则a xi-1bxi第6页e
3、nd通常将a,b区间n等分，h=(b-a)/n,xi=a+ih,称为梯形公式 2、重积分 重积分数值计算可经过单积分组累计算第7页end我们利用梯形法，先将a,b区间m等分，hx=(b-a)/m,xi=a+ihx,i=0,1,m再将c(xi),d(xi)区间n等分,hy(i)=(d(xi)-c(xi)/n,yij=c(xi)+jhy(i),则G(xi)Aabd(x)c(x)其中第8页endM文件dblquad2.m给出二重积分数值计算法。I=dblquad2(f_name,a,b,c_lo,d_hi,m,n)其中f_name为被积函数f(x,y)字符串,其中x为标量,y为向量,c_lo和d_h
4、i是y下上限;a,b为x下上限；m,n 为x和y方向等分数(缺省值100)。第9页end四、使用四、使用MATLAB解 将被积函数及y上下限函数写为M函数文件，再调用dblquad2即可求解。例1trapz 梯形法积分 quad 变步长数值积分int 符号积分 quad8 高精度数值积分dblquad 矩形区域二重积分第10页end 1、梯形积分法 trapz 是最基本数值积分方法,精度低。z=trapz(x,y),返回积分近似值，其中x 表示积分区间离散化向量;y是与x同维数向量,表示被积函数。例2解 clear;x=-1:0.1:1;y=exp(-x.2);trapz(x,y)第11页en
5、d 2、变步长数值积分 z=quad8(fun,a,b,tol)返回积分近似值，其中 fun表示被积函数M函数名,a、b 表示积分下上限,tol为精度,缺省值为1e-30。注1：quad使用方法与quad8相同，但quad8精 度较高，且对假收敛和假奇异积分 含有一定适应性，而quad较差。注2：trapz,quad,quad8都不能用于求广义 积分。第12页end3、重积分 矩形区域二重积分.z=dblquad(fun,a,b,c,d)其中fun 表示被积函数 f M函数名,a,b表示变量 x 下上限，c,d表示变量 y 下上限。dblquad只能求矩形区域二重积分，不如上述M函数dblqu
6、ad2适用面广。第13页end 4、符号积分 int(s)符号表示式s 不定积分；int(s,v)符号表示式s 关于变量v不定 积分；int(s,a,b)符号表示式s 定积分,a,b分别 为下上限；int(s,v,a,b)符号表示式s 关于变量v 从 a 到 b定积分。第14页end例 3(引例)现在来求大蛋糕重量和表面积：解 syms h;r=2-(exp(2*h)+exp(-2*h)/5;vpa(int(pi*r2,h,0,1),5),rn=subs(r,h,1);vpa(int(2*pi*r,h,0,1)+pi*rn2,5)例 4 二分之一径为5m球形水罐充满了水，底部有二分之一径为b=0.1m小孔漏水，问多少时间以后，水面下降至离底部0.5m?解 水从孔漏出速度由以下能量方程决定 g(z+R)=u2/2 ，u是速度,z表示从球心 测量水面高度,g为重力加速度。第15页end考虑在时间dt内水面改变dz，漏水体积为 uAdt=-x2dz其中x为高度z水面半径,A=b2因为R2=z2+x2得dt=在顶部z=R水降到0.5m时，z=0.5-R，从而 t=5m0.5m0第16页
浙ICP备2021020529号-1 | 浙B2-20240490 
