Aufgaben:Aufgabe 5.4: Vergleich von Rechteck- und Hanningfenster: Unterschied zwischen den Versionen

Aus LNTwww
Wechseln zu:Navigation, Suche
Maintenance script (Diskussion | Beiträge)
Add English interlanguage link
 
(10 dazwischenliegende Versionen von 3 Benutzern werden nicht angezeigt)
Zeile 3: Zeile 3:
}}
}}


[[Datei:P_ID1166__Sig_A_5_4_neu.png|250px|right|Beispiel für die Spektralanalyse]]
[[Datei:P_ID1166__Sig_A_5_4_neu.png|250px|right|frame|Beispiele für die Spektralanalyse]]


Gegeben sei der prinzipielle Zeitverlauf eines periodischen Signals:  
Gegeben sei der prinzipielle Zeitverlauf eines periodischen Signals:  
   
   
$$x(t)  =  A_1 \cdot \cos (2 \pi \cdot f_1 \cdot t)  +  A_2 \cdot \cos (2 \pi \cdot f_2 \cdot t) \hspace{0.05cm}.$$
:$$x(t)  =  A_1 \cdot \cos (2 \pi \cdot f_1 \cdot t)  +  A_2 \cdot \cos (2 \pi \cdot f_2 \cdot t) \hspace{0.05cm}.$$


Unbekannt und damit zu schätzen seien dessen Parameter $A_1$, $f_1$, $A_2$ und $f_2$.
Unbekannt und damit zu schätzen sind dessen Parameter  $A_1$,  $f_1$,  $A_2$  und  $f_2$.


Nach Gewichtung des Signals mit der Fensterfunktion $w(t)$ wird das Produkt $y(t) = x(t) \cdot w(t)$ einer [[Signaldarstellung/Diskrete_Fouriertransformation_(DFT)|Diskreten Fouriertransformation]] (DFT) mit den Parametern $N = 512$ und $T_{\rm P}$ unterworfen. Die Zeitdauer $T_{\rm P}$ des zu analysierenden  Signalausschnitts kann vom Benutzer beliebig eingestellt werden.
Nach Gewichtung des Signals mit der Fensterfunktion  $w(t)$  wird das Produkt  $y(t) = x(t) \cdot w(t)$  einer  [[Signaldarstellung/Diskrete_Fouriertransformation_(DFT)|Diskreten Fouriertransformation]]  (DFT) mit den Parametern  $N = 512$  und  $T_{\rm P}$  unterworfen.  Die Zeitdauer  $T_{\rm P}$  des zu analysierenden  Signalausschnitts kann vom Benutzer beliebig eingestellt werden.


