Advertisement
Guest User

Untitled

a guest
Mar 24th, 2017
63
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.49 KB | None | 0 0
  1. # Sampling frequency
  2. fs <- 1000
  3. # Sampling period
  4. ts <- 1/fs
  5. # Length of signal
  6. l <- 1500
  7. # Time vector
  8. t <- seq(0, l - 1)*ts
  9. # Signal
  10. s <- 0.4*sin(2*pi*50*t) + sin(2*pi*120*t)
  11. # Noise
  12. x <- s + 2*rnorm(length(t))
  13. # Plot signal
  14. qplot(t, x) + geom_line()
  15.  
  16. # FFT
  17. y <- fft(x)
  18. # Two-sided spectrum
  19. p2 <- abs(y/l)
  20. # Convert to one-sided spectrum
  21. p1 <- p2[1:(l/2 + 1)]
  22. p1[2:(length(p1) - 1)] <- 2*p1[2:(length(p1) - 1)]
  23. # Define frequency
  24. f <- fs*seq(0, l/2)/l
  25. # Plot spectrum
  26. qplot(f, p1) + geom_line()
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement