Difference between revisions of "Aufgaben:Exercise 5.4Z: On the Hanning Window"
(24 intermediate revisions by 5 users not shown) | |||
Line 1: | Line 1: | ||
− | {{quiz-Header|Buchseite= | + | {{quiz-Header|Buchseite=Signal_Representation/Spectrum_Analysis |
}} | }} | ||
− | [[File:P_ID1168__Sig_Z_5_4.png|250px|right|Hanning | + | [[File:P_ID1168__Sig_Z_5_4.png|250px|right|frame|Characterisation of the Hanning window]] |
− | In | + | In this task, important properties of the frequently used Hanning window are to be derived. The continuous-time representation in the interval from $-T_{\rm P}/2$ to $+T_{\rm P}/2$ is here as follows: |
− | $$w(t)= {\rm cos}^2(\pi \cdot | + | :$$w(t)= {\rm cos}^2(\pi \cdot |
− | {t}/{T_{\rm P}})= 0.5\cdot \ | + | {t}/{T_{\rm P}})= 0.5\cdot \big(1 + {\rm cos}(2\pi \cdot |
− | {t}/{T_{\rm P}}) \ | + | {t}/{T_{\rm P}}) \big ) |
\hspace{0.05cm}.$$ | \hspace{0.05cm}.$$ | ||
− | + | Outside the symmetric time domain of duration $T_{\rm P}$: $w(t) \equiv 0$. | |
− | + | The upper graph shows the discrete-time representation $w(\nu) = w({\nu} \cdot T_{\rm A})$, where $T_{\rm A}$ is smaller than $T_{\rm P}$ by a factor $N = 32$. The definition range of the discrete time variable $ν$ extends from $-16$ to $+15$. | |
− | In | + | In the lower graph, the Fourier transform $W(f)$ of the continuous-time window function $w(t)$ is shown logarithmically. The abscissa is normalised to $f_{\rm A} = 1/T_{\rm P}$. One can see: |
− | * | + | *The equidistant values $W({\mu} \cdot f_{\rm A})$ are zero except for $μ = 0$ and $μ = ±1$. |
− | * | + | *The main lobe thus extends to the frequency range $|f| ≤ 2 · f_{\rm A}$. |
− | *$W(f)$ | + | *$W(f)$ is largest in magnitude outside the main lobe for $f = ±2.5 · f_{\rm A}$. |
− | * | + | *Thus, the minimum distance between the main and the side lobes applies here: |
− | $$ | + | :$$D_{\text{min (main lobe / side lobes) } } = 20 \cdot {\rm lg}\hspace{0.15cm} |
\frac{|W(0)|}{|W(2.5 \cdot f_{\rm A})|} \hspace{0.15cm}{\rm (in}\hspace{0.1cm}{\rm dB)}\hspace{0.05cm}.$$ | \frac{|W(0)|}{|W(2.5 \cdot f_{\rm A})|} \hspace{0.15cm}{\rm (in}\hspace{0.1cm}{\rm dB)}\hspace{0.05cm}.$$ | ||
− | |||
− | |||
− | |||
− | === | + | |
+ | |||
+ | |||
+ | |||
+ | |||
+ | ''Hints:'' | ||
+ | *This task belongs to the chapter [[Signal_Representation/Spectrum_Analysis|Spectral Analysis]]. | ||
+ | *Note that the frequency resolution $f_{\rm A}$ is equal to the reciprocal of the adjustable parameter $T_{\rm P}$. | ||
+ | *In this exercise we use the ${\rm si}$–function: ${\rm si}(x) = \sin(x)/x$. | ||
+ | |||
+ | |||
+ | |||
+ | |||
+ | |||
+ | ===Questions=== | ||
<quiz display=simple> | <quiz display=simple> | ||
− | { | + | {Give the discrete-time coefficients $w(ν)$ of the Hanning window analytically. What are the numerical values for $ν = 0$, $ν = 1$ and $ν = -\hspace{0.05cm}8$? |
|type="{}"} | |type="{}"} | ||
− | $w(ν = 0) \ =$ { 1 1% } | + | $w(ν = 0) \hspace{0.37cm} = \ $ { 1 1% } |
− | $w(ν = 1) \ =$ { 0.99 1% } | + | $w(ν = 1) \hspace{0.37cm} = \ $ { 0.99 1% } |
− | $w(ν = -\hspace{0. | + | $w(ν = -8) \hspace{0.03cm} = \ $ { 0.5 1% } |
− | { | + | {Calculate the spectral function $W(f)$ in general. Which of the following statements are correct? |
|type="[]"} | |type="[]"} | ||
− | - $W(f)$ | + | - $W(f)$ yields complex results for special frequency values. |
− | + $W(f)$ | + | + $W(f)$ is even with respect to $f$ , i.e. $W(-f) = W(+f)$. |
− | + | + | + The spectral value $W(f = 0)=0.5/f_{\rm A}$ and thus real. |
− | { | + | {What are $W(f = ±f_{\rm A})$ and the $\text{6dB}$ bandwidth normalised to $f_{\rm A}$? |
|type="{}"} | |type="{}"} | ||
− | $B_{\rm 6\hspace{0.05cm}dB}\hspace{-0.05cm}/\hspace{-0.05cm}f_{\rm A} \ =$ { 2 1% } | + | $W(±f_{\rm A}) \hspace{0.15cm} = \ $ { 0.25 1% } $\ \cdot \ 1/f_{\rm A}$ |
+ | $B_{\rm 6\hspace{0.05cm}dB}\hspace{-0.05cm}/\hspace{-0.05cm}f_{\rm A} \hspace{0.2cm} = \ $ { 2 1% } | ||
− | { | + | {What is the minimum distance between the main lobe and the side lobes. |
|type="{}"} | |type="{}"} | ||
− | $ | + | $D_{\text{min (main lobe / side lobes) } } \ = \ $ { 32.3 1% } $\ \rm dB$ |
</quiz> | </quiz> | ||
− | === | + | ===Solution=== |
{{ML-Kopf}} | {{ML-Kopf}} | ||
− | '''1 | + | '''(1)''' After trigonometric transformation, the result for the continuous-time window function is: |
− | + | :$$w(t) = {\rm cos}^2(\pi \cdot | |
− | + | {t}/{T_{\rm P}}) = 0.5+ 0.5\cdot {\rm | |
− | + | cos}(2\pi \cdot {t}/{T_{\rm P}})\hspace{0.05cm}.$$ | |
+ | *After time discretisation with $ν = t/T_{\rm A}$ and $T_{\rm P}/T_{\rm A} = N = 32$, one obtains for the discrete-time window: | ||
+ | :$$w(\nu) = w(\nu \cdot T_{\rm A}) = 0.5+ 0.5\cdot {\rm | ||
+ | cos}(2\pi \cdot {\nu}/{N})\hspace{0.8cm} | ||
+ | \Rightarrow \hspace{0.3cm}w(\nu = 0) \hspace{0.15 cm}\underline{= 1},$$ | ||
+ | :$$w(\nu = 1) = 0.5+ 0.5\cdot {\rm cos}( | ||
+ | \frac{\pi}{16})\hspace{0.15 cm}\underline{ = 0.99}, $$ | ||
+ | :$$w(\nu = -8)=0.5+ | ||
+ | 0.5\cdot {\rm cos}( \frac{-\pi}{2}) \hspace{0.15 cm}\underline{= | ||
+ | 0.5}\hspace{0.05cm}.$$ | ||
− | |||
+ | '''(2)''' The correct <u>solutions are 2 and 3:</u>: | ||
+ | *The periodic continuation of $w(t)$ corresponding to the period $T_{\rm P}$ yields a (periodic) signal with a DC and a cosine component. | ||
+ | *From this follows with $f_{\rm A} = 1/T_{\rm P}$: | ||
+ | :$${\rm P}\{w(t)\} = 0.5+0.5\cdot {\rm | ||
+ | cos}(2\pi \cdot f_{\rm A} \cdot t) | ||
+ | \hspace{0.2cm}\circ\!\!-\!\!\!-\!\!\!-\!\!\bullet\, | ||
+ | \hspace{0.2cm}0.5\cdot {\rm \delta}(f) + 0.25\cdot {\rm | ||
+ | \delta}(f \pm f_{\rm A}))\hspace{0.05cm}.$$ | ||
+ | *The time-limited signal $w(t)$ is obtained from ${\rm P}\{w(t)\}$ by multiplication with a rectangle of amplitude $1$ and duration $T_{\rm P}$. | ||
+ | *Its spectrum $W(f)$ is thus obtained from the convolution of the above spectral function with the function $T_{\rm P} · {\rm si}(π \cdot f \cdot T_{\rm P}) = 1/f_{\rm A} · {\rm si}(π \cdot f/f_{\rm A})$: | ||
+ | :$$w(t) | ||
+ | \circ\!\!-\!\!\!-\!\!\!-\!\!\bullet\, | ||
+ | W(f) = \frac{0.5}{f_{\rm A}}\cdot {\rm si}( \frac{\pi f}{f_{\rm | ||
+ | A}})+ \frac{0.25}{f_{\rm A}}\cdot {\rm si}(\pi \cdot | ||
+ | \frac{f-f_{\rm A}}{f_{\rm A}})+ \frac{0.25}{f_{\rm A}}\cdot {\rm | ||
+ | si}(\pi \cdot \frac{f+f_{\rm A}}{f_{\rm A}})\hspace{0.05cm}.$$ | ||
+ | *This spectral function is even and also real for all frequencies $f$ . The spectral value at frequency $f = 0$ gives the window area: | ||
+ | :$$W(f=0) = | ||
+ | \frac{0.5}{f_{\rm A}}= | ||
+ | \int_{-\infty}^{+\infty}w(t)\hspace{0.05cm}{\rm | ||
+ | d}t\hspace{0.05cm}.$$ | ||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | + | '''(3)''' From the result of sub-task '''(2)''' it also follows: | |
+ | :$$W(f = ±f_{\rm A}) = W(0)/2\hspace{0.15cm}\underline{ = 0.25} \cdot 1/{f_{\rm A}}.$$ | ||
+ | *Due to the monotonic course in the range $|f| < f_{\rm A}$ , the magnitude function $|W(f)|$ has dropped to half of the maximum for the first time exactly at $± f_{\rm A}$ . | ||
+ | *Thus, $B_{\rm 6\hspace{0.05cm}dB}\hspace{-0.05cm}/\hspace{-0.05cm}f_{\rm A} \;\underline{=2}$. | ||
+ | |||
+ | |||
− | [ | + | '''(4)''' The largest spectral magnitude outside the main lobe occurs at $f = ±2.5 f_{\rm A}$. With the result of subtask '''(2)''' holds: |
+ | :$$W(f = 2.5 \cdot f_{\rm A}) = \frac{0.5}{f_{\rm A}}\cdot {\rm si}(2.5 \pi ) | ||
+ | +\frac{0.25}{f_{\rm A}}\cdot {\rm si}(1.5 \pi )+\frac{0.25}{f_{\rm A}}\cdot {\rm si}(3.5 \pi )= \frac{0.25}{\pi \cdot f_{\rm A}}\left[ \frac{2}{2.5}-\frac{1}{1.5}-\frac{1}{3.5}\right] \approx -\frac{0.0121}{ f_{\rm A}}\hspace{0.05cm}.$$ | ||
+ | *This gives the minimum distance between the main lobe and the side lobes: | ||
+ | :$$D_{\text{min (main lobe / side lobes) } } = 20 \cdot {\rm lg}\hspace{0.15cm} | ||
+ | \frac{|W(0)|}{|W(2.5 \cdot f_{\rm A})|} = 20 \cdot {\rm lg}\hspace{0.15cm} | ||
+ | \frac{0.5}{0.0121}\hspace{0.15 cm}\underline{\approx 32.3\,\,{\rm dB}}\hspace{0.05cm}.$$ | ||
− | |||
{{ML-Fuß}} | {{ML-Fuß}} | ||
__NOEDITSECTION__ | __NOEDITSECTION__ | ||
− | [[Category: | + | [[Category:Signal Representation: Exercises|^5.4 Spectrum Analysis^]] |
Latest revision as of 04:18, 18 September 2022
In this task, important properties of the frequently used Hanning window are to be derived. The continuous-time representation in the interval from $-T_{\rm P}/2$ to $+T_{\rm P}/2$ is here as follows:
- $$w(t)= {\rm cos}^2(\pi \cdot {t}/{T_{\rm P}})= 0.5\cdot \big(1 + {\rm cos}(2\pi \cdot {t}/{T_{\rm P}}) \big ) \hspace{0.05cm}.$$
Outside the symmetric time domain of duration $T_{\rm P}$: $w(t) \equiv 0$.
The upper graph shows the discrete-time representation $w(\nu) = w({\nu} \cdot T_{\rm A})$, where $T_{\rm A}$ is smaller than $T_{\rm P}$ by a factor $N = 32$. The definition range of the discrete time variable $ν$ extends from $-16$ to $+15$.
In the lower graph, the Fourier transform $W(f)$ of the continuous-time window function $w(t)$ is shown logarithmically. The abscissa is normalised to $f_{\rm A} = 1/T_{\rm P}$. One can see:
- The equidistant values $W({\mu} \cdot f_{\rm A})$ are zero except for $μ = 0$ and $μ = ±1$.
- The main lobe thus extends to the frequency range $|f| ≤ 2 · f_{\rm A}$.
- $W(f)$ is largest in magnitude outside the main lobe for $f = ±2.5 · f_{\rm A}$.
- Thus, the minimum distance between the main and the side lobes applies here:
- $$D_{\text{min (main lobe / side lobes) } } = 20 \cdot {\rm lg}\hspace{0.15cm} \frac{|W(0)|}{|W(2.5 \cdot f_{\rm A})|} \hspace{0.15cm}{\rm (in}\hspace{0.1cm}{\rm dB)}\hspace{0.05cm}.$$
Hints:
- This task belongs to the chapter Spectral Analysis.
- Note that the frequency resolution $f_{\rm A}$ is equal to the reciprocal of the adjustable parameter $T_{\rm P}$.
- In this exercise we use the ${\rm si}$–function: ${\rm si}(x) = \sin(x)/x$.
Questions
Solution
- $$w(t) = {\rm cos}^2(\pi \cdot {t}/{T_{\rm P}}) = 0.5+ 0.5\cdot {\rm cos}(2\pi \cdot {t}/{T_{\rm P}})\hspace{0.05cm}.$$
- After time discretisation with $ν = t/T_{\rm A}$ and $T_{\rm P}/T_{\rm A} = N = 32$, one obtains for the discrete-time window:
- $$w(\nu) = w(\nu \cdot T_{\rm A}) = 0.5+ 0.5\cdot {\rm cos}(2\pi \cdot {\nu}/{N})\hspace{0.8cm} \Rightarrow \hspace{0.3cm}w(\nu = 0) \hspace{0.15 cm}\underline{= 1},$$
- $$w(\nu = 1) = 0.5+ 0.5\cdot {\rm cos}( \frac{\pi}{16})\hspace{0.15 cm}\underline{ = 0.99}, $$
- $$w(\nu = -8)=0.5+ 0.5\cdot {\rm cos}( \frac{-\pi}{2}) \hspace{0.15 cm}\underline{= 0.5}\hspace{0.05cm}.$$
(2) The correct solutions are 2 and 3::
- The periodic continuation of $w(t)$ corresponding to the period $T_{\rm P}$ yields a (periodic) signal with a DC and a cosine component.
- From this follows with $f_{\rm A} = 1/T_{\rm P}$:
- $${\rm P}\{w(t)\} = 0.5+0.5\cdot {\rm cos}(2\pi \cdot f_{\rm A} \cdot t) \hspace{0.2cm}\circ\!\!-\!\!\!-\!\!\!-\!\!\bullet\, \hspace{0.2cm}0.5\cdot {\rm \delta}(f) + 0.25\cdot {\rm \delta}(f \pm f_{\rm A}))\hspace{0.05cm}.$$
- The time-limited signal $w(t)$ is obtained from ${\rm P}\{w(t)\}$ by multiplication with a rectangle of amplitude $1$ and duration $T_{\rm P}$.
- Its spectrum $W(f)$ is thus obtained from the convolution of the above spectral function with the function $T_{\rm P} · {\rm si}(π \cdot f \cdot T_{\rm P}) = 1/f_{\rm A} · {\rm si}(π \cdot f/f_{\rm A})$:
- $$w(t) \circ\!\!-\!\!\!-\!\!\!-\!\!\bullet\, W(f) = \frac{0.5}{f_{\rm A}}\cdot {\rm si}( \frac{\pi f}{f_{\rm A}})+ \frac{0.25}{f_{\rm A}}\cdot {\rm si}(\pi \cdot \frac{f-f_{\rm A}}{f_{\rm A}})+ \frac{0.25}{f_{\rm A}}\cdot {\rm si}(\pi \cdot \frac{f+f_{\rm A}}{f_{\rm A}})\hspace{0.05cm}.$$
- This spectral function is even and also real for all frequencies $f$ . The spectral value at frequency $f = 0$ gives the window area:
- $$W(f=0) = \frac{0.5}{f_{\rm A}}= \int_{-\infty}^{+\infty}w(t)\hspace{0.05cm}{\rm d}t\hspace{0.05cm}.$$
(3) From the result of sub-task (2) it also follows:
- $$W(f = ±f_{\rm A}) = W(0)/2\hspace{0.15cm}\underline{ = 0.25} \cdot 1/{f_{\rm A}}.$$
- Due to the monotonic course in the range $|f| < f_{\rm A}$ , the magnitude function $|W(f)|$ has dropped to half of the maximum for the first time exactly at $± f_{\rm A}$ .
- Thus, $B_{\rm 6\hspace{0.05cm}dB}\hspace{-0.05cm}/\hspace{-0.05cm}f_{\rm A} \;\underline{=2}$.
(4) The largest spectral magnitude outside the main lobe occurs at $f = ±2.5 f_{\rm A}$. With the result of subtask (2) holds:
- $$W(f = 2.5 \cdot f_{\rm A}) = \frac{0.5}{f_{\rm A}}\cdot {\rm si}(2.5 \pi ) +\frac{0.25}{f_{\rm A}}\cdot {\rm si}(1.5 \pi )+\frac{0.25}{f_{\rm A}}\cdot {\rm si}(3.5 \pi )= \frac{0.25}{\pi \cdot f_{\rm A}}\left[ \frac{2}{2.5}-\frac{1}{1.5}-\frac{1}{3.5}\right] \approx -\frac{0.0121}{ f_{\rm A}}\hspace{0.05cm}.$$
- This gives the minimum distance between the main lobe and the side lobes:
- $$D_{\text{min (main lobe / side lobes) } } = 20 \cdot {\rm lg}\hspace{0.15cm} \frac{|W(0)|}{|W(2.5 \cdot f_{\rm A})|} = 20 \cdot {\rm lg}\hspace{0.15cm} \frac{0.5}{0.0121}\hspace{0.15 cm}\underline{\approx 32.3\,\,{\rm dB}}\hspace{0.05cm}.$$