Difference between revisions of "Digital Signal Transmission/Carrier Frequency Systems with Non-Coherent Demodulation"

From LNTwww
 
(56 intermediate revisions by 7 users not shown)
Line 1: Line 1:
 
   
 
   
 
{{Header
 
{{Header
|Untermenü=Verallgemeinerte Beschreibung digitaler Modulationsverfahren
+
|Untermenü=Generalized Description of Digital Modulation Methods
 
|Vorherige Seite=Trägerfrequenzsysteme mit kohärenter Demodulation
 
|Vorherige Seite=Trägerfrequenzsysteme mit kohärenter Demodulation
 
|Nächste Seite=Beschreibungsgrößen digitaler Kanalmodelle
 
|Nächste Seite=Beschreibungsgrößen digitaler Kanalmodelle
 
}}
 
}}
  
== Rayleigh– und Riceverteilung (1) ==
+
== Rayleigh and Rice Distribution ==
 
<br>
 
<br>
Die für eine kohärente Demodulation erforderliche Schätzung des Phasenwinkels aus dem ankommenden Signal ist bei vielen Anwendungen nicht oder nur eingeschränkt möglich. So führt die Bewegung eines Mobilteilnehmers mit hoher Geschwindigkeit zu sehr schnellen zeitlichen Änderungen des Phasenwinkels  <i>&#981;</i>, was dessen ausreichend genaue Bestimmung erschwert oder gar verhindert.<br>
+
The estimation of the phase angle from the incoming signal,&nbsp; which is required for coherent demodulation,&nbsp; is not possible or only possible to a limited extent in many applications.&nbsp; For example,&nbsp; the movement of a mobile subscriber at high speed leads to very rapid temporal changes in the phase angle &nbsp;$\phi$,&nbsp; which makes its sufficiently accurate determination difficult or even impossible.<br>
  
Diese Tatsache führt zu den nichtkohärenten Demodulationsverfahren mit dem Vorteil reduzierter Komplexität, allerdings mit erhöhter Verfälschungswahrscheinlichkeit. Bei der Herleitung der Gleichungen stößt man auf zwei Wahrscheinlichkeitsdichtefunktionen, die hier vorneweg angegeben werden:
+
This fact leads to the&nbsp; '''non-coherent demodulation processes'''&nbsp; with the advantage of reduced complexity,&nbsp; but with increased error probability.&nbsp; In the derivation of the equations one encounters two probability density functions,&nbsp; which are given here in advance:
*Die Rayleighverteilung erhält man für die WDF der Zufallsgröße <i>y</i> mit Realisierung <i>&eta;</i>, die sich aus den beiden gaußverteilten und statistisch unabhängigen Komponenten <i>u</i> und <i>&upsilon;</i> (beide mit der gleichen Streuung <i>&sigma;<sub>n</sub></i>) wie folgt ergibt:
+
 
 +
