标签：function plt 4.1 def 微分 12 数值 return 囚徒
数值微分近似
#囚徒4.0_13_数值微分近似
import numpy as np
import matplotlib.pylab as plt
#求 数值微分 导数
def numerical_diff(f, x):
h = 1e-4 # 0.0001
return (f(x+h) - f(x-h)) / (2*h)
#函数
def function_1(x):
return 0.01*x**2 + 0.1*x
def tangent_line(f, x):
d = numerical_diff(f, x) #数值微分
print(d)
y = f(x) - d*x
return lambda t: d*t + y #输入是t输出是d*t + y 这里是条直线
x = np.arange(0.0, 20.0, 0.1)
y = function_1(x)
plt.xlabel("x")
plt.ylabel("f(x)")
#在这里求取在 x=5这一点的数值微分的直线
tf = tangent_line(function_1, 5)
y2 = tf(x)
plt.plot(x, y)
plt.plot(x, y2)
plt.show()
标签：function,
plt,
4.1,
def,
微分,
12,
数值,
return,
囚徒
From： https://www.cnblogs.com/qt-pyq/p/17840986.html