Die Datenpunkte $P(x|y)$ werden von Tigerjython mit einer Nachkommastelle zurückgegeben. Die berechneten Fourierkoeffizienten haben aber 13 Nachkommastellen. Damit wird unsere Datenreduktion um den Faktor 13 verschlechtert. Abhilfe schafft runden der Fourierkoeffizienten auf zwei oder drei Nachkommastellen.

round(float Zahl, Anzahl Nachkommastellen)

Da unsere Fourierkoeffizienten komplexe Zahlen sind müssen wir wie folgt vorgehen.

c=complex(round(c.real,2),round(c.imag,2))

Sei $n$ die Anzahl Punkte, wie viele Fourierkoeffizienten $c_k$ können, dürfen wir berechnen? Dazu müssen wir uns die Formel zur Berechnung der Fourierkoeffizienten genauer anschauen. $$c_k=\int_0^1 f(t)\cdot e^{-2 \pi ikt}dt \approx \sum_{j=0}^{n-1}f(j\cdot \Delta t)\cdot e^{-2 \pi ikj\cdot \Delta t}\cdot \Delta t$$ Relevant ist dabei die komplexe Exponentialfunktion, welche $2\pi$-periodisch ist. Wir berechnen die $c_k$ über eine endliche Summe. $$c_k=\sum_{j=0}^{n-1}f(j\cdot \Delta t)\cdot e^{-2 \pi ikj\cdot \Delta t}\cdot \Delta t$$ Die komplexe Exponentialfunktion als Funktion von $k$ lautet dann $$f_j(k)= e^{-2 \pi ikj\cdot \Delta t} \qquad j\in\{0,1,...,n-1\}$$ Zudem ist $\Delta t = \frac{1}{n-1}$. Setzen wir dies ein, so erkennen wir $$f_j(k)= e^{-2 \pi i k j\cdot \frac{1}{n-1}}$$ Sobald $k=n-1$ ist, wiederholt sich die Funktion. D.h. wir haben eine Periode von $n-1$. Es gilt: $$f_j(k)=f_j(k+\lambda\cdot (n-1)) $$ Wenn $k_{max}=n-1$ ist, dann würden wir (2n-1) Fourierkoeffizienten berechnen, da unsere $c_k$ von $-(n-1)\leq k \leq (n-1)$ durchlaufen. Dies sind zuviele Fourierkoeffzienten. Es gilt die Regel, dass aus $n$ Datenpunkten höchsten $n$ Fourierkoeffizienten berechnet werden können. Somit liegt der maximale Index $k_{max}$ bei $$k_{max} = \lfloor \frac{n-1}{2} \rfloor$$. Das sind i.d.R. immer noch zu viele Fourierkoeffizienten.

Aufgabe Ermittle für den Datesatz “KSBG.txt”