Aufgaben:Aufgabe 5.5Z: Rechenaufwand für die FFT: Unterschied zwischen den Versionen

Aus LNTwww
Wechseln zu:Navigation, Suche
 
(11 dazwischenliegende Versionen von 2 Benutzern werden nicht angezeigt)
Zeile 3: Zeile 3:
 
}}
 
}}
  
[[Datei:P_ID1179__Sig_Z_5_5.png|right|Zum Rechenaufwand der FFT]]
+
[[Datei:P_ID1179__Sig_Z_5_5.png|right|frame|Letzte Stufe der FFT für  $N=8$]]
Der FFT–Algorithmus realisiert eine Diskrete Fouriertransformation mit dem kleinstmöglichen Rechenaufwand, wenn der Parameter $N$ eine Zweierpotenz ist. Im Einzelnen sind zur Durchführung einer FFT folgende Rechenschritte notwendig:
+
Der FFT–Algorithmus  (Fast Fourier Transform)  realisiert eine  „Diskrete Fouriertransformation”  mit dem kleinstmöglichen Rechenaufwand, wenn der Parameter  $N$  eine Zweierpotenz ist.  
*Die FFT geschieht in ${\rm ld} \ N$ Stufen, wobei in jeder Stufe die genau gleiche Anzahl an Rechenoperationen durchzuführen ist. „ld” steht hier als Abkürzung für den Logarithmus zur Basis 2.
+
 
*Die Grafik zeigt die dritte und letzte Stufe für das Beispiel $N = 8$. Man erkennt, dass in dieser und auch den anderen Stufen jeweils $N/2$ Basisoperationen durchzuführen sind.
+
Im Einzelnen sind zur Durchführung einer FFT folgende Rechenschritte notwendig:
*In jeder Basisoperation, die man häufig auch als '''Butterfly''' bezeichnet, werden aus den beiden komplexen Eingangsgrößen $E_1$ und $E_2$ zwei komplexe Ausgänge berechnet:
+
*Die FFT geschieht in  ${\rm log_2} \ N$  Stufen, wobei in jeder Stufe die genau gleiche Anzahl an Rechenoperationen durchzuführen ist.  
 +
*Die Grafik zeigt die dritte und letzte Stufe für das Beispiel  $N = 8$.  Man erkennt, dass in dieser und auch den anderen Stufen jeweils  $N/2$  Basisoperationen durchzuführen sind.
 +
*In jeder Basisoperation, die man häufig auch als  '''Butterfly'''  bezeichnet, werden aus den beiden komplexen Eingangsgrößen  $E_1$  und  $E_2$  zwei komplexe Ausgänge berechnet:
 
:$$ A_1  = E_1 + E_2 \cdot w^{\hspace{0.04cm}\mu}, $$
 
:$$ A_1  = E_1 + E_2 \cdot w^{\hspace{0.04cm}\mu}, $$
 
:$$ A_2 = E_1 - E_2 \cdot w^{\hspace{0.04cm} \mu}\hspace{0.05cm}.$$
 
:$$ A_2 = E_1 - E_2 \cdot w^{\hspace{0.04cm} \mu}\hspace{0.05cm}.$$
*Hierbei bezeichnet $w =  {\rm e}^{-{\rm j} 2\pi/N}$ den komplexen Drehfaktor. Für das Beispiel $N = 8$ hat dieser folgenden Wert:   $w =  {\rm e}^{-{\rm j} \hspace{0.05cm}\cdot \hspace{0.05cm} \pi/4} = \cos(45^\circ) - {\rm j} \cdot \sin(45^\circ)\hspace{0.05cm}.$
+
*Hierbei bezeichnet  $w =  {\rm e}^{-{\rm j}\hspace{0.05cm}\cdot \hspace{0.05cm} 2\pi/N}$  den komplexen Drehfaktor.  Für   $N = 8$  ergibt sich der Wert  $w =  {\rm e}^{-{\rm j} \hspace{0.05cm}\cdot \hspace{0.05cm} \pi/4} = \cos(45^\circ) - {\rm j} \cdot \sin(45^\circ)\hspace{0.05cm}.$
*Der Exponent $\mu$ dieses komplexen Drehfaktors kann alle ganzzahligen Werte zwischen $0$ und $N/2-1$ annehmen. Für $N = 8$ gilt:
+
*Der Exponent  $\mu$  für diesen komplexen Drehfaktor kann alle ganzzahligen Werte zwischen  $0$  und  $N/2-1$  annehmen.  Für  $N = 8$  gilt:
 
