实验四 教堂顶部曲面面积.ppt

* * 　　本实验主要涉及微积分，通过实验将复习 曲面面积的计算、重积分和Taylor展开等知 识；另外将介绍重积分的数值计算法和取得函 数近似解析表达式的摄动方法。 一、实验目的 二、实际问题 　　某个阿拉伯国家有一座著名的伊斯兰教堂，它以中央大厅的金色巨大拱形圆顶名闻遐迩．因年久失修．国王下令将教堂顶部重新贴金箔装饰．据档案记载，大厅的顶都形状为半球面．其半径为 30 米．考虑到可能的损耗和其他技术因素，实际用量 将会比教堂顶部面积多 1.5 % ．据此，国王的财政大臣拨出了可制造5750m2 有规定厚度金箔的黄金。建筑商人哈桑略通数学．他计算了一下，觉得黄金会有盈余。于是他以较低的承包价得到了这项装饰工程．但在施工前的测量中，工程师发现教堂顶部实际上并非是一个精确的半球面而是半椭球面，其半立轴恰是30m，而长半轴和短半轴分别是30.6m 和29.6m ．这一来哈桑犯了愁．他担心黄金是否还有盈余？甚至可能短缺．最后的结果究竟如何呢？ 三、数学模型 　　取椭球面中心为坐标原点建立直角坐标系，则教堂顶部半椭球面的方程可写为 　　作为被积函数它的原函数容易吗？如不易，怎么办？ 引进变量代换 X=a＊r＊cosθ，y=b＊r＊sinθ 则有 这个积分相当复杂，若记 那么(4.3)中关于r的积分 　　这里μ＝1的情况要对表达式求极限。注意到μ的表达式(4.4)，若将式(4.5)代入(4.3)得到一个极为复杂的积分式。事实上，这是一个无法以初等函数形式来表达的积分，因此必须用近似方法来处理它。考虑到这一积分形式相当复杂，我们宁可直接对式(4.3)来进行处理。 四、数值积分方法 　　对于二重积分，可以如同一元函数定积分那样，将区域划分为小块，然后在每个小区域上对被积函数作近似简化求积，再把所得的值求和即可。 　　考虑矩形区域D:{0≤ s ≤c , 0≤ t ≤d}上的二重积分 　　将D划分为m×n个相等的小{si-1≤s≤si , tj-1≤t≤tj}，其中si和tj分别是s和t方向的分点： si=ik , tj=jh (i =0,1,2, … , m ; j =0,1,2, … ,n) 而k=c/m , h=d/n , 那么小矩形上的积分可写为 记 则 若对这两个单积分都用梯形法，就有 (4.7) 而 (4.8) 这样可求得在D上的积分Ι的近似值 当将分点增加一倍使得 而记 那么对 (4.10) 的两次积分都用Simpson 法,就得到 (4.11) 从而 (4.12) 于是有积分I 的近似值为 奇数 (4.13) 　　为了形象地表现由式(4.8)和式(4.12)给出的计算格式,常用的办法是如下图那样在分割的小矩形顶点上标出相应的函数值和式中的加权数. 1 1 (i,j) (i+1,j+1 ) 1(i-1,j-1) 　　1 １ (格式(4.８)) 16 (i,j) 1 　　　 4 　　1 １ 1 　　　　　　 1 (i-1,j-1) (格式(4.1２)) 　　对于矩形区域的重积分,还存在其他的计算格式; 另外,对于一些规则区域,也有相应的近似计算格式 ． 　　对于问题中的面积,由于积分表达式(4.５)　 在r=1处的奇性(被积函数分母为零),不能直接对此 式用上述数值方法.但只要作一个变换 （4.14） 就可以将(4.5)改写为 (4.15) 现在我们可以用数值积分的方法了.把 二重积分(4.15)的积分区域 分割为m2个小矩形,利用近似格式(4.10),我们得到 如表4.1的计算结果 表4.1 m S m S 2 5621.42 16 5679.83 4 5679.78 24 5679.82 6 5679.89 44 5679.81 10 5679.84 100 5679.81 从所得结果看,取面积值为5679.81(m2)已经 极为精确,若仅要求精确到0.1 (m2) ,那么计算量 可大大减少,此时 ,而 S ≈ 5679.81(m2) 加上技术与损耗等原因,教堂顶部实际使用金箔总面积为 显然,建筑商人哈桑在金箔上将入不敷出,从而招受损失. 建议同学们用近似格式(4.12)格式再算一次,以作 比较. 数值积分的方法给出了近似结果,以下我们再介绍 另一种近似方法. 五.摄动方法 简单地说,摄动方法就是对解析式中的小参 数进行展开,从而求得近似解析解的方法.应用与 积分计算,常常是采取将被积函数(或其部分)