\documentclass[12pt,draftclsnofoot,onecolumn]{IEEEtran}

%\IEEEoverridecommandlockouts
\usepackage{graphicx}
% \usepackage{subcaption}

%\ifCLASSOPTIONcompsoc
%    \usepackage[caption=false, font=normalsize, labelfont=sf, textfont=sf]{subfig}
%\else
%\usepackage[caption=false, font=footnotesize]{subfig}
%\fi
\usepackage{cite}
\usepackage{url}
\usepackage{mdwmath}  
\usepackage{soul}
\usepackage[table]{xcolor}
\usepackage{amsfonts,amsmath,amssymb,amsxtra,amsbsy,floatflt}
\usepackage{indentfirst} 
\usepackage{tabularx}
\usepackage{multirow}
\usepackage{nth}
\usepackage{comment}
\usepackage{epstopdf}
\usepackage{hhline}
\usepackage{siunitx}
\usepackage{booktabs}
\usepackage[bottom]{footmisc}
\usepackage[acronym]{glossaries}
\usepackage{upgreek}
\usepackage{algorithm}
\usepackage{algorithmicx,algpseudocode}
\usepackage{bm}
\usepackage[utf8]{inputenc}
\usepackage{lipsum}
\usepackage{array} 
\usepackage{epsfig}
\usepackage{times}
\usepackage{framed}
\usepackage{enumitem}
\usepackage{float}
\usepackage{etoolbox}
\usepackage{algpseudocode}
\usepackage{lipsum}
\usepackage{soul}
\usepackage{mathtools}
\usepackage{bbm}
\usepackage{makecell}
\usepackage{pifont}
\usepackage[export]{adjustbox} 
\usepackage{latexsym}
\usepackage{multicol}
\usepackage{eurosym}
\usepackage{xspace}
\usepackage{pdfpages}
\usepackage{verbatim}
\usepackage{stfloats}
\usepackage{blindtext}
\input{Definitions.tex}
\newcommand{\rmRIS}{\scriptscriptstyle\mathrm{RIS}}
\newcommand{\rmBS}{\scriptscriptstyle\mathrm{BS}}
\newcommand{\rmT}{\scriptscriptstyle\mathrm{T}}
\newcommand{\rmR}{\scriptscriptstyle\mathrm{R}}
\newcommand{\rmS}{\scriptscriptstyle\mathrm{S}}
\newcommand{\rmd}{\scriptscriptstyle\mathrm{d}}
\newcommand{\rmL}{\scriptscriptstyle\mathrm{L}}
\newcommand{\rmO}{\scriptscriptstyle\mathrm{O}}
\newcommand{\rmG}{\scriptscriptstyle\mathrm{G}}

\input{acronyms_list}

\begin{document}
\section{Changing variable of normal distribution}
%https://andrewcharlesjones.github.io/journal/rayleigh.html
Consider two random variables $X$ and $Y$ that are independent and identically distributed from a zero-mean Gaussian with variance of $1$

\[X, Y \stackrel{iid}{\sim} \mathcal{N}(0, 1)\]


Now, consider the transformation to two new random variables $R$, $\Theta$ such that

\[R = \sqrt{X^2 + Y^2}\]

\[\Theta = \tan^{-1}\big(\frac{Y}{X}\big)\]

Also recall that we can transform from $R$ and $\Theta$ to $X$ and $Y$ by the following relationship:

\[X = r \cos \theta \quad \textrm{and} \quad Y = r \sin \theta.\]

\subsection{Computing the transformation}

Now, in order to compute the joint and marginal distributions of $R$ and $\Theta$, we can use the Jacobian transformation technique. Recall that the Jacobian transformation method defines the relationship between sets of random variables that have a one-to-one mapping between them. In our case, since we know the densities of $X$ and $Y$, we can write the joint density of $R$ and $\Theta$ as

\[f_{R, \Theta}(r, \theta) = f_{X, Y}(g_1^{-1}(r), g_2^{-1}(\theta))|J|\]

where $g_1^{-1}(r) = r \cos \theta$, $g_2^{-1}(\theta) = r \sin \theta$, and $\lvert J\rvert$ is the Jacobian of the transformation. That is:

