Psth

PSTH Computes the peri-stimulus time histogram for neural spike times. The routine plots the trial averaged spike rate as a function of time.
R = PSTH(TIMES, BINSIZE, FS,NTRIALS,TRIALLEN)
R = PSTH(TIMES, BINSIZE, FS,NTRIALS,TRIALLEN ,AXESHANDLE)
TIMES - spike times (samples)
BINSIZE - binwidth (ms)
FS - sampling rate (hz)
NTRIALS - number of trials
TRIALLEN - length of a trial (samples)
R - spike rate (spikes/s)

Can handle large datasets and uses Matlab's built-in histc function for optimized performance.

An example:
%spike times can be specified in continuous time
%here we have 3 trials and a trial length of 1000 samples
t = [10, 250, 900, 1300, 1600, 2405, 2900];

%the same spike times can also be specified per trial
t2 =[10, 250, 900, 300, 600, 405, 900];
r = psth(t,10,1000,3,1000) ;
r2 = psth(t2,10,1000,3,1000);