:$$w^0 = 1,\hspace{0.2cm}w^1 = {1}/{\sqrt{2}}- {\rm j}
 
:$$w^0 = 1,\hspace{0.2cm}w^1 = {1}/{\sqrt{2}}- {\rm j}
 
\cdot{1}/{\sqrt{2}},\hspace{0.2cm}w^2 = - {\rm
 
\cdot{1}/{\sqrt{2}},\hspace{0.2cm}w^2 = - {\rm
Zeile 17: Zeile 19:
 
\cdot{1}/{\sqrt{2}} \hspace{0.05cm}.$$
 
\cdot{1}/{\sqrt{2}} \hspace{0.05cm}.$$
  
Mit dieser Aufgabe sollen die für die FFT erforderliche Anzahl von Rechenoperationen $(O_{\rm FFT})$ ermittelt und mit dem für die DFT angebbaren Wert $O_{\rm DFT} ≈ 8\cdot N^2$ verglichen werden.
+
Mit dieser Aufgabe sollen die für die FFT erforderliche Anzahl von Rechenoperationen  $(\mathcal{O}_{\rm FFT})$  ermittelt und mit dem für die DFT angebbaren Wert  $\mathcal{O}_{\rm DFT} ≈ 8\cdot N^2$  verglichen werden.
  
 
Zu beachten ist:
 
Zu beachten ist:
 +
*Sinnvollerweise werden die Potenzen von  $w$  vor dem eigentlichen Algorithmus berechnet und in einer Lookup–Tabelle abgelegt.  Die hierfür notwendigen Operationen sollen hier unberücksichtigt bleiben.
 +
*Die Bitumkehroperation – eine Umsortierung, die vor der ersten Stufe durchzuführen ist – soll bei dieser Abschätzung ebenfalls nicht berücksichtigt werden.
 +
 +
 +
 +
 +
 +
 +
 +
''Hinweis:''
 +
*Die Aufgabe gehört zum  Kapitel  [[Signaldarstellung/Fast-Fouriertransformation_(FFT)|Fast-Fouriertransformation (FFT)]].
 +
  
*Sinnvollerweise werden die Potenzen von $w$ vor dem eigentlichen Algorithmus berechnet und in einer Lookup–Tabelle abgelegt. Die hierfür notwendigen Operationen sollen deshalb unberücksichtigt bleiben.
 
*Die Bitumkehroperation – eine Umsortierung, die vor der ersten Stufe durchzuführen ist – soll bei dieser Abschätzung ebenfalls nicht berücksichtigt werden.
 
  
''Hinweise:''
 
*Die Aufgabe gehört zum  Kapitel [[Signaldarstellung/Fast-Fouriertransformation_(FFT)|Fast-Fouriertransformation (FFT)]].
 
*Sollte die Eingabe des Zahlenwertes „0” erforderlich sein, so geben Sie bitte „0.” ein.
 
  
  
Zeile 32: Zeile 41:
  
 
<quiz display=simple>
 
<quiz display=simple>
{Wieviele reelle Additionen (<i>A</i><sub>A</sub>) erfordert eine komplexe Addition?
+
{Wieviele reelle Additionen&nbsp; $(A_{\rm A})$&nbsp; erfordert eine komplexe Addition?
 
|type="{}"}
 
|type="{}"}
$A_A$ = { 2 3% }
+
$A_{\rm A} \hspace{0.3cm} = \ $ { 2 }
  
{Wieviele reelle Additionen (<i>A</i><sub>M</sub>) und Multiplikationen (<i>M</i><sub>M</sub>) sind für eine komplexe Multiplikation erforderlich?
+
{Wieviele reelle Additionen&nbsp; $(A_{\rm M})$&nbsp; und Multiplikationen&nbsp; $(M_{\rm M})$&nbsp; sind für eine komplexe Multiplikation erforderlich?
 
|type="{}"}
 
|type="{}"}
$A_M$ = { 2 3% }
+
$A_{\rm M} \hspace{0.3cm} = \  $ { 2 }
 
+
$M_{\rm M} \hspace{0.2cm} = \  $ { 4 }
$M_M$ = { 4 3% }
 
  
{Wieviele komplexe Additionen/Subtraktionen (<i>a</i><sub>B</sub>) erfordert eine einzige Basisoperation (&bdquo;Butterfly&rdquo;)?   
+
{Wieviele komplexe Additionen/Subtraktionen&nbsp; $(a_{\rm B})$&nbsp; erfordert eine einzige Basisoperation (&bdquo;Butterfly&rdquo;)?   
Wieviele komplexe Multiplikationen (<i>m</i><sub>B</sub>) sind pro Basisoperation notwendig?
+
<br>Wieviele komplexe Multiplikationen&nbsp; $(m_{\rm B})$&nbsp; sind pro Basisoperation notwendig?
 
|type="{}"}
 
|type="{}"}
$a_B$ = { 2 3% }
+
$a_{\rm B} \hspace{0.32cm} = \ $  { 2 }
 +
$m_{\rm B} \hspace{0.2cm} = \ $  { 1 }
  
$m_B$ = { 1 3% }
+
{Wieviele Rechenoperationen&nbsp; (Additionen, Subtraktionen, Multiplikationen gleichermaßen)&nbsp; erfordert eine Basisoperation?
 
 
{Wieviele Rechenoperationen (Additionen, Subtraktionen, Multiplikationen gleichermaßen) erfordert eine Basisoperation?
 
 
|type="{}"}
 
|type="{}"}
$O_B$ = { 10 3% }
+
$\mathcal{O}_{\rm B} \ = \ $ { 10 }
  
{Wieviele reelle Operationen erfordert der gesamte FFT&ndash;Algorithmus? Welche Zahlenwerte ergeben sich für <i>N</i> = 16 und <i>N</i> = 1024?
+
{Wieviele reelle Operationen&nbsp; $(\mathcal{O}_{\rm FFT})$&nbsp; erfordert der gesamte FFT&ndash;Algorithmus?&nbsp; Welche Werte ergeben sich für&nbsp; $N = 16$&nbsp; und&nbsp; $N = 1024$?
 
|type="{}"}
 
|type="{}"}
$O_{FFT} (N = 16 )$ = { 320 3% }
+
$N = 16\text{:} \hspace{0.65cm} \mathcal{O}_{\rm FFT} \ = \ $ { 320 3% }
 
+
$N = 1024\text{:} \hspace{0.2cm} \mathcal{O}_{\rm FFT} \ = \ $ { 51200 3% }
$O_{FFT} (N = 1024 )$ = { 51200 3% }
 
  
{Wie groß ist der Zeitgewinn <i>G</i><sub>FFT</sub> = <i>O</i><sub>DFT</sub>/<i>O</i><sub>FFT</sub> für <i>N</i> = 16 und <i>N</i> = 1024?
+
{Wie groß ist der Zeitgewinn&nbsp; $G_{\rm FFT} = \mathcal{O}_{\rm DFT} - \mathcal{O}_{\rm FFT}$&nbsp; für&nbsp; $N = 16$&nbsp; und&nbsp; $N = 1024$?
 
|type="{}"}
 
|type="{}"}
$G_{FFT} (N = 16 )$ = { 6.4 3% }
+
$N = 16\text{:} \hspace{0.65cm} G_{\rm FFT} \ = \ $ { 6.4 3% }
 
+
$N = 1024\text{:} \hspace{0.2cm} G_{\rm FFT} \ = \ $ { 164 3% }
$G_{FFT} (N = 16 )$ = { 164 3% }
 
  
 
</quiz>
 
</quiz>
Zeile 69: Zeile 74:
 
===Musterlösung===
 
===Musterlösung===
 
{{ML-Kopf}}
 
{{ML-Kopf}}
<b>a)</b>&nbsp;&nbsp;Jede komplexe Addition erfordert zwei reelle Additionen:
+
'''(1)'''&nbsp; Jede komplexe Addition erfordert zwei reelle Additionen:
$$(R_1 + {\rm j} \cdot I_1) + (R_2 + {\rm j} \cdot I_2) = (R_1 +
+
:$$(R_1 + {\rm j} \cdot I_1) + (R_2 + {\rm j} \cdot I_2) = (R_1 +
 
R_2) + {\rm j} \cdot (I_1 + I_2)\hspace{0.3cm}\Rightarrow \hspace{0.3cm}
 
R_2) + {\rm j} \cdot (I_1 + I_2)\hspace{0.3cm}\Rightarrow \hspace{0.3cm}
 
\hspace{0.15 cm}\underline{ A_{\rm A} =
 
\hspace{0.15 cm}\underline{ A_{\rm A} =
 
2}\hspace{0.05cm}.$$
 
2}\hspace{0.05cm}.$$
<b>b)</b>&nbsp;&nbsp;Eine jede komplexe Multiplikation benötigt vier reelle Multiplikationen und zwei reelle Additionen:
+
 
