Fig5.6s <- function() { # postscript(file = "Fig5.6.ps", height = 5.46, width = 5.46, horiz = F) par(mgp = c(1.4, 0.4, 0), mar = c(1.5, 1.5, 0, 0)) rs <- c(57, 14, 55, 51, 30, 0, 53, 44, 34, 53, 49, 2) .Random.seed <- rs x <- (1:1024)/1024 y <- sin(2 * pi * x) y[513:1024] <- y[513:1024] + sin(4 * pi * x[513:1024]) y[769:1024] <- y[769:1024] + ((1:256) * sin((10 * pi * (1:256))/256)) / 128 y <- y + (1:1024)/1024 y[201:280] <- y[201:280] + sin(((1:80) * 30 * pi)/80) y <- y + rnorm(1024) * 0.131366 time.scale.plot(dwt(y, wavelet = "d10")) # graphics.off() NULL }