如題，請賜教～謝謝！

http://delphi.ktop.com.tw/topic.php?topic_id=29001 http://delphi.ktop.com.tw/topic.php?topic_id=41840    參考一些前輩們的文章 ---------------------- 台上一分鐘，台下十年功

------
台上一分鐘，台下十年功

Dear Ray24, 以前二篇主題我有看過了~~第一篇有看沒有懂，第二篇是寫好的程式，無法載入我的程式之中，我想要做的是一維陣列的數值傅利葉轉換，請問您還可以提供其它方法嗎？謝謝！ ben

■【問題】[問題]富利業轉換中的虛部運算  http://delphi.ktop.com.tw/topic.php?TOPIC_ID=54689 不知道您是否看過這篇? 有公式應該就可以算了。    吃軟也吃硬 dllee.ktop.com.tw StatPlus 系統資源監測器 @ SoftKing VMASK - ViewMove Automation Software Kernel

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

to dllee: 這篇文章我看了還是不會使用耶(數學不好) 如果我有一個陣列  int A[10]={0,0,0,0.2,0.5,1,0.5,0.2,0,0}; 該如何引用以下之前所post的文章程式呢？  請再指教~~感恩~~      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 的次方 } } }

引言:這篇文章我看了還是不會使用耶(數學不好)

我不清楚您是那裡不會... 那篇我主要是說明 complex 型態的使用。 如果您說您的數學不好... 那我建議您還是別玩 FFT，不然，您就算作出來，您也不知道那是什麼東西， 就算別人拿可用的 Source 給您，您也無法辨認是否正確。 如果您要找原始碼... 目前在 KTOP 上是沒有完整的 FFT 原始碼可以下載， 那篇只有片段，同時，您必需了解 FFT 的演算法才能了解程式， 我也附件了我的作業在 FFT 演算法上的解說，而這是大部分介紹 FFT 的書都會有的。 使用 google 去查 FFT C code, FFT C code http://www.google.com/search?q=FFT C code http://www.google.com/search?q=FFT C++ code 裡面就有 C/C 的 FFT 完整原始碼，同樣，如果您不懂數學， 您可能也不會知道那原始碼各項所需要的輸入參數是什麼、如何傳入、傳回的又是什麼... 如果您沒有教科書，只是老闆要你作一些 FFT 的東西，建議您還是先找本 書看一看，不然，用 google 查 FFT algorithm 或 FFT 演算法 也可以了解 FFT 是如何演算的。 以您想要的把 A[10] 用 FFT 作轉換，基本上，如果是老師出的題目， 就表示他除了要考你 FFT 演算法外，還要考你 FFT 所代表的意義， 您可知道 A[10] 是無法作 FFT 的，您必需把它變成 A[8] (去掉2個元素) 或變成 A[16] (補上6個元素) 成為 2^N 次方的數量才能使用 FFT。 (這只是 FFT 的限制，不是它所代表的意義) 吃軟也吃硬 dllee.ktop.com.tw StatPlus 系統資源監測器 @ SoftKing VMASK - ViewMove Automation Software Kernel

Dear dllee,   非常感謝您詳細的解說，基本上，我只是要使用FFT的函數庫來運用我的結果，我試了您給的網頁上幾個例題，結果RUN出來的答案都與matlab的fft結果不同。我以以下為例~~~ 底下有二個例子~~      嗯~~~我自己再去網頁上找找好了~~感恩~        double x[]={1.0,2.0,3,4,5,6,7,8},y[]={0,0,0,0,0,0,0,0};    int DFT(int dir,int m,double *x1,double *y1) {    long i,k;    double arg;    double cosarg,sinarg;    double *x2=NULL,*y2=NULL;       x2 =(double *) malloc(m*sizeof(double));    y2 =(double *) malloc(m*sizeof(double));    if (x2 == NULL || y2 == NULL)       return(FALSE);       for (i=0;i

如果您只有 8 點，用手算都可以算得出來  請參考: ■我的 href="http://dllee.ktop.com.tw" target="blank">吃軟也吃硬 dllee.ktop.com.tw StatPlus 系統資源監測器 @ SoftKing VMASK - ViewMove Automation Software Kernel

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

嗯~~了解 ~~~ 但是我的數值是數百個點，只是例子寫成八個。 在這裡其它篇有使用到BCB呼叫MATLAB函式庫， 我再試試，成功再POST上來與大家分享~~ 感恩~~

你可以使用SDL VCL for BCB中的VCL就可以作FFT

感謝！~~~~我來試試看~~~ 有問題再向各位請教 ~~~

感謝sppkenny提供這個資訊，我上網download它的元件庫。 http://www.lohninger.com/delfcomp.html http://www.lohninger.com/download_sdlsuite.html 還滿簡單使用的，只是做轉換後，結果有一些些許的不同。不曉得是那裡出問題。或者是定義不同。