\begin{align}
\lvert J\rvert &= \begin{vmatrix} \frac{\delta x}{\delta r} & \frac{\delta y}{\delta r} \\ \frac{\delta x}{\delta \theta} & \frac{\delta y}{\delta \theta} \end{vmatrix}
\end{align}

First, let’s compute the Jacobian of the transformation:

\begin{align}
\det J &= \begin{vmatrix} \frac{\delta x}{\delta r} & \frac{\delta y}{\delta r} \nonumber \\
\frac{\delta x}{\delta \theta} & \frac{\delta y}{\delta \theta} \end{vmatrix} \nonumber  \\
&= \begin{vmatrix} \cos \theta & \sin \theta \\ -r\sin\theta & r\cos\theta \end{vmatrix} \\
&= r\cos^2\theta + r\sin^2\theta \\ &= r
\end{align}

Also, because $X$ and $Y$ are independent, we can simply multiply their individual densities to get their joint density:

\[f_{X, Y}(x, y) = \frac{1}{2\pi} e^{-\frac{x^2 + y^2}{2}}.\]

Writing out the full density of $R$ and $\Theta$, we have

\begin{align}
f_{R, \Theta}(r, \theta) &=  f_{X, Y}(g_1^{-1}(r), g_2^{-1}(\theta))\lvert J\rvert \\
&= \frac{1}{2\pi} e^{-\frac{r^2\cos^2\theta * r^2\sin^2\theta}{2}}r \\
&= \frac{1}{2\pi} e^{\frac{r^2}{2}}r
\end{align}
where the $r$ at the end comes from the Jacobian computed above. Notice that $\theta$ doesn’t appear in this joint distribution, meaning that it’s constant in $\theta$, and $\theta$ is uniformly distributed.

To compute the marginal distributions for $R$ and $\Theta$, we can integrate out $\Theta$ and $R$, respectively.

\[f_R(r) = \int_0^{2\pi} f_{R, \Theta}(r, \theta) d\theta = re^{-\frac{r^2}{2}}\]

\[f_\Theta(\theta) = \int_0^{\infty} f_{R, \Theta}(r, \theta) dr = \frac{1}{2\pi}\]

Notice that $f_{R, \Theta}(r, \theta) = f_R(r)f_\Theta(\theta)$, indicating that $R$ and $\Theta$ are independent.

The distribution $f_R(r)$ above is known as a \textbf{Rayleigh distribution}, a common density whose practical usefulness is discussed later. It’s a special case of the chi distribution, and also has close ties to the chi-squared and gamma distributions.


\section{Changing variable of non-zero mean circular Gaussian distribution}
%https://andrewcharlesjones.github.io/journal/rayleigh.html
Consider two random variables $X$ and $Y$ that are independent and distributed from a non-zero-mean Gaussian with variance of $\frac{\sigma^2}{2}$

\[X \sim \mathcal{N}(\mu_x, \frac{\sigma^2}{2}) \quad \textrm{and} \quad Y \sim \mathcal{N}(\mu_y, \frac{\sigma^2}{2})\]

Their probability density function is

\[f_X(x) = \frac{1}{\sqrt{\pi}\sigma} e^{-\frac{(x-\mu_X)^2}{\sigma^2}} \quad \textrm{and} \quad f_Y(y) = \frac{1}{\sqrt{\pi}\sigma} e^{-\frac{(y-\mu_Y)^2}{\sigma^2}}\]

The pdf of the circular Gaussian distribution is
\[f_{XY}(x,y) = \frac{1}{\pi\sigma^2} \exp{\left(-\frac{(x-\mu_X)^2 + (y-\mu_Y)^2}{\sigma^2}\right)}\]


Now, consider the transformation to two new random variables $R$, $\Theta$ such that

\[R = \sqrt{X^2 + Y^2} \quad \textrm{and} \quad \Theta = \tan^{-1}\big(\frac{Y}{X}\big)\]

Also recall that we can transform from $R$ and $\Theta$ to $X$ and $Y$ by the following relationship:

\[x = r \cos \theta \quad \textrm{and} \quad y = r \sin \theta.\]

