clc
clear all
close all

load('Datenlogger_auf_x-Achse_aus_Matlab_org_Daten')

%messwerte_mean=messwerte_mean(messwerte(:,1));   % Mittelwert der Messwerte abziehen, Verschiebeung auf 0 Achse

MaxLog=max(messwerte_mean)
sqrtlinks=MaxLog/sqrt(2)

messwerte_mean_norm=messwerte_mean/0.42;        %Umrechnungsfakrot Filtec D169 Fernbereich  0.42mV/micro m     
messwerte_mean_norm10=messwerte_mean_norm*10;

%figure('Units','centimeters','Position',[10,10,20,8],'Color','white');
plot(zeitwerte,messwerte_mean_norm10);
xlabel('Zeit t [s]');
ylabel('Amplitude [mm]');
xlim([75 85]);
ylim([-15 15]);
grid on; box off;
hold on;

Maximun_norm=max(messwerte_mean_norm)


pos_index=find(messwerte_mean_norm10>0);
neg_index=find(messwerte_mean_norm10<0);
%Speichere Werte in jeweilige Variable
zeitwerte_pos=zeitwerte(pos_index);
zeitwerte_neg=zeitwerte(neg_index);
pos_werte=messwerte_mean_norm10(pos_index);
neg_werte=abs(messwerte_mean_norm10(neg_index));

[pks_pos_zeit, zeitwerte_peak_pos]=findpeaks(pos_werte, zeitwerte_pos);
[pks_neg_zeit, zeitwerte_peak_neg]=findpeaks(neg_werte, zeitwerte_neg);

% for i=1:length(messwerte_mean_norm10)
%     index_pks_pos(i)=find(messwerte_mean_norm10==pks_pos(i));
%     %zeit_pks=zeitwerte(index_pks_pos(i));
% end

%t_pos=zeitwerte();
%t_neg=zeitwerte(index_peak_neg);

figure(2)
subplot(311)
plot(zeitwerte_peak_pos,pks_pos_zeit)
title('Positive Peaks')
xlim([75 85]);
ylim([0 15]);
grid on
subplot(312)
plot(zeitwerte_peak_neg,pks_neg_zeit)
xlim([75 85]);
ylim([0 15]);
title('Negative Peaks')
grid on
subplot(313)
plot(pks_pos_zeit)
hold all
plot(pks_neg_zeit)
%xlim([75 85]);
%ylim([0 15]);
title('Positive und negative Peaks übereinander')
hold off
grid on

%%
f=zeitwerte*0;
f1=74.93755;
f2=79.10075;
t1=79;
t2=81;
for i=1:length(zeitwerte)
    f(i)=(f2-f1)/(t2-t1)*zeitwerte(i)+f1-(f2-f1)/(t2-t1)*t1;
    
end

figure(3)
plot(f,messwerte_mean_norm10)
ylabel('Amplitude [mm]')
xlabel('Frequenz [Hz]')
title('Verlauf von Stahl aus Experiment über die Frequenz')


%% Neue Frequenzwertplots
f_pos=f(pos_index);
f_neg=f(neg_index);
[pks_pos_freq, f_peak_pos]=findpeaks(pos_werte, f_pos);
[pks_neg_freq, f_peak_neg]=findpeaks(neg_werte, f_neg);

figure(4)
plot(f_peak_pos,pks_pos_freq)