Rwavelet Vignette
Christophe Chesneau and Fabien Navarro
2020-12-11
1 1d examples
1.1 Load and plot a signal from the MakeSignal function
J <- 10; n <- 2^J; t <- (1:n) / n
name <- c('Bumps')
f <- MakeSignal(name, n)
plot(t, f, xlab="t", ylab="f(t)",
type='l', lwd=1.2, main=name)
1.2 Plot a noisy version of the original function
SNR <- 4
set.seed(1)
ssig <- sd(f)
sigma <- ssig / SNR
y <- f + rnorm(n, mean=0, sd=sigma)
plot(t, y, xlab="t", ylab="y", main=paste("Noisy", name))