Model-based fit procedure for power-law-like spectra

$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.