$$(R_1 + {\rm j} \cdot I_1)  (R_2 + {\rm j} \cdot I_2) = (R_1 \cdot
+
 
 +
'''(2)'''&nbsp; Eine jede komplexe Multiplikation benötigt vier reelle Multiplikationen und zwei reelle Additionen:
 +
:$$(R_1 + {\rm j} \cdot I_1)  (R_2 + {\rm j} \cdot I_2) = (R_1 \cdot
 
R_2 - I_1 \cdot I_2) + {\rm j} \cdot (R_1 \cdot I_2 + R_2 \cdot
 
R_2 - I_1 \cdot I_2) + {\rm j} \cdot (R_1 \cdot I_2 + R_2 \cdot
I_1)$$
+
I_1)\hspace{0.3cm}
$$\Rightarrow \hspace{0.3cm}\hspace{0.15 cm}\underline{A_{\rm M} = 2,\hspace{0.3cm}M_{\rm M} =
+
\Rightarrow \hspace{0.3cm}\hspace{0.15 cm}\underline{A_{\rm M} = 2,\hspace{0.3cm}M_{\rm M} =
 
4} \hspace{0.05cm}.$$
 
4} \hspace{0.05cm}.$$
<b>c)</b>&nbsp;&nbsp;Die Basisoperationen lauten mit den komplexen Eingangsgrößen <i>E</i><sub>1</sub>, <i>E</i><sub>2</sub> und <i>w</i><sup><i>&mu;</i></sup>:
+
 
