compute welch method, very neat
time-frequency analysis with wavelet works too, but it needs to be tested