numpy.unwrap#

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

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

这通过改变那些与前驱元素的绝对差大于max(discont, period/2)的元素来解开信号p,将其值改为其period互补值。

对于默认情况,其中period\(2\pi\)discont\(\pi\),这将解开一个弧度相位p,使得相邻差异永远不会大于\(\pi\),方法是为某个整数\(k\)添加\(2k\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.])