harmonic_cursors_demo.m

% Load Example Data
load handel
t=linspace(0,length(y)/Fs,length(y));
%  Plot Time Series
h_fig=figure;
subplot(2,1,1)
plot(t,y);
harmonic_cursors(gca)

NFFT=round(length(y)/2);
[Pxx,F] = PSD(y,NFFT,Fs);

% Plot Power Spectral Density
subplot(2,1,2)
plot(F,Pxx);
set(gca,'XLim',[0 1500])
harmonic_cursors(gca)