nieshuyan1 2018-12-24 15:48 采纳率: 0%
浏览 806

C语言fft变换,滤波,逆变换结果不对


```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);
  • 写回答

1条回答 默认 最新

  • CSDN-Ada助手 CSDN-AI 官方账号 2022-09-20 19:31
    关注
    不知道你这个问题是否已经解决, 如果还没有解决的话:

    如果你已经解决了该问题, 非常希望你能够分享一下解决方案, 写成博客, 将相关链接放在评论区, 以帮助更多的人 ^-^
    评论

报告相同问题?

悬赏问题

  • ¥15 如何在scanpy上做差异基因和通路富集?
  • ¥20 关于#硬件工程#的问题,请各位专家解答!
  • ¥15 关于#matlab#的问题:期望的系统闭环传递函数为G(s)=wn^2/s^2+2¢wn+wn^2阻尼系数¢=0.707,使系统具有较小的超调量
  • ¥15 FLUENT如何实现在堆积颗粒的上表面加载高斯热源
  • ¥30 截图中的mathematics程序转换成matlab
  • ¥15 动力学代码报错,维度不匹配
  • ¥15 Power query添加列问题
  • ¥50 Kubernetes&Fission&Eleasticsearch
  • ¥15 報錯:Person is not mapped,如何解決?
  • ¥15 c++头文件不能识别CDialog