Numpy数值积分的实现

更新时间：2023年02月23日 09:21:21 作者：微小冷

本文主要介绍了Numpy数值积分的实现，文中通过示例代码介绍的非常详细，对大家的学习或者工作具有一定的参考学习价值，需要的朋友们下面随着小编来一起学习学习吧

求和

在Numpy中可以非常方便地进行求和或者连乘操作，对于形如 x 0 , x 1 , ⋯ , x_n的数组而言，其求和 ∑x_i或者连乘 ∏x_i分别通过sum和prod实现。

x = np.arange(10)
print(np.sum(x))    # 返回45
print(np.prod(x))   # 返回0

这两种方法均被内置到了数组方法中，

x += 1
x.sum()     # 返回55
x.prod()    # 返回3628800

有的时候数组中可能会出现坏数据，例如

x = np.arange(10)/np.arange(10)
print(x)
# [nan  1.  1.  1.  1.  1.  1.  1.  1.  1.]

其中x[0]由于是0/0，得到的结果是nan，这种情况下如果直接用sum或者prod就会像下面这样

>>> x.sum()
nan
>>> x.prod()
nan

为了避免这种尴尬的现象发生，numpy中提供了nansum和nanprod，可以将nan排除后再进行操作

>>> np.nansum(x)
9.0
>>> np.nanprod(x)
1.0

累加和累乘

和连加连乘相比，累加累乘的使用频次往往更高，尤其是累加，相当于离散情况下的积分，意义非常重大。

from matplotlib.pyplot as plt
xs = np.arange(100)/10
ys = np.sin(xs)
ys1 = np.cumsum(ys)/10
plt.plot(xs, ys)
plt.plot(xs, ys1)
plt.show()

效果如图所示

在这里插入图片描述

cumprood可以实现累乘操作，即

x = np.arange(1, 10)
print(np.cumprod(x))
# [     1      2      6     24    120    720   5040  40320 362880]

与sum, prod相似，cumprod和cumsum也提供了相应的nancumprod, nancumsum函数，用以处理存在nan的数组。

>>> x = np.arange(10)/np.arange(10)
<stdin>:1: RuntimeWarning: invalid value encountered in true_divide
>>> np.cumsum(x)
array([nan, nan, nan, nan, nan, nan, nan, nan, nan, nan])
>>> np.nancumsum(x)
array([0., 1., 2., 3., 4., 5., 6., 7., 8., 9.])
>>> np.nancumprod(x)
array([1., 1., 1., 1., 1., 1., 1., 1., 1., 1.])

trapz

cumsum操作是比较容易理解的，可以理解为离散化的差分，比如

>>> x = np.arange(5)
>>> y = np.cumsum(x)
>>> print(x)
array([0, 1, 2, 3, 4])
>>> print(y)
array([ 0,  1,  3,  6, 10])

trap为梯形积分求解器，同样对于[0,1,2,3,4]这样的数组，那么稍微对高中知识有些印象，就应该知道[0,1]之间的积分是，此即梯形积分

>>> np.trapz(x)
8.0

接下来对比一下trapz和cumsum作用在 sin ⁡ x \sin x sinx上的效果

from matplotlib.pyplot as plt
xs = np.arange(100)/10
ys = np.sin(xs)
y1 = np.cumsum(ys)/10
y2 = [np.trapz(ys[:i+1], dx=0.1) for i in range(100)]
plt.plot(xs, y1)
plt.plot(xs, y2)
plt.show()

结果如图，可见二者差别极小。

在这里插入图片描述

到此这篇关于Numpy数值积分的实现的文章就介绍到这了,更多相关Numpy数值积分内容请搜索脚本之家以前的文章或继续浏览下面的相关文章希望大家以后多多支持脚本之家！

您可能感兴趣的文章:

对python实时得到鼠标位置的示例讲解
今天小编就为大家分享一篇对python实时得到鼠标位置的示例讲解，具有很好的参考价值，希望对大家有所帮助。一起跟随小编过来看看吧
2018-10-10
python中如何读取xlsx文件
这篇文章主要介绍了python中如何读取xlsx文件问题,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教
2024-02-02
python语音识别的转换方法
语音识别技术，也被称为自动语音识别，目标是以电脑自动将人类的语音内容转换为相应的文字。应用包括语音拨号、语音导航、室内设备控制、语音文档检索、简单的听写数据录入等。本文给大家介绍python语音识别的方法，感兴趣的朋友一起看看吧
2021-10-10
在Keras中CNN联合LSTM进行分类实例
这篇文章主要介绍了在Keras中CNN联合LSTM进行分类实例，具有很好的参考价值，希望对大家有所帮助。一起跟随小编过来看看吧
2020-06-06
PyTorch变分自编码器的构建与应用小结
变分自编码器是一种强大的深度学习模型,用于学习数据的潜在表示并能生成新的数据点,使用PyTorch实现VAE不仅可以加深对生成模型的理解,还可以利用其灵活性进行各种实验,这篇文章主要介绍了PyTorch变分自编码器的构建与应用,需要的朋友可以参考下
2024-07-07
seek引发的python文件读写的问题及解决
这篇文章主要介绍了seek引发的python文件读写的问题及解决,文中通过示例代码介绍的非常详细，对大家的学习或者工作具有一定的参考学习价值,需要的朋友可以参考下
2019-07-07
Python 等分切分数据及规则命名的实例代码
这篇文章主要介绍了Python 等分切分数据及规则命名的实例代码，代码简单易懂，非常不错，具有一定的参考借鉴价值,需要的朋友可以参考下
2019-08-08
Python用tkinter实现自定义记事本的方法详解
这篇文章主要为大家详细介绍了Python用tkinter实现自定义记事本的方法，文中示例代码介绍的非常详细，具有一定的参考价值，感兴趣的小伙伴们可以参考一下，希望能够给你带来帮助
2022-03-03
Python实现计算图像RGB均值方式
这篇文章主要介绍了Python实现计算图像RGB均值方式，具有很好的参考价值，希望对大家有所帮助。一起跟随小编过来看看吧
2020-06-06
Python 常用的安装Module方式汇总
这篇文章主要介绍了Python 常用的安装Module方式汇总,需要的朋友可以参考下
2017-05-05