das Wochenende stand mal ganz im Zeichen der Mathematik und von Python. So wollte ich versuchen, die schnelle Fouriertransformation mal zu programmieren. Der "langsame" Teil ging auch ganz gut. Ein Programm kann schonmal diskret fourier transformieren.

Leider komme ich mit der wirklichen schnellen Variante nicht hin. Diese ist ja rekursiv definiert und sollte durch die Halbierung des Intervalls entsprechend rasant sein. (Meine "langsame" Variante braucht für 1024 Stützwerte immerhin nahezu eine halbe Minute). Meine bisherigen Versuche, die Rekursion hinzubekommen, waren weniger erfolgreich. Grundsätzlich versuche ich mir eine Funktion fft zu definieren, die einen Vektor und die Anzahl der Runden übergeben bekommt. Sollte die Anzahl eins sein, wird kalkuliert und ansonsten wird der Vektor immer wieder halbiert und mit den halbierten Vektoren die Funktion wieder aufgerufen. (Ich habe leider gerade keine Möglichkeit den originalen Code hierhin zu kopieren. Der folgt später.)
Da das nun nicht wie erwartet funzt und das Wochenende sich dem Ende neigt, wollte ich fragen, ob ihr evtl. Pythoncode kennt, in dem ich mal "spicken" kann? Beim Gurgeln war ich nicht sehr erfolgreich (bis auf NumPy).