unwrap函数详解

一般在我们计算一个系统相频特性时，就要用到反正切函数提取相位，计算机中反正切函数规定，在一、二象限中的角度为0～pi，三四象限的角度为0～-pi。但实际得到的结果会发生相位跳变，跳变幅度为2pi，这就叫相位的卷绕。unwrap函数的作用就是解卷绕，使相位在pi处不发生跳变，从而反应出真实的相位变化。
实际在默认的情况下，unwrap在检查到数据前后两点的差距在超过 pi 的时候，就认为有跳变，当然其阈值也能通过函数设置。

函数的完整格式：

unwrap(pha,tol,dim)；

pha 可为数列或矩阵，tol 为判断的阈值，dim 参数可以设置对矩阵的行，还是列进行unwrap操作，
如果对列每一列操作，那么第三个参数可以省略或填1
而要对每一行操作，那么第三个参数填2

unwrap(pha,[],1)

表示对矩阵数据每一列进行相位矫正，使用默认的检测跳变标准 pi。

当一列数据为：

其unwrap后结果为：

当 2.5025与 -2.9863之间超过阈值 pi 后视为发生跳变，将2.5025加上或减去 2*pi，以此类推。
分享一下我当时的思路，最终以纯c代码实现，经测试与matlab中完全无误。
具体实现代码和步骤解析如下链接：https://download.csdn.net/download/cao_jie_xin/11208965