#--------------------------------------------------------
#  Somma di due onde sinusoidali fasate 
#  rispetto a sfasamento.R viene plottata  
#  anche la curva della somma usando la formula,
#  ottenuta, per semplificare i conti, per il
#  solo caso in cui le ampiezze delle due onde sono uguali
#
#   GdA gen 2022
# --------------------------------------------------------


A1  <- A2  <- 1   # solo stessa ampiezza (per semplificare)
                  # il caso in cui si usa la formula della somma
nu <- 200      # Hz
d.phi <- 1/4*pi  # sfasamento: cambiare il valore per vedere l'effetto 

om <- 2*pi*nu
T  <- 1/nu

t.max <- 3*T
ylim <- c(-1,1) * (A1+A2)

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

curve(A1*cos(om*t), 0, t.max, xname='t', n=1001, col='blue', lwd=2,
      ylim=ylim, xlab=' t  (s)', ylab='f(t)',
      main=sprintf("nu = %.0f;   Delta phi = %.0f°", nu, d.phi*180/pi))
Sys.sleep(0.5)
curve(A2*cos(om*t + d.phi), xname='t', n=1001, col='green', lwd=2, add=TRUE)
Sys.sleep(0.5)
curve(A1*cos(om*t)+A2*cos(om*t + d.phi), xname='t', n=1001, col='red',
      lwd=2, add=TRUE)
Sys.sleep(0.5)
curve(2*A1*cos(om*t +d.phi/2)*cos(d.phi/2), xname='t', n=1001, col='yellow',
      lwd=2, lty=2, add=TRUE)
par(mar=old.mar)