From the new infrared (IR) reflectivity and time-domain terahertz spectra combined with available high-frequency dielectric data above the megahertz range in a broad temperature range of 10 to 900 K, a full picture of the soft- and central-mode behavior in the classical relaxor ferroelectric $\mathrm{Pb}(\mathrm{M}{\mathrm{g}}_{1/3}\mathrm{N}{\mathrm{b}}_{2/3}){\mathrm{O}}_{3}$ (PMN) is suggested. A detailed comparison is made with the recent hyper-Raman spectroscopy data [Hehlen et al., Phys. Rev. Lett. 117, 155501 (2016)] and also with other available experiments based on inelastic light and neutron scattering. It is revealed that each type of experiment provides slightly different data. The closest agreement is with the hyper-Raman data: both techniques yield the same number of soft-mode components and the same high-temperature softening towards the temperature ${T}^{*}\ensuremath{\approx}400\phantom{\rule{0.16em}{0ex}}\mathrm{K}$. In addition to evaluation of the IR-terahertz data using fitting with the standard factorized form of the dielectric function, we performed a successful fitting of the same data using the effective medium approach (EMA), originally based on the assumption that the mesoscopic structure of PMN consists of randomly oriented uniaxially anisotropic polar nanodomains (PNDs) with somewhat harder transverse optical polar modes in the direction along the local PND dipole [Hlinka et al., Phys. Rev. Lett. 96, 027601 (2006)]. Evaluation using Bruggeman EMA modeling has been successfully applied in the entire investigated temperature range. These results suggest that the response perpendicular to the local dipole moment, at high temperatures induced by random fields rather than PNDs, undergoes a classical softening from high temperatures with permittivity obeying the Curie-Weiss law, ${\ensuremath{\varepsilon}}_{\ensuremath{\perp}}=C/(T\ensuremath{-}{T}_{C})$, where $C=1.7\ifmmode\times\else\texttimes\fi{}{10}^{5}\phantom{\rule{0.16em}{0ex}}\mathrm{K}$ and ${T}_{C}=380\phantom{\rule{0.16em}{0ex}}\mathrm{K}$, whereas the response parallel to it shows no softening. Below the Burns temperature, $\ensuremath{\sim}620\phantom{\rule{0.16em}{0ex}}\mathrm{K}$, a gigahertz relaxation ascribed to flipping of the PNDs emerges from the soft-mode response, slows down, and broadens, remaining quite strong towards the cryogenic temperatures, where it can be assigned to fluctuations of the PND boundaries.