Fig5.5s <- function() { # postscript(file = "Fig5.5.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)) 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) time.scale.plot(dwt(y, wavelet = "d10")) # graphics.off() NULL }