Tengo que hacer una pregunta complementaria a este post:
https://stackoverflow.com/questions/1449290/algorithmically-get-amplitude-and-phase-of-sine-wave
Mi objetivo es encontrar la amplitud y la fase de alguna señal. Para una mejor comprensión ahora quiero asumir que esta señal es una simple onda sinusoidal.
El código que tengo actualmente es:
t <- seq(0, 100, length.out = 1000)
s <- sin(t)
qplot(t, s, geom = "line")
f <- fft(s)
a2 <- Mod(f/length(s))
qplot(t, a2, geom = "line")Por lo que he entendido, el segundo gráfico visualiza ahora la transformada de Fourier de doble cara, ¿no? Entonces, ¿qué tengo que hacer para terminar con amplitude = 1 ? Al ver el segundo gráfico, ¿tendría que sumar los dos picos?
Y, con respecto a la fase. Sé que tendría que usar algo como p2 = Arg(f) Pero, ¿cómo termino con phase = 0 ?
¡Muchas gracias!
