The sum I am looking for is the following sum as $M \to \infty$:
$$ L(\omega) = \sum_{m = 1}^{M} \frac{\sin\left( N \frac{\omega_m - \omega}{2} \right)}{\sin\left( \frac{\omega_m - \omega}{2} \right)} \cos\left(N \frac{\omega_m - \omega}{2} + \beta_m \right) $$
where
- $\omega_m$ is a random number from a Gaussian distribution having the parameters (mean $\mu$ and variance $\sigma^2$). $$ \omega_m \sim \mathcal{N}(\mu, \sigma^2) $$
- the $\beta_m$ are random numbers drawn from a uniform distribution from $-\pi$ to $+\pi$,
$$ \omega_m \sim \mathcal{U}[-\pi, +\pi] $$$$ \beta_m \sim \mathcal{U}[-\pi, +\pi] $$