numpy.unwrap#

numpy.unwrap(p, discont=None, axis=-1, *, period=6.283185307179586)[source]#

通过对相对于周期的较大增量取补码来解开。

这通过将与前一个元素的绝对差值大于 max(discont, period/2) 的元素更改为其 period 补码值来解开信号 p

对于默认情况下 period\(2\pi\)discont\(\pi\) 的情况,这通过为某个整数 \(k\) 添加 \(2k\pi\) 来解开一个弧度相位 p,使得相邻差值永远不会大于 \(\pi\)

参数:
parray_like

输入数组。

discontfloat, 可选

值之间的最大不连续性,默认值为 period/2。小于 period/2 的值将被视为 period/2。若要产生与默认值不同的效果,discont 应大于 period/2

axisint, 可选

unwrap 将在其上执行操作的轴,默认值为最后一个轴。

periodfloat, 可选

输入包装的范围大小。默认情况下,它是 2 pi

版本 1.21.0 中的新增功能。

返回值:
outndarray

输出数组。

另请参阅

rad2deg, deg2rad

备注

如果 p 中的不连续性小于 period/2,但大于 discont,则不会执行解开,因为取补码只会使不连续性更大。

示例

>>> import numpy as np
>>> phase = np.linspace(0, np.pi, num=5)
>>> phase[3:] += np.pi
>>> phase
array([ 0.        ,  0.78539816,  1.57079633,  5.49778714,  6.28318531]) # may vary
>>> np.unwrap(phase)
array([ 0.        ,  0.78539816,  1.57079633, -0.78539816,  0.        ]) # may vary
>>> np.unwrap([0, 1, 2, -1, 0], period=4)
array([0, 1, 2, 3, 4])
>>> np.unwrap([ 1, 2, 3, 4, 5, 6, 1, 2, 3], period=6)
array([1, 2, 3, 4, 5, 6, 7, 8, 9])
>>> np.unwrap([2, 3, 4, 5, 2, 3, 4, 5], period=4)
array([2, 3, 4, 5, 6, 7, 8, 9])
>>> phase_deg = np.mod(np.linspace(0 ,720, 19), 360) - 180
>>> np.unwrap(phase_deg, period=360)
array([-180., -140., -100.,  -60.,  -20.,   20.,   60.,  100.,  140.,
        180.,  220.,  260.,  300.,  340.,  380.,  420.,  460.,  500.,
        540.])