phaseplot2 {seewave}

Phase-phase 2D plot of a time wave

Description

This functions returns a 2D representation of a time wave against a delayed version of itself.

Usage

phaseplot2(wave, f, channel = 1, tau = 1, type = "l",
xlab = "x(t)",
ylab = paste("x(t+", tau, ")", sep = ""), ...)

Arguments

wave

an R object.

f

sampling frequency of wave (in Hz). Does not need to be specified if embedded in wave.

channel

channel of the R object, by default left channel (1).

tau

the time delay to apply in number of samples.

type

type of plot that should be drawn. See plot for details (by default "l" for lines).

xlab

title of the x axis.

ylab

title of the y axis.

...

other plot parameters.

Details

The principle consists in displaying in a single x-y graph the original time wave with a delayed version of itself. The delay is controlled with the argument tau that needs to be specified in number of samples. The conversion of tau in second is obtained by calculating tau/f, with f the sampling frequency.

Value

Nothing is returned except an x-y plot.

Note

Phase-phase plot can be used to test non-linearity.

Author(s)

Jerome Sueur

References

Kantz H, Schreiber T (2003) Non linear time series analysis. Cambridge University Press.

See Also

phaseplot

Examples

s <- synth(d=0.05, f=44100, cf=440, out="Wave")
n <- noisew(d=0.05, f=44100, out="Wave")
par(mfrow=c(2,1))
phaseplot2(s)
phaseplot2(n)