tomwalters@0: % method of class @signal
tomwalters@0: % 
tomwalters@0: %   INPUT VALUES:
tomwalters@0: %  
tomwalters@0: %   RETURN VALUE:
tomwalters@0: %
tomwalters@0: % 
tomwalters@0: % (c) 2003-2008, University of Cambridge, Medical Research Council 
tomwalters@0: % Maintained by Tom Walters (tcw24@cam.ac.uk), written by Stefan Bleeck (stefan@bleeck.de)
tomwalters@0: % http://www.pdn.cam.ac.uk/cnbh/aim2006
tomwalters@0: % $Date: 2008-06-10 18:00:16 +0100 (Tue, 10 Jun 2008) $
tomwalters@0: % $Revision: 585 $
tomwalters@0: 
tomwalters@0: function rethandle=plot(sig,border,stil,ax)
tomwalters@0: % usage: handle=plot(sig,border,stil)
tomwalters@0: % plots the class signal. This plot works a little bit different from usual.
tomwalters@0: % Parameter can either be: plot(signal,stil) or plot(signal,border) or plot(signal,border,stil)
tomwalters@0: % the border can consist of times and y-value or only the times that should be plotted (in seconeds!)
tomwalters@0: 
tomwalters@0: if nargin<4
tomwalters@0:     ax=gca;
tomwalters@0: end
tomwalters@0: if nargin<3
tomwalters@0:     stil='b-';
tomwalters@0: end
tomwalters@0: 
tomwalters@0: if isfield(vargin{1}
tomwalters@0: 
tomwalters@0: 
tomwalters@0: sr=getsr(sig);
tomwalters@0: durationbin=size(sig.werte,1);
tomwalters@0: duration=durationbin/sr;
tomwalters@0: if nargin<2
tomwalters@0:     t1=1;
tomwalters@0:     t2=durationbin;
tomwalters@0:     border=[t1 t2 min(sig.werte(t1:t2)) max(sig.werte(t1:t2))];
tomwalters@0: end
tomwalters@0: 
tomwalters@0: if isstruct(border) || isempty(border)
tomwalters@0:     options=border;
tomwalters@0:     t1=1;
tomwalters@0:     t2=durationbin;
tomwalters@0:     border=[t1 t2 min(sig.werte(t1:t2)) max(sig.werte(t1:t2))];
tomwalters@0: else
tomwalters@0:     options=[];
tomwalters@0:     if ischar(border)
tomwalters@0:         stil=border;
tomwalters@0:         t1=1;
tomwalters@0:         t2=durationbin;
tomwalters@0:         border=[t1 t2 min(sig.werte(t1:t2)) max(sig.werte(t1:t2))];
tomwalters@0:     else
tomwalters@0:         if isempty(border)
tomwalters@0:             t1=1;
tomwalters@0:             t2=durationbin;
tomwalters@0:             border=[t1 t2 min(sig.werte(t1:t2)) max(sig.werte(t1:t2))];
tomwalters@0:         else
tomwalters@0:             t1=border(1);
tomwalters@0:             t2=border(2);
tomwalters@0:             x1=time2bin(sig,t1);
tomwalters@0:             x2=time2bin(sig,t2);
tomwalters@0:             nr=size(border,2);
tomwalters@0:             if nr==2    % wenn nur die x-Werte angegeben werden
tomwalters@0:                 border=[x1 x2 min(sig.werte(x1+1:x2)) max(sig.werte(x1+1:x2))];
tomwalters@0:             else
tomwalters@0:                 border(1)=x1;
tomwalters@0:                 border(2)=x2;
tomwalters@0:             end
tomwalters@0:         end
tomwalters@0:     end
tomwalters@0: end
tomwalters@0: 
tomwalters@0: if ~isfield(options,'is_log');
tomwalters@0:     is_log=0;
tomwalters@0: else
tomwalters@0:     is_log=options.is_log;
tomwalters@0: end
tomwalters@0: 
tomwalters@0: if isfield(options,'minimum_time');
tomwalters@0:     border(1)=time2bin(sig,options.minimum_time);
tomwalters@0:     minimum_time=options.minimum_time;
tomwalters@0: else
tomwalters@0:     minimum_time=bin2time(sig,border(1));
tomwalters@0: end
tomwalters@0: 
tomwalters@0: if isfield(options,'maximum_time');
tomwalters@0:     border(2)=time2bin(sig,options.maximum_time);
tomwalters@0:     maximum_time=options.maximum_time;
tomwalters@0: else
tomwalters@0:     maximum_time=bin2time(sig,border(2));
tomwalters@0: end
tomwalters@0: 
tomwalters@0: 
tomwalters@0: % if the time scale is reversed (time from left to right)
tomwalters@0: if isfield(options,'time_reversed');
tomwalters@0:     time_reversed=options.time_reversed;
tomwalters@0: else
tomwalters@0:     time_reversed=0;
tomwalters@0: end
tomwalters@0: 
tomwalters@0: start_time=getminimumtime(sig);
tomwalters@0: 
tomwalters@0: min_x_screen=border(1);   % einer wird abgezogen damit wir bei Null beginnen, nicht beim ersten bin, was blöde aussieht
tomwalters@0: max_x_screen=border(2);
tomwalters@0: minshowy=border(3);
tomwalters@0: maxy=border(4);
tomwalters@0: 
tomwalters@0: xvals=getxvalues(sig);
tomwalters@0: % if we are dealing with milliseconds (very often), then the units are in
tomwalters@0: % ms as well
tomwalters@0: if strfind(sig.unit_x,'(ms)')>0
tomwalters@0:     multiplier=1000;
tomwalters@0: else
tomwalters@0:     multiplier=1;
tomwalters@0: end
tomwalters@0: 
tomwalters@0: 
tomwalters@0: xvals=xvals.*multiplier;
tomwalters@0: yvals=getvalues(sig);
tomwalters@0: 
tomwalters@0: % this is the plotting command:
tomwalters@0: handle=plot(ax,xvals,yvals);
tomwalters@0: 
tomwalters@0: if time_reversed
tomwalters@0:     set(ax,'XDir','reverse')   % turn them around, because the higher values shell end on the right
tomwalters@0: else
tomwalters@0:     set(ax,'XDir','normal')   % normale ausrichtung
tomwalters@0: end
tomwalters@0: 
tomwalters@0: if is_log
tomwalters@0:     set(ax,'XScale','log')
tomwalters@0:     t=minimum_time*multiplier;
tomwalters@0:     tix=[t 2*t 4*t 8*t 16*t 32*t 64*t 128*t 256*t 512*t 1024*t];
tomwalters@0:     if ~isempty(sig.x_tick_labels)
tomwalters@0:         ti=sig.x_tick_labels;
tomwalters@0:         set(ax,'XTicklabel',ti);
tomwalters@0:     end
tomwalters@0:     set(ax,'XTick',tix);
tomwalters@0: else % its linear
tomwalters@0:     set(ax,'XScale','linear')
tomwalters@0: end
tomwalters@0: 
tomwalters@0: miny=border(3);
tomwalters@0: maxy=border(4);
tomwalters@0: y=[miny maxy];
tomwalters@0: if miny==maxy
tomwalters@0:     maxy=miny+1;
tomwalters@0:     miny=miny-1;
tomwalters@0: end
tomwalters@0: axis([minimum_time*multiplier maximum_time*multiplier miny*1.05 maxy*1.05]);
tomwalters@0: 
tomwalters@0: 
tomwalters@0: 
tomwalters@0: xlabel(sig.unit_x);
tomwalters@0: ylabel(sig.unit_y);
tomwalters@0: title(sig.name,'Interpreter','none');
tomwalters@0: 
tomwalters@0: if nargout==1
tomwalters@0:     rethandle=handle;
tomwalters@0: end
tomwalters@0: 
tomwalters@0: return