From a somewhat heuristic derivation in a physics problem, I found the identity
$$ \prod_{\ell=1}^{N-1}\left[1\pm x \cos\left(\frac{\ell \pi}{N}\right)\right] = \frac{\lambda_{+}^N - \lambda_{-}^N}{\lambda_+ - \lambda_-}; \qquad \lambda_{\pm} = \frac{1}{2} \left(1 \pm \sqrt{1-x^2}\right) $$
I assume this identity is already known (for example, the $x=1$ case reduces to a version of the well known formula discussed here, here, and many other times on math.stackexchange), but I am unfamiliar with the relevant literature.
Does anyone know/Can anyone construct an analytical proof for this result?