Difference between revisions of "Aufgaben:Exercise 3.10: Rayleigh Fading"

From LNTwww
(Die Seite wurde neu angelegt: „ {{quiz-Header|Buchseite=Stochastische Signaltheorie/Weitere Verteilungen }} right| :Häufig beschreibt man ein Bandpass…“)
 
m (Text replacement - "root mean square" to "standard deviation")
 
(29 intermediate revisions by 4 users not shown)
Line 1: Line 1:
  
{{quiz-Header|Buchseite=Stochastische Signaltheorie/Weitere Verteilungen
+
{{quiz-Header|Buchseite=Theory_of_Stochastic_Signals/Further_Distributions
 
}}
 
}}
  
[[File:P_ID177__Sto_A_3_10.png|right|]]
+
[[File:P_ID177__Sto_A_3_10.png|right|frame|Rayleigh distribution:  real part, imaginary part, magnitude ]]
:H&auml;ufig beschreibt man ein Bandpass&uuml;bertragungssystem im sogenannten <i>&auml;quivalenten Tiefpassbereich</i>. Diese Darstellung f&uuml;hrt zu einem komplexen Signal (vgl. Buch 1, Kapitel 4):
+
One often describes a band-pass transmission system in the so-called&nbsp; "equivalent low-pass domain".&nbsp; This representation leads according to the chapter&nbsp; "Band-pass signals"&nbsp; in the book&nbsp; [[Signal Representation]]&nbsp; to a complex signal:
:$$\it z(t)=x(t)+\rm \rm j \cdot \it y(t).$$
+
:$$z(t)=x(t)+ {\rm j} \cdot y(t).$$
  
:Der Realteil <i>x</i>(<i>t</i>) kennzeichnet hierbei die <i>Inphasekomponente</i> und der Imagin&auml;rteil die <i>Quadraturkomponente</i>.
+
The real part&nbsp; $x(t)$&nbsp; denotes the&nbsp; "inphase component"&nbsp; and the imaginary part&nbsp; $y(t)$&nbsp; the&nbsp; "quadrature component".
  
:Bei einem Mobilfunksystem, bei dem zwischen dem mobilen Teilnehmer und der sog. Basisstation keine Sichtverbindung besteht, gelangt somit das Funksignal ausschlie&szlig;lich auf indirekten Wegen (Brechung, Streuung, Reflexion usw.) zum Empf&auml;nger. In diesem Fall ist folgendes Modell anwendbar:
+
In a mobile radio system where there is no line of sight between the mobile user and the base station,&nbsp; the radio signal thus reaches the receiver exclusively by indirect means&nbsp; (refraction,&nbsp; scattering,&nbsp; reflection, etc.);  
  
:*Der Realteil <i>x</i>(<i>t</i>) und auch der Imagin&auml;rteil <i>y</i>(<i>t</i>) sind jeweils gau&szlig;verteilt und mittelwertfrei.
+
In this case,&nbsp; the following model is applicable:
 +
*The real part&nbsp; $x(t)$&nbsp; and also the imaginary part&nbsp; $y(t)$&nbsp; are both Gaussian distributed and zero mean.
 +
*$x(t)$&nbsp; and&nbsp; $y(t)$&nbsp; each have the same standard deviation&nbsp; $\sigma$&nbsp; and are independent of each other.
 +
*Internal bindings of the signals&nbsp; $x(t)$&nbsp; and&nbsp; $y(t)$&nbsp; due to the Doppler effect shall not be considered here.
  
:*<i>x</i>(<i>t</i>) und <i>y</i>(<i>t</i>) besitzen jeweils die gleiche Streuung <i>&sigma;</i> und sind voneinander unabh&auml;ngig.
 
  
:*Innere Bindungen der Signale <i>x</i>(<i>t</i>) und <i>y</i>(<i>t</i>) aufgrund des Dopplereffekts sollen hier nicht beachtet werden.
+
The complex signal&nbsp; $z(t)$&nbsp; can also be represented by magnitude and phase:
 +
:$$ z(t)= a(t)\cdot {\rm e}^{\rm j \phi(t)}.$$
 +
*Because of symmetry with respect to&nbsp; $x(t)$&nbsp; and&nbsp; $y(t)$&nbsp; the phase&nbsp; $\phi(t)$&nbsp; is uniformly distributed.
 +
*In contrast,&nbsp; the magnitude&nbsp; $a(t) = |z(t)|$&nbsp; is Rayleigh distributed,&nbsp; which has led to the naming&nbsp; "Rayleigh fading".
 +
