#--------------------------------------------------------
#  Battimenti prodotti da due onde sinusoidali 
#  Per sentirli (!):
#   - aprire https://www.szynalski.com/tone-generator/
#     su due tab del browser;
#   - settare le frequenze dei due generatori 
#     (usando le freccie ai lati del valore di frequenza
#      si possono cambiare le frequenze a passi di 1Hz,
#      anche mentre il generatore è attivo);
#   - cliccare su PLAY dei due generatori
#   
#  GdA gen 2022
# --------------------------------------------------------


A1  <- 1    
A2  <- 1
nu1 <- 200  # Hz
nu2 <- 210  # Hz    # figure sul sito ottenute con 202 Hz e 210 Hz

om1 <- 2*pi*nu1
om2 <- 2*pi*nu2
T1  <-1/nu1

t.max <- 200*T1
ylim <- c(-1,1) * (A1+A2)

old.mar = par("mar")
par(mar=c(4.0,4.0,0.1,0.1))   # margini ottimizzati

curve(A1*cos(om1*t), 0, t.max, xname='t', n=10001, col='blue',
      ylim=ylim, xlab=' t  (s)', ylab='f(t)')
Sys.sleep(1)
curve(A2*cos(om2*t), xname='t', n=10001, col='green', add=TRUE)
Sys.sleep(1)
curve(A1*cos(om1*t)+A2*cos(om2*t), xname='t', n=10001, col='red', add=TRUE)

par(mar=old.mar)