tomwalters@0: % generating function for 'aim-mat'
tomwalters@0: % 
tomwalters@0: %   INPUT VALUES:
tomwalters@0: %  
tomwalters@0: %   RETURN VALUE:
tomwalters@0: %
tomwalters@0: % 
bleeck@3: % (c) 2011, University of Southampton
bleeck@3: % Maintained by Stefan Bleeck (bleeck@gmail.com)
bleeck@3: % download of current version is on the soundsoftware site: 
bleeck@3: % http://code.soundsoftware.ac.uk/projects/aimmat
bleeck@3: % documentation and everything is on http://www.acousticscale.org
bleeck@3: 
tomwalters@0: 
tomwalters@0: function displaydualprofile(sai,options,frame_number,ax)
tomwalters@0: if nargin<4
tomwalters@0:     ax=gca;
tomwalters@0: end
tomwalters@0: 
tomwalters@0: 
tomwalters@0: % Test if the frame_number is available
tomwalters@0: if length(sai)<frame_number
tomwalters@0:     % no its not
tomwalters@0:     return
tomwalters@0: end
tomwalters@0: if ~(isfield(sai{frame_number}, 'interval_sum') & isfield(sai{frame_number}, 'frequency_sum'))
tomwalters@0:     return
tomwalters@0: end
tomwalters@0: 
tomwalters@0: % ?????? per Definition ????
tomwalters@0: minimum_time_interval=options.minimum_time_interval;  % in ms
tomwalters@0: maximum_time_interval=options.maximum_time_interval;
tomwalters@0: nr_labels = options.nr_labels;
tomwalters@0: 
tomwalters@0: % Normalize the profiles
tomwalters@0: fq_sum = sai{frame_number}.frequency_sum;
tomwalters@0: int_sum = sai{frame_number}.interval_sum;
tomwalters@0: if (getnrpoints(fq_sum)~=0)
tomwalters@0:     int_sum = int_sum/getnrpoints(fq_sum);
tomwalters@0: end
tomwalters@0: if (getnrpoints(int_sum)~=0)
tomwalters@0:     fq_sum = (fq_sum/getnrpoints(int_sum))*options.scalefactor*1.7;
tomwalters@0: end
tomwalters@0: 
tomwalters@0: cla;
tomwalters@0: %Plot both profiles into one figure
tomwalters@0: % frequency profile
tomwalters@0: fqp = getvalues(fq_sum)';
tomwalters@0: plot(ax,sai{frame_number}.channel_center_fq, fqp,'r');
tomwalters@0: hold on
tomwalters@0: 
tomwalters@0: % time interval profile
tomwalters@0: tip=getvalues(int_sum);
tomwalters@0: tip_x = bin2time(sai{frame_number}.interval_sum, 1:length(tip));            % Get the times
tomwalters@0: tip_x = tip_x((tip_x>=(minimum_time_interval/1000)) & tip_x<=(maximum_time_interval/1000));  
tomwalters@0: tip = tip(time2bin(sai{frame_number}.interval_sum,tip_x(1)):time2bin(sai{frame_number}.interval_sum,tip_x(end)));
tomwalters@0: % tip_x is in ms. Change to Hz
tomwalters@0: tip_x = 1./tip_x;
tomwalters@0: plot(tip_x, tip, 'b');
tomwalters@0: set(ax,'XScale','log');    
tomwalters@0: 
tomwalters@0: % Now lable it !
tomwalters@0: xlabel('Frequency [Hz]');
tomwalters@0: set(ax, 'YAxisLocation','right');
tomwalters@0: ti=50*power(2,[0:nr_labels]);
tomwalters@0: set(ax,'XTick', ti);
tomwalters@0: set(ax, 'XLim',[1000/maximum_time_interval sai{frame_number}.channel_center_fq(end)])
tomwalters@0: set(options.handles.checkbox6, 'Value',0);
tomwalters@0: set(options.handles.checkbox7, 'Value',0);
tomwalters@0: hold off
tomwalters@0: 
tomwalters@0: 
tomwalters@0: return