使用fftw c库获得前4次谐波的幅度和相位

Question

我有360个值的周期性系列。如何使用fftw c库来获得一次，二次，三次和四次谐波的幅度和相位？

如果我这样做N = 360

fftw_plan_r2r_1d(N,input_array,output_array,FFTW_R2HC,FFTW_ESTIMATE);

我可以得到离散傅立叶变换，其中output_array [i]是halfcomplex数组的第i个元素的实部，而output_array [N-i]是第i个值的虚部。

所以，如果我想要一次谐波的幅度或幅度（??），我应该这样做吗？

ampl_1sth=sqrt(output_array[1]*output_array[1]+output_array[N-1]*output_array[N-1])

这是对的吗？那么二次谐波的幅度是相同的，但是使用output_array [2]，依此类推吧？

然后我怎样才能得到第一，第二，第三和第四谐波的相位？

谢谢。