$$ A_1  = E_1 + E_2 \cdot w^{\hspace{0.04cm}\mu}, $$
+
 
$$ A_2 = E_1 - E_2 \cdot w^{\hspace{0.04cm} \mu}\hspace{0.05cm}.$$
+
'''(3)'''&nbsp; Die Basisoperationen lauten mit den komplexen Eingangsgrößen&nbsp; $E_1$,&nbsp; $E_2$&nbsp; und&nbsp; $w^{\hspace{0.04cm}\mu}$:
Dies bedeutet eine komplexe Multiplikation und zwei komplexe Additionen:
+
:$$ A_1  = E_1 + E_2 \cdot w^{\hspace{0.04cm}\mu},\hspace{0.5cm}  A_2 = E_1 - E_2 \cdot w^{\hspace{0.04cm} \mu}\hspace{0.05cm}.$$
$$\hspace{0.15 cm}\underline{a_{\rm B} = 2, \hspace{0.2cm}m_{\rm B} = 1}
+
*Dies bedeutet eine komplexe Multiplikation und zwei komplexe Additionen: &nbsp; $\hspace{0.15 cm}\underline{a_{\rm B} = 2, \hspace{0.2cm}m_{\rm B} = 1}
  \hspace{0.05cm}.$$
+
  \hspace{0.05cm}.$
