Difference between revisions of "Aufgaben:Exercise 5.2Z: DFT of a Triangular Pulse"
m (Oezdemir moved page Aufgabe 5.2Z: DFT eines Dreieckimpulses to Exercise 5.2Z: DFT of a Triangular Pulse) |
|||
Line 3: | Line 3: | ||
}} | }} | ||
− | [[File:P_ID1140__Sig_Z_5_2.png|right|frame| | + | [[File:P_ID1140__Sig_Z_5_2.png|right|frame|Discretisation of a triangular pulse]] |
− | + | 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}$$ | ||
− | + | The signal parameters have the following values: | |
* Amplitude A=4 V, | * Amplitude A=4 V, | ||
− | * | + | * equivalent pulse duration Δt=T=1ms. |
− | + | 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)=A⋅T⋅si2(πfT). | :X(f)=A⋅T⋅si2(πfT). | ||
− | + | 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(\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}.$$ | ||
− | + | 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: | ||
− | '' | + | ''Hints:'' |
− | * | + | *This task belongs to the chapter [[Signal_Representation/Discrete_Fourier_Transform_(DFT)|Discrete Fourier Transformation (DFT)]]. |
− | * | + | *You can check your solutions with the interactive applet [[Applets:Diskrete_Fouriertransformation_und_Inverse|Discrete Fourier Transform and Inverse]] . |
Line 48: | Line 48: | ||
− | === | + | ===Question=== |
<quiz display=simple> | <quiz display=simple> | ||
Line 86: | Line 86: | ||
</quiz> | </quiz> | ||
− | === | + | ===Solution=== |
{{ML-Kopf}} | {{ML-Kopf}} | ||
'''(1)''' Aus der Grafik ergeben sich mit A = 4 \ {\rm V} folgende Werte: | '''(1)''' Aus der Grafik ergeben sich mit A = 4 \ {\rm V} folgende Werte: |
Revision as of 21:14, 20 March 2021
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:
- This task belongs to the chapter Discrete Fourier Transformation (DFT).
- You can check your solutions with the interactive applet Discrete Fourier Transform and Inverse .
Question
Solution
- {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\%.