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

Aus LNTwww
Wechseln zu:Navigation, Suche
 
(22 dazwischenliegende Versionen von 3 Benutzern werden nicht angezeigt)
Zeile 1: Zeile 1:
  
{{quiz-Header|Buchseite=*Buch*/*Kapitel*
+
{{quiz-Header|Buchseite=Signaldarstellung/Spektralanalyse
 
}}
 
}}
  
[[Datei:P_ID1166__Sig_A_5_4_neu.png|250px|right|Beispiel für die Spektralanalyse (Aufgabe A5.4)]]
+
[[Datei:P_ID1166__Sig_A_5_4_neu.png|250px|right|frame|Beispiele für die Spektralanalyse]]
  
Gegeben sei der prinzipielle Zeitverlauf x(t) eines periodischen Signals. Unbekannt sind die Parameter A1, f1, A2 und f2:
+
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)
+
:$$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}.$$
\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  [[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 dem Fenster w(t) wird das Produkt y(t) = x(t) · w(t) einer Diskreten Fouriertransformation (DFT) mit den Parametern N = 512 und TP unterworfen. Die Zeitdauer TP des analysierten 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:
Für die Fensterung stehen folgende Funktionen zur Verfügung, die jeweils für |t| > TP/2 identisch 0 sind:
+
*Das  '''Rechteckfenster''':
das Rechteckfenster:
 
 
   
 
   
$${w} (\nu)  = \left\{ \begin{array}{c} 1 \\
+
:$${w} (\nu)  = \left\{ \begin{array}{c} 1 \\
 
  0 \\  \end{array} \right.\quad
 
  0 \\  \end{array} \right.\quad
 
\begin{array}{*{10}c}    {\rm{f\ddot{u}r}}
 
\begin{array}{*{10}c}    {\rm{f\ddot{u}r}}
Zeile 22: Zeile 24:
 
\end{array}$$
 
\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  '''Hanning–Fenster''':
 
   
 
   
$${w} (\nu)  = \left\{ \begin{array}{c} 0.5 + 0.5 \cdot \cos (2 \pi \cdot \frac{\nu}{N}) \\
+
:$${w} (\nu)  = \left\{ \begin{array}{c} 0.5 + 0.5 \cdot \cos (2 \pi \cdot {\nu}/{N}) \\
 
  0 \\  \end{array} \right.\quad
 
  0 \\  \end{array} \right.\quad
 
\begin{array}{*{10}c}    {\rm{f\ddot{u}r}}
 
\begin{array}{*{10}c}    {\rm{f\ddot{u}r}}
Zeile 35: Zeile 37:
 
\end{array}$$  
 
\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
 
\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
 
\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}.$$
 
{\rm si}(\pi \cdot \frac{f+f_{\rm A}}{f_{\rm A}})\hspace{0.05cm}.$$
  
Beachten Sie, dass die Frequenzauflösung fA gleich dem Kehrwert des einstellbaren Parameters TP ist. W(f) ist die Fouriertransformierte der zeitkontinuierlichen Fensterfunktion w(t), während die oben angegebene Funktion w(ν) die zeitdiskrete Gewichtungsfunktion angibt.
+
$W(f)$  ist hierbei die Fouriertransformierte der zeitkontinuierlichen Fensterfunktion  $w(t)$, während  $w(ν)$  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  $Y(f)$  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})
 
  0.5\,\, {\rm V}\cdot {\rm \delta} (f \pm  1.125\,\,{\rm kHz})
 
  \hspace{0.05cm}.$$
 
  \hspace{0.05cm}.$$
  
In der obigen Grafik sind zwei weitere Spektralfunktionen YB(f) und YC(f) abgebildet, die sich ergeben, wenn ein 1 kHz–Signal mittels DFT analysiert wird und der DFT–Parameter TP = 8.5 ms ungünstig gewählt ist.
+
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.
+
 
Hinweis: Diese Aufgabe bezieht sich auf den Theorieteil von Kapitel 5.4.
+
*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  [[Signaldarstellung/Spektralanalyse|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===
 
===Fragebogen===
Zeile 55: Zeile 70:
 
<quiz display=simple>
 
<quiz display=simple>
  
{Welche der folgenden Aussagen treffen mit Sicherheit zu, wenn die DFT das Ausgangsspektrum YA(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 TP = 4 ms verwendet.
+
- Es wurde der DFT–Parameter&nbsp; $T_{\rm P} = 4\ \text{ms}$&nbsp; verwendet.
+ Das DFT–Spektrum YA(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 TP = 8 ms, wenn das Eingangsspektrum X(f) = YA(f) anliegt? Geben Sie die Gewichte der Diraclinien bei f1 = 1 kHz und f2 = 1.125 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}) =$ { 0.625 3% } V
+
$G(f_1 = 1.000 \ \text{kHz})\ = \ $ { 0.625 3% } &nbsp;$\text{V}$
$G(f_1 = 1.125 \text{kHz}) =$ { 0.5 3% } V
+
$G(f_2 = 1.125 \ \text{kHz})\ = \ $ { 0.5 3% } &nbsp;$\text{V}$
  
{Wir betrachten das 1 kHz–Cosinussignal x(t). Welches Spektrum – YB(f) oder YC(f) – ergibt sich mit dem Rechteck– bzw. dem Hanning–Fenster, wenn der DFT-Parameter TP = 8.5 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="()"}
- YB(f) ergibt sich bei Rechteckfensterung.
+
- $Y_{\rm B}(f)$ ergibt sich bei Rechteckfensterung.
+ YB(f) ergibt sich mit dem Hanning-Fenster.
+
+ $Y_{\rm B}(f)$ ergibt sich mit dem Hanning-Fenster.
  
 
</quiz>
 
</quiz>
Zeile 76: Zeile 91:
 
===Musterlösung===
 
===Musterlösung===
 
{{ML-Kopf}}
 
{{ML-Kopf}}
'''1.''' a)  Bei Verwendung des Hanning–Fensters müssten zunächst 3 Diracfunktionen zu erkennen sein, auch wenn x(t) nur eine Frequenz beinhaltet ⇒ es wurde das Rechteckfenster verwendet.
+
'''(1)'''&nbsp; Richtig sind die <u>Lösungsvorschläge 1 und 4</u>:
Mit TP = 4 ms ergibt sich für die Frequenzauflösung fA = 1/TP = 0.25 kHz. Damit liegt die Frequenz f2 nicht im vorgegebenen Raster und Y(f) würde sich aus sehr vielen Diraclinien zusammensetzen. Das heißt: die dritte Aussage ist falsch.
+
*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.
Wie aus der nachfolgenden Grafik hervorgeht, hat x(t) die Periodendauer T0 = 8 ms. Wählt man den DFT–Parameter gleich TP = 8 ms (oder ein ganzzahliges Vielfaches davon), so stimmt die periodische Fortsetzung P{x(t)} im Intervall |t| ≤ TP/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. Richtig sind somit die Lösungsvorschläge 1 und 4.
+
*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.
 +
[[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|250px|right|Beispielsignal 1 zur Spektralanalyse (ML zu Aufgabe A5.4)]]
+
*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.  
  
b) Wegen TP = 8 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]]
 
   
 
   
$$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}, \\
 
  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
+
  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}}, \\
 
  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
+
  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}}, \\
 
  V}}, \\
  G(f = 1.250\,{\rm kHz}) & = & 0.5\, {\rm V}\cdot 0.25 = 0.125\, {\rm
+
  G(f = 1.250\,{\rm kHz}) & = 0.5\, {\rm V}\cdot 0.25 = 0.125\, {\rm
 
  V}
 
  V}
  \hspace{0.05cm}.$$
+
  \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|250px|right|Beispielsignal 2 zur Spektralanalyse (ML zu Aufgabe A5.4)]]
+
Die Grafik zeigt die Abschwächung der Ränder durch die Gewichtungsfunktion&nbsp; $w(t)$&nbsp; des Hanning–Fensters.
c)  Das Rechteck–Fenster liefert dann ein sehr stark verfälschtes Ergebnis, wenn die Fensterbreite TP (wie hier) nicht an die Frequenz des Cosinussignals angepasst ist. In diesem Fall ist das Hanning–Fenster besser geeignet. Daraus folgt: Richtig ist der zweite Lösungsvorschlag.
+
<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 frequenzdisktrete Signaldarstellung^]]
+
[[Category:Aufgaben zu Signaldarstellung|^5. Zeit- und frequenzdiskrete Signaldarstellung^]]

Aktuelle Version vom 19. Mai 2021, 15:42 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.