二重积分的数值方法
《数值分析课程设计》
报 告
专业:
学号:
学生姓名:
指导教师:
一、题目
数值积分中二重积分探究。
二、理论
数值积分就是用数值方法近似计算定积分。其原理很简单,就是将积分核用插值多项式替代,用多项式的结果近似定积分的值。
一般常用的方法是,将积分区间等分为个子区间,即取步长
,
子区间端点为(k=0,1,…,n),在每个子区间上套用插值积分公式,再将个区间的结果累加起来。
比较常用的有梯形公式,是在每个子区间上用1阶多项式(即直线段)近似并积分的结果:
另外一种在实际应用中很受欢迎的方法是,在每个子区间上用2阶多项式(即抛物线)近似并积分,得到著名的辛普森(Simpson)公式:
其中。
三、方法、算法与程序设计
Ⅰ.辛普森公式求二重积分
考虑二重积分,它是曲面与平面区域R围成的体积,对于矩形区域,可将它写成累次积分
。
若用复合辛普森公式,可分别将,分成N,M等份,步长
,,
先对积分,应用复合辛普森公式,令,,则
从而得
。对每个积分再分别用复合辛普森公式即可求得积分值。MATLAB程序见附录1,MATLAB中自带自适应辛普森公式dblquad(),对于变量区域同样适用。
对于变量区域,