Probability Density of Rayleigh Fading

From LNTwww

Allgemeine Beschreibung des Mobilfunkkanals


Im Folgenden wird zur Vereinfachung der Schreibweise auf den Zusatz „TP” verzichtet. Somit liegt das reelle Signal s(t) = 1 am Eingang des Mobilfunkkanals an und das Ausgangssignal r(t) ist komplexwertig. Zusätzliche Rauschprozesse werden ausgeschlossen.

Das Funksignal s(t) kann den Empfänger über eine Vielzahl von Pfaden erreichen, wobei die einzelnen Signalanteile in unterschiedlicher Weise gedämpft und verschieden lang verzögert werden. Allgemein kann man für das Tiefpass–Empfangssignal ohne Berücksichtigung von thermischem Rauschen schreiben:

\[r(t)= \sum_{k=1}^{K} \alpha_{k}(t) \cdot {\rm e}^{\hspace{0.05cm}{\rm j}\hspace{0.02cm}\cdot \hspace{0.02cm} \phi_{k}(t)} \cdot s(t - \tau_{k}) \hspace{0.05cm}.\]

Hierbei sind folgende Bezeichnungen verwendet:

  • Der zeitabhängige Dämpfungsfaktor auf dem k–ten Pfad ist αk(t).
  • Der zeitabhängige Phasenverlauf auf dem k–ten Pfad ist ϕk(t).
  • Die Laufzeit auf dem k–ten Pfad ist τk.

Die Anzahl K der sich (zumindest geringfügig) unterscheidenden Pfade ist meist sehr groß und für eine direkte Modellierung ungeeignet. Das Modell lässt sich aber entscheidend vereinfachen, wenn man jeweils Pfade mit näherungsweise gleichen Verzögerungen zusammenfasst. Man unterscheidet somit nur noch M Hauptpfade, die durch großräumige Wegeunterschiede und damit merkliche Laufzeitunterschiede gekennzeichnet sind:

\[r(t)= \sum_{m=1}^{M} \hspace{0.1cm} \sum_{n=1}^{N_m} \alpha_{m,\hspace{0.01cm}n}(t) \cdot {\rm e}^{\hspace{0.05cm}{\rm j}\hspace{0.02cm}\cdot \hspace{0.02cm} \phi_{m,\hspace{0.02cm}n}(t)} \cdot s(t - \tau_{m,\hspace{0.02cm}n}) \hspace{0.05cm}.\]

Die beiden bisher angegebenen Gleichungen sind identisch. Eine Vereinfachung ergibt sich, wenn man für jeden Hauptpfad m die Nm Laufzeiten, die sich durch Reflexionen an Feinstrukturen sowie eventuell durch Beugungs– und Brechungserscheinungen geringfügig unterscheiden, durch eine mittlere Laufzeit ersetzt:

\[\tau_{m} = \frac{1}{N_m} \cdot \sum_{n=1}^{N_m} \tau_{m,\hspace{0.02cm}n} \hspace{0.05cm}.\]

Damit erhält man das folgende wichtige Zwischenergebnis:

\[r(t)= \sum_{m=1}^{M} z_m(t) \cdot s(t - \tau_{m}) \hspace{0.5cm} {\rm mit} \hspace{0.5cm} z_m(t) = \sum_{n=1}^{N_m} \alpha_{m,\hspace{0.01cm}n}(t) \cdot {\rm e}^{\hspace{0.05cm}{\rm j}\hspace{0.02cm}\cdot \hspace{0.02cm} \phi_{m,\hspace{0.02cm}n}(t)} \hspace{0.05cm}.\]

Modellierung von nichtfrequenzselektivem Fading (1)


Ausgehend von der soeben hergeleiteten Gleichung

\[r(t)= \sum_{m=1}^{M} z_m(t) \cdot s(t - \tau_{m}) \hspace{0.5cm} {\rm mit} \hspace{0.5cm} z_m(t) = \sum_{n=1}^{N_m} \alpha_{m,\hspace{0.01cm}n}(t) \cdot {\rm e}^{\hspace{0.05cm}{\rm j}\hspace{0.02cm}\cdot \hspace{0.02cm} \phi_{m,\hspace{0.02cm}n}(t)} \hspace{0.05cm}\]

können zwei wichtige Sonderfälle abgeleitet werden:

  • Gibt es mehr als einen Hauptpfad (M ≥ 2), so spricht man von Mehrwegeausbreitung. Wie im Kapitel 2 noch gezeigt werden wird, kommt es dann – je nach Frequenz – zu konstruktiven oder destruktiven Überlagerungen bis hin zu völliger Auslöschung. Für manche Frequenzen erweist sich die Mehrwegeausbreitung als günstig, für andere als extrem ungünstig. Man bezeichnet diesen Effekt auch als frequenzselektives Fading.
  • Bei nur einem Hauptpfad (M = 1, auf den Index „1” verzichten wir in diesem Fall) vereinfacht sich die obige Gleichung wie folgt:
