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