python 二重积分 辛普森法_二重积分的数值方法.doc
二重积分的数值方法
《数值分析课程设计》
报 告
专业:
学号:
学生姓名:
指导教师:
一、题目
数值积分中二重积分探究。
二、理论
数值积分就是用数值方法近似计算定积分。其原理很简单,就是将积分核用插值多项式替代,用多项式的结果近似定积分的值。
一般常用的方法是,将积分区间等分为个子区间,即取步长
,
子区间端点为(k=0,1,…,n),在每个子区间上套用插值积分公式,再将个区间的结果累加起来。
比较常用的有梯形公式,是在每个子区间上用1阶多项式(即直线段)近似并积分的结果:
另外一种在实际应用中很受欢迎的方法是,在每个子区间上用2阶多项式(即抛物线)近似并积分,得到著名的辛普森(Simpson)公式:
其中。
三、方法、算法与程序设计
Ⅰ.辛普森公式求二重积分
考虑二重积分,它是曲面与平面区域R围成的体积,对于矩形区域,可将它写成累次积分
。
若用复合辛普森公式,可分别将,分成N,M等份,步长
,,
先对积分,应用复合辛普森公式,令,,则
从而得
。对每个积分再分别用复合辛普森公式即可求得积分值。MATLAB程序见附录1,MATLAB中自带自适应辛普森公式dblquad(),对于变量区域同样适用。
对于变量区域,写成累次积分的形式:
进行数值计算的表达式为:
上面的表达式中、表示权重,取决于一维积分方法。我们常用复合辛普森公式,先对内积分进行计算,在计算外积分,与矩形区域情况基本一致。
Ⅱ高斯求积公式求二重积分
在高斯求积公式中,若取权函数,区间为,则得公式
。
勒让德多项式是区间上的正交多项式,因此,勒让德多多项式的零点就是求积公式的高斯点。
若取的零点做节点构造求积公式
;
若取的零点构造求积公式
;
当时,求积公式为
同样先用高斯求积公式求内积分,再求外积分,可得二重积分值。
四、算例、应用实例
算例:
计算二重积分。
(1)若区域,试分别用复合辛普森公式(取n=4)及高斯求积公式(取n=4)求积分。
(2)若区域用复合辛普森公式(取n=4)求积分。
解:
(1)
=
对各个积分应用复合辛普森公式。
也可应用MATLAB中的函数进行计算,程序见附录2。
先将区域变换为区域,其中
,等价于,有
。
对于取时的高斯求积公式节点及系数,即
,,,,
,,,
用的高斯积分公式计算积分I,
(2)
,
等分为4等份,对应值为的值,用
节点应用辛普森公式对内积分求积,再用复合辛普森公式对外积分求积,
也可用MATLAB中的函数实现,结果和程序如下(附录3)。
五、参考文献
【1】 数值分析 李庆扬,王朝能,易大义 清华大学出版社
【2】 数值分析课程设计 陈越,童若锋 浙江大学出版社
【3】 MATLAB教程 张志涌 北京航空航天大学出版社
六、附录
附录1:
function q=DblSimpson(f,a,A,b,B,m,n)if(m==1 && n==1)??%辛普森公式q=((B-b)*(A-a)/9)*(subs(sym(f),findsym(sym(f)),{a,b})+...subs(sym(f),findsym(sym(f)),{a,B})+...subs(sym(f),findsym(sym(f)),{A,b})+...subs(sym(f),findsym(sym(f)),{A,B})+...4*subs(sym(f),findsym(sym(f)),{(A-a)/2,b})+...4*subs(sym(f),findsym(sym(f)),{(A-a)/2,B})+...4*subs(sym(f),findsym(sym(f)),{a,(B-b)/2})+... 4*subs(sym(f),findsym(sym(f)),{A,(B-b)/2})+...16*subs(sym(f),findsym(sym(f)),{(A-a)/2,(B-b)/2}));else?? %复合辛普森公式q=0;for i=0:n-1for j=0:m-1x=a+2*i*(A-a)/2/n;y=b+2*j*(B-b)/2/m;?????????
python 二重积分 辛普森法_二重积分的数值方法.doc相关推荐
- python驼峰命名法_蛇形命名,驼峰命名,你的选择?
关于今天要讲的,小编想想还有点小激动.那就是关于变量的命名,这个一直都是容易引起程序员论战的话题.如何命名才更具可读性.易写性以及明义性? 程序命名常使用的是英文单词,当他们被作为变量名时,有全小写. ...
- python 斗地主最优解法_一步步解析Python斗牛游戏的概率
过年回家,都会约上亲朋好友聚聚会,会上经常会打麻将,斗地主,斗牛.在这些游戏中,斗牛是最受欢迎的,因为可以很多人一起玩,而且没有技术含量,都是看运气(专业术语是概率). 斗牛的玩法是: 1.把牌中的J ...
- python交叉验证法_详解python实现交叉验证法与留出法
在机器学习中,我们经常在训练集上训练模型,在测试集上测试模型.最终的目标是希望我们的模型在测试集上有最好的表现. 但是,我们往往只有一个包含m个观测的数据集D,我们既要用它进行训练,又要对它进行测试. ...
- python指数描述法_季节指数法的计算步骤
展开全部 1.收集历年(通常至少有三年)各月或各季的统计资料(观察值).62616964757a686964616fe59b9ee7ad9431333365643662 2.求出各年同月或同季观察值的 ...
- python二重积分_用python求一重积分和二重积分的例子
首先是对一元函数求积分,使用Scipy下的integrate函数: from scipy import integrate def g(x): return (1-x**2)**0.5 #用integ ...
- 八皇后问题python回溯_解决Python基于回溯法子集树模板实现8皇后问题
这篇文章主要介绍了Python基于回溯法子集树模板实现8皇后问题,简单说明了8皇后问题的原理并结合实例形式分析了Python回溯法子集树模板解决8皇后问题的具体实现技巧,需要的朋友可以参考下 本文实例 ...
- python回溯方法的模板_实例讲解Python基于回溯法子集树模板实现图的遍历功能
这篇文章主要介绍了Python基于回溯法子集树模板实现图的遍历功能,结合实例形式分析了Python使用回溯法子集树模板针对图形遍历问题的相关操作技巧与注意事项,需要的朋友可以参考下 本文实例讲述了Py ...
- 数据结构与算法(Python)– 回溯法(Backtracking algorithm)
数据结构与算法(Python)– 回溯法(Backtracking algorithm) 1.回溯法 回溯算法实际上一个类似枚举的搜索尝试过程,主要是在搜索尝试过程中寻找问题的解,当发现已不满足求解条 ...
- 用python玩转数据测试答案_MOOC_用Python玩转数据_测试答案
利息.据测股息.红利所得以每次收入额为应纳税所得额.A:错B:对 利息保障倍数中,试答利息支出:试答A:不包括财务费用中的利息支出B:包括计入固定资产成本的资本化利息C:包括财务费用中的利息支出D:不 ...
最新文章
- libcurl+ncurses 分段range批量下载和进度条显示源码实例
- Ubuntu 14.04 64位上安装wps office软件
- remote: 此仓库的限制大小为: 2048 MB, 您的使用已经超出限额
- java设计模式 观察者模式_理解java设计模式之观察者模式
- form-validation-engine中的正则表达式
- 2021牛客暑期多校训练营1 I-Increasing Subsequence(期望dp+优化)
- oracle外表日期,三步教会您掌握oracle外表(external table)
- LeetCode 116. 填充每个节点的下一个右侧节点指针(递归循环)
- 使用计算机教学的意义,信息技术对教学的意义
- 网络配置管理大大减轻网络管理员的负担
- 怎么解苹果7手机ID锁
- 安卓设备如何ROOT?玩转ROOT,让你的安卓手机更强更好用
- 【windows屏幕扩展】把你多余屏幕利用起来,spacedesk屏幕扩展超低延迟解决方案
- 余涛-Kinect技术在企业级的应用及展望
- fatal error C1001: INTERNAL COMPILER ERROR(compiler file 'msc1.cpp', line 1786)解决方法
- color.cpp:7456: error: (-215) scn == 3 || scn == 4 的解决办法
- 英语语法——句子分类
- JavaScript ES6 特性
- 【matlab】GMSK高斯最小频移键控
- wish平台怎么样?wish跨境电商好做吗?
热门文章
- 自考07172信息安全历年真题 四套
- 北大青鸟:手机软件开发指引
- CSRNet: Dilated Convolutional Neural Networks for Understanding the Highly Congested Scenes
- pacemaker之三节点drbd(单primary)
- IDEA启动报错问题解决
- pat考试1018 锤子剪刀布
- 游戏软件测试点,游戏签到--软件测试点练习
- 小白的jquery学习之路之天猫04 效果 fadeTo animate方法 天猫图片滑动效果
- HCIP-DATACOM-带解析-1-50题(821)
- matlab序列负轴,MATLAB中雷达信号回波图为什么会有负半轴