Proving $ \sum_{n \in \mathbb{Z} } \left[\frac{\sin (n \alpha + \theta) }{ n \alpha + \theta} \right]^2 = \frac{\pi}{\alpha} $

I will prove the stated equality for the case $|\alpha|\leq\pi$. (I believe it is false when $|\alpha|>\pi$.)

The Poisson summation formula tells us that if $f(x)$ is a function and $$ F(q) = \int_{-\infty}^{+\infty}dx\, f(x)\, e^{-i q x} $$ is its continuous Fourier transform (FT), then $$ \sum_{n = -\infty}^{+\infty} f(n)\;=\; \sum_{n = -\infty}^{+\infty} F(2\pi n) $$

For this problem, let $$ f(x) = {\mathrm{sinc}(\alpha x + \theta)}^2\, , $$ where $\mathrm{sinc}(x) = \sin(x)/x$. Then the FT of $f(x)$ is \begin{align*} F(q) &= \int_{-\infty}^{+\infty}dx\; {\mathrm{sinc}(\alpha x + \theta)}^2\; e^{-i q x}\\ &= e^{+i \, q\, \theta/\alpha} \int_{-\infty}^{+\infty}dx\; {\mathrm{sinc}(\alpha x)}^2\; e^{-i q x} \qquad \text{Change of variables: $x\rightarrow x - \theta/\alpha$}\\ &= e^{+i \, q\, \theta/\alpha} \frac{1}{|\alpha|}\int_{-\infty}^{+\infty}dz\; {\mathrm{sinc}(z)}^2\; e^{-i u z} \qquad z = |\alpha|x\, , \; u = q/|\alpha| \end{align*}

The FT of the square of a sinc function is a "tent" function: \begin{equation*} \int_{-\infty}^{+\infty}dz\; {\mathrm{sinc}(z)}^2\; e^{-i u z}\;=\; \begin{cases} \pi \left(1 - \frac{|u|}{2}\right) & |u| < 2\\ 0 & \text{else} \end{cases} \end{equation*} I'm not going to prove this here so as to not clutter up the derivation, but it can be found in various sources, and follows from the convolution theorem for Fourier transforms and the fact that the FT of a sinc function is a "top hat" function.

Putting this all together, we arrive at: \begin{equation} F(q) \;=\; e^{+i \, q\, \theta/\alpha} \; \frac{\pi}{|\alpha|} \; \begin{cases} 1 - \frac{|q|}{2|\alpha|} & |q| < 2|\alpha|\\ 0 & \text{else} \end{cases} \qquad\qquad (1) \end{equation}

The Poisson summation formula tells us that \begin{equation*} \sum_{n = -\infty}^{+\infty} {\left[\frac{\sin(n\alpha + \theta)}{n\alpha + \theta}\right]}^2\;=\; \sum_{n = -\infty}^{+\infty} f(n)\;=\; \sum_{n = -\infty}^{+\infty} F(2\pi n) \end{equation*}

In performing the latter sum over $F(2\pi n)$ using Eq. (1), we see that since $F(q)$ has compact support on $|q| < 2|\alpha|$, we need only sum $n$ over $|n| < |\alpha|/\pi$. When $|\alpha| \leq \pi$, only the $n = 0$ term contributes to the sum, and we are left with: \begin{equation*} \sum_{n = -\infty}^{+\infty} {\left[\frac{\sin(n\alpha + \theta)}{n\alpha + \theta}\right]}^2\;=\; F(0) \;=\; \frac{\pi}{|\alpha|} \qquad\qquad (|\alpha| \leq \pi) \end{equation*}

When $|\alpha| > \pi$, more terms with $n \neq 0$ enter the sum, and the expression becomes more complicated, the $\theta$ dependence does not go away, and I'm not sure there's a nice, closed-form expression. For instance, evaluating the sum in Mathematica for $\alpha = 4$ and $\theta = 0$ yields: \begin{equation*} \sum_{n = -\infty}^{+\infty} {\left[\frac{\sin(4 n)}{4 n}\right]}^2\;=\; \frac{6\pi - \pi^2}{8} \;\neq \; \frac{\pi}{4}\, . \end{equation*}

Edited to add:

The general sum is given by \begin{equation*} \sum_{n = -\infty}^{+\infty} {\left[\frac{\sin(n\alpha + \theta)}{n\alpha + \theta}\right]}^2\;=\; \frac{\pi}{|\alpha|}\left[ 1 \;+\; 2\sum_{n = 1}^{\left\lfloor |\alpha|/\pi \right\rfloor} \cos(2\pi \,n\, \theta\, /\, \alpha) \left(1 - \frac{\pi\, n}{|\alpha|}\right) \right]\, . \end{equation*} If this has a "nice" closed form, I'm unaware of it.