unwrap函数详解

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

函数的完整格式:

unwrap(pha,tol,dim);

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

unwrap(pha,[],1)

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

当一列数据为:
Matlab中unwrap函数内容详解-编程知识网
其unwrap后结果为:
Matlab中unwrap函数内容详解-编程知识网
当 2.5025与 -2.9863之间超过阈值 pi 后视为发生跳变,将2.5025加上或减去 2*pi,以此类推。
分享一下我当时的思路,最终以纯c代码实现,经测试与matlab中完全无误。
具体实现代码和步骤解析如下链接:https://download.csdn.net/download/cao_jie_xin/11208965