\subsection{Computing the transformation}

Now, in order to compute the joint and marginal distributions of $R$ and $\Theta$, we can use the Jacobian transformation technique. Recall that the Jacobian transformation method defines the relationship between sets of random variables that have a one-to-one mapping between them. In our case, since we know the densities of $X$ and $Y$, we can write the joint density of $R$ and $\Theta$ as

\[f_{R, \Theta}(r, \theta) = f_{X, Y}(g_1^{-1}(r), g_2^{-1}(\theta))|J|\]

where $g_1^{-1}(r) = r \cos \theta$, $g_2^{-1}(\theta) = r \sin \theta$, and $\lvert J\rvert$ is the Jacobian of the transformation. That is:

\begin{align}
\lvert J\rvert &= \begin{vmatrix} \frac{\delta x}{\delta r} & \frac{\delta y}{\delta r} \\ \frac{\delta x}{\delta \theta} & \frac{\delta y}{\delta \theta} \end{vmatrix}
\end{align}

First, let’s compute the Jacobian of the transformation:

\begin{align}
\det J &= \begin{vmatrix} \frac{\delta x}{\delta r} & \frac{\delta y}{\delta r} \\
\frac{\delta x}{\delta \theta} & \frac{\delta y}{\delta \theta} \end{vmatrix} 
= \begin{vmatrix} \cos \theta & \sin \theta \\ -r\sin\theta & r\cos\theta \end{vmatrix} \\
&= r\cos^2\theta + r\sin^2\theta = r
\end{align}

Let us also define the equivalences $\mu_X = \mu_R \cos \mu_\Theta$ and $\mu_Y = \mu_R \sin \mu_\Theta$

Writing out the full density of $R$ and $\Theta$, we have



\begin{align}
f_{R, \Theta}(r, \theta) &=  f_{X, Y}(g_1^{-1}(r), g_2^{-1}(\theta))\lvert J\rvert \nonumber \\
&= \frac{r}{\pi\sigma^2} \exp{\left(-\frac{( r \cos \theta-\mu_R \cos \mu_\Theta)^2 + (r \sin \theta-\mu_R \sin \mu_\Theta)^2}{\sigma^2}\right)} \label{eq:full_density_1}
\end{align}


Let us simplify the exponent
\begin{align}
& ( r \cos \theta-\mu_R \cos \mu_\Theta)^2 + (r \sin \theta-\mu_R \sin \mu_\Theta)^2 \nonumber \\
=& r^2 \cos^2 \theta-2\mu_R \cos \mu_\Theta r \cos \theta + \mu_R^2 \cos^2 \mu_\Theta +
 r^2 \sin^2 \theta-2\mu_R \sin \mu_\Theta r \sin \theta + \mu_R^2 \sin^2 \mu_\Theta \nonumber \\
 =& r^2 + \mu_R^2 -2r \mu_R (\cos \mu_\Theta \cos \theta + \sin \mu_\Theta \sin \theta )\nonumber \\
 =& r^2 + \mu_R^2 -2r \mu_R \cos (\mu_\Theta - \theta)\nonumber
\end{align}
and rewrite \eqref{eq:full_density_1} as
\begin{align}
f_{R, \Theta}(r, \theta) &=
\frac{r}{\pi\sigma^2} \exp{\left(-\frac{r^2 + \mu_R^2 -2 r \mu_R \cos (\mu_\Theta - \theta)}{\sigma^2}\right)} \label{eq:full_density_2}
\end{align}


Let us compute the marginal distribution for $\Theta$

\begin{align}
f_\Theta(\theta) &= \int_0^{\infty} f_{R, \Theta}(r, \theta) dr\label{eq:full_density_3}
\end{align}



From the equivalence of $x^2-x2A = (x - A)^2 - A^2$ we rewrite the exponent as
\begin{align}
r^2 -2r \mu_R \cos (\mu_\Theta - \theta)\nonumber &=  (r - \mu_R \cos (\mu_\Theta - \theta))^2 - \mu_R^2 \cos^2 (\mu_\Theta - \theta) \nonumber
\end{align}

