幽人拾遗 2022-08-24 17:45 采纳率: 0%
浏览 75
已结题

fftw库与matlab进行傅里叶变换结果不一致

我用fftw库对一维数组double array[ ] = { 0, 1, 2, 3, 4, 5, 6, 7, 8, 9 }进行FFT变换,但是得到的结果与matlab的结果不一样,两种fft结果前六个数据是一样的,后三个数据就很奇怪。原代码如下:

#include "fftw3.h"  
#include<stdio.h>
#include<iostream>
#include<vector>

using namespace std;
int main()
{
    //****************************ifft********************************
    double array[] = { 0, 1, 2, 3, 4, 5, 6, 7, 8, 9 };
    double* out;
    double* err;
    int i, size = 10;

    fftw_complex* out_cpx;

    fftw_plan fft;
    fftw_plan ifft;
    out_cpx = (fftw_complex*)fftw_malloc(sizeof(fftw_complex) * size);
    out = (double*)malloc(size * sizeof(double));
    err = (double*)malloc(size * sizeof(double));

    fft = fftw_plan_dft_r2c_1d(size, array, out_cpx, FFTW_ESTIMATE);  //Setup fftw plan for fft
    fftw_execute(fft);

    cout << "fft输出结果" << endl;
    for (int k = 0; k < size; k++)
    {
        cout << "(" << out_cpx[k][0] << "," << out_cpx[k][1] << ") ";
    }
    cout << endl;

    ifft = fftw_plan_dft_c2r_1d(size, out_cpx, out, FFTW_ESTIMATE);   //Setup fftw plan for ifft


    fftw_execute(ifft);

    for (i = 0; i < size; i++)
    {
        err[i] = (array[i] - out[i]);
        printf("%f\t%f\n", (array[i]), out[i] / size);//需要做归一化处理
    }



    fftw_destroy_plan(fft);
    fftw_destroy_plan(ifft);
    fftw_free(out_cpx);
    free(err);
    free(out);
    //***************************************ifft*********************

    system("pause");//暂停
    return 0;
}


c++运行结果

img

matlab运行结果

img


matlab代码:

img

  • 写回答

1条回答 默认 最新

  • 幽人拾遗 2022-08-27 21:46
    关注

    已解决

    本回答被题主选为最佳回答 , 对您是否有帮助呢?
    评论

报告相同问题?

问题事件

  • 系统已结题 9月4日
  • 已采纳回答 8月27日
  • 创建了问题 8月24日

悬赏问题

  • ¥15 nslt的可用模型,或者其他可以进行推理的现有模型
  • ¥15 arduino上连sim900a实现连接mqtt服务器
  • ¥15 vncviewer7.0安装后如何正确注册License许可证,激活使用
  • ¥15 phython如何实现以下功能?查找同一用户名的消费金额合并2
  • ¥66 关于人体营养与饮食规划的线性规划模型
  • ¥15 基于深度学习的快递面单识别系统
  • ¥15 Multisim仿真设计地铁到站提醒电路
  • ¥15 怎么用一个500W电源给5台60W的电脑供电
  • ¥15 请推荐一个轻量级规则引擎,配合流程引擎使用,规则引擎负责判断出符合规则的流程引擎模板id
  • ¥15 Excel表只有年月怎么计算年龄