set.seed(42)
n <- 600
t <- seq_len(n)
p1 <- 30
p2 <- 240
xy <- (1 + 0.6 * sin(t*2*pi/p2)) * sin(t*2*pi/p1) + 2 * sin(t*2*pi/p2) +
rnorm(n, sd = 0.5) + t * 0.01
inter_dt <- round(runif(length(xy), min = 0.5, max = 1.5),1)
dt <- cumsum(inter_dt)
dec <- extricate(xy, dt, nimf = 7, sifting = 10,
repl = 1, comb = 40, factor_noise = 10,
speak = TRUE, output_sifting = TRUE)
symmetry(dec$m)
plot_emd(dec, select = c(6,8,9), pdf = FALSE, adapt.axis = TRUE)
Run the code above in your browser using DataLab