{\color{blue}
in the exponent it actually is
\begin{align}
r^2 -2r \mu_R \cos (\mu_\Theta - \theta) + \mu_R^2 \nonumber &=  (r - \mu_R \cos (\mu_\Theta - \theta))^2 - \mu_R^2 \cos^2 (\mu_\Theta - \theta) + \mu_R^2 \nonumber \\
r^2 -2r \mu_R \cos (\mu_\Theta - \theta) + \mu_R^2 \nonumber &=  (r - \mu_R \cos (\mu_\Theta - \theta))^2 + \mu_R^2 \sin^2 (\mu_\Theta - \theta) \nonumber \\
\end{align}
}

And substitute it in \eqref{eq:full_density_2}. We proceed with the execution and simplification of \eqref{eq:full_density_3}


\begin{align}
f_\Theta(\theta) &= \int_0^{\infty} f_{R, \Theta}(r, \theta) dr \nonumber \\
& = \frac{1}{\pi\sigma^2}  \exp{\left(-\frac{\mu_R^2}{\sigma^2}\right)}
\int_0^{\infty} r \exp{\left(-\frac{r^2 -2r \mu_R \cos (\mu_\Theta - \theta)}{\sigma^2}\right)} dr \nonumber \\
& = \frac{1}{\pi\sigma^2}  \exp{\left(-\frac{\mu_R^2}{\sigma^2}\right)}
\int_0^{\infty} r \exp{\left(-\frac{(r - \mu_R \cos (\mu_\Theta - \theta))^2 - \mu_R^2 \cos^2 (\mu_\Theta - \theta)}{\sigma^2}\right)} dr  \nonumber \\
& = \frac{1}{\pi\sigma^2}  \exp{\left(-\frac{\mu_R^2}{\sigma^2}+\frac{  \mu_R^2 \cos^2 (\mu_\Theta - \theta)}{\sigma^2}\right)}
\int_0^{\infty} r \exp{\left(-\frac{(r - \mu_R \cos (\mu_\Theta - \theta))^2}{\sigma^2}\right)} dr \nonumber \\
& = \frac{1}{\pi\sigma^2}  \exp{\left(-\frac{\mu_R^2}{\sigma^2}\left(1-\cos^2 (\mu_\Theta - \theta)\right)\right)}
\int_0^{\infty} r \exp{\left(-\left(\frac{r - \mu_R \cos (\mu_\Theta - \theta)}{\sigma}\right)^2\right)} dr \nonumber
\end{align}

hint From Wolfram 
\[\int e^{-(x-A)^2/B^2} x dx = \frac{1}{2} B \left(B \left(-e^{-\left(\frac{A -  x}{B}\right)^2}\right)- \sqrt{\pi} A \,\text{erf}\left(\frac{A - x}{B}\right) \right)+ constant \]


let us introduce the variable
\begin{align}
u &= \frac{r - \mu_R \cos (\mu_\Theta - \theta)}{\sigma} \nonumber \\
r &= u\sigma + \mu_R \cos (\mu_\Theta - \theta)  \nonumber \\
du = \frac{dr}{\sigma} & \quad \textrm{and} \quad dr = \sigma du \nonumber
\end{align}


The integral becomes

Let introduce $\alpha = \cos (\mu_\Theta - \theta)$


