python 导数 vector_Python。Scipy数值导数

不久前我发现了scipy.misc.derivative。在def derivative(func, x0, dx=1.0, n=1, args=(), order=3)

令人难以置信的是它的速度和精确度。

所以，我决定弄清楚它是怎么工作的。

但是我不完全理解这个函数的代码，它在derivative中使用：

^{pr2}$

据我所知，这个函数在高阶导数的展开式中找到系数，但是它是如何做到的呢？在

代码：def central_diff_weights(Np, ndiv=1):

"""

Return weights for an Np-point central derivative.

Assumes equally-spaced function points.

If weights are in the vector w, then

derivative is w[0] * f(x-ho*dx) + ... + w[-1] * f(x+h0*dx)

Parameters

----------

Np : int

Number of points for the central derivative.

ndiv : int, optional

Number of divisions. Default is 1.

Notes

-----

Can be inaccurate for large number of points.

"""

if Np < ndiv + 1:

raise ValueError("Number of points must be at least the derivative order

+ 1.")

if Np % 2 == 0:

raise ValueError("The number of points must be odd.")

from scipy import linalg

ho = Np >> 1

x = arange(-ho,ho+1.0)

x = x[:,newaxis]

X = x**0.0

for k in range(1,Np):

X = hstack([X,x**k])

w = product(arange(1,ndiv+1),axis=0)*linalg.inv(X)[ndiv]

return w

python 导数 vector_Python。Scipy数值导数相关推荐

干货 | 用Python做图像处理：图像导数实战
导读:数量庞大的图像和视频充斥着我们的生活.我们需要对图片进行检索.分类等操作时,利用人工手段显然是不现实的.于是,计算机视觉相关技术便应运而生,并且得到了快速的发展.本文以时下最流行的Python语 ...
python编程求导数_SciPy函数求导数
17. SciPy求函数的导数在SciPy里提供了很多的方法函数可以实现对某函数进行求导和求积分的操作. SciPy的求导相对简单也容易理解.已知函数$f(x)$求其在$x_0$的导数即$f'(x_ ...
pilt图像处理_干货 | 用Python做图像处理：图像导数实战
导读:数量庞大的图像和视频充斥着我们的生活.我们需要对图片进行检索.分类等操作时,利用人工手段显然是不现实的.于是,计算机视觉相关技术便应运而生,并且得到了快速的发展.本文以时下最流行的Python语 ...
python解决微分方程（数值解法）
Python求解微分方程(数值解法) 对于一些微分方程来说,数值解法对于求解具有很好的帮助,因为难以求得其原方程. 比如方程: 但是我们知道了它的初始条件,这对于我们叠代求解很有帮助,也是必须的. 那 ...
python中的scipy基础知识_python中SciPy是什么？
python中Numpy常用于计算二维数组计算,而python的另一个库SciPy库与Numpy有着密切的关系,是需要通过Numpy为基础,同时也是通过Numpy数据来操控科学计算.常见的是插值运算. ...
Python库之Scipy库的简介、安装详细
目录 Scipy库的简介 Scipy库的安装 Scipy库的简介 Scipy高级科学计算库:和Numpy联系很密切,Scipy一般都是操控Numpy数组来进行科学计算.统计分析,所以可以说是基于Num ...
Python语言学习之数值、小数、空格那些事：python和数值、小数、空格的使用方法之详细攻略
Python语言学习之数值.小数.空格那些事:python和数值.小数.空格的使用方法之详细攻略目录 Python与数值那些事 1.python保留两位小数/保留小数点位数 Python去掉空格的方 ...
python语言整数类型-Python 的内置数值类型
Python 是一种敏捷的.动态类型化的.极富表现力的开源编程语言,可以被自由地安装到多种平台上.Python 代码是被解释的.如果您对编辑.构建和执行循环较为熟悉,则 Python 代码对您来说更简 ...
7.边缘检测：2D运算——回顾、高斯滤波器2D的导数、Sigma对导数的影响_1
目录回顾高斯滤波器2D的导数 Sigma对导数的影响回顾我们要完成我们的边缘检测这个单元,然后它会被用在你以后要做的事情上. 上单元我们讲了边的概念以及它们是如何与梯度和函数导数的大小相关的. ...

python 导数 vector_Python。Scipy数值导数

python 导数 vector_Python。Scipy数值导数相关推荐

最新文章

热门文章