wolffd@0: function [h,out] = datawavstem(wavin, varargin)
wolffd@0: %[h,out] = datawavplot(data)
wolffd@0: 
wolffd@0: 
wolffd@0: [sr,miny,maxy,minx,maxx,unused]=process_options(varargin,...
wolffd@0:     'sr',44800,'miny',-1,'maxy',-1,'minx',-1,'maxx',-1);
wolffd@0: out= wavin;
wolffd@0: 
wolffd@0: m=min(wavin);
wolffd@0: ma=max(wavin);
wolffd@0: add=(ma-m)*0.1;
wolffd@0: if miny ==-1
wolffd@0:     miny=m-add;
wolffd@0: end
wolffd@0: if maxy==-1
wolffd@0:     maxy=ma+add;
wolffd@0: end
wolffd@0: wavlens=round(length(wavin)/sr*100)/100;
wolffd@0: timetbl=([0:length(wavin)-1]./length(wavin))*wavlens;
wolffd@0: h=stem(timetbl,out,'r');
wolffd@0: 
wolffd@0: if minx ==-1
wolffd@0:     minx=0;
wolffd@0: end
wolffd@0: if maxx==-1
wolffd@0:     maxx=wavlens;
wolffd@0: end
wolffd@0: axis([minx maxx miny maxy]);
wolffd@0: set(gca,'FontSize',8);
wolffd@0: h=get(0,'CurrentFigure');
wolffd@0: axis xy;
wolffd@0: xlabel('Time[s]','FontSize',8)
wolffd@0: ylabel('Amplitude','FontSize',8)