Considering the following random vectors:
\begin{align} \textbf{h} &= [h_{1}, h_{2}, \ldots, h_{M}]^{T} \sim \mathcal{CN}\left(\textbf{0}_{M},d\textbf{I}_{M \times M}\right), \\[8pt] \textbf{w} &= [w_{1}, w_{2}, \ldots, w_{M}]^{T} \sim \mathcal{CN}\left(\textbf{0}_{M},\frac{1}{p}\textbf{I}_{M \times M}\right), \\[8pt] \textbf{y} &= [y_{1}, y_{2}, \ldots, y_{M}]^{T} \sim \mathcal{CN}\left(\textbf{0}_{M},\left(d + \frac{1}{p}\right)\textbf{I}_{M \times M}\right), \end{align}
where $\textbf{y} = \textbf{h} + \textbf{w}$ and therefore, $\textbf{y}$ and $\textbf{h}$ are not independent.
I'm trying to find the following expectation:
$$\mathbb{E} \left[ \frac{\textbf{h}^{H} \textbf{y}\textbf{y}^{H} \textbf{h}}{ \| \textbf{y} \|^{4} } \right],$$
where $\| \textbf{y} \|^{4} = (\textbf{y}^{H} \textbf{y}) (\textbf{y}^{H} \textbf{y}$).
In order to find the desired expectation, I'm applying the following approximation:
$$\mathbb{E} \left[ \frac{\textbf{x}}{\textbf{z}} \right] \approx \frac{\mathbb{E}[\textbf{x}]}{\mathbb{E}[\textbf{z}]} - \frac{\text{cov}(\textbf{x},\textbf{z})}{\mathbb{E}[\textbf{z}]^{2}} + \frac{\mathbb{E}[\textbf{x}]}{\mathbb{E}[\textbf{z}]^{3}}\text{var}(\mathbb{E}[\textbf{z}]).$$
However, applying this approximation to the desired expectation is time consuming and prone to errors as it involves expansions with lots of terms .
I have been wondering if there is a more direct/smarter way of finding the desired expectation.
$\textbf{UPDATE 21-04-2018}$: I've created a simulation in order to identify the pdf shape of the ratio inside of the expectation operator and as can be seen below it seems much like the pdf of a Gaussian random variable. Additionally, I've also noticed that the ratio results in real valued terms, the imaginary part is always equal to zero.
Is there another kind of approximation that can be used to find the expectation (one analytical/closed form result and not only the simulated value of the expection) given that the pdf looks like a Gaussian and probably can be approximated as such?