# Clear the Workspace rm (list=ls()) # Generate a 5 second 1 kHz signal at 44.1 kHz sampling f <- 1000 amp <- 26213 tau <- 0.015 omega <- 2*pi*f time <- 5 rate <- 44100 len <- time *rate # Increment time samples at 44.1 kHz t <- seq(rate*time)/rate audio <- amp*sin(omega*t)*(1-exp(-t/tau))*(1-exp(-(time-t)/tau)) # Append zero a 1/2 second before and 1/2 second after t1 <- seq(rate/2)*0 final <- as.integer(c(t1, audio, t1)) # Save the file to Pulse1k.txt write( final, file="Pulse1k.txt", ncolumns = 1 )