1:
2:
3:
4:
5:
6:
7:
8:
9:
10:
11:
12:
13:
14:
15:
16:
17:
18:
19:
20:
21:
22:
23:
24:
25:
26:
27:
28:
29:
30:
31:
32:
33:
34:
35:
36:
37:
38:
39:
40:
41:
42:
43:
44:
45:
46:
47:
48:
49:
50:
51:
52:
53:
54:
55:
56:
57:
58:
59:
60:
61:
62:
63:
64:

function fftwf_plan_dft_r2c_1d(n : Integer; inData : PSingle; outData : PSingle; flags : Longword) : Pointer; cdecl; external 'fftw3.dll';
function fftwf_plan_dft_c2r_1d(n : Integer; inData : PSingle; outData : PSingle; flags : Longword) : Pointer; cdecl; external 'fftw3.dll';
procedure fftwf_destroy_plan(plan : Pointer); cdecl; external 'fftw3.dll';
procedure fftwf_execute(plan : Pointer); cdecl; external 'fftw3.dll';


procedure TForm1.Button1Click(Sender: TObject);
Type
 complex = record
  real, imag: single;
 end;
Var
 FFT,FFT2, IFFT: Pointer;
 x,y: array of single;
 w,v,u: array of complex;
 N, I: Integer;
begin
 N := 8;
 SetLength(x, N);
 SetLength(y, N);
 SetLength(w, N div 2+1); //diese Länge des Ergebnisses ist im FFTW Manual so zu finden
 SetLength(v, N div 2+1);
 SetLength(u, N div 2+1);

 FFT := fftwf_plan_dft_r2c_1d(8, @x[0], @w[0],64); //Aufrufkonvention des FFTW Pakets....
 FFT2 := fftwf_plan_dft_r2c_1d(8, @y[0], @v[0],64);
 IFFT := fftwf_plan_dft_c2r_1d(8, @u[0], @x[0],64);

 x[0] := 0; //das sind die beiden Array mit den reellen Polynomen X[0] ist 10^0, X[1] 10^1 usw...
 x[1] := 0;
 x[2] := 0;
 x[3] := 0;
 x[4] := 0;
 x[5] := 0;
 x[6] := 0;
 x[7] := 0;

 y[0] := 0;
 y[1] := 0;
 y[2] := 0;
 y[3] := 0;
 y[4] := 0;
 y[5] := 0;
 y[6] := 0;
 y[7] := 0;


 fftwf_execute(FFT); //auf beide arrays den FFT anwenden
 fftwf_execute(FFT2);

 for i:=0 to (n div 2) do begin
   u[i].real:=v[i].real*w[i].real-v[i].imag*w[i].imag; //Ergebnis im Komplexen stellenweise multiplizieren
   u[i].imag:=v[i].imag*w[i].real-v[i].real*w[i].imag;
 end;

 fftwf_execute(IFFT);

 for I := 0 to N-1 do x[I] := x[I] / N; //Ergebnis durch Anzahl der Stellen teilen
 for I := 0 to N-1 do Memo2.Lines.Add(FloatToStr(x[I]));  //und mir ausgeben

 fftwf_destroy_plan(FFT);
 fftwf_destroy_plan(FFT2);
 fftwf_destroy_plan(IFFT);
end;