Start with jointly Gaussian PDF (I, Q, I^, Q^)
Covariance matrix of (I, Q, I^, Q^)
where
bn is n-th moment of Doppler spectral power density
For uniform angles of arrives
where (n-1)!! = 1 3 5 .....(n-1)