diff userProgramsTim/plotIFRAN.m @ 38:c2204b18f4a2 tip

End nov big change

--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/userProgramsTim/plotIFRAN.m	Mon Nov 28 13:34:28 2011 +0000
@@ -0,0 +1,54 @@
+function plotIFRAN(data_matrix,start_time,end_time,sfreq,BFlist,plothandle)
+
+%function that produces a plot similar to the plots in Secker-Walker, JASA
+%1990 Fig.2 (Neurograms) from a IFRAN data matrix
+%Tim Jürgens, January 2011
+
+%plot 50ms right from the middle of the signal
+time_axis = [0:1/sfreq:(size(data_matrix,2)-1)/sfreq];
+%start_time = size(data_matrix,2)/2/sfreq-0.025; %start time to plot in s
+%temporal_length = min([0.05 size(data_matrix,2)/sfreq]); %length to plot in s, at least 50ms or the length of the stimulus
+[tmp, start_time_index] = min(abs(start_time-time_axis));
+%end_time = start_time+temporal_length;
+[tmp, end_time_index] = min(abs(end_time-time_axis));
+
+
+%smoothing using a 4-point hamming window
+hamm_window = hamming(4);
+for iCounter = 1:size(data_matrix,1)
+    for jCounter = 2:length(time_axis)-2 %start with 2 and end with 2 samples
+        %less the length of time_axis in order not to get in conflict with the length of
+        %the hamm_window
+        smoothed_data_matrix(iCounter,jCounter-1) = ...
+            data_matrix(iCounter,(jCounter-1):(jCounter+2))*hamm_window./sum(hamm_window);
+    end
+end
+
+smoothed_time_axis = time_axis(2:end-2);
+
+if ~exist('plothandle'), plothandle=figure; end
+
+verticalshift = 1000; %vertical shift of single time series in z-coordinate units
+Tickvector = [];
+TickLabels = [];
+for iCounter = 1:size(data_matrix,1)
+    set(gcf,'Currentaxes',plothandle);
+    
+    %multiply by 1000 to set abscissa to ms units
+    verticalposition = -verticalshift*(iCounter-1);
+    plot(1000.*smoothed_time_axis(start_time_index:end_time_index), ...
+        smoothed_data_matrix(size(data_matrix,1)-iCounter+1,start_time_index:end_time_index)+verticalposition, ...
+        'k','LineWidth',0.5);
+    xlim([1000*start_time 1000*end_time]);
+    if mod(iCounter,5) == 1 %set best frequency as a label every 5 channels
+        Tickvector = [Tickvector verticalposition];
+        TickLabels = [TickLabels; BFlist(size(data_matrix,1)-iCounter+1)];
+    end
+    hold on;
+end
+set(gca,'yTick',Tickvector(end:-1:1),'yTickLabel',num2str(TickLabels(end:-1:1),'%4.0f'));
+ylabel('Best Frequency of Filter (Hz)');
+xlabel('Stimulus Time (ms)');
+set(gca,'YLim',[-(size(data_matrix,1))*verticalshift verticalshift*1.5]);
+box on;
+hold off;
\ No newline at end of file