Python数值计算(33)——simpson 3/8积分公式

1. 背景知识

既然前的Simpson可以通过使用三个点构造二次曲线近似积分,那么,如果点数增加到了4个,然后不就可以构造三次多项式的曲线,实现对目标值的积分吗?

如果采用和上一节介绍的同样的方法,我们可以推导出:

其中

2. 算法实现

通过上述原理可以很容易使用Python实现:

python 复制代码
# Simpson's 3/8 Rule for Numerical Integration
import numpy as np
def simpson38(f, a, b):
    h=(b-a)/3
    x=np.linspace(a,b,4)
    y=np.vectorize(f)(x)
    return 3*h*(y[0]+3*y[1]+3*y[2]+y[3])/8.0

如果使用该方法,计算humps函数在区间[0,1]上的数值积分,其效果如下:

可见误差相比之前有了很大的改善,但还是有27%,还是偏大的。

3. 复合simpson 3/8积分公式

和前面类似,如果将一个大的积分区间,划分为更多个较小的积分区间,在每个积分区间上使用Simpson 3/8积分公式,则可以更加贴近原函数,和复合Simpson 1/3积分公式要求点的个数为偶数类似,这个也要求点的个数为3的倍数。

最终实现效果如下:

相关推荐
大佬,救命!!!6 小时前
整理python快速构建数据可视化前端的Dash库
python·信息可视化·学习笔记·dash·记录成长
孔丘闻言6 小时前
python调用mysql
android·python·mysql
Teletele-Lin6 小时前
Miniconda安装与VSCode搭建远程Python、Jupyter开发环境
vscode·python·jupyter·环境配置·远程开发
伊玛目的门徒6 小时前
告别 OpenAI SDK:如何使用 Python requests 库调用大模型 API(例如百度的ernie-4.5-turbo)
python·openai·requests·大模型调用·ernie-4.5-turbo
sinat_602035367 小时前
模块与包的导入
运维·服务器·开发语言·python
计算机学姐7 小时前
基于Python的旅游数据分析可视化系统【2026最新】
vue.js·后端·python·数据分析·django·flask·旅游
恋雨QAQ8 小时前
python函数和面向对象
开发语言·python
天雪浪子8 小时前
Python入门教程之逻辑运算符
开发语言·python
空灵之海8 小时前
Ubuntu系统安全合规配置
linux·ubuntu·系统安全·1024程序员节
张子夜 iiii8 小时前
实战项目-----在图片 hua.png 中,用红色画出花的外部轮廓,用绿色画出其简化轮廓(ε=周长×0.005),并在同一窗口显示
人工智能·pytorch·python·opencv·计算机视觉