function graphBuoyPSD()
    S = load('buoy.mat');
    z = S.buoy/100; % convert to meters since original data was in cm
    fs = 1.28;
    Nfft = 256;
    [P,f] = myPSD(z,fs,Nfft);
    figure(1)
        plot(f,P)
        xlim([0,0.4])
        texText('f [Hz]','x')
        texText('mean power','y')
        graphToPDF('../buoyPSD.pdf',4,3)
    figure(2)
        dB = 10*log10(P/max(P));
        plot(f,dB)
        texText('f [Hz]','x')
        texText('mean power [dB]','y')
        ylim([-25,1])
        xlim([0,0.4])
        graphToPDF('../buoyPSDlog.pdf',4,3)
end