% spikedatenanalyse
% zeigt alle in der Datei spikedaten_kurz gespeicherten Antworten eines
% Neurons zusammen mit dem Stimulus grafisch an.

load spikedaten_kurz

bin=0.1;         %samplerate 10000Hz entspricht bins von 0.1ms

[zeilen spalten]=size(daten_kurz);

t=bin:bin:bin*spalten;  %Zeitvektor

%Stimulus plotten:
subplot(2,1,2)
plot(t,stimulus_kurz)   %Stimulus fuer alle Durchlaeufe gleich
axis([0 1000 -1 3])     %damit man den Stimulus besser sieht

%Gleiche Skalierung fuer alle trials rausfinden:
mini=min(daten_kurz(:));   %Minimum in allen trials
mini=10*(floor(mini/10));  %Abrunden auf naechstkleineren Zehner
maxi=max(daten_kurz(:));   %Maximum in allen trials
maxi=10*(ceil(maxi/10));   %Aufrunden auf naechstgroesseren Zehner

%Probeweise Schwelle festlegen:
schwelle=2.*ones(1,spalten);

%Alle trials einmal zusammen mit der Schwelle ansehen:
for i=1:1:zeilen
    subplot(2,1,1)
    plot(t,daten_kurz(i,:),'b',t,schwelle,'r')
    axis([0 1000 mini maxi])
    pause
end

