pilot_spec function

Calculate initial power spectral density estimates