*As a further quantity we define the instantaneous power
 +
:$$ p(t)=x^{\rm 2}(t)+y^{\rm 2}( t)=a^{\rm 2}(t).$$
  
:Das komplexe Signal <i>z</i>(<i>t</i>) kann man auch nach Betrag und Phase darstellen:
+
The random variable&nbsp; $p$&nbsp; is&nbsp; [[Theory_of_Stochastic_Signals/Exponentially_Distributed_Random_Variables#One-sided_exponential_distribution|(one-sided) exponentially distributed]].&nbsp; Its PDF is for&nbsp; $p>0$:
:$$ z(t)= a(t)\cdot \rm e^{\rm j \it \phi(\it t)}.$$
+
:$$f_p(p)=\frac{1}{2\sigma^{\rm 2}}\cdot {\rm e}^{ -p/(\sigma^{\rm 2})}.$$
  
:Aufgrund der Symmetrie bez&uuml;glich <i>x</i>(<i>t</i>) und <i>y</i>(<i>t</i>) ist die Phase <i>&#981;</i>(<i>t</i>) gleichverteilt. Dagegen ist der Betrag <i>a</i>(<i>t</i>) = |<i>z</i>(<i>t</i>)| rayleighverteilt, was zu der Namensgebung <i>Rayleighfading</i> gef&uuml;hrt hat.
+
For all negative&nbsp; $p$&ndash;values, of course&nbsp; $f_p(p)= 0$&nbsp; holds,&nbsp; since&nbsp; $p$&nbsp; denotes a power.
  
:Als weitere Gr&ouml;&szlig;e definieren wir noch die Momentanleistung
 
:$$\it p(t)=x^{\rm 2}\it (t)+\it y^{\rm 2}(\it t)=\it a^{\rm 2}(t).$$
 
  
:Die Zufallsgröße <i>p</i> ist hier (einseitig) exponentialverteilt. Deren WDF lautet f&uuml;r <i>p</i> > 0:
 
:$$\it f_p(p)=\frac{\rm 1}{\rm 2\it \sigma^{\rm 2}}\cdot \rm e^{\it -p/(\rm 2\sigma^{\rm 2})}.$$
 
  
:F&uuml;r alle negativen <i>p</i>-Werte gilt natürlich <i>f<sub>p</sub></i>(<i>p</i>) = 0, da <i>p</i> eine Leistung kennzeichnet.
 
  
:Im Folgenden werden die Streuungen der beiden Gau&szlig;schen Zufallsgr&ouml;&szlig;en <i>x</i> und <i>y</i> stets zu <i>&sigma;</i> = 1 vorausgesetzt. Alle Gr&ouml;&szlig;en sind deshalb als normiert zu verstehen.
+
Hints:  
 +
*This exercise belongs to the chapter&nbsp; [[Theory_of_Stochastic_Signals/Further_Distributions|Further Distributions]].
 +
*In particular, reference is made to the section&nbsp; [[Theory_of_Stochastic_Signals/Further_Distributions#Rayleigh_PDF|Rayleigh PDF]].
 +
*The standard deviation of the two Gaussian random variables&nbsp; $x$&nbsp; and&nbsp; $y$ &nbsp;are each&nbsp; $\sigma = 1$.  
 +
*All variables are therefore to be understood as normalized.
  
:<b>Hinweis</b>: Diese Aufgabe bezieht sich auf die Seite Rayleighverteilung im Kapitel 3.7.
+
  
  
  
===Fragebogen===
+
 
 +
===Questions===
  
 
<quiz display=simple>
 
<quiz display=simple>
{Welche der nachfolgenden Aussagen sind stets zutreffend?
+
{Which of the following statements are always true?
 
|type="[]"}
 
|type="[]"}
+ Kleine Momentanleistungen sind wahrscheinlicher als gro&szlig;e.
+
+ Small instantaneous powers are more likely than large ones.
- Die Phase <i>&#981;</i>(<i>t</i>) = &pi;/2 bedeutet auch &bdquo;Imagin&auml;rteil <i>y</i>(<i>t</i>) = 0&rdquo;.
+
- The phase&nbsp; $\phi(t) = \pi/2$&nbsp; also means&nbsp; "imaginary part $y(t) = 0$".
+ Die Phase <i>&#981;</i>(<i>t</i>) = &ndash;&pi;/2 bedeutet auch &bdquo;Realteil <i>x</i>(<i>t</i>) = 0&rdquo;.
+
+ The phase&nbsp; $\phi(t) = -\pi/2$&nbsp; also means&nbsp; "real part $x(t) = 0$".
  
  
{Mit welcher Wahrscheinlichkeit ist die Momentanleistung <i>p</i>(<i>t</i>) > 4?
+
{With what probability is the (normalized) instantaneous power&nbsp; $p(t)$&nbsp; greater than&nbsp; $4$?
 
|type="{}"}
 
|type="{}"}
$Pr(p(t) > 4)$ = { 0.135 3% }
+
${\rm Pr}(p(t) > 4) \ = \ $ { 13.5 3% } $\ \%$
  
  
{Wie gro&szlig; ist die Wahrscheinlichkeit, dass der Betrag <i>a</i>(<i>t</i>) gr&ouml;&szlig;er als 2 ist?
+
{What is the probability that the magnitude&nbsp; $a(t)$&nbsp; is greater than&nbsp; $2$?
 
|type="{}"}
 
|type="{}"}
$Pr(a(t) > 2)$ = { 0.135 3% }
+
${\rm Pr}(a(t) > 2) \ = \ $ { 13.5 3% } $\ \%$
  
  
{Berechnen Sie - ausgehend von <i>f<sub>p</sub></i>(<i>p</i>) - die WDF der Zufallsgr&ouml;&szlig;e <i>a</i>. Welcher WDF-Wert ergibt sich f&uuml;r <i>a</i> = 1?
+
{Calculate,&nbsp; starting from&nbsp; $f_p(p)$&nbsp; the PDF of the random variable&nbsp; $a$.&nbsp; What PDF value results for&nbsp; $a=1$?
 
|type="{}"}
 
|type="{}"}
$f_a(a = 1)$ = { 0.607 3% }
+
$f_a(a = 1)\ = \ $ { 0.607 3% }
  
  
 
</quiz>
 
</quiz>
  
===Musterlösung===
+
===Solution===
 
{{ML-Kopf}}
 
{{ML-Kopf}}
:<b>1.</b>&nbsp;&nbsp;Die erste Aussage trifft aufgrund der Exponentialverteilung für <i>p</i>(<i>t</i>) zu. Ein Phasenwinkel <i>&#981;</i>(<i>t</i>) von &plusmn;90&deg; bedeutet, dass der Realteil <i>x</i>(<i>t</i>) = 0 ist. Bei positivem Imagin&auml;rteil &nbsp;&#8658;&nbsp; <i>y</i>(<i>t</i>) > 0 ist der Phasenwinkel <i>&#981;</i>(<i>t</i>) = +90&deg;, bei negativem Imagin&auml;rteil betr&auml;gt der Phasenwinkel &ndash;90&deg;. Richtig sind also <u>der erste und der dritte Lösungsvorschlag</u>.
+
'''(1)'''&nbsp; Correct are&nbsp; <u>the first and the third proposed solutions</u>:
 +
*The first statement is true because of the exponential distribution for&nbsp; $p(t)$.
 +
* A phase angle&nbsp; $\phi(t) = \pm \pi/2 \ (\pm 90^\circ)$&nbsp; means that the real part&nbsp; $x(t) = 0$.
 +
* For positive imaginary part&nbsp;&#8658;&nbsp; $y(t) > 0$&nbsp; the phase angle&nbsp; $\phi(t) = +90^\circ$,&nbsp; for negative imaginary part the phase angle&nbsp; $\phi(t) = -90^\circ$.
 +
 
 +
 
 +
 
 +
'''(2)'''&nbsp; With&nbsp; $\sigma = 1$&nbsp; holds for the PDF of the instantaneous power:
 +
:$$f_p(p)= {\rm 1}/{\rm 2}\cdot{\rm e}^{-p/\rm 2}.$$
 +
 
 +
*The probability we are looking for is thus:
 +
:$${\rm Pr}(p(t) > 4) = \int_{\rm 4}^{\infty}{1}/{2}\cdot{\rm e}^{-p/\rm 2}\,{\rm d} p={\rm e}^{\rm -2} \rm \hspace{0.15cm}\underline{=13.5\%}.$$
 +
 
  
:<b>2.</b>&nbsp;&nbsp;Mit <i>&sigma;</i> = 1 gilt f&uuml;r die WDF der Momentanleistung:
 
:$$\it f_p(p)= {\rm 1}/{\rm 2}\cdot\rm e^{-\it p/\rm 2}.$$
 
  
:Die gesuchte Wahrscheinlichkeit ist demnach:
+
'''(3)'''&nbsp; Since&nbsp; $p(t) =a^2(t)$&nbsp; holds and moreover &nbsp;$a(t) < 0$&nbsp; is not possible,&nbsp; the event &nbsp;$a(t) > 2$&nbsp; is identical to the event &nbsp;$p(t) > 4$:
:$$\rm Pr(\it p (\it t)> \rm 4) = \int_{\rm 4}^{\infty}\frac{1}{2}\cdot\rm e^{-\it p/\rm 2}\,{\rm d}\it p=\rm e^{\rm -2} \hspace{0.15cm}\underline{=0.135}.$$
 
  
:<b>3.</b>&nbsp;&nbsp;Da <i>p</i> = <i>a</i><sup>2</sup> gilt und <i>a</i> < 0 nicht möglich, ist das Ereignis &bdquo;<i>a</i> > 2&rdquo; identisch mit dem Ereignis &bdquo;<i>p</i> > 4&rdquo;. Es ergibt sich die gleiche Wahrscheinlichkeit <u>0.135</u> wie unter (b) berechnet.
+
*It results in the same probability&nbsp; ${\rm Pr}(a(t) > 2)  \hspace{0.15cm}\underline{=13.5\%}$&nbsp; as calculated in&nbsp; '''(2)'''&nbsp;.
  
:<b>4.</b>&nbsp;&nbsp;Allgemein gilt:
 
:$$\it f_a(a)=\frac{\it f_p(p)}{|g'(p)|}\Big |_{\, \it p=h(a)}.$$
 
  
:Die Transformationskennlinie lautet:
+
'''(4)'''&nbsp; In general:
:$$\it g'(p)=\frac{ \rm d\it a}{ \rm d\it p}=\frac{{\rm1}}{\rm 2 \cdot \it \sqrt{p}}.$$
+
:$$f_a(a)=\frac{f_p(p)}{|g'(p)|}\Bigg |_{\, p=h(a)}.$$
  
:Diese ist stets positiv. Daraus folgt für die Rayleigh&ndash;WDF:
+
*The transformation characteristic is:
:$$\it f_a(a)=\sqrt{\it p}\cdot\rm e^{\it -p/\rm 2}\Big|_{\it p=a^{\rm 2}}=\it a\cdot \rm e^{\, \it -a^{\rm 2}/\rm 2}.$$
+
:$$g\hspace{0.05cm}'(p)=\frac{ {\rm d} a}{ {\rm d} p}=\frac{1}{2 \cdot \sqrt{p}}.$$
  
:F&uuml;r <i>a</i> = 1 ergibt sich somit der Wert e<sup>&ndash;0.5</sup> <u>&asymp; 0.607</u>.
+
*This is always positive.&nbsp; It follows for the Rayleigh PDF:
 +
:$$f_a(a)=\sqrt {p}\cdot {\rm e}^{ -p/\rm 2}\Big|_{ p=a^{\rm 2}}= a\cdot {\rm e}^{\, -a^{\rm 2}/\rm 2}.$$
  
 +
*For&nbsp; $a = 1$&nbsp; this gives the value&nbsp; $f_a(a = 1)= {\rm e}^{-0.5}\hspace{0.15cm}\underline{=0.607} $.
 
{{ML-Fuß}}
 
{{ML-Fuß}}
  
  
  
[[Category:Aufgaben zu Stochastische Signaltheorie|^3.7 Weitere Verteilungen^]]
+
[[Category:Theory of Stochastic Signals: Exercises|^3.7 Further Distributions^]]

Latest revision as of 12:01, 17 February 2022

Rayleigh distribution:  real part, imaginary part, magnitude

One often describes a band-pass transmission system in the so-called  "equivalent low-pass domain".  This representation leads according to the chapter  "Band-pass signals"  in the book  Signal Representation  to a complex signal:

$$z(t)=x(t)+ {\rm j} \cdot y(t).$$

The real part  $x(t)$  denotes the  "inphase component"  and the imaginary part  $y(t)$  the  "quadrature component".

In a mobile radio system where there is no line of sight between the mobile user and the base station,  the radio signal thus reaches the receiver exclusively by indirect means  (refraction,  scattering,  reflection, etc.);

In this case,  the following model is applicable:

  • The real part  $x(t)$  and also the imaginary part  $y(t)$  are both Gaussian distributed and zero mean.
  • $x(t)$  and  $y(t)$  each have the same standard deviation  $\sigma$  and are independent of each other.
  • Internal bindings of the signals  $x(t)$  and  $y(t)$  due to the Doppler effect shall not be considered here.


The complex signal  $z(t)$  can also be represented by magnitude and phase:

$$ z(t)= a(t)\cdot {\rm e}^{\rm j \phi(t)}.$$
  • Because of symmetry with respect to  $x(t)$  and  $y(t)$  the phase  $\phi(t)$  is uniformly distributed.
  • In contrast,  the magnitude  $a(t) = |z(t)|$  is Rayleigh distributed,  which has led to the naming  "Rayleigh fading".
  • As a further quantity we define the instantaneous power
$$ p(t)=x^{\rm 2}(t)+y^{\rm 2}( t)=a^{\rm 2}(t).$$

The random variable  $p$  is  (one-sided) exponentially distributed.  Its PDF is for  $p>0$:

$$f_p(p)=\frac{1}{2\sigma^{\rm 2}}\cdot {\rm e}^{ -p/(\sigma^{\rm 2})}.$$

For all negative  $p$–values, of course  $f_p(p)= 0$  holds,  since  $p$  denotes a power.



Hints:

  • This exercise belongs to the chapter  Further Distributions.
  • In particular, reference is made to the section  Rayleigh PDF.
  • The standard deviation of the two Gaussian random variables  $x$  and  $y$  are each  $\sigma = 1$.
  • All variables are therefore to be understood as normalized.




Questions

1

Which of the following statements are always true?

Small instantaneous powers are more likely than large ones.
The phase  $\phi(t) = \pi/2$  also means  "imaginary part $y(t) = 0$".
The phase  $\phi(t) = -\pi/2$  also means  "real part $x(t) = 0$".

2

With what probability is the (normalized) instantaneous power  $p(t)$  greater than  $4$?

${\rm Pr}(p(t) > 4) \ = \ $

$\ \%$

3

What is the probability that the magnitude  $a(t)$  is greater than  $2$?

${\rm Pr}(a(t) > 2) \ = \ $

$\ \%$

4

Calculate,  starting from  $f_p(p)$  the PDF of the random variable  $a$.  What PDF value results for  $a=1$?

$f_a(a = 1)\ = \ $


Solution

(1)  Correct are  the first and the third proposed solutions:

  • The first statement is true because of the exponential distribution for  $p(t)$.
  • A phase angle  $\phi(t) = \pm \pi/2 \ (\pm 90^\circ)$  means that the real part  $x(t) = 0$.
  • For positive imaginary part ⇒  $y(t) > 0$  the phase angle  $\phi(t) = +90^\circ$,  for negative imaginary part the phase angle  $\phi(t) = -90^\circ$.


(2)  With  $\sigma = 1$  holds for the PDF of the instantaneous power:

$$f_p(p)= {\rm 1}/{\rm 2}\cdot{\rm e}^{-p/\rm 2}.$$
  • The probability we are looking for is thus:
$${\rm Pr}(p(t) > 4) = \int_{\rm 4}^{\infty}{1}/{2}\cdot{\rm e}^{-p/\rm 2}\,{\rm d} p={\rm e}^{\rm -2} \rm \hspace{0.15cm}\underline{=13.5\%}.$$


(3)  Since  $p(t) =a^2(t)$  holds and moreover  $a(t) < 0$  is not possible,  the event  $a(t) > 2$  is identical to the event  $p(t) > 4$:

  • It results in the same probability  ${\rm Pr}(a(t) > 2) \hspace{0.15cm}\underline{=13.5\%}$  as calculated in  (2) .


(4)  In general:

$$f_a(a)=\frac{f_p(p)}{|g'(p)|}\Bigg |_{\, p=h(a)}.$$
  • The transformation characteristic is:
$$g\hspace{0.05cm}'(p)=\frac{ {\rm d} a}{ {\rm d} p}=\frac{1}{2 \cdot \sqrt{p}}.$$
  • This is always positive.  It follows for the Rayleigh PDF:
$$f_a(a)=\sqrt {p}\cdot {\rm e}^{ -p/\rm 2}\Big|_{ p=a^{\rm 2}}= a\cdot {\rm e}^{\, -a^{\rm 2}/\rm 2}.$$
  • For  $a = 1$  this gives the value  $f_a(a = 1)= {\rm e}^{-0.5}\hspace{0.15cm}\underline{=0.607} $.