*The&nbsp; [[Theory_of_Stochastic_Signals/Further_Distributions#Rayleigh_PDF|'''Rayleigh distribution''']]&nbsp; is obtained for the power density function&nbsp; $\rm (PDF)$&nbsp; of the random variable&nbsp; $y$&nbsp; with realization&nbsp; $\eta$,&nbsp; which is obtained from the two Gaussian distributed and statistically independent components&nbsp; $u$&nbsp; and&nbsp; $v$&nbsp; $($both with the same standard deviation&nbsp; $\sigma_n)$&nbsp; as follows:
  
 
::<math>y = \sqrt{u^2 + v^2} \hspace{0.3cm} \Rightarrow \hspace{0.3cm} p_y (\eta) ={\eta}/{\sigma_n^2}
 
::<math>y = \sqrt{u^2 + v^2} \hspace{0.3cm} \Rightarrow \hspace{0.3cm} p_y (\eta) ={\eta}/{\sigma_n^2}
Line 17: Line 18:
 
  \hspace{0.05cm}.</math>
 
  \hspace{0.05cm}.</math>
  
*Die Riceverteilung erhält man unter sonst gleichen Randbedingungen für den Fall, dass bei einer der Komponenten (entweder <i>u</i> oder <i>&upsilon;</i>) noch eine Konstante <i>C</i> addiert wird:
+
*The&nbsp; [[Theory_of_Stochastic_Signals/Further_Distributions#Rice_PDF|'''Rice distribution''']]&nbsp; is obtained under the same boundary conditions for the case where a constant &nbsp; $C$&nbsp; is added to one of the components $($either to&nbsp;$u$&nbsp; or &nbsp;$v)$:&nbsp;
  
 
::<math>y = \sqrt{(u+C)^2 + v^2} \hspace{0.3cm} \Rightarrow \hspace{0.3cm} p_y (\eta) = {\eta}/{\sigma_n^2}
 
::<math>y = \sqrt{(u+C)^2 + v^2} \hspace{0.3cm} \Rightarrow \hspace{0.3cm} p_y (\eta) = {\eta}/{\sigma_n^2}
Line 23: Line 24:
 
  \hspace{0.05cm}.</math>
 
  \hspace{0.05cm}.</math>
  
Für die Riceverteilung benötigt man die <i>modifizierte Besselfunktion nullter Ordnung</i>, deren Definition und Reihenentwicklung wie folgt lauten:
+
:The Rice distribution uses the&nbsp; "modified zero-order Bessel function",&nbsp; whose definition and series expansion are as follows:
 +
::<math>{\rm I }_0 (x) = \frac{1}{ 2\pi} \cdot  \int_{0}^{2\pi} {\rm e }^{-x \hspace{0.03cm}\cdot \hspace{0.03cm}\cos(\alpha)} \,{\rm d} \alpha \hspace{0.2cm} \approx \hspace{0.2cm}
 +
\sum_{k = 0}^{\infty} \frac{(x/2)^{2k}}{k! \cdot \Gamma (k+1)}
 +
\hspace{0.05cm}.</math>
 +
The graph shows the probability density functions of Rayleigh and Rice distributions.&nbsp; To this representation it is to be noted:
 +
[[File:P ID2082 Dig T 4 5 S1 version1_ret.png|right|frame|Rayleigh resp. Rice PDF for $&sigma;_n = 0.5$|class=fit]]
 +
*The Rice distribution is determined by the two parameters&nbsp; $C$&nbsp; and&nbsp; $\sigma_n$.&nbsp; With&nbsp; $C = 0$,&nbsp; the Rice PDF is identical to the Rayleigh PDF.<br>
 +
 
 +
*The Rayleigh PDF with larger&nbsp; $\sigma_n$&nbsp; is of the same shape as shown for&nbsp; $\sigma_n = 0.5$,&nbsp; but wider and lower in the ratio of standard deviation.<br>
  
:<math>{\rm I }_0 (x) = \frac{1}{ 2\pi} \cdot  \int_{0}^{2\pi} {\rm e }^{-x \hspace{0.03cm}\cdot \hspace{0.03cm}\cos(\alpha)} \,{\rm d} \alpha \hspace{0.2cm} \approx \hspace{0.2cm}
+
*$\sigma_n$&nbsp; indicates the standard deviation of the two Gaussian distributed random variables&nbsp; $u$&nbsp; and&nbsp; $v$&nbsp; and not the standard deviation of the Rayleigh distributed random variable&nbsp; $y$.&nbsp; Rather,&nbsp; the following applies to the latter:
  \sum_{k = 0}^{\infty} \frac{(x/2)^{2k}}{k! \cdot \Gamma (k+1)}
+
 
 +
::<math>\sigma_y = \sigma_n  \cdot  \sqrt{2 - {\pi}/{2 }} \hspace{0.2cm} \approx \hspace{0.2cm} 0.655 \cdot \sigma_n 
 +
\hspace{0.05cm}.</math>
 +
 
 +
*The Rayleigh PDF is extremely asymmetric,&nbsp; recognizable by the large value for the&nbsp;  [[Theory_of_Stochastic_Signals/Expected_Values_and_Moments#Some_common_central_moments|"third order central moment"]] &nbsp; &rArr; &nbsp; "Charlier's skewness":
 +
:$$\mu_3/\sigma_y \approx 0.27.$$
 +
 
 +
*The Rice PDF is more symmetrical the larger the quotient&nbsp; $C/\sigma_n$.&nbsp; For&nbsp; $C/\sigma_n \ge 4$:&nbsp; $\mu_3 \approx 0$ &nbsp; &rArr; &nbsp; symmetrical PDF.
 +
 
 +
*The larger&nbsp; $C/\sigma_n$&nbsp; is,&nbsp; the more the Rice PDF $($with &nbsp;$C$, &nbsp;$\sigma_n)$&nbsp; approaches a&nbsp; [[Theory_of_Stochastic_Signals/Gaussian_Distributed_Random_Variables#Probability_density_function_.E2.80.93_Cumulative_density_function|"Gaussian PDF"]]&nbsp; with mean&nbsp; $C$&nbsp; and standard deviation&nbsp; $\sigma_n$:
 +
 
 +
:$$p_y (\eta) \approx \frac{1}{\sqrt{2\pi} \cdot \sigma_n}
 +
  \cdot {\rm exp } \left [ - \frac{(\eta - C)^2}{2 \sigma_n^2}\right ]$$
 +
:$$ \Rightarrow \hspace{0.3cm} m_y = C\hspace{0.05cm},\hspace{0.2cm}\sigma_y = \sigma_n
 +
\hspace{0.05cm}.$$
 +
 
 +
== Non-coherent demodulation of&nbsp; "on–off keying"&nbsp; (OOK)==
 +
<br>
 +
We consider&nbsp;  [[Digital_Signal_Transmission/Carrier_Frequency_Systems_with_Coherent_Demodulation#On.E2.80.93off_keying_.282.E2.80.93ASK.29|"on&ndash;off keying"]]&nbsp; $(\rm OOK$&nbsp; or&nbsp; $\rm 2&ndash;ASK)$&nbsp; in the equivalent low-pass range.
 +
[[File:EN_Dig_A_1_9_neu.png|right|frame|Coherent and non-coherent demodulation of&nbsp; "on-off keying" |class=fit]]
 +
 +
*In the case of coherent demodulation&nbsp; (left graph),&nbsp; the signal space constellation of the received signal&nbsp; $\boldsymbol{r}(t)$&nbsp; is the same as that of the transmitted signal&nbsp; $\boldsymbol{s}(t)$&nbsp; and consists of two points.
 +
 
 +
*The decision boundary&nbsp; $G$&nbsp; lies in the middle between these points&nbsp; $\boldsymbol{r}_0$&nbsp; and&nbsp; $\boldsymbol{r}_1$.
 +
 +
*The arrows mark the rough direction of noise vectors that may cause transmission errors.
 +
 
 +
 
 +
On the other hand, for non-coherent demodulation (right graph):
 +
*The point&nbsp; $\boldsymbol{r}_1 = \boldsymbol{s}_1 = 0$&nbsp; is still preserved.&nbsp; But&nbsp; $\boldsymbol{r}_0 = \boldsymbol{s}_0 \cdot {\rm e}^{\hspace{0.02cm}{\rm j}\hspace{0.03cm}\phi}$&nbsp; can lie on any point of the circle around&nbsp; $\boldsymbol{s}_0$,&nbsp; since&nbsp; $\phi$&nbsp; is unknown.<br>
 +
 
 +
*The decision process in consideration of the AWGN noise happens in two dimensions&nbsp; $($indicated by arrows in the right graph$)$.<br>
 +
 
 +
*The decision area&nbsp; $I_1$&nbsp; is a circle whose radius&nbsp; $G$&nbsp; is an optimizable parameter. The decision area&nbsp; $I_0$ &nbsp;lies outside the circle.<br>
 +
 
 +
 
 +
Thus,&nbsp; the structure of the optimal OOK receiver&nbsp; (in the equivalent low&ndash;pass range)&nbsp; is fixed &nbsp; &rArr; &nbsp;  see second graph:
 +
[[File:EN_Dig_T_4_5_S2b_neu.png|right|frame|Receiver for non-coherent OOK demodulation&nbsp; (complex signals are labeled blue) |class=fit]]
 +
 +
*The input signal&nbsp; $\boldsymbol{r}(t) = \boldsymbol{s}(t) \cdot {\rm e}^{\hspace{0.02cm}{\rm j}\hspace{0.03cm}\phi} + \boldsymbol{n}(t)$&nbsp; is generally complex because of the phase angle&nbsp; $\phi$&nbsp; and because of the complex noise term&nbsp; $\boldsymbol{n}(t)$.&nbsp;
 +
 
 +
*Consequently required is now the correlation between the complex received signal&nbsp;  $\boldsymbol{r}(t)$&nbsp; and a complex basis function&nbsp; $\boldsymbol{\xi}_1(t)$.<br>
 +
 
 +
*The result is the&nbsp; (complex)&nbsp; detection value&nbsp; $\boldsymbol{r}$,&nbsp; from which the magnitude&nbsp; $y = |\boldsymbol{r}(t)|$&nbsp; is formed as a real decision input variable.<br>
 +
 
 +
*If the decision value&nbsp; $y \gt G$,&nbsp;  then&nbsp; $m_0$&nbsp; is output as the estimated value,&nbsp; otherwise&nbsp; $m_1$.
 +
 
 +
*Thus,&nbsp; the symbol error probability for equally probable symbols <br>&rArr; &nbsp; ${\rm Pr}(\boldsymbol{m}_0) = {\rm Pr}(\boldsymbol{m}_1) =1/2$&nbsp; is:
 +
 
 +
::<math>p_{\rm S} = {\rm Pr}({\cal{E}}) =  {1}/{ 2} \cdot \int_{0}^{G} p_{y\hspace{0.05cm}\vert \hspace{0.05cm}m} (\eta \hspace{0.05cm}|\hspace{0.05cm} m_0) \,{\rm d} \eta
 +
  + {1}/{ 2} \cdot \int_{G}^{\infty} p_{y\hspace{0.05cm}\vert \hspace{0.05cm}m} (\eta \hspace{0.05cm}\vert \hspace{0.05cm} m_1) \,{\rm d} \eta 
 +
\hspace{0.05cm}.</math>
 +
 
 +
*However,&nbsp; due to the Rice PDF &nbsp; $p_{y\hspace{0.05cm}|\hspace{0.05cm}m} (\eta\hspace{0.05cm}|\hspace{0.05cm}m_0)$ &nbsp; and the Rayleigh PDF &nbsp; $p_{y\hspace{0.05cm}|\hspace{0.05cm}m} (\eta\hspace{0.05cm}|\hspace{0.05cm}m_1)$, &nbsp; this probability can only be determined numerically.
 +
 
 +
* The optimal decision limit&nbsp; $G$&nbsp; has to be determined beforehand as the solution of the following equation:
 +
::<math>p_{y\hspace{0.05cm}\vert \hspace{0.05cm}m} (G \hspace{0.05cm}|\hspace{0.05cm}m_0) = p_{y\hspace{0.05cm}|\hspace{0.05cm}m} (G \hspace{0.05cm}|\hspace{0.05cm}m_1)
 
  \hspace{0.05cm}.</math>
 
  \hspace{0.05cm}.</math>
  
Die Grafik zeigt Rayleigh&ndash; und Rice&ndash;Wahrscheinlichkeitsdichtefunktionen. Die Bildbeschreibung folgt auf der nächsten Seite.<br>
+
{{GraueBox|TEXT=
 +
[[File:P ID3148 Dig T 4 5 S2c version1.png|right|frame|Density functions for "OOK, non-coherent"]]
 +
$\text{Example 1:}$&nbsp; The graph shows the result of this equation for&nbsp; $\sigma_n = 0.5$&nbsp; and&nbsp; $C = 2$, where the (red) Rice PDF is approximated by a Gaussian PDF with mean&nbsp; $C$&nbsp; and standard deviation&nbsp; $\sigma_n$.&nbsp;
 +
 
 +
One can see from this sketch:
 +
*The optimal decision value &nbsp;$($here: &nbsp; $G \approx 1.25)$&nbsp; results from the intersection of the two PDF curves.<br>
 +
 
 +
*The symbol error probability&nbsp; $p_{\rm S}$&nbsp; is the sum of the two colored areas.&nbsp; In the example:&nbsp; $p_{\rm S} \approx 5\%$ results.
 +
 
  
[[File:P ID2082 Dig T 4 5 S1 version1.png|Rayleigh- und Rice-WDF|class=fit]]<br>
+
&rArr; &nbsp; You can determine the error probability for other values of&nbsp; $C$&nbsp; and&nbsp; $\sigma_n$&nbsp; as well as the optimal decision value&nbsp; $G$&nbsp; using the&nbsp; HTML5/JavaScript applet&nbsp;  [[Applets:Coherent_and_Non-Coherent_On-Off_Keying|"Coherent and Non-coherent On-Off Keying"]].&nbsp; }}
  
== Rayleigh– und Riceverteilung (2) ==
+
== Non-coherent demodulation of binary FSK (2&ndash;FSK)==
 
<br>
 
<br>
Die Grafik am Seitenende zeigt nochmals die Wahrscheinlichkeitsdichtefunktionen von Rayleigh&ndash; und Riceverteilung. Zu dieser Darstellung ist anzumerken:
+
As already shown in the&nbsp; [[Digital_Signal_Transmission/Carrier_Frequency_Systems_with_Coherent_Demodulation#Binary_frequency_shift_keying_.282.E2.80.93FSK.29|"last chapter"]],&nbsp; "Binary Frequency Shift Keying"&nbsp; $\rm (2&ndash;FSK)$&nbsp; in the equivalent low-pass range can be described by the basis functions
*Die [http://en.lntwww.de/Stochastische_Signaltheorie/Weitere_Verteilungen#Riceverteilung Riceverteilung] ist durch die beiden Parameter <i>C</i> und <i>&sigma;<sub>n</sub></i> bestimmt. Mit <i>C</i> = 0 ist die Rice&ndash;WDF identisch mit der [http://en.lntwww.de/Stochastische_Signaltheorie/Weitere_Verteilungen#Rayleighverteilung Rayleigh&ndash;WDF].<br>
+
::<math>\xi_1(t) = \sqrt{1/T} \cdot {\rm e}^{\hspace{0.05cm}+{\rm j} \hspace{0.03cm}\cdot \hspace{0.03cm} \pi \hspace{0.03cm}\cdot \hspace{0.03cm} h \hspace{0.03cm}\cdot \hspace{0.03cm}t/T}\hspace{0.05cm},\hspace{0.2cm} 0 \le t \le T\hspace{0.05cm},</math>
 +
::<math> \xi_2(t) = \sqrt{1/T} \cdot {\rm e}^{\hspace{0.05cm}-{\rm j} \hspace{0.03cm}\cdot \hspace{0.03cm} \pi \hspace{0.03cm}\cdot \hspace{0.03cm} h \hspace{0.03cm}\cdot \hspace{0.03cm}t/T}\hspace{0.05cm},\hspace{0.2cm} 0 \le t \le T \hspace{0.05cm}</math>
 +
 
 +
To achieve orthogonality between these two complex basis functions,&nbsp; the&nbsp; [[Digital_Signal_Transmission/Carrier_Frequency_Systems_with_Coherent_Demodulation#Binary_frequency_shift_keying_.282.E2.80.93FSK.29|"modulation index"]]&nbsp; $h$&nbsp; must be chosen to be integer:
 +
 
 +
::<math>< \hspace{-0.05cm}\xi_1(t) \hspace{0.1cm}  \cdot  \hspace{0.1cm} \xi_2(t) \hspace{-0.05cm}> \hspace{0.2cm}= 0
 +
  \hspace{0.3cm} \Rightarrow \hspace{0.3cm} h = 2 \cdot \Delta f_{\rm A} \cdot T\hspace{0.05cm}= 1, 2, 3, \text{...}</math>
 +
 
 +
The diagram shows the structure for non-coherent orthogonal demodulation of binary FSK.<br>
 +
 
 +
[[File:EN_Dig_T_4_5_S3a_neu.png|right|frame|Non-coherent demodulation of the binary FSK|class=fit]]
 +
 
 +
In the noise-free case &nbsp; &#8658; &nbsp;  $n(t) \equiv 0$&nbsp; applies to the outputs of the two correlators:
  
*Die Rayleigh&ndash;WDF mit größerem <i>&sigma;<sub>n</sub></i> ist formgleich mit der gezeichneten Kurve (<i>&sigma;<sub>n</sub></i> = 0.5), jedoch im Verhältnis der Streuungen breiter und niedriger.<br>
+
::<math>r_1 = \hspace{0.2cm} < \hspace{-0.05cm}r(t) \hspace{0.1cm}  \cdot  \hspace{0.1cm} \xi_1(t) \hspace{-0.05cm}> \hspace{0.2cm}= 0\hspace{0.05cm},  \hspace{0.4cm} {\rm if}\hspace{0.15cm} m = m_1\hspace{0.05cm},</math>
 +
::<math> r_2 = \hspace{0.2cm} <  \hspace{-0.05cm}r(t) \hspace{0.1cm}  \cdot  \hspace{0.1cm} \xi_2(t) \hspace{-0.05cm}> \hspace{0.2cm}= 0\hspace{0.05cm}, \hspace{0.4cm} {\rm if}\hspace{0.15cm} m = m_0\hspace{0.05cm}.</math>
  
*<i>&sigma;<sub>n</sub></i> gibt die Streuungen der beiden gaußverteilten Zufallsgrößen <i>u</i> und <i>&upsilon;</i> an (beide haben gleiche Streuung) und nicht die Streuung der rayleighverteilten Zufallsgröße <i>y</i>. Für diese gilt vielmehr:
+
After respective magnitude formation &nbsp; &#8658; &nbsp; $y_1 = |r_1|, \ \ y_2 = |r_2|$&nbsp; the following decision rule is applicable:
 +
::<math>\hat{m} =
 +
\left\{ \begin{array}{c} m_0 \\
 +
m_1  \end{array} \right.\quad
 +
\begin{array}{*{1}c} {\rm if}\hspace{0.15cm} y_1 > y_2 \hspace{0.05cm},
 +
\\  {\rm if}\hspace{0.15cm} y_1 < y_2  \hspace{0.05cm}.\\ \end{array}</math>
  
::<math>\sigma_y = \sigma_n  \cdot  \sqrt{2 - {\pi}/{2 }} \hspace{0.2cm} \approx \hspace{0.2cm} 0.655 \cdot \sigma_n  
+
For a simpler realization of the decision,&nbsp; the difference&nbsp; $y_1 - y_2$&nbsp; can also be evaluated with the decision boundary&nbsp; $G = 0$.&nbsp; <br>
 +
 
 +
== Error probability with non-coherent 2&ndash;FSK demodulation==
 +
<br>
 +
In the following,&nbsp; the error probability is calculated under the assumption that&nbsp; $m = m_0$&nbsp; was sent.
 +
 
 +
*Under the further assumption of equally probable binary messages&nbsp; $m_0$&nbsp; and&nbsp; $m_1$,&nbsp; the absolute error probability is exactly the same:
 +
:$${\rm Pr}({\cal{E}}) =  {\rm Pr}({\cal{E}}\hspace{0.05cm} | \hspace{0.05cm}m_0)
 +
\hspace{0.05cm}.$$
 +
 
 +
*With&nbsp; $m = m_0$&nbsp; we get for the complex correlation output values&nbsp; $r_i$&nbsp; and their magnitudes&nbsp; $y_i$:<br>
 +
::<math>r_1 = \sqrt{E} \cdot {\rm e}^{{\rm j}\phi} + n_1\hspace{0.3cm} \Rightarrow \hspace{0.3cm}y_1 = |r_1|\hspace{0.15cm}{\rm is}\hspace{0.15cm}{\rm Rice}\hspace{0.15cm}{\rm distributed}
 +
\hspace{0.05cm},</math>
 +
::<math> r_2 =  n_2\hspace{0.3cm} \Rightarrow \hspace{0.3cm}y_2 = |r_2|\hspace{0.15cm}{\rm is}\hspace{0.15cm}{\rm Rayleigh}\hspace{0.15cm}{\rm distributed}
 +
\hspace{0.05cm}.</math>
 +
 
 +
:Here,&nbsp; $E$&nbsp; due to&nbsp; $M = 2$&nbsp; represents the&nbsp; "average symbol energy"&nbsp; $(E_{\rm S})$&nbsp; and the&nbsp; "average bit energy"&nbsp;  $(E_{\rm B})$&nbsp; equally.&nbsp; $n_1$&nbsp; and&nbsp; $n_2$&nbsp; are uncorrelated complex noise variables with mean zero and variance&nbsp; $2 \cdot \sigma_n^2$.
 +
 
 +
*Thus,&nbsp; the&nbsp; [[Theory_of_Stochastic_Signals/Two-Dimensional_Random_Variables#Joint_probability_density_function|"joint probability density function"]]&nbsp; is:
 +
::<math>p_{y_1,\hspace{0.03cm} y_2 \hspace{0.03cm}| \hspace{0.03cm}m} ( \eta_1, \eta_2 \hspace{0.05cm}| \hspace{0.05cm}m_0) =
 +
  p_{y_1 \hspace{0.01cm}| \hspace{0.03cm}m} ( \eta_1 \hspace{0.05cm}| \hspace{0.05cm}m_0) \cdot
 +
  p_{y_2 \hspace{0.03cm}| \hspace{0.03cm}m} ( \eta_2 \hspace{0.05cm}| \hspace{0.05cm}m_0)
 +
\hspace{0.05cm},</math>
 +
::$$p_{y_1 \hspace{0.01cm}| \hspace{0.03cm}m} ( \eta_1 \hspace{0.05cm}| \hspace{0.05cm}m_0) = {\eta_1}/{\sigma_n^2}
 +
\cdot {\rm e }^{ - ({\eta_1^2 + E})/({2 \sigma_n^2}) }\cdot {\rm I }_0 \left [{\eta_1 \cdot  \sqrt{E}}/{ \sigma_n^2}\right ]
 +
\hspace{0.05cm},$$
 +
::$$ p_{y_2 \hspace{0.01cm}| \hspace{0.03cm}m} ( \eta_2 \hspace{0.05cm}| \hspace{0.05cm}m_0) = {\eta_2}/{\sigma_n^2}
 +
\cdot {\rm e }^{  - \eta_2^2 /({2 \sigma_n^2}) } \hspace{0.05cm}.$$
 +
 
 +
*The error probability is generally obtained as follows:
 +
 
 +
::<math>{\rm Pr}({\cal{E}}) = \int_{0}^{\infty} \int_{\eta_1}^{\infty}
 +
  p_{y_1,\hspace{0.03cm} y_2 \hspace{0.03cm}| \hspace{0.03cm}m} ( \eta_1, \eta_2 \hspace{0.05cm}| \hspace{0.05cm}m_0) \,\,{\rm d} \eta_2\,\,{\rm d} \eta_1 = \int_{0}^{\infty} p_{y_1 \hspace{0.01cm}| \hspace{0.03cm}m} ( \eta_1 \hspace{0.05cm}| \hspace{0.05cm}m_0) \cdot \int_{\eta_1}^{\infty}
 +
  p_{y_2 \hspace{0.03cm}| \hspace{0.03cm}m} ( \eta_2 \hspace{0.05cm}| \hspace{0.05cm}m_0)\,\,{\rm d} \eta_2\,\,{\rm d} \eta_1
 +
\hspace{0.05cm}.</math>
 +
 
 +
{{BlaueBox|TEXT=  
 +
$\text{Conclusion:}$&nbsp; For the&nbsp; '''error probability with non-coherent demodulation of the binary FSK''',&nbsp; one obtains after some mathematical transformations the surprisingly simple result
 +
 
 +
::<math>p_{\rm S} = {\rm Pr}({\cal{E} }) = {1}/{2} \cdot {\rm e}^{-E_{\rm S}/(2 N_0)}
 
  \hspace{0.05cm}.</math>
 
  \hspace{0.05cm}.</math>
  
*Die Rayleighverteilung ist extrem unsymmetrisch, erkennbar am (relativ) großen Wert für das [http://en.lntwww.de/Stochastische_Signaltheorie/Erwartungswerte_und_Momente#Zentralmomente Zentralmoment 3. Ordnung:] <i>&mu;</i><sub>3</sub>/<i>&sigma;<sub>y</sub></i><sup>3</sup> &asymp; 0.27.<br>
+
*For comparison,&nbsp; the result for&nbsp; '''coherent demodulation'''&nbsp; is here given again:
 +
::<math>p_{\rm S} = {\rm Pr}({\cal{E} }) = {\rm Q}(\sqrt{ E_{\rm S}/N_0})
 +
\hspace{0.05cm}.</math>}}
 +
 
 +
 
 +
{{BlaueBox|TEXT= 
 +
$\text{Derivation:}$&nbsp; This result is now to be derived in some calculation steps.&nbsp; We start from the following equations:
 +
 
 +
:$${\rm Pr}({\cal{E} }) = \int_{0}^{\infty} p_{y_1 \hspace{0.01cm} \vert\hspace{0.03cm}m} ( \eta_1 \hspace{0.05cm} \vert \hspace{0.05cm}m_0) \cdot \int_{\eta_1}^{\infty}
 +
  p_{y_2 \hspace{0.03cm}\vert \hspace{0.03cm}m} ( \eta_2 \hspace{0.05cm}\vert \hspace{0.05cm}m_0)\,\,{\rm d} \eta_2\,\,{\rm d} \eta_1
 +
\hspace{0.05cm},\hspace{0.5cm}\text{with}$$
 +
::$$p_{y_1 \hspace{0.01cm}\vert \hspace{0.03cm}m} ( \eta_1 \hspace{0.05cm}\vert \hspace{0.05cm}m_0) = {\eta_1}/{\sigma_n^2} \cdot {\rm e }^{  - ({\eta_1^2 + E})/({2 \sigma_n^2}) }\cdot {\rm I }_0 \left [{\eta_1 \cdot  \sqrt{E} }/{ \sigma_n^2}\right ] \hspace{0.05cm},$$
 +
::$$p_{y_\hspace{0.01cm}\vert \hspace{0.03cm}m} ( \eta_2 \hspace{0.05cm}\vert \hspace{0.05cm}m_0) = {\eta_2}/{\sigma_n^2}  \cdot {\rm e }^{  - \eta_2^2 /({2 \sigma_n^2}) } \hspace{0.05cm}.$$
 +
 
 +
<b>(1)</b> &nbsp; The inner integral gives the probability that the Rayleigh distributed random variable&nbsp; $\eta_2$&nbsp; is larger than&nbsp; $\eta_1$&nbsp; &ndash; see solution to&nbsp; [[Aufgaben:Exercise_4.17Z:_Rayleigh_and_Rice_Distribution|"Exercise 4.17Z"]]:
 +
 
 +
:$$\int_{\eta_1}^{\infty}  p_{y_2 \hspace{0.03cm}\vert\hspace{0.03cm}m} ( \eta_2 \hspace{0.05cm}\vert \hspace{0.05cm}m_0)\,\,{\rm d} \eta_2 =
 +
{\rm e }^{  - \eta_1^2 /({2 \sigma_n^2}) } \hspace{0.3cm}
 +
\Rightarrow \hspace{0.3cm} {\rm Pr}({\cal{E} }) = \int_{0}^{\infty}{\eta_1}/{\sigma_n^2}
 +
\cdot {\rm e }^{  - ({2\eta_1^2 + E})/({2 \sigma_n^2}) }\cdot {\rm I }_0
 +
\left [ {\eta_1 \cdot \sqrt{E} }/{ \sigma_n^2}\right ]\,\,{\rm d} \eta_1
 +
\hspace{0.05cm}.$$
 +
 
 +
<b>(2)</b> &nbsp; With the&nbsp; (arbitrary)&nbsp; substitutions&nbsp; $C_0^2  = E/4$&nbsp; and&nbsp; $\sigma_0^2 = \sigma_n^2/2$,&nbsp; we obtain:
  
*Die Riceverteilung ist um so symmetrischer, je größer das Verhältnis <i>C</i>/<i>&sigma;<sub>n</sub></i> von deterministischer und stochastischer Komponente ist. Für <i>C</i>/<i>&sigma;<sub>n</sub></i> &#8805; 4 ist <i>&mu;</i><sub>3</sub> nahezu 0.<br>
+
:$${\rm Pr}({\cal{E} }) = \int_{0}^{\infty} \frac{\eta_1}{2 \cdot \sigma_0^2}
 +
\cdot {\rm exp } \left [ - \frac{2 \eta_1^2 + 4 C_0^2}{4 \sigma_0^2}\right ] \cdot {\rm I }_0
 +
\left [ \frac{\eta_1 \cdot  2C_0}{ 2 \sigma_0^2}\right ]\,\,{\rm d} \eta_1 = \int_{0}^{\infty} \frac{\eta_1}{2 \cdot \sigma_0^2}
 +
\cdot {\rm exp } \left [ - \frac{\eta_1^2 + 2 C_0^2}{2 \sigma_0^2}\right ] \cdot {\rm I }_0
 +
\left [ \frac{\eta_1 \cdot  C_0}{  \sigma_0^2}\right ]\,\,{\rm d} \eta_1
 +
\hspace{0.05cm}.$$
  
*Weiterhin ist zu erkennen, dass sich die Riceverteilung (mit den Parametern <i>C</i> und <i>&sigma;<sub>n</sub></i>) immer mehr einer Gaußverteilung mit Mittelwert <i>C</i> und Streuung <i>&sigma;<sub>n</sub></i> annähert, je größer der Quotient <i>C</i>/<i>&sigma;<sub>n</sub></i> ist:
+
<b>(3)</b> &nbsp; By shifting fractions in front of the integral,&nbsp; we succeed that the integrand again describes a&nbsp; [[Theory_of_Stochastic_Signals/Further_Distributions#Rice_PDF|"Rice distribution"]]:&nbsp;
  
::<math>p_y (\eta) \approx \frac{1}{\sqrt{2\pi} \cdot \sigma_n}
+
::<math>{\rm Pr}({\cal{E} }) = {1}/{2} \cdot {\rm exp } \left [ - \frac{ C_0^2}{2 \sigma_0^2}\right ] \cdot \int\limits_{0}^{\infty} \frac{\eta_1}{ \sigma_0^2}
  \cdot {\rm exp } \left [ - \frac{(\eta - C)^2}{2 \sigma_n^2}\right ]
+
  \cdot {\rm exp } \left [ - \frac{\eta_1^2 +  C_0^2}{2 \sigma_0^2}\right ] \cdot {\rm I }_0
  \hspace{0.3cm} \Rightarrow \hspace{0.3cm} m_y = C\hspace{0.05cm},\hspace{0.2cm}\sigma_y = \sigma_n
+
  \left [ \frac{\eta_1 \cdot  C_0}{ \sigma_0^2}\right ]\,\,{\rm d} \eta_1
 
  \hspace{0.05cm}.</math>
 
  \hspace{0.05cm}.</math>
  
[[File:P ID2086 Dig T 4 5 S1 version1.png|Rayleigh– und Rice–WDF|class=fit]]<br>
+
<b>(4)</b> &nbsp; The integrand now describes the Rice PDF.&nbsp; The integral over the entire definition area from&nbsp; $0$&nbsp; to &nbsp;$+\infty$&nbsp; yields the value one,&nbsp; as for any PDF,&nbsp; so that holds:
 +
 
 +
::<math>{\rm Pr}({\cal{E} }) = {1}/{2} \cdot {\rm exp } \left [ - \frac{ C_0^2}{2 \sigma_0^2}\right ]\hspace{0.05cm}.</math>
  
== Nichtkohärente Demodulation von On–Off–Keying (1) ==
+
<b>(5)</b> &nbsp; Finally, with&nbsp; $C_0^2  = E/4$&nbsp; and&nbsp; $\sigma_0^2 = \sigma_n^2/2$&nbsp; and the generally valid relation&nbsp; $\sigma_n^2 = N_0$,&nbsp; we obtain:
 +
:$${\rm Pr}({\cal{E} })
 +
= {1}/{2} \cdot {\rm exp } \left [ - \frac{ C_0^2}{2 \sigma_0^2}\right ] = {1}/{2} \cdot {\rm exp } \left [ - \frac{ E_{\rm S}/4}{N_{\rm 0}/2}\right ] \hspace{0.3cm}\Rightarrow \hspace{0.3cm} p_{\rm S} = {\rm Pr}({\cal{E} }) = {1}/{2} \cdot {\rm e}^{-E_{\rm S}/(2 N_0)}\hspace{0.05cm}.$$
 +
&nbsp; &nbsp; &nbsp; &nbsp; $E_{\rm S}$&nbsp; here indicates the average signal energy per symbol, which is equal to the normalization energy&nbsp; $E$&nbsp; for FSK.
 +
<div align="right">$\text{q.e.d.}$</div>}}<br>
 +
 
 +
{{GraueBox|TEXT=
 +
[[File:EN_Dig_T_4_5_S3b.png|right|frame|FSK symbol error probability for coherent and non-coherent demodulation|class=fit]] 
 +
$\text{Example 2:}$&nbsp; The graph compares the error probability curves of both demodulation methods as a function of the AWGN quotient&nbsp; $E_{\rm S}/N_0$.&nbsp;
 +
<br><br><br><br>
 +
One can see:
 +
#The non-coherent BFSK&nbsp; (red curve)&nbsp; requires a &nbsp; $0.8 \ \rm dB$ &nbsp; larger&nbsp;$E_{\rm S}/N_0$&nbsp; compared to the coherent BFSK&nbsp; (blue curve)&nbsp; at&nbsp; $p_{\rm S}= 10^{-5}$.&nbsp;<br><br> 
 +
#At &nbsp; $p_{\rm S}= 10^{-3}$ &nbsp; the distance is even&nbsp; $1.3 \ \rm dB$.<br><br>
 +
#In contrast,&nbsp; the distance between the coherent binary FSK from the coherent BPSK is always&nbsp; $1.3 \ \rm dB$&nbsp; regardless of the error probability.
 +
<br clear =all>}}
 +
 
 +
 
 +
== Non-coherent demodulation of multi-level FSK==
 
<br>
 
<br>
Wir betrachten [http://en.lntwww.de/Digitalsignal%C3%BCbertragung/Tr%C3%A4gerfrequenzsysteme_mit_koh%C3%A4renter_Demodulation#On.E2.80.93Off.E2.80.93Keying_bzw._2.E2.80.93ASK On&ndash;Off&ndash;Keying] im äquivalenten Tiefpassbereich. Bei kohärenter Demodulation (linke Grafik) ist die Signalraumkonstellation des Empfangssignals gleich der des Sendesignals und besteht aus zwei Punkten. Die Entscheidungsgrenze <i>G</i> liegt in der Mitte zwischen diesen Punkten <b><i>r</i><sub>0</sub></b> und <b><i>r</i><sub>1</sub></b>. Die Pfeile markieren Rauschvektoren, die eventuell zu Übertragungsfehlern führen.<br>
+
[[File:P ID2089 Dig T 4 5 S4a version1.png|right|frame|Orthogonal $M$&ndash;level FSK for &nbsp;$M= 3$]]
 +
 
 +
We now consider the message set&nbsp; $\{m_1, m_2,\hspace{0.05cm}\text{ ...} \hspace{0.05cm}, m_{M}\}$&nbsp; and denote&nbsp; $M$&nbsp; as the&nbsp; "level number".
 +
#As in the case of binary FSK,&nbsp; a prerequisite for the application of the modulation process&nbsp; "Frequency Shift Keying"&nbsp; and at the same time of a non-coherent demodulator is an integer modulation index&nbsp; $h$.<br><br>
 +
#In this case,&nbsp; the&nbsp; $M$&ndash;level FSK is orthogonal and a signal space constellation results as shown in the adjacent diagram for the special case&nbsp; $M = 3$.&nbsp; <br><br>
 +
 
 +
The non-coherent demodulator is sketched below.  
 +
*Compared to the&nbsp; [[Digital_Signal_Transmission/Carrier_Frequency_Systems_with_Non-Coherent_Demodulation#Non-coherent_demodulation_of_binary_FSK_.282.E2.80.93FSK.29|"receiver structure for binary FSK"]],&nbsp; this receiver differs only by&nbsp; $M$&nbsp; branches instead of only two,&nbsp; which provide the comparison values&nbsp; $y_1$,&nbsp; $y_2$, ... , $y_M$.<br>
 +
 
 +
*To calculate the error probability,&nbsp; we assume that&nbsp; $m_1$&nbsp; was sent. This means that the decision is correct if the largest detection output value is&nbsp; $y_1$:&nbsp;
 +
 
 +
::<math>{\rm Pr}({\cal{C}}) \hspace{-0.1cm}  =  \hspace{-0.1cm} {\rm Pr} \big [ (y_2 < y_1) \cap (y_3 < y_1) \cap ... \cap (y_{M} < y_1)
 +
  \hspace{0.05cm}| \hspace{0.05cm} m = m_1\big ] =  {\rm Pr} \left [ \hspace{0.1cm} \bigcap\limits_{k = 2}^M (y_k < y_1) 
 +
  \hspace{0.05cm}| \hspace{0.05cm}m = m_1\right ]
 +
  \hspace{0.01cm}.</math>
 +
[[File:EN_Dig_T_4_5_S4b.png|right|frame|Non-coherent receiver structure for &nbsp;$M$&ndash;level FSK|class=fit]]
 +
<br clear=all>
 +
{{BlaueBox|TEXT= 
 +
$\text{Conclusions:}$&nbsp;
 +
* The &nbsp; '''error probability of M&ndash;level FSK in non-coherent demodulation''' &nbsp; is equal to&nbsp; ${\rm Pr}({\cal{E} }) = 1 - {\rm Pr}({\cal{C} })$.
 +
 +
*In the following,&nbsp; it is shown that&nbsp; ${\rm Pr}({\cal{C} })$&nbsp; can be represented as follows:
 +
 
 +
::<math>{\rm Pr}({\cal{C} }) = \sum_{i = 0}^{M-1} (-1)^i \cdot {M-1 \choose i }  \cdot \frac{1}{i+1} \cdot {\rm exp } \left [ - \frac{i \cdot E_{\rm S} }
 +
  {(i+1) \cdot N_0}\right ]  
 +
  \hspace{0.05cm}\hspace{0.5cm} \text{with}\hspace{0.5cm}  E_{\rm S} = E_{\rm B} \cdot {\rm log_2}(M)\hspace{0.05cm}.</math>
 +
 
 +
*In the special case&nbsp; $M = 2$,&nbsp; of course, the result obtained in the&nbsp; [[Digital_Signal_Transmission/Carrier_Frequency_Systems_with_Coherent_Demodulation|"last section"]]&nbsp; is again obtained:
 +
 
 +
::<math>p_{\rm S} = {\rm Pr}({\cal{E} }) = {1}/{2} \cdot {\rm e}^{-E_{\rm S}/(2 N_0)}
 +
\hspace{0.05cm}.</math>}}
 +
 
 +
 
 +
{{BlaueBox|TEXT= 
 +
$\text{Derivation:}$&nbsp; This result shall now be derived in some calculation steps.
 +
*We assume that&nbsp; $m_1$&nbsp; has been sent.&nbsp; This means that the decision is correct if the largest detection output value is&nbsp; $y_1$.&nbsp;
 +
 
 +
*Otherwise,&nbsp; there are certain analogies to the derivation of the&nbsp; [[Digital_Signal_Transmission/Carrier_Frequency_Systems_with_Non-Coherent_Demodulation#Error_probability_with_non-coherent_2.E2.80.93FSK_demodulation|"BFSK error probability"]].<br>
 +
 
 +
 
 +
<b>(1)</b> &nbsp; Using the conditional probability density&nbsp; $p_{y_1 \hspace{0.01cm} \vert\hspace{0.03cm}m} ( \eta_1 \hspace{0.05cm} \vert \hspace{0.05cm}m_1)$,&nbsp; we obtain:
 +
 
 +
::<math>{\rm Pr}({\cal{C} }) =  \int_{0}^{\infty} {\rm Pr} \left [ \hspace{0.1cm} \bigcap\limits_{k = 2}^M (y_k < y_1) 
 +
  \hspace{0.05cm}\vert\hspace{0.05cm}y_1 = \eta_1, m = m_1\right ] \cdot p_{y_1 \hspace{0.01cm} \vert\hspace{0.03cm}m} ( \eta_1 \hspace{0.05cm} \vert \hspace{0.05cm}m_1) \,\,{\rm d} \eta_1
 +
  \hspace{0.05cm}.</math>
 +
 
 +
<b>(2)</b> &nbsp; The decision values&nbsp; $y_2$,&nbsp; $y_3$,&nbsp; ... , $y_M$&nbsp; are statistically independent for given&nbsp; $y_1$.&nbsp; Therefore:
 +
 
 +
::<math>{\rm Pr}({\cal{C} }) =  \int_{0}^{\infty}  \left \{ {\rm Pr} \big [ (y_2 < y_1) 
 +
  \hspace{0.05cm}\vert \hspace{0.05cm}y_1 = \eta_1, m = m_1\big ] \right \}^{M-1} \cdot p_{y_1 \hspace{0.01cm} \vert\hspace{0.03cm}m} ( \eta_1 \hspace{0.05cm} \vert \hspace{0.05cm}m_1) \,\,{\rm d} \eta_1
 +
  \hspace{0.05cm}.</math>
 +
 
 +
<b>(3)</b> &nbsp; The selected value&nbsp; $y_2$&nbsp; conditioned on&nbsp; $m_1$&nbsp; has a Rayleigh distribution with parameter&nbsp; $\sigma_n^2$:
 +
 
 +
::<math>{\rm Pr} \big [ (y_2 < y_1) 
 +
  \hspace{0.05cm} \vert \hspace{0.05cm}y_1 = \eta_1, m = m_1\big ] \hspace{-0.1cm} = \hspace{-0.1cm} \int_{0}^{\eta_1}  p_{y_2 \hspace{0.01cm} \vert\hspace{0.03cm}m} ( \eta_2 \hspace{0.05cm} \vert \hspace{0.05cm}m_1)\,\,{\rm d} \eta_2=  1 - {\rm exp } \big [ - {\eta_1^2 }/({2 \sigma_n^2})\big ] = 1 - a  \hspace{0.2cm}{\rm(abbreviation)}
 +
  \hspace{0.05cm}.</math>
  
[[File:P ID2083 Dig T 4 5 S2a version1.png|Kohärente und nichtkohärente Demodulation von On-Off-Keying|class=fit]]<br>
+
<b>(4)</b> &nbsp; Now we are looking for the expression &nbsp; $(1 -a)^{M-1}$,&nbsp; for which with the abbreviation from&nbsp; <b>(3)</b>&nbsp; holds:
  
Dagegen gilt bei nichtkohärenter Demodulation:
+
::<math> (1-a)^{M-1} \hspace{-0.1cm}  =  \hspace{-0.1cm} \sum_{i = 0}^{M-1} (-1)^i \cdot {M-1 \choose i } \cdot (-1)^i  \cdot a^i = \sum_{i = 0}^{M-1} (-1)^i \cdot {M-1 \choose i }  \cdot {\rm exp } \left [ - \frac{i \cdot \eta_1^2 }{2 \sigma_n^2}\right ]
*Der Punkt <i><b>r</b></i><sub>1</sub> = <i><b>s</b></i><sub>1</sub> = 0 bleibt weiter erhalten.<br>
+
  \hspace{0.05cm}.</math>
  
*Dagegen kann <i><b>r</b></i><sub>0</sub> = <i><b>s</b></i><sub>0</sub> &middot; exp (j<i>&#981;</i>) auf jeden Punkt des Kreises um <i><b>s</b></i><sub>0</sub> liegen, da <i>&#981;</i> unbekannt ist.<br>
+
<b>(5)</b> &nbsp; $y_1$&nbsp; has a&nbsp; [[Digital_Signal_Transmission/Carrier_Frequency_Systems_with_Non-Coherent_Demodulation#Rayleigh_and_Rice_Distribution|"Rice distribution"]] conditioned on&nbsp; $m=m_1$.&nbsp; Thus,&nbsp; the probability of a correct decision can be expressed in the following form:
  
*Der Entscheidungsprozess unter Berücksichtigung des AWGN&ndash;Rauschens ist nun 2&ndash;dimensional zu interpretieren, wie es durch die Pfeile in der rechten Grafik angedeutet ist.<br>
+
::<math>{\rm Pr}({\cal{C} }) \hspace{-0.1cm}  =  \hspace{-0.1cm} \sum_{i = 0}^{M-1} (-1)^i \cdot {M-1 \choose i }  \cdot \int_{0}^{\infty} {\rm exp } \left [ - \frac{i \cdot \eta_1^2 }{2 \sigma_n^2}\right ] \cdot \frac{\eta_1}{ \sigma_n^2}\cdot 
 +
  {\rm exp } \left [ - \frac{\eta_1^2 +  E_{\rm S} }{2 \sigma_n^2}\right ] \cdot {\rm I }_0
 +
\left [ \frac{\eta_1 \cdot  \sqrt{E_{\rm S} } }{  \sigma_n^2}\right ] \,\,{\rm d} \eta_1
 +
  \hspace{0.05cm}.</math>
  
*Das Entscheidungsgebiet <i>I</i><sub>1</sub> ist nun ein Kreis, dessen Radius <i>G</i> ein optimierbarer Parameter ist. Das Entscheidungsgebiet <i>I</i><sub>0</sub> liegt außerhalb dieses Kreises.<br><br>
+
<b>(6)</b> &nbsp; Substitutions succeed in shaping the integrand according to the Rice distribution.&nbsp; Since each probability density integrates to one,&nbsp; we obtain:
  
Damit liegt die Strukur des optimalen OOK&ndash;Empfängers (im äquivalenten Tiefpassbereich) fest.<br>
+
::<math>{\rm Pr}({\cal{C} }) = \sum_{i = 0}^{M-1} (-1)^i \cdot {M-1 \choose i }  \cdot \frac{1}{i+1} \cdot {\rm exp } \left [ - \frac{i \cdot E_{\rm S} }
 +
  {(i+1) \cdot N_0}\right ]
 +
  \hspace{0.05cm}\hspace{0.5cm} \text{with}\hspace{0.5cm}  E_{\rm S} = E_{\rm B} \cdot {\rm log_2}(M)\hspace{0.05cm}.</math>
 +
 
 +
<b>(7)</b> &nbsp; The special case&nbsp; $M = 2$&nbsp; leads to exactly the same result as calculated for the binary FSK:
 +
 
 +
::<math>{\rm Pr}({\cal{C} }) = (-1)^0 \cdot {2-1 \choose 0 }  \cdot \frac{1}{0+1} \cdot {\rm exp } \left [ - \frac{0 \cdot E_{\rm S} }
 +
  {(i+1) \cdot N_0}\right ] + (-1)^1 \cdot {2-1 \choose 1 }  \cdot \frac{1}{1+1} \cdot {\rm exp } \left [ - \frac{1 \cdot E_{\rm S} }
 +
  {(i+1) \cdot N_0}\right ] </math>
 +
::<math>  \Rightarrow \hspace{0.5cm} {\rm Pr}({\cal{C} }) = 1 - {1}/{2} \cdot {\rm e }^{-E_{\rm S}/(2N_0)}
 +
\hspace{0.5cm} \Rightarrow \hspace{0.5cm} {\rm Pr}({\cal{E} }) = 1 - {\rm Pr}({\cal{C} })  = {1}/{2} \cdot {\rm e }^{-E_{\rm S}/(2N_0)}
 +
  \hspace{0.05cm}.</math>
 +
<div align="right">$\text{q.e.d.}$</div>}}
 +
 
 +
== Exercises for the chapter==
 +
<br>
 +
[[Aufgaben:Exercise_4.17:_Non-Coherent_On-Off_Keying|Exercise 4.17: Non-Coherent On-Off Keying]]
  
[[File:P ID3147 Dig T 4 5 S2b version1.png|Empfänger für nichtkohärente OOK-Demodulation|class=fit]]<br>
+
[[Aufgaben:Exercise_4.17Z:_Rayleigh_and_Rice_Distribution|Exercise 4.17Z: Rayleigh and Rice Distribution]]
  
 +
[[Aufgaben:Exercise_4.18:_Non-Coherent_FSK_Demodulation|Exercise 4.18: Non-Coherent FSK Demodulation]]
  
 +
[[Aufgaben:Exercise_4.18Z:_BER_of_Coherent_and_Non-Coherent_FSK|Exercise 4.18Z: BER of Coherent and Non-Coherent FSK]]
  
 +
[[Aufgaben:Exercise_4.19:_Orthogonal_Multilevel_FSK|Exercise 4.19: Orthogonal Multilevel FSK]]
  
 
{{Display}}
 
{{Display}}

Latest revision as of 13:43, 10 October 2022

Rayleigh and Rice Distribution


The estimation of the phase angle from the incoming signal,  which is required for coherent demodulation,  is not possible or only possible to a limited extent in many applications.  For example,  the movement of a mobile subscriber at high speed leads to very rapid temporal changes in the phase angle  $\phi$,  which makes its sufficiently accurate determination difficult or even impossible.

This fact leads to the  non-coherent demodulation processes  with the advantage of reduced complexity,  but with increased error probability.  In the derivation of the equations one encounters two probability density functions,  which are given here in advance:

  • The  Rayleigh distribution  is obtained for the power density function  $\rm (PDF)$  of the random variable  $y$  with realization  $\eta$,  which is obtained from the two Gaussian distributed and statistically independent components  $u$  and  $v$  $($both with the same standard deviation  $\sigma_n)$  as follows:
\[y = \sqrt{u^2 + v^2} \hspace{0.3cm} \Rightarrow \hspace{0.3cm} p_y (\eta) ={\eta}/{\sigma_n^2} \cdot {\rm exp } \left [ - {\eta^2}/{ (2\sigma_n^2)}\right ] \hspace{0.05cm}.\]
  • The  Rice distribution  is obtained under the same boundary conditions for the case where a constant   $C$  is added to one of the components $($either to $u$  or  $v)$: 
\[y = \sqrt{(u+C)^2 + v^2} \hspace{0.3cm} \Rightarrow \hspace{0.3cm} p_y (\eta) = {\eta}/{\sigma_n^2} \cdot {\rm exp } \left [ - ({\eta^2 + C^2})/(2 \sigma_n^2) \right ] \cdot {\rm I }_0 \left [{\eta \cdot C}/{ \sigma_n^2}\right ] \hspace{0.05cm}.\]
The Rice distribution uses the  "modified zero-order Bessel function",  whose definition and series expansion are as follows:
\[{\rm I }_0 (x) = \frac{1}{ 2\pi} \cdot \int_{0}^{2\pi} {\rm e }^{-x \hspace{0.03cm}\cdot \hspace{0.03cm}\cos(\alpha)} \,{\rm d} \alpha \hspace{0.2cm} \approx \hspace{0.2cm} \sum_{k = 0}^{\infty} \frac{(x/2)^{2k}}{k! \cdot \Gamma (k+1)} \hspace{0.05cm}.\]

The graph shows the probability density functions of Rayleigh and Rice distributions.  To this representation it is to be noted:

Rayleigh resp. Rice PDF for $σ_n = 0.5$
  • The Rice distribution is determined by the two parameters  $C$  and  $\sigma_n$.  With  $C = 0$,  the Rice PDF is identical to the Rayleigh PDF.
  • The Rayleigh PDF with larger  $\sigma_n$  is of the same shape as shown for  $\sigma_n = 0.5$,  but wider and lower in the ratio of standard deviation.
  • $\sigma_n$  indicates the standard deviation of the two Gaussian distributed random variables  $u$  and  $v$  and not the standard deviation of the Rayleigh distributed random variable  $y$.  Rather,  the following applies to the latter:
\[\sigma_y = \sigma_n \cdot \sqrt{2 - {\pi}/{2 }} \hspace{0.2cm} \approx \hspace{0.2cm} 0.655 \cdot \sigma_n \hspace{0.05cm}.\]
  • The Rayleigh PDF is extremely asymmetric,  recognizable by the large value for the  "third order central moment"   ⇒   "Charlier's skewness":
$$\mu_3/\sigma_y \approx 0.27.$$
  • The Rice PDF is more symmetrical the larger the quotient  $C/\sigma_n$.  For  $C/\sigma_n \ge 4$:  $\mu_3 \approx 0$   ⇒   symmetrical PDF.
  • The larger  $C/\sigma_n$  is,  the more the Rice PDF $($with  $C$,  $\sigma_n)$  approaches a  "Gaussian PDF"  with mean  $C$  and standard deviation  $\sigma_n$:
$$p_y (\eta) \approx \frac{1}{\sqrt{2\pi} \cdot \sigma_n} \cdot {\rm exp } \left [ - \frac{(\eta - C)^2}{2 \sigma_n^2}\right ]$$
$$ \Rightarrow \hspace{0.3cm} m_y = C\hspace{0.05cm},\hspace{0.2cm}\sigma_y = \sigma_n \hspace{0.05cm}.$$

Non-coherent demodulation of  "on–off keying"  (OOK)


We consider  "on–off keying"  $(\rm OOK$  or  $\rm 2–ASK)$  in the equivalent low-pass range.

Coherent and non-coherent demodulation of  "on-off keying"
  • In the case of coherent demodulation  (left graph),  the signal space constellation of the received signal  $\boldsymbol{r}(t)$  is the same as that of the transmitted signal  $\boldsymbol{s}(t)$  and consists of two points.
  • The decision boundary  $G$  lies in the middle between these points  $\boldsymbol{r}_0$  and  $\boldsymbol{r}_1$.
  • The arrows mark the rough direction of noise vectors that may cause transmission errors.


On the other hand, for non-coherent demodulation (right graph):

  • The point  $\boldsymbol{r}_1 = \boldsymbol{s}_1 = 0$  is still preserved.  But  $\boldsymbol{r}_0 = \boldsymbol{s}_0 \cdot {\rm e}^{\hspace{0.02cm}{\rm j}\hspace{0.03cm}\phi}$  can lie on any point of the circle around  $\boldsymbol{s}_0$,  since  $\phi$  is unknown.
  • The decision process in consideration of the AWGN noise happens in two dimensions  $($indicated by arrows in the right graph$)$.
  • The decision area  $I_1$  is a circle whose radius  $G$  is an optimizable parameter. The decision area  $I_0$  lies outside the circle.


Thus,  the structure of the optimal OOK receiver  (in the equivalent low–pass range)  is fixed   ⇒   see second graph:

Receiver for non-coherent OOK demodulation  (complex signals are labeled blue)
  • The input signal  $\boldsymbol{r}(t) = \boldsymbol{s}(t) \cdot {\rm e}^{\hspace{0.02cm}{\rm j}\hspace{0.03cm}\phi} + \boldsymbol{n}(t)$  is generally complex because of the phase angle  $\phi$  and because of the complex noise term  $\boldsymbol{n}(t)$. 
  • Consequently required is now the correlation between the complex received signal  $\boldsymbol{r}(t)$  and a complex basis function  $\boldsymbol{\xi}_1(t)$.
  • The result is the  (complex)  detection value  $\boldsymbol{r}$,  from which the magnitude  $y = |\boldsymbol{r}(t)|$  is formed as a real decision input variable.
  • If the decision value  $y \gt G$,  then  $m_0$  is output as the estimated value,  otherwise  $m_1$.
  • Thus,  the symbol error probability for equally probable symbols
    ⇒   ${\rm Pr}(\boldsymbol{m}_0) = {\rm Pr}(\boldsymbol{m}_1) =1/2$  is:
\[p_{\rm S} = {\rm Pr}({\cal{E}}) = {1}/{ 2} \cdot \int_{0}^{G} p_{y\hspace{0.05cm}\vert \hspace{0.05cm}m} (\eta \hspace{0.05cm}|\hspace{0.05cm} m_0) \,{\rm d} \eta + {1}/{ 2} \cdot \int_{G}^{\infty} p_{y\hspace{0.05cm}\vert \hspace{0.05cm}m} (\eta \hspace{0.05cm}\vert \hspace{0.05cm} m_1) \,{\rm d} \eta \hspace{0.05cm}.\]
  • However,  due to the Rice PDF   $p_{y\hspace{0.05cm}|\hspace{0.05cm}m} (\eta\hspace{0.05cm}|\hspace{0.05cm}m_0)$   and the Rayleigh PDF   $p_{y\hspace{0.05cm}|\hspace{0.05cm}m} (\eta\hspace{0.05cm}|\hspace{0.05cm}m_1)$,   this probability can only be determined numerically.
  • The optimal decision limit  $G$  has to be determined beforehand as the solution of the following equation:
\[p_{y\hspace{0.05cm}\vert \hspace{0.05cm}m} (G \hspace{0.05cm}|\hspace{0.05cm}m_0) = p_{y\hspace{0.05cm}|\hspace{0.05cm}m} (G \hspace{0.05cm}|\hspace{0.05cm}m_1) \hspace{0.05cm}.\]
Density functions for "OOK, non-coherent"

$\text{Example 1:}$  The graph shows the result of this equation for  $\sigma_n = 0.5$  and  $C = 2$, where the (red) Rice PDF is approximated by a Gaussian PDF with mean  $C$  and standard deviation  $\sigma_n$. 

One can see from this sketch:

  • The optimal decision value  $($here:   $G \approx 1.25)$  results from the intersection of the two PDF curves.
  • The symbol error probability  $p_{\rm S}$  is the sum of the two colored areas.  In the example:  $p_{\rm S} \approx 5\%$ results.


⇒   You can determine the error probability for other values of  $C$  and  $\sigma_n$  as well as the optimal decision value  $G$  using the  HTML5/JavaScript applet  "Coherent and Non-coherent On-Off Keying"

Non-coherent demodulation of binary FSK (2–FSK)


As already shown in the  "last chapter",  "Binary Frequency Shift Keying"  $\rm (2–FSK)$  in the equivalent low-pass range can be described by the basis functions

\[\xi_1(t) = \sqrt{1/T} \cdot {\rm e}^{\hspace{0.05cm}+{\rm j} \hspace{0.03cm}\cdot \hspace{0.03cm} \pi \hspace{0.03cm}\cdot \hspace{0.03cm} h \hspace{0.03cm}\cdot \hspace{0.03cm}t/T}\hspace{0.05cm},\hspace{0.2cm} 0 \le t \le T\hspace{0.05cm},\]
\[ \xi_2(t) = \sqrt{1/T} \cdot {\rm e}^{\hspace{0.05cm}-{\rm j} \hspace{0.03cm}\cdot \hspace{0.03cm} \pi \hspace{0.03cm}\cdot \hspace{0.03cm} h \hspace{0.03cm}\cdot \hspace{0.03cm}t/T}\hspace{0.05cm},\hspace{0.2cm} 0 \le t \le T \hspace{0.05cm}\]

To achieve orthogonality between these two complex basis functions,  the  "modulation index"  $h$  must be chosen to be integer:

\[< \hspace{-0.05cm}\xi_1(t) \hspace{0.1cm} \cdot \hspace{0.1cm} \xi_2(t) \hspace{-0.05cm}> \hspace{0.2cm}= 0 \hspace{0.3cm} \Rightarrow \hspace{0.3cm} h = 2 \cdot \Delta f_{\rm A} \cdot T\hspace{0.05cm}= 1, 2, 3, \text{...}\]

The diagram shows the structure for non-coherent orthogonal demodulation of binary FSK.

Non-coherent demodulation of the binary FSK

In the noise-free case   ⇒   $n(t) \equiv 0$  applies to the outputs of the two correlators:

\[r_1 = \hspace{0.2cm} < \hspace{-0.05cm}r(t) \hspace{0.1cm} \cdot \hspace{0.1cm} \xi_1(t) \hspace{-0.05cm}> \hspace{0.2cm}= 0\hspace{0.05cm}, \hspace{0.4cm} {\rm if}\hspace{0.15cm} m = m_1\hspace{0.05cm},\]
\[ r_2 = \hspace{0.2cm} < \hspace{-0.05cm}r(t) \hspace{0.1cm} \cdot \hspace{0.1cm} \xi_2(t) \hspace{-0.05cm}> \hspace{0.2cm}= 0\hspace{0.05cm}, \hspace{0.4cm} {\rm if}\hspace{0.15cm} m = m_0\hspace{0.05cm}.\]

After respective magnitude formation   ⇒   $y_1 = |r_1|, \ \ y_2 = |r_2|$  the following decision rule is applicable:

\[\hat{m} = \left\{ \begin{array}{c} m_0 \\ m_1 \end{array} \right.\quad \begin{array}{*{1}c} {\rm if}\hspace{0.15cm} y_1 > y_2 \hspace{0.05cm}, \\ {\rm if}\hspace{0.15cm} y_1 < y_2 \hspace{0.05cm}.\\ \end{array}\]

For a simpler realization of the decision,  the difference  $y_1 - y_2$  can also be evaluated with the decision boundary  $G = 0$. 

Error probability with non-coherent 2–FSK demodulation


In the following,  the error probability is calculated under the assumption that  $m = m_0$  was sent.

  • Under the further assumption of equally probable binary messages  $m_0$  and  $m_1$,  the absolute error probability is exactly the same:
$${\rm Pr}({\cal{E}}) = {\rm Pr}({\cal{E}}\hspace{0.05cm} | \hspace{0.05cm}m_0) \hspace{0.05cm}.$$
  • With  $m = m_0$  we get for the complex correlation output values  $r_i$  and their magnitudes  $y_i$:
\[r_1 = \sqrt{E} \cdot {\rm e}^{{\rm j}\phi} + n_1\hspace{0.3cm} \Rightarrow \hspace{0.3cm}y_1 = |r_1|\hspace{0.15cm}{\rm is}\hspace{0.15cm}{\rm Rice}\hspace{0.15cm}{\rm distributed} \hspace{0.05cm},\]
\[ r_2 = n_2\hspace{0.3cm} \Rightarrow \hspace{0.3cm}y_2 = |r_2|\hspace{0.15cm}{\rm is}\hspace{0.15cm}{\rm Rayleigh}\hspace{0.15cm}{\rm distributed} \hspace{0.05cm}.\]
Here,  $E$  due to  $M = 2$  represents the  "average symbol energy"  $(E_{\rm S})$  and the  "average bit energy"  $(E_{\rm B})$  equally.  $n_1$  and  $n_2$  are uncorrelated complex noise variables with mean zero and variance  $2 \cdot \sigma_n^2$.
\[p_{y_1,\hspace{0.03cm} y_2 \hspace{0.03cm}| \hspace{0.03cm}m} ( \eta_1, \eta_2 \hspace{0.05cm}| \hspace{0.05cm}m_0) = p_{y_1 \hspace{0.01cm}| \hspace{0.03cm}m} ( \eta_1 \hspace{0.05cm}| \hspace{0.05cm}m_0) \cdot p_{y_2 \hspace{0.03cm}| \hspace{0.03cm}m} ( \eta_2 \hspace{0.05cm}| \hspace{0.05cm}m_0) \hspace{0.05cm},\]
$$p_{y_1 \hspace{0.01cm}| \hspace{0.03cm}m} ( \eta_1 \hspace{0.05cm}| \hspace{0.05cm}m_0) = {\eta_1}/{\sigma_n^2} \cdot {\rm e }^{ - ({\eta_1^2 + E})/({2 \sigma_n^2}) }\cdot {\rm I }_0 \left [{\eta_1 \cdot \sqrt{E}}/{ \sigma_n^2}\right ] \hspace{0.05cm},$$
$$ p_{y_2 \hspace{0.01cm}| \hspace{0.03cm}m} ( \eta_2 \hspace{0.05cm}| \hspace{0.05cm}m_0) = {\eta_2}/{\sigma_n^2} \cdot {\rm e }^{ - \eta_2^2 /({2 \sigma_n^2}) } \hspace{0.05cm}.$$
  • The error probability is generally obtained as follows:
\[{\rm Pr}({\cal{E}}) = \int_{0}^{\infty} \int_{\eta_1}^{\infty} p_{y_1,\hspace{0.03cm} y_2 \hspace{0.03cm}| \hspace{0.03cm}m} ( \eta_1, \eta_2 \hspace{0.05cm}| \hspace{0.05cm}m_0) \,\,{\rm d} \eta_2\,\,{\rm d} \eta_1 = \int_{0}^{\infty} p_{y_1 \hspace{0.01cm}| \hspace{0.03cm}m} ( \eta_1 \hspace{0.05cm}| \hspace{0.05cm}m_0) \cdot \int_{\eta_1}^{\infty} p_{y_2 \hspace{0.03cm}| \hspace{0.03cm}m} ( \eta_2 \hspace{0.05cm}| \hspace{0.05cm}m_0)\,\,{\rm d} \eta_2\,\,{\rm d} \eta_1 \hspace{0.05cm}.\]

$\text{Conclusion:}$  For the  error probability with non-coherent demodulation of the binary FSK,  one obtains after some mathematical transformations the surprisingly simple result

\[p_{\rm S} = {\rm Pr}({\cal{E} }) = {1}/{2} \cdot {\rm e}^{-E_{\rm S}/(2 N_0)} \hspace{0.05cm}.\]
  • For comparison,  the result for  coherent demodulation  is here given again:
\[p_{\rm S} = {\rm Pr}({\cal{E} }) = {\rm Q}(\sqrt{ E_{\rm S}/N_0}) \hspace{0.05cm}.\]


$\text{Derivation:}$  This result is now to be derived in some calculation steps.  We start from the following equations:

$${\rm Pr}({\cal{E} }) = \int_{0}^{\infty} p_{y_1 \hspace{0.01cm} \vert\hspace{0.03cm}m} ( \eta_1 \hspace{0.05cm} \vert \hspace{0.05cm}m_0) \cdot \int_{\eta_1}^{\infty} p_{y_2 \hspace{0.03cm}\vert \hspace{0.03cm}m} ( \eta_2 \hspace{0.05cm}\vert \hspace{0.05cm}m_0)\,\,{\rm d} \eta_2\,\,{\rm d} \eta_1 \hspace{0.05cm},\hspace{0.5cm}\text{with}$$
$$p_{y_1 \hspace{0.01cm}\vert \hspace{0.03cm}m} ( \eta_1 \hspace{0.05cm}\vert \hspace{0.05cm}m_0) = {\eta_1}/{\sigma_n^2} \cdot {\rm e }^{ - ({\eta_1^2 + E})/({2 \sigma_n^2}) }\cdot {\rm I }_0 \left [{\eta_1 \cdot \sqrt{E} }/{ \sigma_n^2}\right ] \hspace{0.05cm},$$
$$p_{y_\hspace{0.01cm}\vert \hspace{0.03cm}m} ( \eta_2 \hspace{0.05cm}\vert \hspace{0.05cm}m_0) = {\eta_2}/{\sigma_n^2} \cdot {\rm e }^{ - \eta_2^2 /({2 \sigma_n^2}) } \hspace{0.05cm}.$$

(1)   The inner integral gives the probability that the Rayleigh distributed random variable  $\eta_2$  is larger than  $\eta_1$  – see solution to  "Exercise 4.17Z":

$$\int_{\eta_1}^{\infty} p_{y_2 \hspace{0.03cm}\vert\hspace{0.03cm}m} ( \eta_2 \hspace{0.05cm}\vert \hspace{0.05cm}m_0)\,\,{\rm d} \eta_2 = {\rm e }^{ - \eta_1^2 /({2 \sigma_n^2}) } \hspace{0.3cm} \Rightarrow \hspace{0.3cm} {\rm Pr}({\cal{E} }) = \int_{0}^{\infty}{\eta_1}/{\sigma_n^2} \cdot {\rm e }^{ - ({2\eta_1^2 + E})/({2 \sigma_n^2}) }\cdot {\rm I }_0 \left [ {\eta_1 \cdot \sqrt{E} }/{ \sigma_n^2}\right ]\,\,{\rm d} \eta_1 \hspace{0.05cm}.$$

(2)   With the  (arbitrary)  substitutions  $C_0^2 = E/4$  and  $\sigma_0^2 = \sigma_n^2/2$,  we obtain:

$${\rm Pr}({\cal{E} }) = \int_{0}^{\infty} \frac{\eta_1}{2 \cdot \sigma_0^2} \cdot {\rm exp } \left [ - \frac{2 \eta_1^2 + 4 C_0^2}{4 \sigma_0^2}\right ] \cdot {\rm I }_0 \left [ \frac{\eta_1 \cdot 2C_0}{ 2 \sigma_0^2}\right ]\,\,{\rm d} \eta_1 = \int_{0}^{\infty} \frac{\eta_1}{2 \cdot \sigma_0^2} \cdot {\rm exp } \left [ - \frac{\eta_1^2 + 2 C_0^2}{2 \sigma_0^2}\right ] \cdot {\rm I }_0 \left [ \frac{\eta_1 \cdot C_0}{ \sigma_0^2}\right ]\,\,{\rm d} \eta_1 \hspace{0.05cm}.$$

(3)   By shifting fractions in front of the integral,  we succeed that the integrand again describes a  "Rice distribution"

\[{\rm Pr}({\cal{E} }) = {1}/{2} \cdot {\rm exp } \left [ - \frac{ C_0^2}{2 \sigma_0^2}\right ] \cdot \int\limits_{0}^{\infty} \frac{\eta_1}{ \sigma_0^2} \cdot {\rm exp } \left [ - \frac{\eta_1^2 + C_0^2}{2 \sigma_0^2}\right ] \cdot {\rm I }_0 \left [ \frac{\eta_1 \cdot C_0}{ \sigma_0^2}\right ]\,\,{\rm d} \eta_1 \hspace{0.05cm}.\]

(4)   The integrand now describes the Rice PDF.  The integral over the entire definition area from  $0$  to  $+\infty$  yields the value one,  as for any PDF,  so that holds:

\[{\rm Pr}({\cal{E} }) = {1}/{2} \cdot {\rm exp } \left [ - \frac{ C_0^2}{2 \sigma_0^2}\right ]\hspace{0.05cm}.\]

(5)   Finally, with  $C_0^2 = E/4$  and  $\sigma_0^2 = \sigma_n^2/2$  and the generally valid relation  $\sigma_n^2 = N_0$,  we obtain:

$${\rm Pr}({\cal{E} }) = {1}/{2} \cdot {\rm exp } \left [ - \frac{ C_0^2}{2 \sigma_0^2}\right ] = {1}/{2} \cdot {\rm exp } \left [ - \frac{ E_{\rm S}/4}{N_{\rm 0}/2}\right ] \hspace{0.3cm}\Rightarrow \hspace{0.3cm} p_{\rm S} = {\rm Pr}({\cal{E} }) = {1}/{2} \cdot {\rm e}^{-E_{\rm S}/(2 N_0)}\hspace{0.05cm}.$$

        $E_{\rm S}$  here indicates the average signal energy per symbol, which is equal to the normalization energy  $E$  for FSK.

$\text{q.e.d.}$


FSK symbol error probability for coherent and non-coherent demodulation

$\text{Example 2:}$  The graph compares the error probability curves of both demodulation methods as a function of the AWGN quotient  $E_{\rm S}/N_0$. 



One can see:

  1. The non-coherent BFSK  (red curve)  requires a   $0.8 \ \rm dB$   larger $E_{\rm S}/N_0$  compared to the coherent BFSK  (blue curve)  at  $p_{\rm S}= 10^{-5}$. 

  2. At   $p_{\rm S}= 10^{-3}$   the distance is even  $1.3 \ \rm dB$.

  3. In contrast,  the distance between the coherent binary FSK from the coherent BPSK is always  $1.3 \ \rm dB$  regardless of the error probability.



Non-coherent demodulation of multi-level FSK


Orthogonal $M$–level FSK for  $M= 3$

We now consider the message set  $\{m_1, m_2,\hspace{0.05cm}\text{ ...} \hspace{0.05cm}, m_{M}\}$  and denote  $M$  as the  "level number".

  1. As in the case of binary FSK,  a prerequisite for the application of the modulation process  "Frequency Shift Keying"  and at the same time of a non-coherent demodulator is an integer modulation index  $h$.

  2. In this case,  the  $M$–level FSK is orthogonal and a signal space constellation results as shown in the adjacent diagram for the special case  $M = 3$. 

The non-coherent demodulator is sketched below.

  • Compared to the  "receiver structure for binary FSK",  this receiver differs only by  $M$  branches instead of only two,  which provide the comparison values  $y_1$,  $y_2$, ... , $y_M$.
  • To calculate the error probability,  we assume that  $m_1$  was sent. This means that the decision is correct if the largest detection output value is  $y_1$: 
\[{\rm Pr}({\cal{C}}) \hspace{-0.1cm} = \hspace{-0.1cm} {\rm Pr} \big [ (y_2 < y_1) \cap (y_3 < y_1) \cap ... \cap (y_{M} < y_1) \hspace{0.05cm}| \hspace{0.05cm} m = m_1\big ] = {\rm Pr} \left [ \hspace{0.1cm} \bigcap\limits_{k = 2}^M (y_k < y_1) \hspace{0.05cm}| \hspace{0.05cm}m = m_1\right ] \hspace{0.01cm}.\]
Non-coherent receiver structure for  $M$–level FSK


$\text{Conclusions:}$ 

  • The   error probability of M–level FSK in non-coherent demodulation   is equal to  ${\rm Pr}({\cal{E} }) = 1 - {\rm Pr}({\cal{C} })$.
  • In the following,  it is shown that  ${\rm Pr}({\cal{C} })$  can be represented as follows:
\[{\rm Pr}({\cal{C} }) = \sum_{i = 0}^{M-1} (-1)^i \cdot {M-1 \choose i } \cdot \frac{1}{i+1} \cdot {\rm exp } \left [ - \frac{i \cdot E_{\rm S} } {(i+1) \cdot N_0}\right ] \hspace{0.05cm}\hspace{0.5cm} \text{with}\hspace{0.5cm} E_{\rm S} = E_{\rm B} \cdot {\rm log_2}(M)\hspace{0.05cm}.\]
  • In the special case  $M = 2$,  of course, the result obtained in the  "last section"  is again obtained:
\[p_{\rm S} = {\rm Pr}({\cal{E} }) = {1}/{2} \cdot {\rm e}^{-E_{\rm S}/(2 N_0)} \hspace{0.05cm}.\]


$\text{Derivation:}$  This result shall now be derived in some calculation steps.

  • We assume that  $m_1$  has been sent.  This means that the decision is correct if the largest detection output value is  $y_1$. 


(1)   Using the conditional probability density  $p_{y_1 \hspace{0.01cm} \vert\hspace{0.03cm}m} ( \eta_1 \hspace{0.05cm} \vert \hspace{0.05cm}m_1)$,  we obtain:

\[{\rm Pr}({\cal{C} }) = \int_{0}^{\infty} {\rm Pr} \left [ \hspace{0.1cm} \bigcap\limits_{k = 2}^M (y_k < y_1) \hspace{0.05cm}\vert\hspace{0.05cm}y_1 = \eta_1, m = m_1\right ] \cdot p_{y_1 \hspace{0.01cm} \vert\hspace{0.03cm}m} ( \eta_1 \hspace{0.05cm} \vert \hspace{0.05cm}m_1) \,\,{\rm d} \eta_1 \hspace{0.05cm}.\]

(2)   The decision values  $y_2$,  $y_3$,  ... , $y_M$  are statistically independent for given  $y_1$.  Therefore:

\[{\rm Pr}({\cal{C} }) = \int_{0}^{\infty} \left \{ {\rm Pr} \big [ (y_2 < y_1) \hspace{0.05cm}\vert \hspace{0.05cm}y_1 = \eta_1, m = m_1\big ] \right \}^{M-1} \cdot p_{y_1 \hspace{0.01cm} \vert\hspace{0.03cm}m} ( \eta_1 \hspace{0.05cm} \vert \hspace{0.05cm}m_1) \,\,{\rm d} \eta_1 \hspace{0.05cm}.\]

(3)   The selected value  $y_2$  conditioned on  $m_1$  has a Rayleigh distribution with parameter  $\sigma_n^2$:

\[{\rm Pr} \big [ (y_2 < y_1) \hspace{0.05cm} \vert \hspace{0.05cm}y_1 = \eta_1, m = m_1\big ] \hspace{-0.1cm} = \hspace{-0.1cm} \int_{0}^{\eta_1} p_{y_2 \hspace{0.01cm} \vert\hspace{0.03cm}m} ( \eta_2 \hspace{0.05cm} \vert \hspace{0.05cm}m_1)\,\,{\rm d} \eta_2= 1 - {\rm exp } \big [ - {\eta_1^2 }/({2 \sigma_n^2})\big ] = 1 - a \hspace{0.2cm}{\rm(abbreviation)} \hspace{0.05cm}.\]

(4)   Now we are looking for the expression   $(1 -a)^{M-1}$,  for which with the abbreviation from  (3)  holds:

\[ (1-a)^{M-1} \hspace{-0.1cm} = \hspace{-0.1cm} \sum_{i = 0}^{M-1} (-1)^i \cdot {M-1 \choose i } \cdot (-1)^i \cdot a^i = \sum_{i = 0}^{M-1} (-1)^i \cdot {M-1 \choose i } \cdot {\rm exp } \left [ - \frac{i \cdot \eta_1^2 }{2 \sigma_n^2}\right ] \hspace{0.05cm}.\]

(5)   $y_1$  has a  "Rice distribution" conditioned on  $m=m_1$.  Thus,  the probability of a correct decision can be expressed in the following form:

\[{\rm Pr}({\cal{C} }) \hspace{-0.1cm} = \hspace{-0.1cm} \sum_{i = 0}^{M-1} (-1)^i \cdot {M-1 \choose i } \cdot \int_{0}^{\infty} {\rm exp } \left [ - \frac{i \cdot \eta_1^2 }{2 \sigma_n^2}\right ] \cdot \frac{\eta_1}{ \sigma_n^2}\cdot {\rm exp } \left [ - \frac{\eta_1^2 + E_{\rm S} }{2 \sigma_n^2}\right ] \cdot {\rm I }_0 \left [ \frac{\eta_1 \cdot \sqrt{E_{\rm S} } }{ \sigma_n^2}\right ] \,\,{\rm d} \eta_1 \hspace{0.05cm}.\]

(6)   Substitutions succeed in shaping the integrand according to the Rice distribution.  Since each probability density integrates to one,  we obtain:

\[{\rm Pr}({\cal{C} }) = \sum_{i = 0}^{M-1} (-1)^i \cdot {M-1 \choose i } \cdot \frac{1}{i+1} \cdot {\rm exp } \left [ - \frac{i \cdot E_{\rm S} } {(i+1) \cdot N_0}\right ] \hspace{0.05cm}\hspace{0.5cm} \text{with}\hspace{0.5cm} E_{\rm S} = E_{\rm B} \cdot {\rm log_2}(M)\hspace{0.05cm}.\]

(7)   The special case  $M = 2$  leads to exactly the same result as calculated for the binary FSK:

\[{\rm Pr}({\cal{C} }) = (-1)^0 \cdot {2-1 \choose 0 } \cdot \frac{1}{0+1} \cdot {\rm exp } \left [ - \frac{0 \cdot E_{\rm S} } {(i+1) \cdot N_0}\right ] + (-1)^1 \cdot {2-1 \choose 1 } \cdot \frac{1}{1+1} \cdot {\rm exp } \left [ - \frac{1 \cdot E_{\rm S} } {(i+1) \cdot N_0}\right ] \]
\[ \Rightarrow \hspace{0.5cm} {\rm Pr}({\cal{C} }) = 1 - {1}/{2} \cdot {\rm e }^{-E_{\rm S}/(2N_0)} \hspace{0.5cm} \Rightarrow \hspace{0.5cm} {\rm Pr}({\cal{E} }) = 1 - {\rm Pr}({\cal{C} }) = {1}/{2} \cdot {\rm e }^{-E_{\rm S}/(2N_0)} \hspace{0.05cm}.\]
$\text{q.e.d.}$

Exercises for the chapter


Exercise 4.17: Non-Coherent On-Off Keying

Exercise 4.17Z: Rayleigh and Rice Distribution

Exercise 4.18: Non-Coherent FSK Demodulation

Exercise 4.18Z: BER of Coherent and Non-Coherent FSK

Exercise 4.19: Orthogonal Multilevel FSK