# program to produce white = independent (= uncorrelated) noise secs = 3 sr = 8000 N = 8192 # fft len y = runif(secs*sr,min = -.8, max=.8) # uniformly dist random numbers #y = rnorm(secs*sr,sd=.8) # uniformly dist normal(0,sd) random numbers plot(y) u = Wave(round(2^(bits-5)*y), samp.rate = sr, bit=16) play(u,"play") Y = fft(y[1:N]) plot(Mod(Y[1:N/2])) # observe modulus of FFT