数值积分——梯形公式和Simpson公式

最新推荐文章于 2023-12-11 22:14:51 发布

Amiyai

最新推荐文章于 2023-12-11 22:14:51 发布

阅读量3.1w

点赞数 19

分类专栏：数值分析文章标签： python 算法

本文链接：https://blog.csdn.net/weixin_45812669/article/details/116719598

版权

数值分析专栏收录该内容

3 篇文章

订阅专栏

本文介绍了数值积分中的梯形公式和Simpson公式，并给出了具体的Python实现。包括基本的梯形公式、Simpson公式及其复化版本，通过实例展示了如何使用这些方法进行数值积分计算。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

数值积分——梯形公式和Simpson公式

梯形公式：

$\displaystyle\int^b_a f(x)dx \approx\frac{b-a}{2}[f(a)+f(b)]$

def Trapezoid(a,b):
    y = (a - b) * ((1/2) * f(a) + (1/2) * f(b))
    return y

Simpson 公式：

$\displaystyle\int^b_a f(x)dx \approx(b-a)[\frac{1}{6}f(a)+\frac{4}{6}f(\frac{a+b}{2}) +\frac{1}{6}f(b)]$

def Simpson(a,b):
    y = (a-b)*((1/6)*f(a) + (4/6)*f((a+b)/2) + (1/6) * f(b))
    return y

测试一下：

def f(x):
    f = np.exp(-x)
    return f

print(Trapezoid(0.5,0.1))
print(Simpson(0.5,0.1))

结果为：

0.3022736155497186
0.298309397365031

复化梯形公式：

$T_n =\displaystyle\int^b_a f(x)dx\approx\frac{h}{2}[f(a)+f(b)+2\sum\limits_{k=1}^{n-1}f(x_k)] \\ h=\frac{b-a}{n},x_k = a + kh,(k =0,1,\dots,n)$

def complex_Trapezoid(n,a,b):
    h = (b-a) / n
    temp_s= 0
    for k in range(1,n):
        x_k = a + k*h
        temp_s += f(x_k)
    I = h/2 * (f(a) + f(b) + 2*temp_s)
    return I

复化Simpson公式：

$\displaystyle\int^b_a f(x)dx\approx\frac{h}{3}[f(a)+f(b)+2\sum\limits_{k=1}^{m-1}f(x_{2k})+4\sum\limits_{k=1}^{m}f(x_{2k-1})] \\ n=2m,h=\frac{b-a}{n},x_k = a + kh,(k =0,1,\dots,n)$

def complex_simpson(n,a,b):
    '''
    n:划分[a,b]区间成n等分
    a: 积分下限
    b: 积分上限
    '''
    h = (b-a) / n        
    temps_odd = 0
    temps_even = 0
    for k in range(1,n):
        if k % 2 == 0:
            x_k = a + k*h
            temps_even += 2 * f(x_k)
        else:
            x_k = a + k*h
            temps_odd += 4 * f(x_k)
    I = (h/3) * (f(a) + f(b) + temps_even + temps_odd)
    return I

测试一下：

def f(x):
    if x == 0:
        return 1
    else:
        f = np.sin(x) / x
        return f

print(complex_Trapezoid(8,0,1))
print(complex_simpson(8,0,1))

结果为：

0.9456908635827014
0.9460833108884719