\small
\begin{align}
f_\Theta(\theta) &= \int_0^{\infty} f_{R, \Theta}(r, \theta) dr \nonumber \\
& = \frac{1}{\pi\sigma^2}  \exp{\left(-\frac{\mu_R^2}{\sigma^2}\left(1-\cos^2 (\mu_\Theta - \theta)\right)\right)}
\int_0^{\infty} r \exp{\left(-\left(\frac{r - \mu_R \cos (\mu_\Theta - \theta)}{\sigma}\right)^2\right)} dr \nonumber \\
& = \frac{1}{\pi\sigma^2}  \exp{\left(-\frac{\mu_R^2}{\sigma^2}\left(1-\alpha^2\right)\right)}
\int_0^{\infty} r \exp{\left(-\left(\frac{r - \mu_R \alpha}{\sigma}\right)^2\right)} dr \nonumber \\
& = \frac{1}{\pi\sigma^2}  \exp{\left(-\frac{\mu_R^2}{\sigma^2}\left(1-\alpha^2\right)\right)}
\int_{\frac{- \mu_R \alpha}{\sigma}}^{\infty} \sigma(u\sigma + \mu_R \alpha) \exp{\left(-u^2\right)} du \nonumber \\
& = \frac{1}{\pi\sigma^2}  \exp{\left(-\frac{\mu_R^2}{\sigma^2}\left(1-\alpha^2\right)\right)}
\left[\int_{\frac{- \mu_R \alpha}{\sigma}}^{\infty} u\sigma^2 \exp{\left(-u^2\right)} du + \int_{\frac{- \mu_R \alpha}{\sigma}}^{\infty}\sigma\mu_R \alpha \exp{\left(-u^2\right)} du 
\right]\nonumber\\
& = \frac{1}{\pi}  \exp{\left(-\frac{\mu_R^2}{\sigma^2}\left(1-\alpha^2\right)\right)}
\int_{\frac{- \mu_R \alpha}{\sigma}}^{\infty} u \exp{\left(-u^2\right)} du  +  \frac{\mu_R \alpha }{\pi\sigma}  \exp{\left(-\frac{\mu_R^2}{\sigma^2}\left(1-\alpha^2\right)\right)}\int_{\frac{- \mu_R \alpha}{\sigma}}^{\infty}\exp{\left(-u^2\right)} du \nonumber\\
& =\frac{1}{\pi} \exp{\left(-\frac{\mu_R^2}{\sigma^2}\left(1-\alpha^2\right)\right)} \left[
\int_{\frac{- \mu_R \alpha}{\sigma}}^{\infty} u \exp{\left(-u^2\right)} du  +  \frac{\mu_R \alpha }{\sigma}  \int_{\frac{- \mu_R \alpha}{\sigma}}^{\infty}\exp{\left(-u^2\right)} du \right] \nonumber\\
& =\frac{1}{\pi} \exp{\left(-\frac{\mu_R^2}{\sigma^2}\left(1-\alpha^2\right)\right)} \left[
\left[-\frac{1}{2}\exp{\left(-u^2\right)}\right]_{\frac{- \mu_R \alpha}{\sigma}}^{\infty}  +  \frac{\mu_R \alpha \sqrt{\pi}}{2\sigma}  \frac{2}{\sqrt{\pi}}\int_{\frac{- \mu_R \alpha}{\sigma}}^{\infty}\exp{\left(-u^2\right)} du \right] \nonumber\\
& =\frac{1}{\pi} \exp{\left(-\frac{\mu_R^2}{\sigma^2}\left(1-\alpha^2\right)\right)} \left[
\frac{1}{2}\exp{\left(-\frac{\mu_R^2 \alpha^2}{\sigma^2}\right)}  +  \frac{\mu_R \alpha \sqrt{\pi}}{2\sigma} \left[\text{erf}(\infty) - \text{erf}\left(\frac{- \mu_R \alpha}{\sigma}\right)\right] \right] \nonumber\\
& =\frac{1}{2\pi} \exp{\left(-\frac{\mu_R^2}{\sigma^2}\left(1-\alpha^2\right)\right)} \left[
\exp{\left(-\frac{\mu_R^2 \alpha^2}{\sigma^2}\right)}  +  \frac{\mu_R \alpha \sqrt{\pi}}{\sigma} \left[1 + \text{erf}\left(\frac{\mu_R \alpha}{\sigma}\right)\right] \right] \nonumber\\
\end{align}
\normalsize



Let us define the SNR $\gamma = \frac{\mu_R^2}{\sigma^2}$ and re substitute $\alpha$

