数学家生日蛋糕.pptVIP

实验四 ;一、引例：数学家的生日蛋糕;解 设高为H，半径 r, 比重为k 若蛋糕是单层圆盘的，则蛋糕的重量和表面积分别为： W=k?Hr2 S = 2?Hr+?r2 若蛋糕是双层的，每层高H/2 ，下层半径r1, 上层半径r2，则 W=k?H(r12+ r22)/2 S = ?H(r1+r2)+?r22 如果蛋糕是n层的，每层高H/n，半径分别r1, …, rn, 则;若蛋糕边缘是曲线r = r(h), 0hH，各层半径近似为ri = r((i-1/2)H/n), i=1, …, n, 那么 当n??，;二、数学理论复习：积分 函数f(x)在区间[a,b]上的积分定义为;三、数值积分：梯形法和重积分 1、梯形法 ;通常将[a,b]区间n等分，h=(b-a)/n, xi=a+ih, ;我们利用梯形法，先将[a,b]区间m等分，hx=(b-a)/m, xi=a+ihx, i=0,1,…,m;M文件dblquad2.m给出二重积分数值计算法。 I=dblquad2(‘f_name’,a,b,‘c_lo’, ‘d_hi’, m,n) 其中‘f_name’为被积函数f(x,y)字符串,其中x为标量, y为向量 ,‘c_lo’和‘d_hi’是y的下上限 ; a , b为x的下上限；m, n 为x和y方向的等分数(缺省值100)。;四、使用MATLAB; 1、梯形积分法 trapz 是最基本的数值积分方法, 精度低。 z=trapz(x,y), 返回积分的近似值，其中x 表示积分区间的离散化向量; y是与x同维数的向量,表示被积函数 。; 2、变步长数值积分 z=quad8(‘fun’,a,b,tol) 返回积分的近似值， 其中 fun表示被积函数的M函数名, a、b 表示积分下上限, tol为精度, 缺省值为1e-30。;3、重积分; 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的定积分。;例 3(引例) 现在来求大蛋糕重量和表面积：解 ? syms h; r=2-(exp(2*h)+exp(-2*h))/5; ?vpa(int(pi*r^2,h,0,1),5) , rn=subs(r,h,1); ? vpa(int(2*pi*r,h,0,1)+pi*rn^2,5);考虑在时间dt内水面变化dz，漏水的体积为 uAdt= - ?x2dz 其中x为高度z水面的半径, A=?b2 由于R2=z2+x2 得dt =