最近有人要我帮求函数的二重积分，于是我勉强为其难的答应了，过程有点艰苦但是还是得坚持。
对于函数求矩形积分我们可以采用插值型求积公式既是牛顿科特斯公式，既是将积分区间 [a,b] 划分为n等分，步长h=(b-a)/n,选取Xk=a+K*h构造出来的牛顿—科特斯公式，牛顿-科特斯公式的系数的n越高一般代数精度越高，但是在n>8的时候不具有稳定性，而辛普森公式是n=4的牛顿-科特斯公式，用辛普森公式可以简单又快速计算出积分结果但是误差大，所以为了提高精度采用复合求积法既是复合辛普森公式。
而所知道的复合辛普森求积公式是求一重积分F（X）的，而我们要求的是二重积分F（X，Y）。F(x,y)=∬▒f（x,y）dxdy,而复合辛普森公式是
![Sn=h/6*[∫▒〖f(a)+4∑_(k=0)^(n-1)▒〖f(x_(k+1/2) )+〗〗 2∑_(k=1)^(n-1)▒〖f(x_k )+f(b)]〗](https://i-blog.csdnimg.cn/blog_migrate/5712932de6931a8632615881b7b9af27.png)
这里求的是x的积分，对于求二重积分F（x,y）而言相当于对x求一次偏导，根据积分的可加性,复合辛普森公式可以变为：
![Sn=h/6*[∫▒〖f(a)+∫▒∑_(k=0)^(n-1)▒f(x_(k+1/2) ) +∫▒∑_(k=1)^(n-1)▒〖f(x_k )+∫▒f(b) 〗〗]](https://i-blog.csdnimg.cn/blog_migrate/bac18574e94f42d933836e9d05d54bc3.png)
对f(x,y)的一次积分有
因为我们要求f(x,y)的二重积分∬▒f（x,y）dxdy，我们可以将①式的每一项修改为
![Vn=h/6*[∬▒〖f(a)+∬▒〖∑_(k=0)^(n-1)▒f(x_(k+1/2) ) +∬▒∑_(k=1)^(n-1)▒〖f(x_k )+∬▒f(b) 〗〗〗]--②](https://i-blog.csdnimg.cn/blog_migrate/f387e3bb832d55269c70bc58c7fbf041.png)
有此结构我们可以对②式进行修改，对每一项求x的偏积分，进而就是对每一项求一次复合辛普森，因此②式可变为
最低0.47元/天 解锁文章
6512
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