<b>d)</b>&nbsp;&nbsp;Im Gegensatz zu den Computern in der Anfangszeit nimmt heute eine Multiplikation keine wesentlich größere Rechenzeit in Anspruch als eine Addition bzw. Subtraktion. Mit den Ergebnissen aus a), b) und c) erhält man für die Gesamtzahl der Rechenoperationen:
+
 
$$ O_{\rm B} = a_{\rm B}\cdot A_{\rm A} + a_{\rm B}\cdot (A_{\rm M}
+
 
 +
'''(4)'''&nbsp; Im Gegensatz zu den ersten Computern nimmt heute eine Multiplikation keine wesentlich größere Rechenzeit in Anspruch als eine Addition bzw. Subtraktion.  
 +
 
 +
*Mit den Ergebnissen der Teilaufgaben&nbsp; '''(1)''',&nbsp; '''(2)'''&nbsp; und&nbsp; '''(3)'''&nbsp; erhält man für die Gesamtzahl der Rechenoperationen:
 +
:$$ \mathcal{O}_{\rm B} = a_{\rm B}\cdot A_{\rm A} + a_{\rm B}\cdot (A_{\rm M}
 
+M_{\rm M} ) = 2 \cdot 2 + 1 \cdot 6\hspace{0.15 cm}\underline{ = 10}
 
+M_{\rm M} ) = 2 \cdot 2 + 1 \cdot 6\hspace{0.15 cm}\underline{ = 10}
 
  \hspace{0.05cm}.$$
 
  \hspace{0.05cm}.$$
<b>e)</b>&nbsp;&nbsp;Insgesamt gibt es ld <i>N</i> Stufen, in denen jeweils <i>N</i>/2 Basisoperationen auszuführen sind:
+
 
$$O_{\rm FFT} = {\rm ld}\hspace{0.1cm}N \cdot \frac{N}{2}\cdot
+
 
O_{\rm B} = 5 \cdot N \cdot {\rm ld}\hspace{0.1cm}N$$
+
'''(5)'''&nbsp; Insgesamt gibt es&nbsp; ${\rm log_2} \ N$ Stufen, in denen jeweils&nbsp; $N/2$&nbsp; Basisoperationen auszuführen sind:
$$\Rightarrow \hspace{0.3cm}O_{\rm FFT}\hspace{0.1cm}(N=16)  =  5\cdot 16 \cdot
+
:$$\mathcal{O}_{\rm FFT} = {\rm log_2}\hspace{0.1cm}N \cdot \frac{N}{2}\cdot
4 \hspace{0.15 cm}\underline{= 320}, $$ $$\hspace{0.3cm}O_{\rm FFT}\hspace{0.1cm}(N=1024)  =  5\cdot 1024
+
\mathcal{O}_{\rm B} = 5 \cdot N \cdot {\rm log_2}\hspace{0.1cm}N$$
 +
:$$\Rightarrow \hspace{0.3cm}\mathcal{O}_{\rm FFT}\hspace{0.1cm}(N=16)  =  5\cdot 16 \cdot
 +
4 \hspace{0.15 cm}\underline{= 320}, \hspace{0.5cm}\mathcal{O}_{\rm FFT}\hspace{0.1cm}(N=1024)  =  5\cdot 1024
 
\cdot 10 \hspace{0.15 cm}\underline{= 51200}
 
