Unwrap lower than the default value (pi)

21 次查看（过去 30 天）
Pranjal Agrawal
Pranjal Agrawal 2023-8-7
评论： Bryan Wingert 2024-7-16
Hi! Can anyone please elaborate on how to eliminate jumps that are less than pi?I have a function which varies from (0,pi/2). I want to unwrap the values in between that range. Can anyone please suggest on how to use this.
Thanks
  2 个评论
Stephen23
Stephen23 2023-8-7
multiply by four, UNWRAP, divide by four
Bryan Wingert
Bryan Wingert 2024-7-16
This is a pretty clever hack of the unwrap function. It scales jumps from the threshold (pi/2 in this case) up to 2*pi, which guarantees that it gets unwrapped, since the unwrap function looks for jumps of 2*pi. (The threshold represents the bidirectional bound, e.g. [-pi,pi] for the default.) More specifically, the scaling factor is approximately 2*pi/threshold. Approximately because numerical error can throw it off sometimes.

请先登录，再进行评论。

请先登录，再回答此问题。

回答（1 个）

Sanjana
Sanjana 2023-8-21
Hi Pranjal,
I understand that you are looking for a way to eliminate jumps that are less than pi. As per the official documentation of the,”unwrap” function, if a jump threshold less than pi is specified, then unwrap function uses the default jump threshold pi.
For eliminating jumps that are less than pi, a custom unwrapping algorithm can be implemented as follows,
% Let P be the phase vector
%Intialize unwrapped phase vector
unwrapped = zeros(size(P));
unwrapped (1) = P(1);
for i = 2:length(P)
diff = P(i)-P(i-1);%Calculating the jump
%check if jump is greater than threshold T
if diff > T
unwrapped(i) = P(i) - 2*pi;
elseif diff < -T
unwrapped(i) = P(i) + 2*pi;
else
unwrapped(i) = P(i);
end
end
Please refer to this following link, for further information,
https://www.mathworks.com/help/matlab/ref/unwrap.html
Hope this helps!
Regards,
Sanjana
  1 个评论
Bryan Wingert
Bryan Wingert 2024-7-16
This only unwraps the point next to the jump, and leaves the rest of the vector untouched. The unwrap function cumulatively shifts the rest of the vector. Use Stephen23's hack above.

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Spectral Measurements 的更多信息

标签

产品

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by