Numpy ndimensional数组的二阶导数

本文将介绍如何在Numpy中计算ndimensional数组的二阶导数。我们将从一维数组开始，并介绍如何在n维数组中计算二阶导数。

阅读更多：Numpy 教程

一维数组

在一维数组中，我们可以使用内置函数numpy.diff()来计算一阶导数。而计算二阶导数则可以通过对一阶导数再次使用numpy.diff()函数来实现。下面是一个示例代码：

import numpy as np

# 生成一维数组
arr = np.array([1, 2, 4, 7, 11, 16])

# 计算一阶导数
deriv1 = np.diff(arr)
print("一阶导数：", deriv1)

# 计算二阶导数
deriv2 = np.diff(deriv1)
print("二阶导数：", deriv2)

输出结果如下：

一阶导数： [1 2 3 4 5]
二阶导数： [1 1 1 1]

多维数组

在多维数组中，我们需要使用numpy.gradient()函数来计算一阶导数。和在一维数组中一样，计算二阶导数则可以通过再次使用numpy.gradient()函数来实现。下面是一个示例代码，其中包含了一个2×3的二维数组和一个2x3x3的三维数组：

import numpy as np

# 生成二维数组
arr_2d = np.array([[1, 2, 4], [7, 11, 16]])

# 计算一阶导数
deriv1_2d = np.gradient(arr_2d)
print("一阶导数：", deriv1_2d)

# 计算二阶导数
deriv2_2d_x = np.gradient(deriv1_2d[0])
deriv2_2d_y = np.gradient(deriv1_2d[1])
print("二阶导数x：", deriv2_2d_x)
print("二阶导数y：", deriv2_2d_y)

# 生成三维数组
arr_3d = np.array([[[1, 2, 4], [7, 11, 16], [2, 5, 9]], [[3, 6, 10], [8, 14, 20], [4, 9, 15]]])

# 计算一阶导数
deriv1_3d = np.gradient(arr_3d)
print("一阶导数：", deriv1_3d)

# 计算二阶导数
deriv2_3d_x = np.gradient(deriv1_3d[0])
deriv2_3d_y = np.gradient(deriv1_3d[1])
deriv2_3d_z = np.gradient(deriv1_3d[2])
print("二阶导数x：", deriv2_3d_x)
print("二阶导数y：", deriv2_3d_y)
print("二阶导数z：", deriv2_3d_z)

输出结果如下：

一阶导数： [array([[6., 9., 12.],
       [6., 9., 12.]]), array([[1., 2., 2.],
       [3., 3., 4.]])]
二阶导数x： [array([0., 0., 0.]), array([2., 1., 0.])]
二阶导数y： [array([0., 0., 0.]), array([0., 0., 0.])]
一阶导数： [array([[ 6.,  9., 12.],
       [ 1.,  3.,  4.],
       [-5., -6., -7.]]), array([[ 1.,  3.,  5.],
       [ 7., 10., 11.],
       [ 1.,  4.,  6.]]), array([[-1., -3., -5.],
       [ 1.,  3.,  5.],
       [ 2.,  4.,  6.]])]
二阶导数x：

[array([[0., 0., 0.],
        [0., 0., 0.],
        [0., 0., 0.]]), array([[ 0., -3., -3.],
        [ 3.,  3.,  3.],
        [ 3.,  3.,  3.]]), array([[2., 1., 0.],
        [2., 1., 0.],
        [2., 1., 0.]])]

二阶导数y：

[array([[0., 0., 0.],
        [0., 0., 0.],
        [0., 0., 0.]]), array([[0., 0., 0.],
        [0., 0., 0.],
        [0., 0., 0.]]), array([[ 0., -2., -2.],
        [ 0., -2., -2.],
        [ 0., -2., -2.]])]

二阶导数z：

[array([[0., 0., 0.],
        [0., 0., 0.],
        [0., 0., 0.]]), array([[0., 0., 0.],
        [0., 0., 0.],
        [0., 0., 0.]]), array([[ 0.,  0.,  0.],
        [ 0.,  0.,  0.],
        [-1., -1., -1.]])]