f.cpp

来自「matlab可以方便地实现付利叶变换」· C++ 代码 · 共 40 行

#include <complex>
#include <fftw3.h>
#include <math.h>
#include <iostream>




#define N 10
using namespace std;


//extern plan_dft_1d()；

void main(int argc, char * argv[]){

    fftw_complex in[N], out[N];
    fftw_plan p;

//一维dft，in 输入，out输出，FFTW_FORWARD 表示 exp 上指数是负号
// out = F in
//where out and in are two vectors of the same length n, and F is a n-by-n matrix with the (j,k) element
// F jk = exp(-i 2 pi j k /n). i 是虚数单位. 

    p= plan_dft_1d(N,in,out,FFTW_FORWARD,FFTW_MEASURE);
    for(int i=0;i <N;i ++) {
        in[i][0]=i;
        in[i][1]=0.0;
    }
   

 /*   fftw_execute(p);

    for(int i=0;i <N;i ++){
        cout<<out[i][0]<<" "<<out[i][1]<<endl;
    }*/

  }