  $1/f^\alpha$ noises are ubiquitous and affect many measurements. These noises
are both a nuisance and a peculiarity of several physical systems; in
dielectrics, glasses and networked liquids it is very common to study this
noise to gather useful information. Sometimes it happens that the noise has a
power-law shape only in a certain frequency range, and contains other important
features, that are however difficult to study because simple fits often fail.
Here I propose a model-based fit procedure that performs well on spectra
obtained in a molecular dynamics simulation.
