A universal method of extraction of the complex dielectric function ϵ(ω)=ϵ1(ω)+iϵ2(ω) from experimentally accessible optical quantities is developed. The central idea is that ϵ2(ω) is parameterized independently at each node of a properly chosen anchor frequency mesh, while ϵ1(ω) is dynamically coupled to ϵ2(ω) by the Kramers-Kronig (KK) transformation. This approach can be regarded as a limiting case of the multi-oscillator fitting of spectra, when the number of oscillators is of the order of the number of experimental points. In the case of the normal-incidence reflectivity from a semi-infinite isotropic sample the new method gives essentially the same result as the conventional KK transformation of reflectivity. In contrast to the conventional approaches, the proposed technique is applicable, without readaptation, to virtually all types of linear-response optical measurements, or arbitrary combinations of measurements, such as reflectivity, transmission, ellipsometry {\it etc.}, done on different types of samples, including thin films and anisotropic crystals.