# Distribution of the Rayleigh quotient

For a research project I need to find the expected value of the generalized Rayleigh quotient:
Here A and B are positive definite deterministic p x p covariance matrices, and w follows a multivariate distribution with circular altitude lines (say, multivariate standard normal). The dimension p is greater than 100.

This problem is easy to solve using simulation; however, I was wondering if somebody might know how this problem could be solved (or approximated) analytically. My first idea was that possibly by the Lindeberg or Lyapunov central limit theorem both the numerator and the denominator are approximately normal distributed, which gives us a ratio of two (correlated) normal random variables, but simulation shows that that is not the case.

Quadratic forms in elliptic distributions and their moments are treated in Mathai, Provost and Hayakawa, Bilinear forms and zonal polynomials (1995), but not to the same extend as in the normal case. As elliptical distributions are usually defined in terms of their characteristic function $e^{it\mu}\xi(t'\Sigma t)$, this function $\xi$ will appear in the solution if one chooses the mgf-approach. Yet, it has never been calculated, afaik.