首页/文章/ 详情

数值积分|牛顿-柯特斯公式

7月前浏览4978

牛顿-柯特斯(Newton-Cotes )公式也叫插值型求积公式。已知    的值。以这n+1个点进行拉格朗日插值,得到n次多项式,再对该n次的多项式求积分。

将积分区间    等分,    

 

则n次拉格朗日插值多项式为:

 

其中

 

那么

 
 

 
 

由     可得

 

这就是牛顿-柯特斯公式。其中,    称为柯特斯系数。

由    式可知,柯特斯系数    与被积函数以及积分区间都无关,只要给出积分区间的等分数n,就可以算出柯特斯系数    。例如,当n=2时

 

对应的牛顿-柯特斯公式为:

 

此即为辛普森(Simpson)公式。

为了便于应用,将柯特斯系数列出,可以快速写出牛顿-柯特斯公式。

牛顿-柯特斯公式的缺点:对于次数较高的多项式而有很大误差(龙格现象),一般取低阶公式计算。

[算例]用牛顿--柯特斯公式计算积分

 

   

 

   

 

   

 

精确值为

 



★★★★往期相关★★★★

数值积分|二元函数的高斯积分

数值积分|高斯积分

数值积分|泰勒(Taylor)公式求积分

Python实现辛普森公式

建模| 利用辛普森公式计算湖泊面积

数值积分| 辛普森公式

来源:数值分析与有限元编程
python
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2024-04-02
最近编辑:7月前
太白金星
本科 慢慢来
获赞 5粉丝 10文章 322课程 0
点赞
收藏
作者推荐

有限元方法为何偏爱高斯积分

有限元中用高斯积分(Gauss)较多,有些参考书会着重介绍高斯积分,譬如王勖成老师的《有限单元法》中会有专门一节。有限元方法为何偏爱高斯积分?个人总结有以下两个原因:一) 积分复杂采用等参坐标构造的等参单元刚度矩阵几乎得不到显式表达式。高斯积分很方便得到单元刚度矩阵。二) 克服剪切锁定等问题在实际问题中如深梁、厚板、厚壳等考虑剪切行为的单元中会出现所谓的剪切自锁,也包括一些体积自锁、膜自锁等单元中的应用。这些出现所谓自锁的单元,用少于完全积分所要求的积分点个数时,可以发现其对自锁能起到缓解作用。减缩积分的概念也容易理解,即采用少于精确积分点个数的积分。比如平面四边形等参数单元求节点位移时采用4个积分点,求单元应力时采用一个积分点。在二维、三维单元中,插值函数会有不同方向多项式的组合,也包括被积函数最终形式可能不是多项式,但仍然采用高斯积分。减缩积分带来的零能模式或沙漏,数值上可以理解为矩阵的奇异性,这就涉及到减缩积分点的布置,譬如在某一个方向采用完全积分,另外一个方向采用减缩积分,或者所有方向统一采用减缩积分。当然还有另外一条技术路径避免剪切自锁,以假设应变为代表,强制对产生自锁的应变进行假设,然后可以仍然采用完全积分进行计算。引入的假设应变相当于添加了约束,可以用约束变分原理来实现。相对于减缩积分避免一些自锁问题,假设应变+完全积分有数值性能相当好单元,譬如在板壳中,MITC这类假设应变壳单元是性能极优异的单元之一。★★★ 往期相关 ★★★数值积分|高斯积分弹性厚板的剪切锁定(shear locking)来源:数值分析与有限元编程

未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习 福利任务 兑换礼品
下载APP
联系我们
帮助与反馈