katkost@3: function plot_peaks(Xmag,Xphase,ploc,t)
katkost@3: figure(1);
katkost@3: subplot(2,1,1)
katkost@3: plot(Xmag);
katkost@3: hold on;
katkost@3: plot(ploc,Xmag(ploc),'r*');
katkost@3: plot([1,length(Xmag)],[t,t],'g');
katkost@3: hold off;
katkost@3: subplot(2,1,2);
katkost@3: plot(Xphase);
katkost@3: hold on;
katkost@3: plot(ploc,Xphase(ploc),'r*');
katkost@3: hold off;
katkost@3: pause;