function AP_Zaehler_neu=spikezaehler(Daten,Schwelle)

[zeilen spalten]=size(Daten); % zeilen=10 mit je Spalten=10000


% Idee-Sammlung: 

% 1: finde erstmal alle Punkte, die Ÿberschwellig sind
%    da "find" spaltenweise durchsucht, transponiere ich die Daten,
%    damit ich den "Mess_Runde"-Index in aufsteigender Reihenfolge habe

[Mess_Punkt  Mess_Runde]=find(Daten'>Schwelle);


% produziere eine 2-Spalten-Matrix mit RundenNummer in Spalte eins und MessPunkt
% Nummer in Splate zwei.
MUP=[Mess_Runde Mess_Punkt];

AP_Zaehler=[];
Counter=0;

% Idee: PrŸfe den Spalteneintrag, der in Daten auf ein Ÿberschwelliges
% Element folgt auf Unterschwelligkeit


% fŸr i=1 bis zeilenanzahl MUP:
% gehe in Daten in die Zeile deren Index an Position 1-MUP und die Spalte+1,
% deren Index an Position 2-MUP steht, lese diesen Wert aus Daten aus und schaue, 
% ob er unterschwellig ist. Ist er das, setze ich Counter um 1 hoch.
for i=1:length(Mess_Runde)
   if Daten(MUP(i,1),MUP(i,2)+1)<Schwelle
      Counter=Counter+1;
      AP_Zaehler( MUP(i,1),1)=Counter; % schreibe den aktuellen Wert von Counter an die
   end                                 % Indexposition von AP_Zaehler, die durch 
                                       % die erste Spalte von MUP
                                       % vorgegeben wird
end  

% AP_Zaehler enthŠlt an Position 1 die Spikeanzahl aus Zeile 1 der Daten
% jedoch an Position 2 die Spikeanzahl aus Zeile 2 + Zeile1

% darum wird ein neuer Zaehler kreiert, bei dem der Wert der vorhergehenden
% Zeile jeweils abgezogen wird
AP_Zaehler_neu(1)=AP_Zaehler(1); 
for i=1:9
  AP_Zaehler_neu(i+1)= AP_Zaehler(i+1)- AP_Zaehler(i);
end
AP_Zaehler_neu'
MW=mean(AP_Zaehler_neu) %Ausgabe Mittelwert
STD=std(AP_Zaehler_neu) % Ausgabe Standardabweichung 
Counter                 % Ausgabe Spike Zahl



    

    










    