\cdot 10 \hspace{0.15 cm}\underline{= 51200}
 
  \hspace{0.05cm}.$$
 
  \hspace{0.05cm}.$$
<b>f)</b>&nbsp;&nbsp; Der Rechenzeitgewinn der FFT gegenüber der herkömmlichen DFT ergibt sich zu:
+
 
$$G_{\rm FFT} = \frac{O_{\rm DFT}}{O_{\rm FFT}} = \frac{8 \cdot N^2} {5 \cdot
+
 
N \cdot {\rm ld}\hspace{0.1cm}N }= 1.6 \cdot \frac{N}{ {\rm
+
'''(6)'''&nbsp; Der Rechenzeitgewinn der FFT gegenüber der herkömmlichen DFT ergibt sich zu:
ld}\hspace{0.1cm}N}$$
+
:$$G_{\rm FFT} = \frac{\mathcal{O}_{\rm DFT}}{\mathcal{O}_{\rm FFT}} = \frac{8 \cdot N^2} {5 \cdot
$$\Rightarrow \hspace{0.3cm}G_{\rm FFT} \hspace{0.1cm}(N=16)  =  1.6 \cdot
+
N \cdot {\rm log_2}\hspace{0.1cm}N }= 1.6 \cdot \frac{N}{ {\rm
\frac{16}{ 4} \hspace{0.15 cm}\underline{= 6.4},$$ $$ \hspace{0.3cm}G_{\rm FFT} \hspace{0.1cm}(N=1024)  =  1.6 \cdot\frac{1024}{ 10}\hspace{0.15 cm}\underline{ \approx 164}
+
log_2}\hspace{0.1cm}N}$$
 +
:$$\Rightarrow \hspace{0.3cm}G_{\rm FFT} \hspace{0.1cm}(N=16)  =  1.6 \cdot
 +
\frac{16}{ 4} \hspace{0.15 cm}\underline{= 6.4}, \hspace{0.5cm}G_{\rm FFT} \hspace{0.1cm}(N=1024)  =  1.6 \cdot\frac{1024}{ 10}\hspace{0.15 cm}\underline{ \approx 164}
 
  \hspace{0.05cm}.$$
 
  \hspace{0.05cm}.$$
  

Aktuelle Version vom 22. Mai 2021, 09:21 Uhr

Letzte Stufe der FFT für  $N=8$

Der FFT–Algorithmus  (Fast Fourier Transform)  realisiert eine  „Diskrete Fouriertransformation”  mit dem kleinstmöglichen Rechenaufwand, wenn der Parameter  $N$  eine Zweierpotenz ist.

Im Einzelnen sind zur Durchführung einer FFT folgende Rechenschritte notwendig:

  • Die FFT geschieht in  ${\rm log_2} \ N$  Stufen, wobei in jeder Stufe die genau gleiche Anzahl an Rechenoperationen durchzuführen ist.
  • Die Grafik zeigt die dritte und letzte Stufe für das Beispiel  $N = 8$.  Man erkennt, dass in dieser und auch den anderen Stufen jeweils  $N/2$  Basisoperationen durchzuführen sind.
  • In jeder Basisoperation, die man häufig auch als  Butterfly  bezeichnet, werden aus den beiden komplexen Eingangsgrößen  $E_1$  und  $E_2$  zwei komplexe Ausgänge berechnet:
$$ A_1 = E_1 + E_2 \cdot w^{\hspace{0.04cm}\mu}, $$
$$ A_2 = E_1 - E_2 \cdot w^{\hspace{0.04cm} \mu}\hspace{0.05cm}.$$
  • Hierbei bezeichnet  $w = {\rm e}^{-{\rm j}\hspace{0.05cm}\cdot \hspace{0.05cm} 2\pi/N}$  den komplexen Drehfaktor.  Für  $N = 8$  ergibt sich der Wert  $w = {\rm e}^{-{\rm j} \hspace{0.05cm}\cdot \hspace{0.05cm} \pi/4} = \cos(45^\circ) - {\rm j} \cdot \sin(45^\circ)\hspace{0.05cm}.$
  • Der Exponent  $\mu$  für diesen komplexen Drehfaktor kann alle ganzzahligen Werte zwischen  $0$  und  $N/2-1$  annehmen.  Für  $N = 8$  gilt:
$$w^0 = 1,\hspace{0.2cm}w^1 = {1}/{\sqrt{2}}- {\rm j} \cdot{1}/{\sqrt{2}},\hspace{0.2cm}w^2 = - {\rm j},\hspace{0.2cm}w^3 = -{1}/{\sqrt{2}}- {\rm j} \cdot{1}/{\sqrt{2}} \hspace{0.05cm}.$$

Mit dieser Aufgabe sollen die für die FFT erforderliche Anzahl von Rechenoperationen  $(\mathcal{O}_{\rm FFT})$  ermittelt und mit dem für die DFT angebbaren Wert  $\mathcal{O}_{\rm DFT} ≈ 8\cdot N^2$  verglichen werden.

Zu beachten ist:

  • Sinnvollerweise werden die Potenzen von  $w$  vor dem eigentlichen Algorithmus berechnet und in einer Lookup–Tabelle abgelegt.  Die hierfür notwendigen Operationen sollen hier unberücksichtigt bleiben.
  • Die Bitumkehroperation – eine Umsortierung, die vor der ersten Stufe durchzuführen ist – soll bei dieser Abschätzung ebenfalls nicht berücksichtigt werden.




Hinweis:



Fragebogen

1

Wieviele reelle Additionen  $(A_{\rm A})$  erfordert eine komplexe Addition?

$A_{\rm A} \hspace{0.3cm} = \ $

2

Wieviele reelle Additionen  $(A_{\rm M})$  und Multiplikationen  $(M_{\rm M})$  sind für eine komplexe Multiplikation erforderlich?

$A_{\rm M} \hspace{0.3cm} = \ $

$M_{\rm M} \hspace{0.2cm} = \ $

3

Wieviele komplexe Additionen/Subtraktionen  $(a_{\rm B})$  erfordert eine einzige Basisoperation („Butterfly”)?
Wieviele komplexe Multiplikationen  $(m_{\rm B})$  sind pro Basisoperation notwendig?

$a_{\rm B} \hspace{0.32cm} = \ $

$m_{\rm B} \hspace{0.2cm} = \ $

4

Wieviele Rechenoperationen  (Additionen, Subtraktionen, Multiplikationen gleichermaßen)  erfordert eine Basisoperation?

$\mathcal{O}_{\rm B} \ = \ $

5

Wieviele reelle Operationen  $(\mathcal{O}_{\rm FFT})$  erfordert der gesamte FFT–Algorithmus?  Welche Werte ergeben sich für  $N = 16$  und  $N = 1024$?

$N = 16\text{:} \hspace{0.65cm} \mathcal{O}_{\rm FFT} \ = \ $

$N = 1024\text{:} \hspace{0.2cm} \mathcal{O}_{\rm FFT} \ = \ $

6

Wie groß ist der Zeitgewinn  $G_{\rm FFT} = \mathcal{O}_{\rm DFT} - \mathcal{O}_{\rm FFT}$  für  $N = 16$  und  $N = 1024$?

$N = 16\text{:} \hspace{0.65cm} G_{\rm FFT} \ = \ $

$N = 1024\text{:} \hspace{0.2cm} G_{\rm FFT} \ = \ $


Musterlösung

(1)  Jede komplexe Addition erfordert zwei reelle Additionen:

$$(R_1 + {\rm j} \cdot I_1) + (R_2 + {\rm j} \cdot I_2) = (R_1 + R_2) + {\rm j} \cdot (I_1 + I_2)\hspace{0.3cm}\Rightarrow \hspace{0.3cm} \hspace{0.15 cm}\underline{ A_{\rm A} = 2}\hspace{0.05cm}.$$