\small
\begin{align}
f_\Theta(\theta) & =\nonumber\\
&= \frac{1}{2\pi} e^{-\gamma\left(1-\cos^2 (\mu_\Theta - \theta)\right)} \left[
e^{-\gamma\cos^2 (\mu_\Theta - \theta)}  +  \sqrt{\gamma}\sqrt{\pi}\cos (\mu_\Theta - \theta)  \left[1 + \text{erf}\left(\sqrt{\gamma}\cos (\mu_\Theta - \theta)\right)\right] \right] \nonumber \\
f_\Theta(\theta) & =\nonumber\\
&= 
\frac{1}{2\pi} e^{-\gamma}  +  \frac{1}{2\pi} e^{-\gamma\sin^2 (\mu_\Theta - \theta)}\sqrt{\gamma}\sqrt{\pi}\cos (\mu_\Theta - \theta)  \left[1 + \text{erf}\left(\sqrt{\gamma}\cos (\mu_\Theta - \theta)\right)\right] \label{eq:marginal_theta}\\
\end{align}
\normalsize

For large values of $\gamma$, i.e., $\mu_R \gg \sigma$, we can approximate \eqref{eq:marginal_theta} as follows
\begin{align}
f_\Theta(\theta) & =\nonumber\\
& = \frac{1}{2\pi} \underbrace{e^{-\gamma}}
_{\approx 0}  +  \frac{1}{2\pi} \underbrace{e^{-\gamma \sin^2 (\mu_\Theta - \theta)}}_{\approx e^{-\gamma(\mu_\Theta - \theta)^2}}\sqrt{\gamma}\sqrt{\pi}\underbrace{\cos (\mu_\Theta - \theta)}_{\approx 1}  \left[1 + \underbrace{\text{erf}\left(\sqrt{\gamma}\cos (\mu_\Theta - \theta)\right)}_{\approx 1}\right]  \nonumber \\
& \approx \frac{\sqrt{\gamma}}{\sqrt{\pi}} e^{-\gamma(\mu_\Theta - \theta)^2}  = \frac{1}{\sqrt{2\pi}\sigma_{\Theta}} e^{-\frac{(\mu_\Theta - \theta)^2}{2\sigma_{\Theta}}} \label{eq:gaussian_approx}
\end{align}
with $\sigma_{\Theta}=\frac{1}{\sqrt{2\gamma}}$

We will have small variations of $\theta$ with respect $\mu_{\theta}$



Matlab verified!!!



\section{For appendix}
%https://andrewcharlesjones.github.io/journal/rayleigh.html
Let us consider a noisy channel sample $h = u + jv \in \mathbb{C}$, with $u$ and $v$ being two statistically independent \glspl{rv} distributed as
$u \sim \mathcal{N}\left(\mu_u, \sigma^2/2\right)$ and $v \sim \mathcal{N}\left(\mu_v, \sigma^2/2\right)$.
$h$ is then distributed as a circular Gaussian \gls{rv}, i.e, $h \sim \mathcal{CN}(\mu_{h},\sigma^2)$, with \gls{pdf} described as
\begin{align}
f_{h}(h)\!\! =\!\! \frac{1}{\pi\sigma^2} e^{-\frac{|h-\mu_{h}|^2}{\sigma^2}} = f_{u,v}(u,v) = \frac{1}{\pi\sigma^2} e^{-\frac{(u-\mu_u)^2 + (v-\mu_v)^2}{\sigma^2}},
\end{align}
with $|\mu_h|^2 = \mu_u^2 + \mu_v^2$ being the signal power, and $\sigma^2$ being the noise power.

Let us introduce two new random variables $r$, and $\theta$ such that
$r = |h| = \sqrt{u^2 + v^2}$ and $\theta = \angle h = \tan^{-1}\left(v/u\right)$. Let us recall also the relationship $u = r \cos \theta$ and $v = r \sin \theta$.

The joint distribution of $r$ and $\theta$ is
\begin{align}
f_{r, \theta}(r, \theta) = f_{u, v}(g_1^{-1}(r), g_2^{-1}(\theta))|J| \label{eq:joint_distr_1}
\end{align}
where $g_1^{-1}(r) = r \cos \theta$, $g_2^{-1}(\theta) = r \sin \theta$, and $\lvert J\rvert$ is the Jacobian of the transformation, which is
\begin{align}
\lvert J\rvert &= \begin{vmatrix} \frac{d\, u}{d\, r} & \frac{d\, v}{d\, r} \\
\frac{d\, u}{d\, \theta} & \frac{d\, v}{d\, \theta} \end{vmatrix} 
= \begin{vmatrix} \cos \theta & \sin \theta \\ -r\sin\theta & r\cos\theta \end{vmatrix} = r\cos^2\theta + r\sin^2\theta = r.
\end{align}

