scipy模块计算导数方法(central_diff

scipy中计算导数有两种方式:

central_diff_weights
derivative

其中第一种方式在scipy帮助中，没有写很清楚，这里重点讲一下。

就举一个例子: 计算下列函数在x=1x=1x=1处的2阶导数
f(x)=3x+x3f(x) = 3^{x} + x^{3}f(x)=3x+x3
利用求导公式，我们很容易得到这个值：9.620846882437746

要利用第一种方法，需要有若干个(N)在求导点附近的函数值，并且需要均匀。N还需要满足两个条件

N为奇数
N大于导数的阶数；比如要计算1阶导数，N>1；计算2阶导数， N>2

我们取N=5, 计算在x=1x=1x=1的附近的5个函数值，左右对称
f(0.8),f(0.9),f(1),f(1.1),f(1.2)f(0.8) ,f(0.9), f(1), f(1.1), f(1.2)f(0.8),f(0.9),f(1),f(1.1),f(1.2)

central_diff_weights(5,2)，返回的是计算2阶导数的各个值的权重，把这5个值和以上5个函数值相乘，并且除以间隔值0.1两次，即为导数

from scipy.misc import central_diff_weights
from scipy.misc import derivative
import numpy as npdef f(x):return 3**x + x**3x = np.r_[0.8:1.2:5j]
y = np.vectorize(f)(x)w = central_diff_weights(5, 2)
print(np.sum(y*w) / 0.1 / 0.1)

计算的值为：9.620841015472953 与理论值很接近了

下面用另一种方法derivative 计算导数，这种方法比较直观

derivative(f, x0=1, dx=0.0001, n=2)

计算的值为：9.620846874724975

scipy模块计算导数方法(central_diff_weights)相关推荐

python卷积计算_Python使用scipy模块实现一维卷积运算示例
本文实例讲述了python使用scipy模块实现一维卷积运算.分享给大家供大家参考,具体如下: 一介绍 signal模块包含大量滤波函数.b样条插值算法等等.下面的代码演示了一维信号的卷积运算. 二 ...
python中Scipy模块求取积分
python中Scipy模块求取积分的方法: SciPy下实现求函数的积分的函数的基本使用,积分,高等数学里有大量的讲述,基本意思就是求曲线下面积之和. 其中rn可认为是偏差,一般可以忽略不计,wi可 ...
python模块下载过程_常用的python模块及安装方法
bsddb3:BerkeleyDB的连接组件 Cheetah-1.0:我比较喜欢这个版本的cheetah cherrypy:一个WEB framework ctypes:用来调用动态链接库 DBUti ...
常用的python模块及安装方法
adodb:我们领导推荐的数据库连接组件 bsddb3:BerkeleyDB的连接组件 Cheetah-1.0:我比较喜欢这个版本的cheetah cherrypy:一个WEB framework c ...
【Scipy高级计算】(2) 常微分方程、洛伦兹吸引子，附python完整代码
大家好,在上一篇博文中,我介绍了如何使用Scipy库计算定积分和二重积分,感兴趣的可以看一下:https://blog.csdn.net/dgvv4/article/details/124226759 ...
Python之常用模块及安装方法
adodb:我们领导推荐的数据库连接组件 bsddb3:BerkeleyDB的连接组件 Cheetah-1.0:我比较喜欢这个版本的cheetah cherrypy:一个WEB framework c ...
浅尝则止 - SciPy科学计算
21. 浅尝则止 - SciPy科学计算 SciPy以NumPy为基础,提供了众多数学.科学.工程计算用的模块,包括但不限于:线性代数.常微分方程求解.信号处理.图像处理.稀疏矩阵处理. 本章适合那些 ...
python二元函数图像在线绘制_Python实现的拟合二元一次函数功能示例【基于scipy模块】...
Python实现的拟合二元一次函数功能示例[基于scipy模块] 发布时间:2020-08-31 00:46:24 来源:脚本之家阅读:121 本文实例讲述了Python实现的拟合二元一次函数功能. ...
python 三维凸包_浅尝则止 - SciPy科学计算 in Python
本文节选自作者的<Python编程基础及应用>视频教程.Python编程基础及应用_哔哩哔哩 (゜-゜)つロ干杯~-bilibiliwww.bilibili.com 21. 浅尝则止 ...

scipy模块计算导数方法(central_diff_weights)

scipy模块计算导数方法(central_diff_weights)相关推荐

最新文章

热门文章