python 离散数值点 求导

在 Python 中，可以使用 NumPy 库中的 diff() 函数来计算离散数值点的导数。diff() 函数可以计算两个相邻元素的差值，并返回一个数组。

例如，如果有一个包含离散数值点的数组 x，以及对应的函数值数组 y，可以使用以下代码计算 y 对 x 的导数：

import numpy as np

# 定义离散数值点 x 和对应的函数值数组 y
x = np.array([1, 2, 3, 4, 5])
y = np.array([1, 4, 9, 16, 25])

# 计算 y 对 x 的导数
dy_dx = np.diff(y) / np.diff(x)

# 打印结果
print(dy_dx)

输出结果为：

[3. 5. 7. 9.]

这个结果表示 y 对 x 的导数在 x=1 处为 3，x=2 处为 5，x=3 处为 7，x=4 处为 9。注意，由于 diff() 函数计算的是相邻元素的差值，因此返回的导数数组比原始数组少一个元素。