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

Analytical Signal and its Spectral Function

From LNTwww

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  fT  is much larger than the bandwidth of the bandpass signal  x(t) .

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

Analytical Signal in the Frequency Domain
X+(f)=[1+sign(f)]X(f)={2X(f)forf>0,0forf<0.

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.


Example of a Spectrum of an Analytical Signal

Example 1: 

The graphic

  • at left shows the (complex) spectrum  X(f)  of the bandpass signal
x(t)=4Vcos(2πfut)+6Vsin(2πfot).
  • and on the right the (complex) spectrum of the analytical signal  x+(t).


Allgemeingültige Berechnungsvorschrift im Zeitbereich


To Derive the Analytical Signal

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+g(f)  and an odd part  X+u(f) :

X+(f)=X+g(f)+X+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+g(f)  von  X+(f)  führt nach der Fouriertransformation zu einem reellen Zeitsignal, der ungerade Anteil  X+u(f)  zu einem imaginären.
  • Es ist offensichtlich, dass  X+g(f)  gleich dem tatsächlichen Fourierspektrum  X(f)  und damit der Realteil von  x+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)+jy(t).
  • Nach den allgemein gültigen Gesetzen der Fouriertransformation entsprechend dem  Zuordnungssatz  gilt somit für die Spektralfunktion des Imaginärteils:
jY(f)=X+u(f)=sign(f)X(f)Y(f)=sign(f)jX(f).
  • Transformiert man diese Gleichung in den Zeitbereich, so wird aus der Multiplikation die  Faltungsoperation, und man erhält:
y(t)=1πtx(t)=1π+x(τ)tτdτ.

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.

Definition:  Für die  Hilberttransformierte  H{x(t)}  einer Zeitfunktion  x(t)  gilt:

y(t)=H{x(t)}=1π+x(τ)tτdτ.
  • 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)=jsign(f)X(f).


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)+jH{x(t)}.
  • Die Hilberttransformierte  H{x(t)}  verschwindet nur für den Fall  x(t)=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)=Re{x+(t)}.

Beispiel 2:  Das Prinzip der Hilbert–Transformation wird durch die folgende Grafik nochmals verdeutlicht:

  • Nach der linken Darstellung  (A)  kommt man vom physikalischen Signal  x(t)  zum analytischen Signal  x+(t), indem man einen Imaginärteil  jy(t)  hinzufügt.
  • Hierbei ist  y(t)=H{x(t)}  eine reelle Zeitfunktion, die sich am einfachsten im Spektralbereich durch die Multiplikation des Spektrums  X(f)  mit  jsign(f)  angeben lässt.
Zur Verdeutlichung der Hilbert–Transformierten

Die rechte Darstellung  (B)  ist äquivalent zu  (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)=jy(t)  ist.


Zeigerdiagrammdarstellung der harmonischen Schwingung


Die Spektralfunktion  X(f)  einer harmonischen Schwingung  x(t)=Acos(2πfTtφ)  besteht bekanntlich aus zwei Diracfunktionen bei den Frequenzen

  • +fT  mit dem komplexen Gewicht  A/2ejφ,
  • fT  mit dem komplexen Gewicht  A/2e+jφ.


Somit lautet das Spektrum des analytischen Signals  (also ohne die Diracfunktion bei der Frequenz  f=fT):

X+(f)=Aejφδ(ffT).

Die dazugehörige Zeitfunktion erhält man durch Anwendung des  Verschiebungssatzes:

x+(t)=Aej(2πfTtφ).

Diese Gleichung beschreibt einen mit konstanter Winkelgeschwindigkeit  ωT=2πfT  drehenden Zeiger.

Beispiel 3:  Aus Darstellungsgründen ist in der folgenden Grafik das Koordinatensystem entgegen der üblichen Darstellung um  90  nach links gedreht (Realteil nach oben, Imaginärteil nach links).

Zeigerdiagramm einer harmonischen Schwingung

Anhand dieser Grafik sind folgende Aussagen möglich:

  • Zum Startzeitpunkt  t=0  liegt der Zeiger der Länge  A  (Signalamplitude) mit dem Winkel  φ  in der komplexen Ebene. Im gezeichneten Beispiel gilt  φ=45.
  • Für Zeiten  t>0  dreht der Zeiger mit konstanter Winkelgeschwindigkeit (Kreisfrequenz)  ω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  T0, 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)=Ii=1Aiejφiδ(ffi).

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.

Zeigerdiagramm eines Verbundes aus drei Schwingungen

Im rechten Bild ist die dazugehörige Zeitfunktion angedeutet. Diese lautet allgemein:

x+(t)=Ii=1Aiej(ωitφi).

Zu dieser Grafik anzumerken:

  • Die Skizze zeigt die Ausgangslage der Zeiger zum Startzeitpunkt  t=0  entsprechend den Amplituden  Ai  und den Phasenlagen  φ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)=[1cos(60)1jsin(60)]+2cos(0)+1cos(180)=1.500j0.866.
  • Für Zeiten  t>0  drehen die drei Zeiger mit unterschiedlichen Winkelgeschwindigkeiten  ωi=2πfi. 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.3Z: Hilbert Transformator

Exercise 4.4Z: Vector Diagram for DSB-AM