\[r(t)= z(t) \cdot s(t - \tau) \hspace{0.05cm}.\]
Die Verzögerung τ bewirkt hier eine für alle Frequenzen konstante Laufzeit, die nicht weiter betrachtet wird. Es gibt nun keine Überlagerungen von Signalanteilen mit merklichen Laufzeitunterschieden und damit auch keine Frequenzabhängigkeit des Gesamtsignals. Man spricht deshalb von nichtfrequenzselektivem Fading oder Flat–Fading. Für dieses gilt:
\[r(t)= z(t) \cdot s(t) \hspace{0.5cm} {\rm mit} \hspace{0.5cm} z(t) = \sum_{n=1}^{N} \alpha_{n}(t) \cdot {\rm e}^{\hspace{0.05cm}{\rm j}\hspace{0.02cm}\cdot \hspace{0.02cm} \phi_{n}(t)} \hspace{0.05cm}. \]

Die Grafik zeigt das Modell zur Erzeugung von nichtfrequenzselektivem Fading. Man spricht auch von Rayleigh–Fading.

Rayleigh–Fading–Kanalmodell

Die Bildbeschreibung folgt auf der nächsten Seite.

Modellierung von nichtfrequenzselektivem Fading (2)


Wir betrachten die multiplikative Verfälschung z(t) entsprechend dem Rayleigh–Modell genauer. Für den komplexen Koeffizienten gilt entsprechend der letzten Seite:

\[z(t) = \sum_{n=1}^{N} \alpha_{n}(t) \cdot {\rm e}^{\hspace{0.05cm}{\rm j}\hspace{0.02cm}\cdot \hspace{0.02cm} \phi_{n}(t) }= \sum_{n=1}^{N} \alpha_{n}(t) \cdot \cos( \phi_{n}( t)) + {\rm j}\cdot \sum_{n=1}^{N} \alpha_{n}(t) \cdot \sin( \phi_{n}( t)) \hspace{0.05cm}. \]

Das Empfangssignal r(t) ergibt sich, wenn man das Sendesignal s(t) mit der Zeitfunktion z(t) multipliziert. Es sei nochmals daran erinnert, dass sich alle Signale bzw. Zeitfunktionen s(t), z(t) und r(t) auf den äquivalenten Tiefpassbereich beziehen.

Rayleigh–Fading–Kanalmodell

Zu obiger Gleichung und der Grafik ist anzumerken:

  • αn(t) und ϕn(t) hängen von den Umgebungsbedingungen ab. ϕn(t) erfasst die verschiedenen Laufzeiten auf den N Pfaden und den Dopplereffekt aufgrund der Bewegung.
  • Die Zeitfunktion z(t) ist eine komplexe Größe, deren Real– und Imaginärteil wir im Folgenden wieder mit x(t) und y(t) bezeichnen.
  • Eine deterministische Beschreibung der Zufallsgröße z(t) = x(t) + j · y(t) ist nicht möglich, vielmehr müssen diese Funktionen durch stochastische Prozesse modelliert werden.
  • Ist die Anzahl N der (leicht) unterschiedlichen Laufzeiten hinreichend groß, so ergeben sich nach dem zentralen Grenzwertsatz Gaußsche Zufallsgrößen x(t) und y(t).
  • x(t) und y(t) sind jeweils mittelwertfrei und besitzen die gleiche Varianz σ2:
\[{\rm E}[x(t)] = {\rm E}[y(t)] = 0\hspace{0.05cm}, \hspace{0.2cm}{\rm E}[x^2(t)] = {\rm E}[y^2(t)] = \sigma^2 \hspace{0.05cm}.\]
  • Zu berücksichtigen ist die Orthogonität von Realteil und Imaginärteil (jeweils Cosinus und Sinus des gleichen Arguments); damit sind sie auch unkorreliert. Nur bei Gaußschen Zufallsgrößen folgt daraus weiter die statistische Unabhängigkeit von x(t) und y(t).
  • Aufgrund des Dopplereffekts gibt es allerdings statistische Bindungen innerhalb des Realteils x(t) und innerhalb des Imaginärteils y(t). Diese werden im Modell durch zwei Digitale Filter erzeugt.

Beispielhafte Signalverläufe bei Rayleigh–Fading (1)


Die Grafiken auf dieser und der nächsten Seite zeigen jeweils durch Simulation gewonnene Signalverläufe von 100 ms Dauer und die dazugehörigen Dichtefunktionen. Es handelt sich um Bildschirmabzüge des Windows–Programms „Mobilfunkkanal” aus dem Praktikum [Söd01][1].

