c - 使用 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] 是半复数数组的第 i 个元素的实部，而 output_array[Ni] 是第 i 个值的虚部。

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

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

那是对的吗？那么二次谐波的幅度将是相同的，但使用 output_array[2] 等等，对吗？

那么我怎样才能得到一、二、三、四次谐波的相位呢？

谢谢你。