We can rewrite \eqref{eq:joint_distr_1} as
\begin{align}
f_{r, \theta}(r, \theta)
&= \frac{r}{\pi\sigma^2}e^{-\frac{r^2 + \mu_r^2 -2 r \mu_r \cos (\mu_\theta - \theta)}{\sigma^2}}, \label{eq:joint_distr_2}
\end{align}
where $\mu_r = \sqrt{\mu_u^2 + \mu_v^2}$, and $\mu_{\theta} = \tan^{-1}(\mu_v/\mu_u)$.

The marginal distribution for $\theta$ is
\begin{align}
f_\theta(\theta) &= \int_0^{\infty} f_{r, \theta}(r, \theta) dr = \frac{1}{\pi\sigma^2}  e^{-\frac{\mu_r^2}{\sigma^2}\sin^2 (\mu_\theta - \theta)}
\int_0^{\infty} r \, e^{-\left(\frac{r - \mu_r \cos (\mu_\theta - \theta)}{\sigma}\right)^2} dr. \label{eq:marginal_distr_1}
\end{align}
The integral is solved by substitution with $t = \frac{r - \mu_R \cos (\mu_\Theta - \theta)}{\sigma}$. From which we have $r = t\sigma + \mu_R \cos (\mu_\theta - \theta)$, $dt = dr/\sigma$, and $dr = \sigma dt$. We can rewrite \eqref{eq:marginal_distr_1} as
\begin{align}
f_\theta(\theta) &= \frac{1}{\pi\sigma^2}  e^{-\frac{\mu_r^2}{\sigma^2}\sin^2 (\mu_\theta - \theta)}
\left(\int_{t|_{r=0}}^{\infty} t\sigma^2 e^{-t^2} dt + \int_{t|_{r=0}}^{\infty}\sigma\mu_r \cos (\mu_\theta - \theta) e^{-t^2} dt \right), \label{eq:marginal_distr_3}
\end{align}
from which we obtain the following marginal distribution
\begin{align}
f_\theta(\theta) & =
\frac{1}{2\pi} e^{-\gamma}  +  \frac{1}{2\pi} e^{-\gamma\sin^2 (\mu_\theta - \theta)}\sqrt{\gamma}\sqrt{\pi}\cos (\mu_\theta - \theta)  \left[1 + \text{erf}\left(\sqrt{\gamma}\cos (\mu_\theta - \theta)\right)\right], \label{eq:marginal_distr_4}
\end{align}
where $\gamma = \frac{\mu_R^2}{\sigma^2}$ is the \gls{snr}, and $\text{erf}(\cdot)$ is the error function.

For large values of $\gamma$, i.e., $\mu_R \gg \sigma$, we can approximate \eqref{eq:marginal_theta} as
\begin{align}
f_\theta(\theta) & = \frac{1}{2\pi} \underbrace{e^{-\gamma}}
_{\approx 0}  +  \frac{1}{2\pi} \underbrace{e^{-\gamma \sin^2 (\mu_\theta - \theta)}}_{\approx e^{-\gamma(\mu_\theta - \theta)^2}}\sqrt{\gamma}\sqrt{\pi}\underbrace{\cos (\mu_\theta - \theta)}_{\approx 1}  \left[1 + \underbrace{\text{erf}\left(\sqrt{\gamma}\cos (\mu_\theta - \theta)\right)}_{\approx 1}\right]  \nonumber \\
& \approx \frac{\sqrt{\gamma}}{\sqrt{\pi}} e^{-\gamma(\mu_\theta - \theta)^2}  = \frac{1}{\sqrt{2\pi}\sigma_{\theta}} e^{-\frac{(\mu_\theta - \theta)^2}{2\sigma_{\theta}}} \label{eq:gaussian_approx}
\end{align}
with $\sigma_{\theta}=\frac{1}{\sqrt{2\gamma}}$, which proves the validity of \eqref{eq:theta_distribution}.
\end{document}