R/Anton.R

In neuroconductor-devel/waveslim: Basic Wavelet Routines for One-, Two-, and Three-Dimensional Signal Processing

AntonB <- function() {

  a0 <- c(0,
          0.02674875741081,
          -0.01686411844287,
          -0.07822326652899,
          0.26686411844288,
          0.60294901823636,
          0.26686411844287, 
          -0.07822326652899,
          -0.01686411844287,
          0.02674875741081,
          0,
          0)

  a1 <- c(0,
          0,
          0,
          0.04563588155712,
          -0.02877176311425,
          -0.29563588155712,
          0.55754352622850,
          -0.29563588155713,
          -0.02877176311425,
          0.04563588155712,
          0,
          0)

  s0 <- c(0,
          0,
          0,
          -0.04563588155712,
          -0.02877176311425,
          0.29563588155712,
          0.55754352622850,
          0.29563588155713,
          -0.02877176311425,
          -0.04563588155712,
          0,
          0)

  s1 <- c(0,
          0.02674875741081,
          0.01686411844287,
          -0.07822326652899,
          -0.26686411844288,
          0.60294901823636,
          -0.26686411844287,
          -0.07822326652899,
          0.01686411844287,
          0.02674875741081,
          0,
          0)

  s0 <- 2 * s0
  s1 <- 2 * s1

  aa0 <- c(0,
           0,
           0.02674875741081,
           -0.01686411844287,
           -0.07822326652899,
           0.26686411844288,
           0.60294901823636,
           0.26686411844287,
           -0.07822326652899,
           -0.01686411844287,
           0.02674875741081,
           0)

  aa1 <- c(0, 
           0,
           0,
           0,
           0.04563588155712,
           -0.02877176311425,
           -0.29563588155712,
           0.55754352622850,
           -0.29563588155713,
           -0.02877176311425,
           0.04563588155712,
           0)

  ss0 <- c(0,
           0,
           -0.04563588155712,
           -0.02877176311425,
           0.29563588155712,
           0.55754352622850,
           0.29563588155713,
           -0.02877176311425,
           -0.04563588155712,
           0,
           0,
           0)

  ss1 <- c(0.02674875741081,
           0.01686411844287,
           -0.07822326652899,
           -0.26686411844288,
           0.60294901823636,
           -0.26686411844287,
           -0.07822326652899,
           0.01686411844287,
           0.02674875741081,
           0,
           0,
           0)

  ss0 <- 2 * ss0
  ss1 <- 2 * ss1

  list(af = list(cbind(a0, a1), cbind(aa0, aa1)),
       sf = list(cbind(s0, s1), cbind(ss0, ss1)))
}