The billiard motion inside an ellipsoid $Q\subset\mathbb{R}^{n+1}$ is completely integrable. Its phase space is a symplectic manifold of dimension $2n$, which is mostly foliated with Liouville tori of dimension n. The motion on each Liouville torus becomes just a parallel translation with some frequency $\omega$ that varies with the torus. Further, any billiard trajectory inside Q is tangent to n caustics $Q_{\lambda_1},\ldots,Q_{\lambda_n}$, so the caustic parameters $\lambda=(\lambda_1,\ldots,\lambda_n)$ are integrals of the billiard map. The frequency map $\lambda\mapsto\omega$ is a key tool for understanding the structure of periodic billiard trajectories. In principle, it is well-defined only for nonsingular values of the caustic parameters. We present two conjectures, fully supported by numerical experiments. We obtain, from one of the conjectures, some lower bounds on the periods. These bounds depend only on the type of the n caustics. We describe the geometric meaning, domain, and range of $\omega$. The map $\omega$ can be continuously extended to singular values of the caustic parameters, although it becomes “exponentially sharp” at some of them. Finally, we study triaxial ellipsoids of $\mathbb{R}^3$. We numerically compute the bifurcation curves in the parameter space on which the Liouville tori with a fixed frequency disappear. We determine which ellipsoids have more periodic trajectories. We check that the previous lower bounds on the periods are optimal, by displaying periodic trajectories with periods four, five, and six whose caustics have the right types. We also give some new insights for ellipses of $\mathbb{R}^2$.