Difference between revisions of "Signal Representation/Analytical Signal and its Spectral Function"
Line 183: | Line 183: | ||
[[Aufgaben:Exercise 4.3Z: Hilbert Transformator|Exercise 4.3Z: Hilbert Transformator]] | [[Aufgaben:Exercise 4.3Z: Hilbert Transformator|Exercise 4.3Z: Hilbert Transformator]] | ||
− | [[Aufgaben:Exercise 4. | + | [[Aufgaben:Exercise 4.4: Vector Diagram for DSB-AM|Exercise 4.4: Vector Diagram for DSB-AM]] |
[[Aufgaben:Exercise 4.4Z: Vector Diagram for DSB-AM|Exercise 4.4Z: Vector Diagram for DSB-AM]] | [[Aufgaben:Exercise 4.4Z: Vector Diagram for DSB-AM|Exercise 4.4Z: Vector Diagram for DSB-AM]] |
Revision as of 20:45, 24 November 2020
Contents
Definition in the Frequency Domain
We consider a real bandpass-like signal $x(t)$ with the corresponding bandpass spectrum $X(f)$, which has an even real and an odd imaginary part with respect to the frequency zero point. It is assumed that the carrier frequency $f_{\rm T}$ is much larger than the bandwidth of the bandpass signal $x(t)$ .
$\text{Definition:}$ The time function belonging to the physical signal $x(t)$ analytical signal $x_+(t)$ is that time function, whose spectrum fulfills the following property
- $$X_+(f)=\big[1+{\rm sign}(f)\big] \cdot X(f) = \left\{ {2 \cdot X(f) \; \hspace{0.2cm}\rm for\hspace{0.2cm} {\it f} > 0, \atop {\,\,\,\, \rm 0 \; \hspace{0.9cm}\rm for\hspace{0.2cm} {\it f} < 0.} }\right.$$
The so called „signum function” is for positive values of $f$ equal to $+1$ and for negative $f$-values equal to $-1$.
- The (double sided) limit value returns $\sign(0) = 0$.
- The index "+" should make clear that $X_+(f)$ has only parts at positive frequencies.
From the graphic you can see the calculation rule for $X_+(f)$:
The actual bandpass spectrum $X(f)$ will
- doubled at the positive frequencies, and
- set to zero at the negative frequencies.
$\text{Example 1:}$
The graphic
- at left shows the (complex) spectrum $X(f)$ of the bandpass signal
- $$x(t) = 4\hspace{0.05cm}{\rm V} \cdot {\cos} ( 2 \pi f_{\rm u} \hspace{0.03cm}t) + 6\hspace{0.05cm}{\rm V} \cdot {\sin} ( 2 \pi f_{\rm o} \hspace{0.03cm}t).$$
- and on the right the (complex) spectrum of the analytical signal $x_{+}(t)$.
Allgemeingültige Berechnungsvorschrift im Zeitbereich
Now we will take a closer look at the spectrum $X_+(f)$ of the analytical signal and divide it into a with respect to $f = 0$ even part $X_{\rm +g}(f)$ and an odd part $X_{\rm +u}(f)$ :
- $$X_+(f) = X_{\rm +g}(f) + X_{\rm +u}(f).$$
All these spectra are generally complex.
If one considers the nbsp; Mapping Theorem of the Fourier transform, then the following statements are possible on the basis of the graphic:
- Der gerade Anteil $X_{\rm +g}(f)$ von $X_{+}(f)$ führt nach der Fouriertransformation zu einem reellen Zeitsignal, der ungerade Anteil $X_{\rm +u}(f)$ zu einem imaginären.
- Es ist offensichtlich, dass $X_{\rm +g}(f)$ gleich dem tatsächlichen Fourierspektrum $X(f)$ und damit der Realteil von $x_{\rm +g}(t)$ gleich dem vorgegebenen Signal $x(t)$ mit Bandpasseigenschaften ist.
- Bezeichnen wir den Imaginärteil mit $y(t)$, so lautet das analytische Signal:
- $$x_+(t)= x(t) + {\rm j} \cdot y(t) .$$
- Nach den allgemein gültigen Gesetzen der Fouriertransformation entsprechend dem Zuordnungssatz gilt somit für die Spektralfunktion des Imaginärteils:
- $${\rm j} \cdot Y(f) = X_{\rm +u}(f)= {\rm sign}(f) \cdot X(f) \hspace{0.3cm}\Rightarrow\hspace{0.3cm}Y(f) = \frac{{\rm sign}(f)}{ {\rm j}}\cdot X(f).$$
- Transformiert man diese Gleichung in den Zeitbereich, so wird aus der Multiplikation die Faltungsoperation, und man erhält:
- $$y(t) = \frac{1}{ {\rm \pi} t} \hspace{0.05cm}\star \hspace{0.05cm}x(t) = \frac{1}{ {\rm \pi}} \cdot \hspace{0.03cm}\int_{-\infty}^{+\infty}\frac{x(\tau)}{ {t - \tau}}\hspace{0.15cm} {\rm d}\tau.$$
Darstellung mit der Hilberttransformation
An dieser Stelle ist es erforderlich, kurz auf eine weitere Spektraltransformation einzugehen, die im Buch Lineare zeitinvariante Systeme noch eingehend behandelt wird.
$\text{Definition:}$ Für die Hilberttransformierte $ {\rm H}\left\{x(t)\right\}$ einer Zeitfunktion $x(t)$ gilt:
- $$y(t) = {\rm H}\left\{x(t)\right\} = \frac{1}{ {\rm \pi} } \cdot \hspace{0.03cm}\int_{-\infty}^{+\infty}\frac{x(\tau)}{ {t - \tau} }\hspace{0.15cm} {\rm d}\tau.$$
- Dieses bestimmte Integral ist nicht auf einfache, herkömmliche Art lösbar, sondern muss mit Hilfe des Cauchy–Hauptwertsatzes ausgewertet werden.
- Entsprechend gilt im Frequenzbereich:
- $$Y(f) = - {\rm j} \cdot {\rm sign}(f) \cdot X(f) \hspace{0.05cm} .$$
Das Ergebnis der letzten Seite lässt sich mit dieser Definition wie folgt zusammenfassen:
- Man erhält aus dem realen, physikalischen Bandpass–Signal $x(t)$ das analytische Signal $x_+(t)$, indem man zu $x(t)$ einen Imaginärteil entsprechend der Hilberttransformierten hinzufügt:
- $$x_+(t) = x(t)+{\rm j} \cdot {\rm H}\left\{x(t)\right\} .$$
- Die Hilberttransformierte $\text{H}\{x(t)\}$ verschwindet nur für den Fall $x(t) = \rm const.$ ⇒ Gleichsignal Bei allen anderen Signalformen ist das analytische Signal $x_+(t)$ somit stets komplex.
- Aus dem analytischen Signal $x_+(t)$ kann das reale Bandpass–Signal in einfacher Weise durch Realteilbildung ermittelt werden:
- $$x(t) = {\rm Re}\left\{x_+(t)\right\} .$$
$\text{Beispiel 2:}$ Das Prinzip der Hilbert–Transformation wird durch die folgende Grafik nochmals verdeutlicht:
- Nach der linken Darstellung $\rm (A)$ kommt man vom physikalischen Signal $x(t)$ zum analytischen Signal $x_+(t)$, indem man einen Imaginärteil ${\rm j} \cdot y(t)$ hinzufügt.
- Hierbei ist $y(t) = {\rm H}\left\{x(t)\right\}$ eine reelle Zeitfunktion, die sich am einfachsten im Spektralbereich durch die Multiplikation des Spektrums $X(f)$ mit $- {\rm j} \cdot \sign(f)$ angeben lässt.
Die rechte Darstellung $\rm (B)$ ist äquivalent zu $\rm (A)$:
- Nun gilt $x_+(t) = x(t) + z(t)$ mit der rein imaginären Funktion $z(t)$.
- Ein Vergleich der beiden Bilder zeigt, dass tatsächlich $z(t) = {\rm j} \cdot y(t)$ ist.
Zeigerdiagrammdarstellung der harmonischen Schwingung
Die Spektralfunktion $X(f)$ einer harmonischen Schwingung $x(t) = A \cdot \text{cos}(2\pi f_{\rm T}t - \varphi)$ besteht bekanntlich aus zwei Diracfunktionen bei den Frequenzen
- $+f_{\rm T}$ mit dem komplexen Gewicht $A/2 \cdot \text{e}^{-\text{j}\hspace{0.05cm}\varphi}$,
- $-f_{\rm T}$ mit dem komplexen Gewicht $A/2 \cdot \text{e}^{+\text{j}\hspace{0.05cm}\varphi}$.
Somit lautet das Spektrum des analytischen Signals $($also ohne die Diracfunktion bei der Frequenz $f =-f_{\rm T})$:
- $$X_+(f) = A \cdot {\rm e}^{-{\rm j} \hspace{0.05cm}\varphi}\cdot\delta (f - f_{\rm T}) .$$
Die dazugehörige Zeitfunktion erhält man durch Anwendung des Verschiebungssatzes:
- $$x_+(t) = A \cdot {\rm e}^{\hspace{0.05cm} {\rm j}\hspace{0.05cm}\cdot\hspace{0.05cm}( 2 \pi f_{\rm T} t \hspace{0.05cm}-\hspace{0.05cm} \varphi)}.$$
Diese Gleichung beschreibt einen mit konstanter Winkelgeschwindigkeit $\omega_{\rm T} = 2\pi f_{\rm T}$ drehenden Zeiger.
$\text{Beispiel 3:}$ Aus Darstellungsgründen ist in der folgenden Grafik das Koordinatensystem entgegen der üblichen Darstellung um $90^\circ$ nach links gedreht (Realteil nach oben, Imaginärteil nach links).
Anhand dieser Grafik sind folgende Aussagen möglich:
- Zum Startzeitpunkt $t = 0$ liegt der Zeiger der Länge $A$ (Signalamplitude) mit dem Winkel $-\varphi$ in der komplexen Ebene. Im gezeichneten Beispiel gilt $\varphi = 45^\circ$.
- Für Zeiten $t > 0$ dreht der Zeiger mit konstanter Winkelgeschwindigkeit (Kreisfrequenz) $\omega_{\rm T}$ in mathematisch positiver Richtung, das heißt entgegen dem Uhrzeigersinn.
- Die Spitze des Zeigers liegt somit stets auf einem Kreis mit Radius $A$ und benötigt für eine Umdrehung genau die Zeit $T_0$, also die Periodendauer der harmonischen Schwingung $x(t)$.
- Die Projektion des analytischen Signals $x_+(t)$ auf die reelle Achse, durch rote Punkte markiert, liefert die Augenblickswerte von $x(t)$.
Zeigerdiagramm einer Summe harmonischer Schwingungen
Für die weitere Beschreibung gehen wir für das analytische Signal von folgendem Spektrum aus:
- $$X_+(f) = \sum_{i=1}^{I}A_i \cdot {\rm e}^{-{\rm j}\hspace{0.05cm}\cdot\hspace{0.05cm} \varphi_i}\cdot\delta (f - f_{i}) .$$
Das linke Bild zeigt ein solches Spektrum für das Beispiel $I = 3$. Wählt man $I$ relativ groß und den Abstand zwischen benachbarten Spektrallinien entsprechend klein, so können mit obiger Gleichung auch (frequenz–) kontinuierliche Spektralfunktionen $X_+(f)$ angenähert werden.
Im rechten Bild ist die dazugehörige Zeitfunktion angedeutet. Diese lautet allgemein:
- $$x_+(t) = \sum_{i=1}^{I}A_i \cdot {\rm e}^{ {\rm j}\hspace{0.05cm}\cdot\hspace{0.05cm}(\omega_i \hspace{0.05cm}\cdot\hspace{0.05cm} t \hspace{0.05cm}-\hspace{0.05cm} \varphi_i)}.$$
Zu dieser Grafik anzumerken:
- Die Skizze zeigt die Ausgangslage der Zeiger zum Startzeitpunkt $t = 0$ entsprechend den Amplituden $A_i$ und den Phasenlagen $\varphi_i$.
- Die Spitze des resultierenden Zeigerverbundes ist durch das violette Kreuz markiert. Man erhält durch vektorielle Addition der drei Einzelzeiger für den Zeitpunkt $t = 0$:
- $$x_+(t= 0) = \big [1 \cdot \cos(60^\circ) - 1 \cdot {\rm j} \cdot \sin(60^\circ) \big ]+ 2 \cdot \cos(0^\circ)+1 \cdot \cos(180^\circ) = 1.500 - {\rm j} \cdot 0.866.$$
- Für Zeiten $t > 0$ drehen die drei Zeiger mit unterschiedlichen Winkelgeschwindigkeiten $\omega_i = 2\pi f_i$. Der rote Zeiger dreht schneller als der grüne, aber langsamer als der blaue Zeiger.
- Da alle Zeiger entgegen dem Uhrzeigersinn drehen, wird sich auch der resultierende Zeiger $x_+(t)$ tendenziell in diese Richtung bewegen. Zum Zeitpunkt $t = 1\,µ\text {s}$ liegt die Spitze des resultierenen Zeigers für die gegebenen Parameterwerte bei
- $$ \begin{align*}x_+(t = 1 {\rm \hspace{0.05cm}µ s}) & = 1 \cdot {\rm e}^{-{\rm j}\hspace{0.05cm}\cdot \hspace{0.05cm}60^\circ}\cdot {\rm e}^{{\rm j}\hspace{0.05cm}2 \pi \hspace{0.05cm}\cdot \hspace{0.1cm}40 \hspace{0.05cm} \cdot \hspace{0.1cm} 0.001} + 2\cdot {\rm e}^{{\rm j}\hspace{0.05cm}2 \pi \hspace{0.05cm}\cdot \hspace{0.1cm}50 \hspace{0.05cm} \cdot \hspace{0.1cm} 0.001}- 1\cdot {\rm e}^{{\rm j}\hspace{0.05cm}2 \pi \hspace{0.05cm}\cdot \hspace{0.1cm}60 \hspace{0.05cm} \cdot \hspace{0.1cm} 0.001} = \\ & = 1 \cdot {\rm e}^{-{\rm j}\hspace{0.05cm}\cdot \hspace{0.05cm}45.6^\circ} + 2\cdot {\rm e}^{{\rm j}\hspace{0.05cm}\cdot \hspace{0.05cm}18^\circ}- 1\cdot {\rm e}^{{\rm j}\hspace{0.05cm}\cdot \hspace{0.05cm}21.6^\circ} \approx 1.673- {\rm j} \cdot 0.464.\end{align*}$$
- Die resultierende Zeigerspitze liegt nun aber nicht wie bei einer einzigen Schwingung auf einem Kreis, sondern es entsteht eine komplizierte geometrische Figur.
Das interaktive Applet Physikalisches Signal & Analytisches Signal verdeutlicht $x_+(t)$ für die Summe dreier harmonischer Schwingungen.
Aufgaben zum Kapitel
Exercise 4.3: Vector Diagram Representation
Exercise 4.3Z: Hilbert Transformator
Exercise 4.4: Vector Diagram for DSB-AM
Exercise 4.4Z: Vector Diagram for DSB-AM