2
$\begingroup$

Given that following two random variables $\textbf{x} \sim \mathcal{CN}(\textbf{0}_{M},\sigma_{x}^{2}\textbf{I}_{M})$ and $\textbf{y} \sim \mathcal{CN}(\textbf{0}_{M},\sigma_{y}^{2}\textbf{I}_{M})$ are independent, what would be the expectation

$$\mathbb{E} \left[ \left| \frac{\textbf{x}^{H} \textbf{y} }{\| \textbf{x} \|^2} \right|^2 \right], $$ where $\mathcal{CN}(.,.)$ is the complex normal random variable.

$\endgroup$
4
  • 1
    $\begingroup$ I think you can use the trace trick for that, $|x^H y|^2 = x^H y y^H x = \mathrm{tr} \, x^H y y^H x = \mathrm{tr}\, x x^H \, y y^H$ and the trace commutes with the expectation... $\endgroup$
    – user114668
    Feb 16, 2019 at 17:36
  • $\begingroup$ @student, thanks! do you think we can apply that trick by Marsaglia? $\endgroup$ Feb 16, 2019 at 17:38
  • 1
    $\begingroup$ In general, Marsaglia’s trick is only useful to get the $x$-linear expectation $E[x \, f(y)]$ for normal and correlated $x$ and $y$. $\endgroup$
    – user114668
    Feb 16, 2019 at 18:09
  • $\begingroup$ @student, with the trace trick I find the first term equal to 1 and the second equal to $\frac{M\sigma_{y}^2}{(M-1)\sigma_{x}^2}$. The $M$ term should not be in the numerator. Do you think there is something wrong? $\endgroup$ Feb 16, 2019 at 20:34

2 Answers 2

5
$\begingroup$

Using the representation in this answer, $$Z= \frac{|\textbf{x}^{H} \textbf{y} |^2}{ |\textbf{x} |^4} =\frac{\sigma_y^2}{\sigma_x^2} \frac{\xi_{2M}\xi_{2}}{(\xi_{2}+\xi_{2M-2})^2},$$ and integrating over the independent chi-squared variables $\xi_2$, $\xi_{2M}$, and $\xi_{2M-2}$ I find the expectation value $$\Rightarrow\mathbb{E}(Z)=\frac{1}{M-1}(\sigma_y/\sigma_x)^2.$$ More generally, the $p$-th moment is finite for $M>p$, given by $$\mathbb{E}(Z^p)=\frac{(\sigma_y/\sigma_x)^2}{{M-1}\choose{p}}.$$

This answer looks simple enough, I wonder which distribution $P(Z)$ has reciprocals of binomial coefficients as its moments?

$\endgroup$
1
  • 1
    $\begingroup$ If you look for a random variable $Z_N$ s.t. $\mathbb{E}(Z_N^p) = {N\choose p}^{-1}$ for $p\leq N$ and undefined for $p \geq N+1$, you can have it with the ratio of an independent $ \boldsymbol{e}\sim $Exp(1) random variable and $ \gamma_{N + 1}\sim$Gamma($N + 1$) random variable. You can check that $ \mathbb{E}((\boldsymbol{e}/\gamma_{N + 1})^p ) = 1/{N\choose p} $. $\endgroup$
    – Synia
    Feb 20, 2019 at 20:34
4
$\begingroup$

To add to the very intriguing accepted answer, here is a quick way to check without doing the difficult bits:

Use the trace trick to factor the mixed expression, $$|x^H y|^2 = x^H y y^H x = \mathrm{tr} \, x^H y y^H x = \mathrm{tr}\, x x^H \, y y^H$$ and move the trace outside to get the expectation $$\mathrm{E}\left[\frac{|x^H y|^2}{\Vert x \Vert^4}\right] = \mathrm{tr} \, E\left[\frac{xx^H}{\Vert x \Vert^4}\right] \Sigma_y\;.$$ Moving the expectation back outside produces $\mathrm{tr} \, xx^H \, \Sigma_y = \sigma_y^2 \, x^H x = \sigma_y^2 \, \Vert x \Vert^2$, since $\Sigma_y = \sigma_y^2 \, I$. The remaining expectation is that of an inverse-chi-squared distribution with $2M$ degrees of freedom, $$E\left[\frac{1}{\Vert x \Vert^2}\right] = \frac{2}{\sigma_x^2} \, \frac{1}{2M-2}$$ so that the final result is indeed $$\frac{\sigma_y^2}{\sigma_x^2} \, \frac{1}{M-1}\;.$$

(This might also be helpful in other situations, since normality is not used until the very end.)

$\endgroup$

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.