Амплитуда и фаза FFT

Question

У меня есть функция звука, которую я пытаюсь разбить на синусы/косинусы. Поэтому я прибегал к быстрой трансформации Фурье. Используя функцию fft(y,inverse=FALSE), я смог преобразовать временную область звука в частотный doman. Результат сложный. Я прочитал, что для преобразования этого вывода, который находится в мнимой форме, и отсеивает необходимую информацию, которая является амплитудой и фазой A(v)cos(2*pi*v+P), для получения амплитуды необходимо использовать выходной сигнал abs(); однако мне трудно найти функцию R, которая получает фазу. В MATLAB функция angle() возвращает фазы БПФ. Какова соответствующая функция в R, чтобы найти фазу?

Update

Спасибо за предложения, ребята; все еще испытывают проблему. Я запускаю FFT на простой функции, чтобы проверить, работает ли она. Моя функция: y=cos(2*pi*(seq(0,10,by=.01)*(1/5)+7.5). Таким образом, частота составляет 1/5 с фазовым сдвигом 7,5.

y=cos(2*pi*(seq(0,10,by=.01)*(1/5)+7.5) 
fty=(y,inverse=F) 
plot(abs(fty),xlim=c(0,10),type="l") 
angle=atan2(Im(fty), Re(fty)) 
> angle[3] 
[1] 1.222766 

Когда я сюжет серии, амплитуда достигает максимума при значении частоты 3 и функции угла (который должен дать мне мою фазу на частоте, при которой амплитуда которой пики при) дает мне фазу 1.2 , Что я делаю не так?

Answer 1

Вы можете обнаружить, что функция atan2 делает то, что вы хотите. Вы дадите ему воображаемые и реальные части ценности, так как прототип - atan2(y, x). Таким образом, вы могли бы сделать:

angle = atan2(Im(value), Re(value));