```C语言fft变换,截取了幅值最高点附近5hz的幅频,其他的地方置0,逆变换后,得到的曲线幅值从低到高,到达正常幅值后又从高到低,正常应该是幅值相等的曲线,请大神帮助。
fft(buff_f,FFT_N);//为4096个点的数据
//计算幅值,频率
float Amp[FFT_N];//幅值
float fre[FFT_N];//频率
float FRE= 1000.0*(float)(s-1)/(time[s-1]- time[0]);//采样频率
for (int i = 0; i < FFT_N; i++)
{
Amp[i] = 2*sqrt(buff_f[i].real*buff_f[i].real+buff_f[i].imag*buff_f[i].imag)/FFT_N;
fre[i]=i*FRE/(float)FFT_N;//%频率序列
}
//找到幅值最大的点,就是基频,由于fft是对称的,只取一半
int max_i=0;//存储最大值时的序列点
for(i=1;i<FFT_N/2;i++)
{
if(Amp[max_i]>=Amp[i-1])
{
basefre_amp[0]=Amp[max_i];//幅值
basefre_amp[1]=fre[max_i];//频率
}
if(Amp[max_i]<Amp[i-1])
{
basefre_amp[0]=Amp[i-1];//幅值
basefre_amp[1]=fre[i-1];//频率
max_i=i-1;
}
}
float dfre=2;
float c=FRE/(float)FFT_N;
for(i=0;i<FFT_N/2;i++)//滤波,将附近的几个点保留,其他滤掉
{
if(max_i*c-dfre>i*c||max_i*c+dfre<i*c)
{
buff_f[i].imag=0;
buff_f[i].real=0;
}
}
for(i=FFT_N/2;i<FFT_N;i++)//滤波,将附近的几个点保留,其他滤掉
{
buff_f[i].imag=buff_f[FFT_N-i-1].imag;
buff_f[i].real=buff_f[FFT_N-i-1].real;
}
//逆变换
IFFT(buff_f,FFT_N);