opencv使用cvDFT 傅立叶转换的例子

 

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


    unsigned int i;
    CvMat* s = cvCreateMat(1,32,CV_32FC1);
    CvMat* d = cvCreateMat(1,32,CV_32FC1);
    CvMat* s2 = cvCreateMat(1,32,CV_32FC1);
    for(i=0;i<32;i++)
        cvSetReal2D(s,0,i,i);

    //显示一下输入信号
    for(i=0;i<32;i++)
        printf("%6.2f ",cvGetReal2D(s,0,i));
    printf(" ");

    //DFT 离散傅立叶变换
    cvDFT(s,d,CV_DXT_FORWARD);     //CV_DXT_FORWARD 代表了正变换：空域－〉频域

    printf("The result of DFT: ");
    for(i=0;i<32;i++)
        printf("%6.2f ",cvGetReal2D(d,0,i));

    //DFT 离散傅立叶逆变换
    cvDFT(d,s2,CV_DXT_INVERSE); //逆变换
    printf(" The result of IDFT: ");
    for(i=0;i<32;i++)
        printf("%6.2f ",cvGetReal2D(s2,0,i));
    printf(" ");

    system("pause");


    return 0;



}

结果：

 

   0.00      1.00      2.00      3.00      4.00      5.00      6.00      7.00      8.00      9.00
  10.00     11.00     12.00     13.00     14.00     15.00     16.00     17.00     18.00     19.00
  20.00     21.00     22.00     23.00     24.00     25.00     26.00     27.00     28.00     29.00
  30.00     31.00
The result of DFT:
496.00    - 16.00    162.45    - 16.00     80.44    - 16.00     52.74    - 16.00     38.63    - 16.00
  29.93    - 16.00     23.95    - 16.00     19.50    - 16.00     16.00    - 16.00     13.13    - 16.00
  10.69    - 16.00      8.55    - 16.00      6.63    - 16.00      4.85    - 16.00      3.18    - 16.00
   1.58    - 16.00
The result of IDFT:
   0.00     32.00     64.00     96.00    128.00    160.00    192.00    224.00    256.00    288.00
320.00    352.00    384.00    416.00    448.00    480.00    512.00    544.00    576.00    608.00
640.00    672.00    704.00    736.00    768.00    800.00    832.00    864.00    896.00    928.00
960.00    992.00
Press any key to  continue  . . .