python计算定积分

聪明乖巧的小狮子

已于 2022-03-21 16:59:37 修改

阅读量1.7w

点赞数 11

分类专栏： python 计算方法文章标签： python 矩阵

于 2022-03-16 15:08:25 首次发布

本文链接：https://blog.csdn.net/m0_52295416/article/details/123526780

版权

python 同时被 2 个专栏收录

16 篇文章 2 订阅

订阅专栏

计算方法

6 篇文章 2 订阅

订阅专栏

这篇博客介绍了如何使用Python的sympy库和numpy库进行数值积分。首先，通过sympy的integrate函数直接求解了从0到1的ex的定积分。其次，采用数值积分的方法，将积分区间分成100000等份，计算每个小矩形的面积并累加，从而近似得到积分值。两种方法都展示了Python在数学计算上的应用。

摘要由CSDN通过智能技术生成

$求\begin{matrix} \int_{0}^{1} e^{x}\, dx\end{matrix}$
方法一：

from sympy import *
import math
x = symbols('x')
print(integrate(math.e**x, (x, 0, 1)))

方法二：
首先把积分区间分成N小段，然后计算每段间隔对应的小矩形面积（底为dx，高为对应的函数值），接着将它们累加求和，得到的积分值。

import numpy as np
#进行积分时，现将积分区间分成N段，N在不超限度时越大结果越精确
N=100000
x=[]
for k in range(0,N):
    x.append(1/N*k)  #将积分区间分成n段放入x数组内
dx=[]  #dx数组代表微分值
y=[]   #y代表N个离散化的x的对应的函数值
#循环得到完整的dx，y数组
for k in range(0,N-1):
    y.append(np.exp(x[k]))
    dx.append(x[k+1]-x[k])
#将x、y、dx转化成numpy的矩阵格式
x=np.matrix(x)
y=np.matrix(y)
dx=np.matrix(dx)
#对应位置的dx，y相乘，累加起来就是积分值
I_num=np.sum(np.multiply(y,dx))