function poor_mans_power_spectrum(rawdata)

samplingrate=10000;
n=length(rawdata);

p = abs(fft(rawdata))/(n/2);
p = p(1:n/2).^2;
freq = [0:n/2-1]/(n/samplingrate);
semilogy(freq,p,'k','linewidth',3);