%Destille
%Kopieren der relevanten Daten aus Gesamtvektor
%Relevant heisst Y-wert größer als Variable limit
limit=8;
n=1;
count=500;   %Zähler für Indexanzeige, zeigt jeden 500.Wert
minwert(1)=0;
mincount=1;
toggle=0;
for i=1:length(Track1)
    
    if Track1(i)>limit
        T(n)=Track1(i);
        n=n+1;
        if toggle==1
            toggle=0;
        end
        %Rückwärtszähler für Indexanzeige, Ausgabe von i wenn 0 erreicht
        count=count-1;
        if count==0
            i
            count=500;
        end
    end
    %Finden der Integratinsintervalle (Minima für T finden)
    if Track1(i)<-25 && toggle==0 
        minwert(mincount)=n;
        mincount= mincount+1;
        toggle=1;
    end
end
%Drei Einzelkurven erzeugen, Grenzen sind die gefundenen Minima für T
k=1;
for m=1:minwert(1)
    TA(m)=T(m);
end
k=1;
for m=(minwert(1)+1):minwert(2)
    TB(k)=T(m);
    k=k+1;
end
k=1;
for m= (minwert(2)+1): (minwert(3)-1)
    TC(k)=T(m);
    k=k+1;
end
%Berechnen der Kurvenwerte für TA, TB und TC
%Intervallschrittweite = 1/Abtastfrequenz
Intvall=1/(5.12e+004);
%Werte für ersten Peak bzw. Bohrung
TAmax=max(TA);
TAintegr=sum(TA)*Intvall;
TAmean=mean(TA);
%Werte für zweiten Peak bzw. Bohrung
TBmax=max(TB);
TBintegr=sum(TB)*Intvall;
TBmean=mean(TB);
%Werte für dritten Peak bzw. Bohrung
TCmax=max(TC);
TCintegr=sum(TC)*Intvall;
TCmean=mean(TC);
%Erzeugen der Ergebnistabelle als Datei
%Reihenfolge der Werte max mean integr.
TAergeb=[TAmax TAmean TAintegr];
TBergeb=[TBmax TBmean TBintegr];
TCergeb=[TCmax TCmean TCintegr];
%
%Datenausgabe im Commandfenster
fprintf('       max         mean    integral \n');
fprintf('TA   %8.2f  %8.2f  %8.2f \n',TAergeb);
fprintf('TB   %8.2f  %8.2f  %8.2f \n',TBergeb);
fprintf('TC   %8.2f  %8.2f  %8.2f \n',TCergeb);