Revision 1def5532c3cefeba955129117e8c23abb7643cf8 authored by Max's Office PC on 20 November 2019, 03:30:51 UTC, committed by Max's Office PC on 20 November 2019, 03:30:51 UTC
plotSpikeDensityfunction.m
function sdf = plotSpikeDensityfunction(spikeTimes, kernel_sd)%% Spike Density Function
%% Get Data
ephysSampleRate = 50000;
spikeTimes = round(spikeTimes * ephysSampleRate);
binned = zeros(spikeTimes(end),1);
for i = 1:length(spikeTimes)
binned(spikeTimes(i)) = 1;
end
%% Kernel
edges = -5 * kernel_sd : 1/ephysSampleRate : 5 * kernel_sd;
kernel = normpdf(edges, 0, kernel_sd);
sdf = conv(binned, kernel);

Computing file changes ...