通过SciPy计算导数(Calculating derivative by SciPy)

通过SciPy计算导数(Calculating derivative by SciPy)

考虑一下这个python脚本

import numpy as np import matplotlib.pyplot as plt x = np.linspace(0,1,100) def y(x,z): return x**2 + z**-1 - 1 for z in np.arange(1,50,1): plt.plot(x, y(x,z))

在固定的x，我该如何绘制

plt.plot(z, y(x,z)) plt.plot(z, y(x,z).derivative) #dy/dz

Consider this python script

import numpy as np import matplotlib.pyplot as plt x = np.linspace(0,1,100) def y(x,z): return x**2 + z**-1 - 1 for z in np.arange(1,50,1): plt.plot(x, y(x,z))

At a fixed x, how can I draw

plt.plot(z, y(x,z)) plt.plot(z, y(x,z).derivative) #dy/dz

最满意答案

您必须使用derivative函数：

import numpy as np import matplotlib.pyplot as plt from scipy.misc import derivative x = np.linspace(0,1,100) def y(z, x): return x**2 + z**(-1) - 1 z = np.arange(1.0, 50.0) dydz = [derivative(y, zi, args=(x, )) for zi in z] plt.plot(z, dydz) plt.show()

截图：

You must use the derivative function:

import numpy as np import matplotlib.pyplot as plt from scipy.misc import derivative x = np.linspace(0,1,100) def y(z, x): return x**2 + z**(-1) - 1 z = np.arange(1.0, 50.0) dydz = [derivative(y, zi, args=(x, )) for zi in z] plt.plot(z, dydz) plt.show()

Screenshot: