Unwrap spectrum phase
unwrap the phase spectrum so it does not wrap around
LocalUnwrap(p, cutoff = cutoff)
p |
phase spectrum |
cutoff |
cut off angle = pi |
Unwrapped spectrum
Algorithm minimizes the incremental phase variation by constraining it to the range [-pi,pi]
Jonathan M. Lees<jonathan.lees.edu>
x <- 1:512 amp <- sin(1*2*pi*x/16) + sin(2*2*pi*x/16) + sin(3*2*pi*x/16) spc <- fft(amp) plot(Mod(spc), type='l') angle <- Arg(spc) plot(angle, type='l') unang <- LocalUnwrap(angle, cutoff =pi ) plot(unang, type='l')
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.