Realteil, Imaginärteil und Phasenverlauf bei Rayleigh-Fading

Die Darstellungen lassen sich wie folgt interpretieren:

  • Der Realteil ist gaußverteilt (siehe rechte obere Grafik), wie auch aus dem Zeitsignalverlauf x(t) hervorgeht. Rot eingezeichnet ist die Gaußsche WDF fx(x) und blau das durch Simulation über 10.000 Abtastwerte gewonnene Histogramm.
  • Im Programm eingestellt war für diese Darstellung eine maximale Dopplerfrequenz von 100 Hz. Deshalb gibt es statistische Bindungen innerhalb der Signale x(t) und y(t). Näheres hierzu finden Sie im Kapitel 1.3.
  • Die WDF fy(y) des Imaginärteils ist identisch mit fx(x). Die Varianz beträgt jeweils σ2 = 0.5. Zwischen x(t) und y(t) bestehen keine statistischen Bindungen; die Signale sind orthogonal.
  • Die Phase ϕ(t) ist gleichverteilt zwischen ±π. Wie aus den Sprungstellen im Phasenverlauf zu erkennen, kann ϕ(t) durchaus größere Werte annehmen. Alle Bereiche (2k±1)π wurden aber bei der Histogrammerstellung auf den Wertebereich –π ... +π projiziert (k ganzzahlig).
  • Die gleichverteilte Phase wird anhand der (hier nicht dargestellten) 2D–WDF verständlich. Diese ist rotationssymmetrisch und dementsprechend gibt es auch keine Vorzugsrichtung:
\[f_{x,\hspace{0.02cm}y}(x, y) = \frac{1}{2\pi \cdot \sigma^2} \cdot {\rm exp} \left [ - \frac{x^2 + y^2}{2\sigma^2} \right ] .\]

Beispielhafte Signalverläufe bei Rayleigh–Fading (2)


Die Grafik zeigt oben nochmals Real– und Imaginärteil von z(t) = x(t) + j · y(t). Darunter gezeichnet sind Verlauf und WDF von Betrag a(t) = |z(t)| und Betragsquadrat p(t) = a2(t) = |z(t)|2.

Realteil, Imaginärteil, Betrag und Betragsquadrat bei Rayleigh-Fading

Aus diesen Darstellungen geht hervor:

  • Der Betrag besitzt eine Rayleigh–WDF  ⇒  Name „Rayleigh–Fading”:
\[f_a(a) = \left\{ \begin{array}{c} a/\sigma^2 \cdot {\rm exp} [ -a^2/(2\sigma^2)] \\ 0 \end{array} \right.\quad \begin{array}{*{1}c} {\rm f\ddot{u}r}\hspace{0.15cm} a \ge 0 \\ {\rm f\ddot{u}r}\hspace{0.15cm} a < 0 \\ \end{array} \hspace{0.05cm}.\]
  • Für die Momente erster bzw. zweiter Ordnung und die Varianz des Betrags a(t) = |z(t)| gilt:
\[{\rm E}[a] = \sigma \cdot \sqrt {{\pi}/{2}}\hspace{0.05cm},\hspace{0.2cm}{\rm E}[a^2] = 2 \cdot \sigma^2 \hspace{0.3cm} \Rightarrow \hspace{0.3cm} {\rm Var}[a] = \sigma_a^2 = \sigma^2 \cdot \left ( 2 - {\pi}/{2}\right ) \hspace{0.05cm}. \]
  • Die WDF des Betragsquadrats p(t) ergibt sich durch nichtlineare Transformation der WDF fa(a) und führt zu einer Exponentialverteilung:
\[f_p(p) = \left\{ \begin{array}{c} 1/(2\sigma^2) \cdot {\rm exp} [ -p/(2\sigma^2)] \\ 0 \end{array} \right.\quad \begin{array}{*{1}c} {\rm f\ddot{u}r}\hspace{0.15cm} p \ge 0 \\ {\rm f\ddot{u}r}\hspace{0.15cm} p < 0 \\ \end{array} \hspace{0.05cm}.\]

Weitere Informationen zum Rayleigh–Fading finden Sie in Aufgabe A1.3 und Aufgabe Z1.3.

Aufgaben


A1.3 Rayleigh–Fading

Zusatzaufgaben:1.3 Nochmals Rayleigh–Fading?

Quellenverzeichnis

  1. Söder, G.: Simulation digitaler Übertragungssysteme. Anleitung zum gleichnamigen Praktikum. Lehrstuhl für Nachrichtentechnik, Technische Universität München, 2001.