phunwrap: Phase unwrapping

In jrevenaugh/TSAUMN: Time Series Analysis Tools Used in ESCI5201 at the University of Minnesota

Description

Unwrap a vector of phase values using Itoh's algorithm.

Usage

phunwrap(p)

Arguments

p	vector of phase values

Value

unwrapped phase values

Examples

s <- rep( 0, 124 )
s[50] <- 1
S <- fft( s )
p <- Arg( S )
up <- phunwrap( p[1:63] )
f <- seq( 0, 0.5, length.out = 63 )
plot( f, up, xlab = "Frequency", ylab = "Phase (rad)", type = "l", main = "Estimated (black) and Theory (red)" )
lines( f, -2 * pi * f * 49, col = "red" )

jrevenaugh/TSAUMN documentation built on Nov. 8, 2019, 2:20 p.m.