一般FFT分析,都是求一个时间历程的频率成分,
用FFTW中的
fftw_plan_dft_r2c_1d
输入是一个实数列向量in[N],
输出是一个复数空间二维矩阵out[N][2],第1列是实数、第2列是复数。
求幅值时,要对复数向量的每个,取向量模,乘以2,再除以点数N。
注意:
1. 如果采样时间不到 1 s, 就设法补成 1 s。
一种方法是补零,(参考 彭真明 老师博文 link2 )
还有一种方法,如果已知频率大概范围,折算出周期后,用正周期的时间历程估计。
2. 如果有直流项(平均值大于零),需要把直流项滤掉。
fftw_complex *out = NULL; |
fftw_plan p; |
in = (double *)fftw_malloc(sizeof(double) * len); |
out = (fftw_complex *)fftw_malloc(sizeof(fftw_complex) * len); |
p = fftw_plan_dft_r2c_1d(len, in, out, FFTW_ESTIMATE); |
fftw_execute(p); |
for (int i = 0; i < len; i++) |
{ |
amp = sqrt(out[i][0]*out[i][0] + out[i][1]*out[i][1]); |
amp = 2.0* amp/len; |
freq = i*1.0/t1; |
} |