离散数据求导

一、先导入函数库

import numpy as np
import matplotlib.pyplot as plt
from scipy import interpolate as inter

二、调用inter.splrep（）自动求解函数关系，k为函数的最高次幂

xx = inter.splrep(T, Y, s=0, k=5)       

二、xx为求解的函数关系，xxx为新的因变量取值集合，直接调用inter.splev求解，der为导数参数，默认为0即求原函数，为1即求一阶导，依次类推

xxx = np.linspace(0, 4, 4097)
y0o = inter.splev(xxx, xx)
yone = inter.splev(xxx, xx, der = 1)
ytwo = inter.splev(xxx, xx, der = 2)
ythree = inter.splev(xxx, xx, der = 3)

plt.plot(xxx, y0o, color="b", linewidth=0.5, label = "优化后轨迹")
plt.show()
plt.plot(xxx, yone, color="b", linewidth=0.5, label = "优化后轨迹")
plt.show()
plt.plot(xxx, ytwo, color="b", linewidth=0.5, label = "优化后轨迹")
plt.show()
plt.plot(xxx, ythree, color="b", linewidth=0.5, label = "优化后轨迹")
plt.show()