Loading [MathJax]/jax/output/HTML-CSS/fonts/TeX/fontdata.js

Difference between revisions of "Aufgaben:Exercise 5.2Z: DFT of a Triangular Pulse"

From LNTwww
Line 3: Line 3:
 
}}
 
}}
  
[[File:P_ID1140__Sig_Z_5_2.png|right|frame|Diskretisierung eines Dreieckimpulses]]
+
[[File:P_ID1140__Sig_Z_5_2.png|right|frame|Discretisation of a triangular pulse]]
Betrachtet wird der skizzierte Dreieckimpuls
+
Consider the sketched triangular momentum
 
:$$x(t) = \left\{ \begin{array}{l} A \cdot \left( 1 - {|t|}/{T} \right ) \\  
 
:$$x(t) = \left\{ \begin{array}{l} A \cdot \left( 1 - {|t|}/{T} \right ) \\  
 
\hspace{0.25cm} 0 \\  \end{array} \right.\quad  
 
\hspace{0.25cm} 0 \\  \end{array} \right.\quad  
Line 12: Line 12:
 
{\left|\hspace{0.005cm} t \hspace{0.05cm} \right| > T.}  \\
 
{\left|\hspace{0.005cm} t \hspace{0.05cm} \right| > T.}  \\
 
\end{array}$$
 
\end{array}$$
Die Signalparameter haben folgende Werte:
+
The signal parameters have the following values:
  
 
* Amplitude  A=4 V,
 
* Amplitude  A=4 V,
  
* äquivalente Impulsdauer  Δt=T=1ms.
+
* equivalent pulse duration  Δt=T=1ms.
  
  
Das Spektrum  X(f)  erhält man durch Anwendung des  [[Signal_Representation/Fourier_Transform_and_Its_Inverse#Das_erste_Fourierintegral|ersten Fourierintegrals]]:
+
The spectrum  X(f)  is obtained by applying  [[Signal_Representation/Fourier_Transform_and_Its_Inverse#The_First_Fourier_Integral|the first Fourier Integral]]:
 
:X(f)=ATsi2(πfT).
 
:X(f)=ATsi2(πfT).
  
Die Spektralfunktion soll nun durch eine ''Diskrete Fouriertransformation''  (DFT) mit  N=8  angenähert werden, wobei die  N  Koeffizienten für den Zeitbereich   ⇒   d(0), ... , d(7)  der Grafik entnommen werden können.
+
The spectral function is now to be approximated by a ''Discrete Fourier Transform''  (DFT) with  N=8 , where the   N  coefficients for the time domain   ⇒   d(0), ... , d(7)  can be taken from the graph.
  
Die dazugehörigen Spektralkoeffizienten  D(0), ... ,  D(7)  sind zu ermitteln, wobei für die Indizes  μ=0, ... , N–1  gilt:
+
The corresponding spectral coefficients  D(0), ... ,  D(7)  are to be determined, whereby for the indices  \mu = 0, ... , N–1  applies:
 
:$$D(\mu) = \frac{1}{N} \cdot \sum_{\nu = 0 }^{N-1}
 
:$$D(\mu) = \frac{1}{N} \cdot \sum_{\nu = 0 }^{N-1}
 
   d(\nu)\cdot  {w}^{\hspace{0.05cm}\nu \hspace{0.05cm} \cdot \hspace{0.05cm}\mu} \hspace{0.05cm}.$$
 
   d(\nu)\cdot  {w}^{\hspace{0.05cm}\nu \hspace{0.05cm} \cdot \hspace{0.05cm}\mu} \hspace{0.05cm}.$$
  
Bezeichnet man den Abstand zweier Abtastwerte im Zeitbereich mit  T_{\rm A}  und den entsprechenden Frequenzabstand zweier Linien mit  f_{\rm A}, so gilt folgender Zusammenhang:
+
If we denote the distance between two samples in the time domain by  T_{\rm A}  and the corresponding frequency distance between two lines by   f_{\rm A}, the following relationship applies:
 
:N \cdot f_{\rm A} \cdot T_{\rm A} = 1 \hspace{0.05cm}.
 
:N \cdot f_{\rm A} \cdot T_{\rm A} = 1 \hspace{0.05cm}.
  
Line 37: Line 37:
  
  
''Hinweise:''  
+
''Hints:''  
*Die Aufgabe gehört zum  Kapitel  [[Signal_Representation/Discrete_Fourier_Transform_(DFT)|Diskrete Fouriertransformation (DFT)]].
+
*This task belongs to the chapter  [[Signal_Representation/Discrete_Fourier_Transform_(DFT)|Discrete Fourier Transformation (DFT)]].
 
   
 
   
*Ihre Lösungen können Sie mit dem interaktiven Applet  [[Applets:Diskrete_Fouriertransformation_und_Inverse|Diskrete Fouriertransformation und Inverse]]  überrprüfen.
+
*You can check your solutions with the interactive applet  [[Applets:Diskrete_Fouriertransformation_und_Inverse|Discrete Fourier Transform and Inverse]] .
  
  
Line 48: Line 48:
  
  
===Fragebogen===
+
===Question===
  
 
<quiz display=simple>
 
<quiz display=simple>
Line 86: Line 86:
 
</quiz>
 
</quiz>
  
===Musterlösung===
+
===Solution===
 
{{ML-Kopf}}
 
{{ML-Kopf}}
 
'''(1)'''&nbsp;  Aus der Grafik ergeben sich mit&nbsp; A = 4 \ {\rm V}&nbsp; folgende Werte:
 
'''(1)'''&nbsp;  Aus der Grafik ergeben sich mit&nbsp; A = 4 \ {\rm V}&nbsp; folgende Werte:

Revision as of 21:14, 20 March 2021

Discretisation of a triangular pulse

Consider the sketched triangular momentum

x(t) = \left\{ \begin{array}{l} A \cdot \left( 1 - {|t|}/{T} \right ) \\ \hspace{0.25cm} 0 \\ \end{array} \right.\quad \begin{array}{*{10}c} {\rm{f\ddot{u}r}} \\ {\rm{f\ddot{u}r}} \\ \end{array}\begin{array}{*{10}c} {\left| \hspace{0.005cm} t\hspace{0.05cm} \right| \le T,} \\ {\left|\hspace{0.005cm} t \hspace{0.05cm} \right| > T.} \\ \end{array}

The signal parameters have the following values:

  • Amplitude  A = 4 \ \text{V},
  • equivalent pulse duration  \Delta t = T = 1 \, \text{ms}.


The spectrum  X(f)  is obtained by applying  the first Fourier Integral:

X(f) = A \cdot T \cdot {\rm si}^2(\pi f T)\hspace{0.05cm}.

The spectral function is now to be approximated by a Discrete Fourier Transform  (DFT) with  N = 8 , where the   N  coefficients for the time domain   ⇒   d(0), ... , d(7)  can be taken from the graph.

The corresponding spectral coefficients  D(0), ... ,  D(7)  are to be determined, whereby for the indices  \mu = 0, ... , N–1  applies:

D(\mu) = \frac{1}{N} \cdot \sum_{\nu = 0 }^{N-1} d(\nu)\cdot {w}^{\hspace{0.05cm}\nu \hspace{0.05cm} \cdot \hspace{0.05cm}\mu} \hspace{0.05cm}.

If we denote the distance between two samples in the time domain by  T_{\rm A}  and the corresponding frequency distance between two lines by   f_{\rm A}, the following relationship applies:

N \cdot f_{\rm A} \cdot T_{\rm A} = 1 \hspace{0.05cm}.




Hints:




Question

1

Geben Sie die Zeitkoeffizienten an. Wie groß sind  d(0)d(3)  und  d(6)?

d(0)\ = \

 \text{V}
d(3)\ = \

 \text{V}
d(6)\ = \

 \text{V}

2

Wie groß ist der Abstand  T_{\rm A}  zweier Zeitabtastwerte?

T_{\rm A}\ = \

 \text{ms}

3

Wie groß ist der Abstand  f_{\rm A}  zweier DFT–Frequenzabtastwerte?

f_{\rm A}\ = \

 \text{kHz}

4

Berechnen Sie den Koeffizienten  D(0)  und interpretieren Sie das Ergebnis.

D(0)\ = \

 \text{V}

5

Berechnen Sie den Koeffizienten  D(2)  und interpretieren Sie das Ergebnis, auch im Hinblick auf die Koeffizienten  D(4)  und  D(6).

D(2)\ = \

 \text{V}

6

Berechnen und interpretieren Sie den DFT–Koeffizienten  D(7).

D(7)\ = \

 \text{V}


Solution

(1)  Aus der Grafik ergeben sich mit  A = 4 \ {\rm V}  folgende Werte:

{d(0) = 4\,{\rm V}, \hspace{0.1cm}d(1) = d(7) = 3\,{\rm V}, \hspace{0.1cm} \hspace{0.1cm}d(2) = d(6) = 2\,{\rm V}, \hspace{0.1cm} \hspace{0.1cm}d(3) = d(5) = 1\,{\rm V}, \hspace{0.1cm} \hspace{0.1cm}d(4) = 0}\hspace{0.05cm}.
\Rightarrow \hspace{0.15 cm}\underline{d(0) = 4\,{\rm V}, \hspace{0.1cm}d(3) = 1\,{\rm V}, \hspace{0.1cm}d(6) = 2\,{\rm V}. \hspace{0.1cm}} \hspace{0.05cm}


(2)  Entsprechend der Grafik gilt  T_{\rm A} = T/4.

  • Mit  T = 1 \ \text{ms}  erhält man somit  \underline{T_{\rm A} = 0.25 \ \text{ms}}.


(3)  Für die Abstände der Abtastwerte im Zeit– und Frequenzbereich gilt:

N \cdot f_{\rm A} \cdot T_{\rm A} = 1 \hspace{0.3cm}\Rightarrow \hspace{0.3cm}f_{\rm A}= \frac{1}{ 8 \cdot 0.25\, {\rm ms}}\hspace{0.15 cm}\underline{= 0.5\, {\rm kHz}}\hspace{0.05cm}.


(4)  Mit  N = 8  und  \mu = 0  folgt aus der DFT–Gleichung:

D(0) = \frac{1}{8}\cdot \sum_{\nu = 0 }^{7} d(\nu) = \frac{1\,{\rm V}}{8}\cdot (4+3+2+1+0+1+2+3)\hspace{0.15 cm}\underline{= 2 \,{\rm V}}\hspace{0.05cm}.
  • Der DFT–Wert D(0) beschreibt den Spektralwert bei  f = 0, wobei folgender Zusammenhang gilt:
X(f=0) = \frac{D(0)}{f_{\rm A}}= \frac{ 2\,{\rm V}}{0.5\,{\rm kHz}}= 4 \cdot 10^{-3}\,{\rm V/Hz}\hspace{0.05cm}.
  • Dieser Wert stimmt mit dem theoretischen Wert   (A \cdot T)  überein.


(5)  Mit  N = 8  und  \mu = 2  erhält man:

D(2) = \frac{1}{8}\cdot \sum_{\nu = 0 }^{7} d(\nu)\cdot {\rm e}^{-{\rm j} \hspace{0.05cm}\cdot \hspace{0.05cm} (\pi /2) \hspace{0.05cm}\cdot \hspace{0.05cm} \nu} = \frac{1}{8}\cdot \sum_{\nu = 0 }^{7} d(\nu)\cdot (-{\rm j})^{\nu} \hspace{0.3cm} \Rightarrow \hspace{0.3cm} = \frac{1\,{\rm V}}{8}\cdot (4-3\cdot{\rm j}-2+{\rm j}-{\rm j}-2+3\cdot{\rm j})\hspace{0.15 cm}\underline{= 0}\hspace{0.05cm}.

Dieses Ergebnis hätte man auch ohne Rechnung vorhersagen können:

  • Die DFT-Koeffizienten  D(\mu)  sind gleichzeitig die Fourierkoeffizienten der im Abstand  T_{\rm P} = 2T  periodifizierten Funktion  x_{\rm Per}(t). Diese ist in der Grafik auf der Angabenseite gestrichelt eingezeichnet.
  • Aufgrund von Symmetrieeigenschaften sind aber alle geradzahligen Fourierkoeffizienten der Funktion  x_{\rm Per}(t)  gleich Null:   ⇒   D(4)\hspace{0.15cm}\underline{=0},   D(6)\hspace{0.15cm}\underline{=0}.


(6)  Der Koeffizient  D(7)  beschreibt die periodifizierte Spektralfunktion bei der Frequenz  f = 7 \cdot f_{\rm A}. Aufgrund der Periodizität und von Symmetrieeigenschaft gilt:

D(7) = D(-1) = D^{\star}(1) \hspace{0.05cm}.

Vorzugsweise berechnen wir diesen DFT–Koeffizienten:

D(1) = \frac{1}{8}\cdot \sum_{\nu = 0 }^{7} d(\nu)\cdot {\rm e}^{-{\rm j} \hspace{0.05cm}\cdot \hspace{0.05cm} (\pi /4) \hspace{0.05cm}\cdot \hspace{0.05cm} \nu} = \frac{1\,{\rm V}}{8}\cdot \left(4 +3\cdot \frac{1 - {\rm j}}{\sqrt{2}}-2\cdot {\rm j}+ \frac{-1 - {\rm j}}{\sqrt{2}}-{\rm j}+ \frac{-1 + {\rm j}}{\sqrt{2}}-{\rm j}+ 2\cdot {\rm j}+3\cdot \frac{1 - {\rm j}}{\sqrt{2}}\right)
\Rightarrow \; \; D(1) = \frac{2 + \sqrt{2}}{4} \approx 0.854{\rm V}\hspace{0.05cm}.

Da  D(1)  rein reell ist, gilt  D(7) = D(1) \; \underline{= 0.854 \ {\rm V}}.

Daraus ergeben sich für die zugehörigen Werte der kontinuierlichen Spektralfunktion:

X(f=-f_{\rm A}) = X(f=+f_{\rm A}) =\frac{D(1)}{f_{\rm A}}= 1.708 \cdot 10^{-3}\,{\rm V/Hz}\hspace{0.05cm}.
  • Wegen der impliziten periodischen Fortsetzung durch die DFT stimmt der so berechnete Wert mit dem tatsächlichen Wert  (4 \cdot A \cdot T/\pi^2 = 1.621 · 10^{-3}\text{ V/Hz}) nicht exakt überein.
  • Der relative Fehler beträgt ca.  5.3\%.