I assume $u_0,u_1,\beta$ are all real. For $k \in \{0,1\}$, define $\alpha_k := \sqrt{e^{j\pi u_k}}$ and define:
$$p_k(z) := \sum_{n=0}^{N-1} e^{j(n-\frac{N-1}{2})\pi u_k}z^n = \sum_{n=0}^{N-1} \alpha_k^{2n-(N-1)}z^n$$
So, $P(z) = p_0(z) - \beta p_1(z)$. Further, define:
$$r_k(z) := z^{-(N-1)}p_k(z^2) = \sum_{n=0}^{N-1} \alpha_k^{2n-(N-1)}z^{2n-(N-1)} = \frac{(\alpha_k z)^N - (\alpha_k z)^{-N}}{(\alpha_k z) - (\alpha_k z)^{-1}} = (\alpha_k z)^{1-N} \frac{(\alpha_k z)^{2N} - 1}{(\alpha_k z)^2 - 1}$$
So, $r_k(z)$ is a meromorphic function with $2N-2$ simple zeros on the unit circle, and $N-1$ poles at 0. Also, the zeros of $r_k(z)$ are square roots of the zeros of $p_k(z)$, which implies $p_k(z)$ has all its zeros on the unit circle.
The crucial feature of $r_k(z)$ is then the following property. Since $|\alpha_k| = 1$ and $z^* = z^{-1}$ on the unit circle, we have that $r_k(z) = (r_k(z))^*$ on the unit circle. That is, $r_k(z)$ is real-valued on the unit circle.
With this, consider the following. If $u_0$ and $u_1$ are not too far apart, then the zeros of $r_0(z)$ and $r_1(z)$ will alternate/interlace on the unit circle. (What "not too far apart" means here will depend on $N$. It's probably something like $|u_0-u_1| \leq \frac{2}{N}$, but don't quote me on it.) In this case, the roots of $p_0(z^2)$ and $p_1(z^2)$ then alternate as well. Finally, this implies that the roots of $p_0(z)$ and $p_1(z)$ also alternate. The rest of the argument is then a sort of Hermite-Kakeya-Obreschkoff (HKO) theorem for the unit circle.
Consider the function:
$$h(z) := \frac{p_0(z)}{p_1(z)}$$
Since $h(z^2) = \frac{r_0(z)}{r_1(z)}$ is real-valued on the unit circle, $h(z)$ therefore is as well. A standard zero counting argument (like in the proof of HKO) then implies $c_0 p_0(z) + c_1 p_1(z)$ has all its zeros on the unit circle for any real $c_0,c_1$. Therefore, $P(z) = p_0(z) - \beta p_1(z)$ has all its zeros on the unit circle.
Since HKO gives an equivalent condition for root location, this further says that if the zeros of $p_0$ and $p_1$ do not alternate (i.e., $|u_0-u_1|$ is too large), then there will be some real $\beta$ such that $P(z)$ does not have zeros only on the unit circle. However, if you are restricting $\beta$ (e.g., to $\beta > 0$), you may still have that the zeros of $P(z)$ are on the unit circle. You would have to tinker with these arguments to get something in that case.