Für die Fensterung stehen folgende Funktionen zur Verfügung, die jeweils für $|t| > T_{\rm P}/2$ identisch 0 sind:
Für die Fensterung stehen zwei Funktionen zur Verfügung, die für  $|t| > T_{\rm P}/2$  jeweils Null sind:
*Das Rechteckfenster:
*Das  '''Rechteckfenster''':
   
   
$${w} (\nu)  = \left\{ \begin{array}{c} 1 \\
:$${w} (\nu)  = \left\{ \begin{array}{c} 1 \\0 \\  \end{array} \right.\quad\begin{array}{*{10}c}    {\rm{f\ddot{u}r}}\\    \\ \end{array}\begin{array}{*{20}c}-N/2 \le \nu < N/2 \hspace{0.05cm}, \\{\rm sonst} \hspace{0.05cm}, \\\end{array}$$
0 \\  \end{array} \right.\quad
\begin{array}{*{10}c}    {\rm{f\ddot{u}r}}
\\    \\ \end{array}\begin{array}{*{20}c}
-N/2 \le \nu < N/2 \hspace{0.05cm}, \\
{\rm sonst} \hspace{0.05cm}, \\
\end{array}$$
   
   
$$W(f) ={1}/{f_{\rm A}}\cdot {\rm si}(\pi \cdot
:$$W(f) ={1}/{f_{\rm A}}\cdot {\rm si}(\pi \cdot{f}/{f_{\rm A}})\hspace{0.05cm},$$
{f}/{f_{\rm A}})\hspace{0.05cm},$$


*das Hanning–Fenster:
*das&nbsp; '''Hanning–Fenster''':
   
   
$${w} (\nu)  = \left\{ \begin{array}{c} 0.5 + 0.5 \cdot \cos (2 \pi \cdot {\nu}/{N}) \\
:$${w} (\nu)  = \left\{ \begin{array}{c} 0.5 + 0.5 \cdot \cos (2 \pi \cdot {\nu}/{N}) \\0 \\  \end{array} \right.\quad\begin{array}{*{10}c}    {\rm{f\ddot{u}r}}\\    \\ \end{array}\begin{array}{*{20}c}-N/2 \le \nu < N/2 \hspace{0.05cm}, \\{\rm sonst} \hspace{0.05cm}, \\\end{array}$$  
0 \\  \end{array} \right.\quad
\begin{array}{*{10}c}    {\rm{f\ddot{u}r}}
\\    \\ \end{array}\begin{array}{*{20}c}
-N/2 \le \nu < N/2 \hspace{0.05cm}, \\
{\rm sonst} \hspace{0.05cm}, \\
\end{array}$$  


$$W(f) ={0.5}/{f_{\rm A}}\cdot {\rm si}(\pi \cdot
:$$W(f) ={0.5}/{f_{\rm A}}\cdot {\rm si}(\pi \cdot\frac{f}{f_{\rm A}})+ {0.25}/{f_{\rm A}}\cdot {\rm si}(\pi\cdot \frac{f-f_{\rm A}}{f_{\rm A}})+ {0.5}/{f_{\rm A}}\cdot{\rm si}(\pi \cdot \frac{f+f_{\rm A}}{f_{\rm A}})\hspace{0.05cm}.$$
\frac{f}{f_{\rm A}})+ {0.25}/{f_{\rm A}}\cdot {\rm si}(\pi
\cdot \frac{f-f_{\rm A}}{f_{\rm A}})+ {0.5}/{f_{\rm A}}\cdot
{\rm si}(\pi \cdot \frac{f+f_{\rm A}}{f_{\rm A}})\hspace{0.05cm}.$$


Beachten Sie, dass die Frequenzauflösung $f_{\rm A}$ gleich dem Kehrwert des einstellbaren Parameters $T_{\rm P}$ ist. $W(f)$ ist die Fouriertransformierte der zeitkontinuierlichen Fensterfunktion $w(t)$, während die oben angegebene Funktion $w(ν)$ die zeitdiskrete Gewichtungsfunktion angibt.
$W(f)$&nbsp; ist hierbei die Fouriertransformierte der zeitkontinuierlichen Fensterfunktion&nbsp; $w(t)$, während&nbsp; $w(ν)$&nbsp; die zeitdiskrete Gewichtungsfunktion angibt.


Im Laufe der Aufgabe wird auf verschiedene Spektralfunktionen $Y(f)$ Bezug genommen, zum Beispiel auf
In der Aufgabe wird auf verschiedene Spektralfunktionen&nbsp; $Y(f)$&nbsp; Bezug genommen, zum Beispiel auf
   
   
$$Y_{\rm A}(f) = 1\, {\rm V}\cdot {\rm \delta} (f \pm  1\,\,{\rm kHz})+
:$$Y_{\rm A}(f) = 1\, {\rm V}\cdot {\rm \delta} (f \pm  1\,\,{\rm kHz})+0.5\,\, {\rm V}\cdot {\rm \delta} (f \pm  1.125\,\,{\rm kHz})\hspace{0.05cm}.$$
0.5\,\, {\rm V}\cdot {\rm \delta} (f \pm  1.125\,\,{\rm kHz})
 
\hspace{0.05cm}.$$
In der obigen Grafik sind zwei weitere Spektralfunktionen&nbsp; $Y_{\rm B}(f)$&nbsp; und&nbsp; $Y_{\rm C}(f)$&nbsp; abgebildet, die sich ergeben, wenn ein&nbsp; $1 \ \text{kHz}$–Signal&nbsp; mittels DFT analysiert wird und der DFT–Parameter&nbsp; $T_{\rm P} = 8.5 \ \text{ms}$&nbsp; ungünstig gewählt ist.
 
*Für eines der Bilder ist das Rechteckfenster zugrunde gelegt, für das andere das Hanning–Fenster.
*Nicht angegeben wird, welche Grafik zu welchem Fenster gehört.
 
 


In der obigen Grafik sind zwei weitere Spektralfunktionen $Y_{\rm B}(f)$ und $Y_{\rm C}(f)$ abgebildet, die sich ergeben, wenn ein $1 \ \text{kHz}$–Signal mittels DFT analysiert wird und der DFT–Parameter $T_{\rm P} = 8.5 \ \text{ms}$ ungünstig gewählt ist.


Für eines der Bilder ist das Rechteckfenster zugrundegelegt, für das andere das Hanning–Fenster. Nicht angegeben wird, welche Spektralfunktion zu welchem Fenster gehört.


''Hinweise:''  
''Hinweise:''  
*Die Aufgabe gehört zum  Kapitel [[Signaldarstellung/Spektralanalyse|Spektralanalyse]].
*Die Aufgabe gehört zum  Kapitel&nbsp; [[Signaldarstellung/Spektralanalyse|Spektralanalyse]].
*Sollte die Eingabe des Zahlenwertes &bdquo;0&rdquo; erforderlich sein, so geben Sie bitte &bdquo;0.&rdquo; ein.
*Beachten Sie, dass die Frequenzauflösung&nbsp; $f_{\rm A}$&nbsp; gleich dem Kehrwert des einstellbaren Parameters&nbsp; $T_{\rm P}$&nbsp; ist.
*Unglücklichweise kollidieren die Indizes von&nbsp; $f_{\rm A}$&nbsp; und&nbsp; $Y_{\rm A}(f)$.&nbsp; Es ist offensichtlich, dass diese nicht in Zusammenhang stehen.&nbsp; Nur zur Sicherheit weisen wir darauf hin.  




Zeile 63: Zeile 52:
<quiz display=simple>
<quiz display=simple>


{Welche der folgenden Aussagen treffen mit Sicherheit zu, wenn die DFT das Ausgangsspektrum $Y_{\rm A}(f)$ anzeigt?
{Welche der folgenden Aussagen treffen mit Sicherheit zu, wenn die DFT das Ausgangsspektrum&nbsp; $Y_{\rm A}(f)$&nbsp; anzeigt?
|type="[]"}
|type="[]"}
+ Zur Gewichtung wurde das Rechteckfenster verwendet.
+ Zur Gewichtung wurde das Rechteckfenster verwendet.
- Zur Gewichtung wurde das Hanning–Fenster verwendet.
- Zur Gewichtung wurde das Hanning–Fenster verwendet.
- Es wurde der DFT–Parameter $T_{\rm P} = 4\ \text{ms}$ verwendet.
- Es wurde der DFT–Parameter&nbsp; $T_{\rm P} = 4\ \text{ms}$&nbsp; verwendet.
+ Das DFT–Spektrum $Y_{\rm A}(f)$ ist identisch mit dem Spektrum $X(f)$.
+ Das DFT–Spektrum $Y_{\rm A}(f)$ ist identisch mit dem tatsächlichen Spektrum $X(f)$.


{Wie lautet $Y(f)$ bei Verwendung des Hanning–Fensters und  $T_{\rm P} = 8 \ \text{ms}$, wenn das Eingangsspektrum $X(f) = Y_{\rm A}(f)$ anliegt? Geben Sie die Gewichte der Diraclinien bei $f_1= 1\ \text{kHz}$ und $f_2 = 1.125\ \text{kHz}$ an.
{Wie lautet&nbsp; $Y(f)$&nbsp; bei Verwendung des Hanning–Fensters und&nbsp; $T_{\rm P} = 8 \ \text{ms}$, wenn das Eingangsspektrum&nbsp; $X(f) = Y_{\rm A}(f)$&nbsp; anliegt? <br>Geben Sie die Gewichte der Diraclinien bei&nbsp; $f_1= 1\ \text{kHz}$&nbsp; und&nbsp; $f_2 = 1.125\ \text{kHz}$&nbsp; an.
|type="{}"}
|type="{}"}
$G(f_1 = 1 \ \text{kHz}) &nbsp;=$ { 0.625 3% } &nbsp;$\text{V}$
$G(f_1 = 1.000 \ \text{kHz})\ = \ $ { 0.625 3% } &nbsp;$\text{V}$
$G(f_1 = 1.125 \ \text{kHz}) &nbsp;=$ { 0.5 3% } &nbsp;$\text{V}$
$G(f_2 = 1.125 \ \text{kHz})\ = \ $ { 0.5 3% } &nbsp;$\text{V}$


{Wir betrachten das $1\ \text{kHz}$–Cosinussignal $x(t)$. Welches Spektrum - $Y_{\rm B}(f)$ oder $Y_{\rm C}(f)$ – ergibt sich mit dem Rechteck– bzw. dem Hanning–Fenster, wenn der DFT-Parameter  $T_{\rm P} = 8.5 \ \text{ms}$ ungünstig gewählt ist?
{Wir betrachten das&nbsp; $1\ \text{kHz}$–Cosinussignal&nbsp; $x(t)$.&nbsp; Welches Spektrum -&nbsp; $Y_{\rm B}(f)$&nbsp; oder&nbsp; $Y_{\rm C}(f)$&nbsp; – ergibt sich mit dem Rechteck– bzw. dem Hanning–Fenster, wenn der DFT-Parameter&nbsp; $T_{\rm P} = 8.5 \ \text{ms}$&nbsp; ungünstig gewählt ist?
|type="[]"}
|type="()"}
- $Y_{\rm B}(f)$ergibt sich bei Rechteckfensterung.
- $Y_{\rm B}(f)$ ergibt sich bei Rechteckfensterung.
+ $Y_{\rm B}(f)$ ergibt sich mit dem Hanning-Fenster.
+ $Y_{\rm B}(f)$ ergibt sich mit dem Hanning-Fenster.


Zeile 84: Zeile 73:
===Musterlösung===
===Musterlösung===
{{ML-Kopf}}
{{ML-Kopf}}
'''1.''' Richtig sind die <u>Lösungsvorschläge 1 und 4</u>:
'''(1)'''&nbsp; Richtig sind die <u>Lösungsvorschläge 1 und 4</u>:
*Bei Verwendung des Hanning–Fensters müssten selbst dann drei Diracfunktionen zu erkennen sein, auch wenn $x(t)$ nur eine Frequenz beinhaltet &nbsp; ⇒ &nbsp; es wurde das Rechteckfenster verwendet.
*Bei Verwendung des Hanning–Fensters müssten selbst dann drei Diracfunktionen zu erkennen sein, wenn&nbsp; $x(t)$&nbsp; nur eine Frequenz beinhaltet &nbsp; ⇒ &nbsp; es wurde das Rechteckfenster verwendet.
*Mit $T_{\rm P} = 4 \ \text{ms}$ ergibt sich für die Frequenzauflösung $f_{\rm A}= 1/T_{\rm P} = 0.25 \ \text{kHz}$ Damit liegt die Frequenz $f_2$ nicht im vorgegebenen Raster und $Y(f)$ würde sich aus sehr vielen Diraclinien zusammensetzen. Das heißt: die dritte Aussage ist falsch.
*Mit&nbsp; $T_{\rm P} = 4 \ \text{ms}$&nbsp; ergibt sich für die Frequenzauflösung&nbsp; $f_{\rm A}= 1/T_{\rm P} = 0.25 \ \text{kHz}$.&nbsp; Damit liegt die Frequenz&nbsp; $f_2$&nbsp; nicht im vorgegebenen Raster und&nbsp; $Y(f)$&nbsp; würde sich aus sehr vielen Diraclinien zusammensetzen.&nbsp; Das heißt: &nbsp; die dritte Aussage ist falsch.
*Wie aus der nachfolgenden Grafik hervorgeht, hat $x(t)$ die Periodendauer $T_{\rm 0} = 8 \ \text{ms}$. Wählt man den DFT–Parameter gleich $T_{\rm P} = 4 \ \text{ms}$ (oder ein ganzzahliges Vielfaches davon), so stimmt die periodische Fortsetzung ${\rm P}\{ x(t)\} $ im Intervall $|t| \leq T_{\rm P}/2$ mit $x(t)$ überein, so dass sich die Gewichtungsfunktion $w(t)$ nicht störend auswirkt: Das DFT–Spektrum $Y(f)$ stimmt somit mit dem tatsächlichen Spektrum überein.
[[Datei:P_ID1167__Sig_A_5_4a.png|right|frame|Signal&nbsp; $y(t)$&nbsp; nach Rechteck&ndash;Fensterung ]]


[[Datei:P_ID1167__Sig_A_5_4a.png|Beispielsignal 1 zur Spektralanalyse]]
*Wie aus der Grafik hervorgeht, hat&nbsp; $x(t)$&nbsp; die Periodendauer&nbsp; $T_{\rm 0} = 8 \ \text{ms}$. Wählt man den DFT–Parameter gleich&nbsp; $T_{\rm P} = 4 \ \text{ms}$&nbsp; (oder ein ganzzahliges Vielfaches davon), so stimmt die periodische Fortsetzung&nbsp; ${\rm P}\{ x(t)\} $&nbsp; im Intervall&nbsp; $|t| \leq T_{\rm P}/2$&nbsp; mit&nbsp; $x(t)$&nbsp; überein, so dass sich die Gewichtungsfunktion&nbsp; $w(t)$&nbsp; nicht störend auswirkt: &nbsp;
*Das DFT–Spektrum&nbsp; $Y(f)$&nbsp; stimmt somit mit dem tatsächlichen Spektrum überein.  




'''2.''' Wegen $T_{\rm 0} = 8 \ \text{ms}$ setzt sich das Hanning–Spektrum $W(f)$ aus drei Diracfunktionen bei positiven Frequenzen und drei dazu achsensymmetrischen Diracs bei negativen Frequenzen zusammen. Für die positiven Frequenzen lautet die Spektralfunktion:
 
'''(2)'''&nbsp; Wegen $T_{\rm 0} = 8 \ \text{ms}$&nbsp; setzt sich das Hanning–Spektrum&nbsp; $W(f)$&nbsp;
*aus drei Diracfunktionen bei positiven Frequenzen  
*und drei dazu achsensymmetrischen Diracs bei negativen Frequenzen  
 
 
zusammen.&nbsp; Für die positiven Frequenzen lautet die Spektralfunktion:
   
   
$$W(f) =0.5\cdot {\rm \delta}(f) + 0.25\cdot {\rm \delta}(f-f_{\rm A})+ 0.25\cdot {\rm \delta}(f+f_{\rm A})\hspace{0.05cm}.$$
:$$W(f) =0.5\cdot {\rm \delta}(f) + 0.25\cdot {\rm \delta}(f-f_{\rm A})+ 0.25\cdot {\rm \delta}(f+f_{\rm A})\hspace{0.05cm}.$$


Das Ausgangsspektrum ergibt sich aus der Faltung zwischen $X(f)$ und $W(f)$. Bei positiven Frequenzen ergeben sich nun vier Diracs mit folgenden Gewichten:
Das Ausgangsspektrum ergibt sich aus der Faltung zwischen&nbsp; $X(f)$&nbsp; und&nbsp; $W(f)$.&nbsp; Bei positiven Frequenzen ergeben sich nun vier Diracs mit folgenden Gewichten:
[[Datei:P_ID1169__Sig_A_5_4b.png|right|frame|Signal&nbsp; $y(t)$&nbsp; nach Hanning&ndash;Fensterung]]
   
   
$$\begin{align*} G(f = 0.875\,{\rm kHz}) & = 1\, {\rm V}\cdot 0.25 = 0.250\, {\rm
:$$\begin{align*} G(f = 0.875\,{\rm kHz}) & = 1\, {\rm V}\cdot 0.25 = 0.250\, {\rm V}, \\G(f = f_1 = 1.000\,{\rm kHz}) & = 1\, {\rm V}\cdot 0.5 + 0.5\, {\rm V}\cdot 0.25 \hspace{0.15 cm}\underline{ = 0.625\, {\rm V}}, \\G(f = f_2 = 1.125\,{\rm kHz}) & = 1\, {\rm V}\cdot 0.25 + 0.5\, {\rm V}\cdot 0.5  \hspace{0.15 cm}\underline{= 0.500\, {\rm V}}, \\G(f = 1.250\,{\rm kHz}) & = 0.5\, {\rm V}\cdot 0.25 = 0.125\, {\rm V}\hspace{0.05cm}.\end{align*}$$
V}, \\
G(f = f_1 = 1.000\,{\rm kHz}) & = 1\, {\rm V}\cdot 0.5 + 0.5\, {\rm V}\cdot 0.25 \hspace{0.15 cm}\underline{ = 0.625\, {\rm
V}}, \\
G(f = f_2 = 1.125\,{\rm kHz}) & = 1\, {\rm V}\cdot 0.25 + 0.5\, {\rm V}\cdot 0.5  \hspace{0.15 cm}\underline{= 0.500\, {\rm
V}}, \\
G(f = 1.250\,{\rm kHz}) & = 0.5\, {\rm V}\cdot 0.25 = 0.125\, {\rm
V}
\hspace{0.05cm}.\end{align*}$$
 
Die folgende Grafik zeigt die Abschwächung der Ränder durch die Gewichtungsfunktion $w(t)$ des Hanning–Fensters.
 
[[Datei:P_ID1169__Sig_A_5_4b.png|Beispielsignal 2 zur Spektralanalyse]]


'''3.''' Das Rechteck–Fenster liefert dann ein sehr stark verfälschtes Ergebnis, wenn die Fensterbreite $T_{\rm P}$ (wie hier) nicht an die Frequenz des Cosinussignals angepasst ist. In diesem Fall ist das Hanning–Fenster besser geeignet. Daraus folgt: Richtig ist der <u>zweite Lösungsvorschlag</u>.
Die Grafik zeigt die Abschwächung der Ränder durch die Gewichtungsfunktion&nbsp; $w(t)$&nbsp; des Hanning–Fensters.
<br clear=all>
'''(3)'''&nbsp; Richtig ist der <u>zweite Lösungsvorschlag</u>:
*Das Rechteck–Fenster liefert dann ein sehr stark verfälschtes Ergebnis, wenn die Fensterbreite&nbsp; $T_{\rm P}$&nbsp; (wie hier) nicht an die Frequenz des Cosinussignals angepasst ist.  
*In diesem Fall ist das Hanning–Fenster besser geeignet.&nbsp; In diesem Fall ergibt sich das gemessene Spektrum&nbsp; $Y_{\rm B}(f)$.
* Aus dem Spektrum&nbsp; $Y_{\rm C}(f)$&nbsp; ist die gesuchte $1\ \rm kHz$&ndash;Linie schlechter zu erkennen.&nbsp; Dieses Spektrum&nbsp; $Y_{\rm C}(f)$&nbsp; ergibt sich nach Rechteck&ndash;Fensterung.  
{{ML-Fuß}}
{{ML-Fuß}}


__NOEDITSECTION__
__NOEDITSECTION__
[[Category:Aufgaben zu Signaldarstellung|^5. Zeit- und frequenzdiskrete Signaldarstellung^]]
[[Category:Aufgaben zu Signaldarstellung|^5. Zeit- und frequenzdiskrete Signaldarstellung^]]
[[en:Aufgaben:Exercise_5.4:_Comparison_of_Rectangular_and_Hanning_Window]]

Aktuelle Version vom 16. März 2026, 14:38 Uhr

Beispiele für die Spektralanalyse

Gegeben sei der prinzipielle Zeitverlauf eines periodischen Signals:

$$x(t) = A_1 \cdot \cos (2 \pi \cdot f_1 \cdot t) + A_2 \cdot \cos (2 \pi \cdot f_2 \cdot t) \hspace{0.05cm}.$$

Unbekannt und damit zu schätzen sind dessen Parameter  $A_1$,  $f_1$,  $A_2$  und  $f_2$.

Nach Gewichtung des Signals mit der Fensterfunktion  $w(t)$  wird das Produkt  $y(t) = x(t) \cdot w(t)$  einer  Diskreten Fouriertransformation  (DFT) mit den Parametern  $N = 512$  und  $T_{\rm P}$  unterworfen.  Die Zeitdauer  $T_{\rm P}$  des zu analysierenden Signalausschnitts kann vom Benutzer beliebig eingestellt werden.

Für die Fensterung stehen zwei Funktionen zur Verfügung, die für  $|t| > T_{\rm P}/2$  jeweils Null sind:

  • Das  Rechteckfenster:
$${w} (\nu) = \left\{ \begin{array}{c} 1 \\0 \\ \end{array} \right.\quad\begin{array}{*{10}c} {\rm{f\ddot{u}r}}\\ \\ \end{array}\begin{array}{*{20}c}-N/2 \le \nu < N/2 \hspace{0.05cm}, \\{\rm sonst} \hspace{0.05cm}, \\\end{array}$$
$$W(f) ={1}/{f_{\rm A}}\cdot {\rm si}(\pi \cdot{f}/{f_{\rm A}})\hspace{0.05cm},$$
  • das  Hanning–Fenster:
$${w} (\nu) = \left\{ \begin{array}{c} 0.5 + 0.5 \cdot \cos (2 \pi \cdot {\nu}/{N}) \\0 \\ \end{array} \right.\quad\begin{array}{*{10}c} {\rm{f\ddot{u}r}}\\ \\ \end{array}\begin{array}{*{20}c}-N/2 \le \nu < N/2 \hspace{0.05cm}, \\{\rm sonst} \hspace{0.05cm}, \\\end{array}$$
$$W(f) ={0.5}/{f_{\rm A}}\cdot {\rm si}(\pi \cdot\frac{f}{f_{\rm A}})+ {0.25}/{f_{\rm A}}\cdot {\rm si}(\pi\cdot \frac{f-f_{\rm A}}{f_{\rm A}})+ {0.5}/{f_{\rm A}}\cdot{\rm si}(\pi \cdot \frac{f+f_{\rm A}}{f_{\rm A}})\hspace{0.05cm}.$$

$W(f)$  ist hierbei die Fouriertransformierte der zeitkontinuierlichen Fensterfunktion  $w(t)$, während  $w(ν)$  die zeitdiskrete Gewichtungsfunktion angibt.

In der Aufgabe wird auf verschiedene Spektralfunktionen  $Y(f)$  Bezug genommen, zum Beispiel auf

$$Y_{\rm A}(f) = 1\, {\rm V}\cdot {\rm \delta} (f \pm 1\,\,{\rm kHz})+0.5\,\, {\rm V}\cdot {\rm \delta} (f \pm 1.125\,\,{\rm kHz})\hspace{0.05cm}.$$

In der obigen Grafik sind zwei weitere Spektralfunktionen  $Y_{\rm B}(f)$  und  $Y_{\rm C}(f)$  abgebildet, die sich ergeben, wenn ein  $1 \ \text{kHz}$–Signal  mittels DFT analysiert wird und der DFT–Parameter  $T_{\rm P} = 8.5 \ \text{ms}$  ungünstig gewählt ist.

  • Für eines der Bilder ist das Rechteckfenster zugrunde gelegt, für das andere das Hanning–Fenster.
  • Nicht angegeben wird, welche Grafik zu welchem Fenster gehört.



Hinweise:

  • Die Aufgabe gehört zum Kapitel  Spektralanalyse.
  • Beachten Sie, dass die Frequenzauflösung  $f_{\rm A}$  gleich dem Kehrwert des einstellbaren Parameters  $T_{\rm P}$  ist.
  • Unglücklichweise kollidieren die Indizes von  $f_{\rm A}$  und  $Y_{\rm A}(f)$.  Es ist offensichtlich, dass diese nicht in Zusammenhang stehen.  Nur zur Sicherheit weisen wir darauf hin.


Fragebogen

1 Welche der folgenden Aussagen treffen mit Sicherheit zu, wenn die DFT das Ausgangsspektrum  $Y_{\rm A}(f)$  anzeigt?

Zur Gewichtung wurde das Rechteckfenster verwendet.
Zur Gewichtung wurde das Hanning–Fenster verwendet.
Es wurde der DFT–Parameter  $T_{\rm P} = 4\ \text{ms}$  verwendet.
Das DFT–Spektrum $Y_{\rm A}(f)$ ist identisch mit dem tatsächlichen Spektrum $X(f)$.

2 Wie lautet  $Y(f)$  bei Verwendung des Hanning–Fensters und  $T_{\rm P} = 8 \ \text{ms}$, wenn das Eingangsspektrum  $X(f) = Y_{\rm A}(f)$  anliegt?
Geben Sie die Gewichte der Diraclinien bei  $f_1= 1\ \text{kHz}$  und  $f_2 = 1.125\ \text{kHz}$  an.

$G(f_1 = 1.000 \ \text{kHz})\ = \ $  $\text{V}$
$G(f_2 = 1.125 \ \text{kHz})\ = \ $  $\text{V}$

3 Wir betrachten das  $1\ \text{kHz}$–Cosinussignal  $x(t)$.  Welches Spektrum -  $Y_{\rm B}(f)$  oder  $Y_{\rm C}(f)$  – ergibt sich mit dem Rechteck– bzw. dem Hanning–Fenster, wenn der DFT-Parameter  $T_{\rm P} = 8.5 \ \text{ms}$  ungünstig gewählt ist?

$Y_{\rm B}(f)$ ergibt sich bei Rechteckfensterung.
$Y_{\rm B}(f)$ ergibt sich mit dem Hanning-Fenster.


Musterlösung

(1)  Richtig sind die Lösungsvorschläge 1 und 4:

  • Bei Verwendung des Hanning–Fensters müssten selbst dann drei Diracfunktionen zu erkennen sein, wenn  $x(t)$  nur eine Frequenz beinhaltet   ⇒   es wurde das Rechteckfenster verwendet.
  • Mit  $T_{\rm P} = 4 \ \text{ms}$  ergibt sich für die Frequenzauflösung  $f_{\rm A}= 1/T_{\rm P} = 0.25 \ \text{kHz}$.  Damit liegt die Frequenz  $f_2$  nicht im vorgegebenen Raster und  $Y(f)$  würde sich aus sehr vielen Diraclinien zusammensetzen.  Das heißt:   die dritte Aussage ist falsch.
Signal  $y(t)$  nach Rechteck–Fensterung
  • Wie aus der Grafik hervorgeht, hat  $x(t)$  die Periodendauer  $T_{\rm 0} = 8 \ \text{ms}$. Wählt man den DFT–Parameter gleich  $T_{\rm P} = 4 \ \text{ms}$  (oder ein ganzzahliges Vielfaches davon), so stimmt die periodische Fortsetzung  ${\rm P}\{ x(t)\} $  im Intervall  $|t| \leq T_{\rm P}/2$  mit  $x(t)$  überein, so dass sich die Gewichtungsfunktion  $w(t)$  nicht störend auswirkt:  
  • Das DFT–Spektrum  $Y(f)$  stimmt somit mit dem tatsächlichen Spektrum überein.


(2)  Wegen $T_{\rm 0} = 8 \ \text{ms}$  setzt sich das Hanning–Spektrum  $W(f)$ 

  • aus drei Diracfunktionen bei positiven Frequenzen
  • und drei dazu achsensymmetrischen Diracs bei negativen Frequenzen


zusammen.  Für die positiven Frequenzen lautet die Spektralfunktion:

$$W(f) =0.5\cdot {\rm \delta}(f) + 0.25\cdot {\rm \delta}(f-f_{\rm A})+ 0.25\cdot {\rm \delta}(f+f_{\rm A})\hspace{0.05cm}.$$

Das Ausgangsspektrum ergibt sich aus der Faltung zwischen  $X(f)$  und  $W(f)$.  Bei positiven Frequenzen ergeben sich nun vier Diracs mit folgenden Gewichten:

Signal  $y(t)$  nach Hanning–Fensterung
$$\begin{align*} G(f = 0.875\,{\rm kHz}) & = 1\, {\rm V}\cdot 0.25 = 0.250\, {\rm V}, \\G(f = f_1 = 1.000\,{\rm kHz}) & = 1\, {\rm V}\cdot 0.5 + 0.5\, {\rm V}\cdot 0.25 \hspace{0.15 cm}\underline{ = 0.625\, {\rm V}}, \\G(f = f_2 = 1.125\,{\rm kHz}) & = 1\, {\rm V}\cdot 0.25 + 0.5\, {\rm V}\cdot 0.5 \hspace{0.15 cm}\underline{= 0.500\, {\rm V}}, \\G(f = 1.250\,{\rm kHz}) & = 0.5\, {\rm V}\cdot 0.25 = 0.125\, {\rm V}\hspace{0.05cm}.\end{align*}$$

Die Grafik zeigt die Abschwächung der Ränder durch die Gewichtungsfunktion  $w(t)$  des Hanning–Fensters.
(3)  Richtig ist der zweite Lösungsvorschlag:

  • Das Rechteck–Fenster liefert dann ein sehr stark verfälschtes Ergebnis, wenn die Fensterbreite  $T_{\rm P}$  (wie hier) nicht an die Frequenz des Cosinussignals angepasst ist.
  • In diesem Fall ist das Hanning–Fenster besser geeignet.  In diesem Fall ergibt sich das gemessene Spektrum  $Y_{\rm B}(f)$.
  • Aus dem Spektrum  $Y_{\rm C}(f)$  ist die gesuchte $1\ \rm kHz$–Linie schlechter zu erkennen.  Dieses Spektrum  $Y_{\rm C}(f)$  ergibt sich nach Rechteck–Fensterung.