1

首先我为我糟糕的英语道歉......

所以这是我的问题。我正在用一个简单的输入信号测试 FFTW3 库，一个连续的。然后我计算 FFT 并得到好的结果：只是频率为 0 的信号，其他一切都为 0。

然后我想用反向 FFT 取回我的输入，但它不起作用。这是我的代码：

fftw_complex* imgIn;
fftw_complex* imgIn2;
fftw_complex* imgOut;
fftw_plan plan;

int taille = 100;
int i;

 //Allocation des entrées et sorties
imgIn = fftw_malloc(sizeof(fftw_complex)*taille);
imgIn2 = fftw_malloc(sizeof(fftw_complex)*taille);
imgOut = fftw_malloc(sizeof(fftw_complex)*taille);

//Remplissage des données d'entrées pour le calcul de la FFT
for(i = 0 ; i < taille ; i++){
    imgIn[i][0] = 1.0;
    imgIn[i][1] = 0.0;
}

//Plan d'execution
plan = fftw_plan_dft_2d(taille/10, taille/10, imgIn, imgOut, FFTW_FORWARD, FFTW_ESTIMATE);

//Execute la FFT
fftw_execute(plan);

//Inverse
plan = fftw_plan_dft_2d(taille/10, taille/10, imgOut, imgIn2, FFTW_BACKWARD, FFTW_ESTIMATE);

for(i = 0 ; i < taille ; i++){
    printf("%d : %g\n%d : %g\n", i, imgIn2[i][0], i, imgIn2[i][1]);
}

如您所见，我只是尝试执行正常的 FFT，然后对其进行反转。问题是我的输出 imgIn2 只是充满了 0，而不是 1 和 0...

那么我的代码有什么问题？

谢谢 ：）

4

1 回答 1

2

您的代码不执行第二个计划。

于 2013-02-26T00:14:30.890 回答