教堂顶部曲面面积的计算方法.docVIP

实验六 教堂顶部曲面面积的计算方法 一、实验目的 本实验主要涉及微机分,通过实验将复习曲面面积的计算,重积分和Taylor展开等知识;另外将介绍重积分的数值计算方法和取得函数近似解析式的摄动方法. 二、实际问题 某个阿拉伯国家有一座著名的伊斯兰教堂,它以中央大厅的金色巨大拱形圆顶名震遐尔.因年久失修,国王下令将教堂顶部重新贴金箔装饰.据档案记载,大厅的顶部形状为半球形,其半径为30m.考虑到可能的损耗和其它技术原因,实际用量将会比教堂顶部面积多1.5%.据此,国王的财政大臣拨出了可制造5,750m SKIPIF 1 0 有规定厚度金箔的黄金.建筑商人哈商略通数学,他计算一下,觉得黄金会有盈余.于是,他以较低的承包价得到了这项装饰工程.但在施工前的测量中,工程师发现教堂顶部实际上并非是一个精确的半球面而是半椭球面,其半立轴恰好30m,而半长轴和半短轴分别是30.6m和29.6m.这一来哈桑犯了愁,他担心黄金是否会有盈余?甚至可能短缺.最后的结果究竟如何呢? 三、数学模型 取球面中心为坐标原点建立直角坐标系,则教堂顶部半椭球面方程可表示为: SKIPIF 1 0 (6.1) 其中R=30,a=30.6,b=29.6,而其表面积为: SKIPIF 1 0 (6.2) 这里积分区域 SKIPIF 1 0 为: SKIPIF 1 0 SKIPIF 1 0 通过简单的计算容易得到: SKIPIF 1 0 引进变量代换： SKIPIF 1 0 则有 SKIPIF 1 0 (6.3) 这个积分相当复杂,是一个无法用初等函数形式表示的积分. 四、数值积分方法 对于二重积分,可以如同一元函数定积分那样,将区域划分为小块,然后在每个区域上对被积函数作近似简化求积,再把所得的值求和即可. 考虑矩形区域 SKIPIF 1 0 上的二重积分 SKIPIF 1 0 (6.4) 将 SKIPIF 1 0 划分作 SKIPIF 1 0 个相等的小矩形{ SKIPIF 1 0 },其中 SKIPIF 1 0 和 SKIPIF 1 0 分别是 SKIPIF 1 0 和 SKIPIF 1 0 方向的分点: SKIPIF 1 0 而 SKIPIF 1 0 ,那么小矩形上的积分可分写为: SKIPIF 1 0 (6.5) 记 SKIPIF 1 0 则 SKIPIF 1 0 若对这两个单积分都用梯形法,就有 SKIPIF 1 0 而 SKIPIF 1 0 (6.6) 这样便可求得在 SKIPIF 1 0 上的积分 SKIPIF 1 0 的近似值 SKIPIF 1 0 (6.7) 当将分点增加一倍使得 SKIPIF 1 0 而记 SKIPIF 1 0 那么对 SKIPIF 1 0 （6.8） SKIPIF 1 0 的两次积分都用Simpson 法,就得到 SKIPIF 1 0 （6.9） 从而 SKIPIF 1 0 SKIPIF 1 0 SKIPIF 1 0 （6.10） 于是有积分I的近似值为 SKIPIF 1 0 , SKIPIF 1 0 奇数 (6.11) 为了形象地表现由式(6.6)和式(6.10)给出的计算格式,常用的办法是如图2.1那样在分割的小矩形顶点上标出相应的函数值和式中的加权数. 1 1 (i,j) 1 4 1 (i+1,j+1 ) 16