double FD_R[255][255]={0},FD_I[255][255]={0};
  for(int u=0;u<255;u  )
   for(int v=0;v<255;v  )
    for(int x=0;x<255;i  )
     for(int y=0;y<255;j  )
     {
         shift_weight=pow(-1,(x y));
         FD_R[u]  =>> >-=>>>
做 >-=>> > =>>
        

都使用 C++ 來寫程式了，就應該使用 C++ 的物件來協助我們處理。 以下我以前寫的 FFT 片段

...
   complex Fodd_Wn,Wn,wn;
   complex **ppcpxFFTData;
...
   for(j=0; j( 1.0 , 0.0);
         wn = complex( wr, -wi);
         for(m=1; m <= stepN2 ; m  )
         {
            for(i=m; i <= N ; i =stepN)
            {
               p = i-1;
               q = p   stepN2;
               Fodd_Wn = ppcpxFFTData[q][j] * Wn;     // Fodd * Wn
               ppcpxFFTData[q][j] = ppcpxFFTData[p][j] - Fodd_Wn;
               ppcpxFFTData[p][j] = ppcpxFFTData[p][j]   Fodd_Wn;
            }
            Wn = Wn * wn;              // 計算 Wn 的次方
         }
      }
   }
...

使用 C 的物件，在使用上，會與書上的數學式相近，也不用去考慮虛實運算的問題，寫作上會比較簡單。 視動科技 VMASK - ViewMove Automation Software Kernel 發表人 - dllee 於 2004/08/13 16:45:24

------
http://www.ViewMove.com

 for(m=1; m <= stepN2 ; m  )
         {
            for(i=m; i <= N ; i =stepN)
            {
               p = i-1;
               q = p   stepN2;

不好意思！想問dllee大大就是關於以上迴圈跟程式碼的執行運算不是很懂，不曉得可以詳細解說？謝謝你的回覆

[quote] 都使用 C++ 來寫程式了，就應該使用 C++ 的物件來協助我們處理。 以下我以前寫的 FFT 片段

...
   complex Fodd_Wn,Wn,wn;
   complex **ppcpxFFTData;
...
   for(j=0; j( 1.0 , 0.0);
         wn = complex( wr, -wi);
         for(m=1; m <= stepN2 ; m  )
         {
            for(i=m; i <= N ; i =stepN)
            {
               p = i-1;
               q = p   stepN2;
               Fodd_Wn = ppcpxFFTData[q][j] * Wn;     // Fodd * Wn
               ppcpxFFTData[q][j] = ppcpxFFTData[p][j] - Fodd_Wn;
               ppcpxFFTData[p][j] = ppcpxFFTData[p][j]   Fodd_Wn;
            }
            Wn = Wn * wn;              // 計算 Wn 的次方
         }
      }
   }
...

使用 C 的物件，在使用上，會與書上的數學式相近，也不用去考慮虛實運算的問題，寫作上會比較簡單。 ==================================================================== 因為沒用過complex物件 希望D大能就我上述所舉之DFT例子 該怎麼應用它呢?

to roger-c: 我的程式是由書本來的，請參考我的作業節錄: http://delphi.ktop.com.tw/topic.php?TOPIC_ID=54769 對於 FFT 算是完整的整理了。    to plihui: 如果您了解數學，就應該了解 complex 的計算及使用， 在程式中使用 

complex A,B;

只是表示宣告 A,B 兩個 complex 就像

int a,b;

是一樣的，接著您可以給它數值，如:

A = complex( 1.0 , 0.0);
B = complex( 0.0 , 1.0);

或是作 - * / 的四則運算，如:

A = A*B   B;

與用整數 int 或是其他的變數型態是沒兩樣的。 視動科技 VMASK - ViewMove Automation Software Kernel 發表人 - dllee 於 2004/08/14 11:01:37

------
http://www.ViewMove.com

謝謝dllee大大的提供我再好好的研究囉

to dllee 請教一下,我想參考你所提到的那本書 可否告知是那一本?

那書是老師自己整理的講義， 能不買書，我就不買了  < >視動科技 < >>

------
http://www.ViewMove.com

各位高手：不能提供完整的FFT程式碼嗎?