powered by
kzsv(y, kza, kz, q, k=3, m=round(0.05*q), tol=1.0e-5)
x <- c(rep(0,4000),rep(0.5,2000),rep(0,4000)) noise <- rnorm(n = 10000, sd = 1.0) # normally-distributed random variates v <- x + noise kzsv(v, q=1000, k=3)
Run the code above in your browser using DataLab