(2)  Eine jede komplexe Multiplikation benötigt vier reelle Multiplikationen und zwei reelle Additionen:

$$(R_1 + {\rm j} \cdot I_1) (R_2 + {\rm j} \cdot I_2) = (R_1 \cdot R_2 - I_1 \cdot I_2) + {\rm j} \cdot (R_1 \cdot I_2 + R_2 \cdot I_1)\hspace{0.3cm} \Rightarrow \hspace{0.3cm}\hspace{0.15 cm}\underline{A_{\rm M} = 2,\hspace{0.3cm}M_{\rm M} = 4} \hspace{0.05cm}.$$


(3)  Die Basisoperationen lauten mit den komplexen Eingangsgrößen  $E_1$,  $E_2$  und  $w^{\hspace{0.04cm}\mu}$:

$$ A_1 = E_1 + E_2 \cdot w^{\hspace{0.04cm}\mu},\hspace{0.5cm} A_2 = E_1 - E_2 \cdot w^{\hspace{0.04cm} \mu}\hspace{0.05cm}.$$
  • Dies bedeutet eine komplexe Multiplikation und zwei komplexe Additionen:   $\hspace{0.15 cm}\underline{a_{\rm B} = 2, \hspace{0.2cm}m_{\rm B} = 1} \hspace{0.05cm}.$


(4)  Im Gegensatz zu den ersten Computern nimmt heute eine Multiplikation keine wesentlich größere Rechenzeit in Anspruch als eine Addition bzw. Subtraktion.

  • Mit den Ergebnissen der Teilaufgaben  (1)(2)  und  (3)  erhält man für die Gesamtzahl der Rechenoperationen:
$$ \mathcal{O}_{\rm B} = a_{\rm B}\cdot A_{\rm A} + a_{\rm B}\cdot (A_{\rm M} +M_{\rm M} ) = 2 \cdot 2 + 1 \cdot 6\hspace{0.15 cm}\underline{ = 10} \hspace{0.05cm}.$$


(5)  Insgesamt gibt es  ${\rm log_2} \ N$ Stufen, in denen jeweils  $N/2$  Basisoperationen auszuführen sind:

$$\mathcal{O}_{\rm FFT} = {\rm log_2}\hspace{0.1cm}N \cdot \frac{N}{2}\cdot \mathcal{O}_{\rm B} = 5 \cdot N \cdot {\rm log_2}\hspace{0.1cm}N$$
$$\Rightarrow \hspace{0.3cm}\mathcal{O}_{\rm FFT}\hspace{0.1cm}(N=16) = 5\cdot 16 \cdot 4 \hspace{0.15 cm}\underline{= 320}, \hspace{0.5cm}\mathcal{O}_{\rm FFT}\hspace{0.1cm}(N=1024) = 5\cdot 1024 \cdot 10 \hspace{0.15 cm}\underline{= 51200} \hspace{0.05cm}.$$


(6)  Der Rechenzeitgewinn der FFT gegenüber der herkömmlichen DFT ergibt sich zu:

$$G_{\rm FFT} = \frac{\mathcal{O}_{\rm DFT}}{\mathcal{O}_{\rm FFT}} = \frac{8 \cdot N^2} {5 \cdot N \cdot {\rm log_2}\hspace{0.1cm}N }= 1.6 \cdot \frac{N}{ {\rm log_2}\hspace{0.1cm}N}$$
$$\Rightarrow \hspace{0.3cm}G_{\rm FFT} \hspace{0.1cm}(N=16) = 1.6 \cdot \frac{16}{ 4} \hspace{0.15 cm}\underline{= 6.4}, \hspace{0.5cm}G_{\rm FFT} \hspace{0.1cm}(N=1024) = 1.6 \cdot\frac{1024}{ 10}\hspace{0.15 cm}\underline{ \approx 164} \hspace{0.05cm}.$$