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
输出数组。
注释
如果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.])