clc; close all;  clear all;

%------------------------Messdaten laden-----------------------
[filename] = uigetfile({'*.mat'});
load(filename)
%------------------------Werksangaben Porsche-----------------------
c_w=0.36;                 % CW Wert; Beachte zusätzliche Fahrzeuganbauten 
A=2.81;                    % [m^2] Querschnittsfläche bzw. Strinfläche
a=3.27;                   % Achsenübersetzung

%------------------------Schleife Bestimmung eps und i -----------------------
eps1=1.2;
eps2=0.23;
eps3=0.15;
i1=4.91;
i2=2.81;
i3=1.84; 
anzahlWerte = size(G, 1);

for k = 1:anzahlWerte
    if ( G(k) == 1)
        it(k) = i1*a;
    end 
    if ( G(k) == 2)
        it(k) = i2*a;
    end 
    if ( G(k) == 3)
        it(k) = i3*a;
    end
end 
for k = 1:anzahlWerte
    if ( G(k) == 1)
        epst(k) = eps1*a;
    end 
    if ( G(k) == 2)
        epst(k) = eps2*a;
    end 
    if ( G(k) == 3)
        epst(k) = eps3*a;
    end
end 
i = it.';
eps = epst.';
%-------------------------- Konstanten und abgeschätzte Werte ------------------------ 

m=2350+2*87*(1+ eps);       % [kg] Fahrzeugmasse mit Einbauten plus 2 Studis 
mue_R=0.024;                % [ ] Rollwiderstand HTW-Prüfgelände
r_dyn=0.358;                % [m] dyn. Reifenhalbdurchmesser 
                            % Berechnung online uter: http://www.flat4.de/calc_tyr.htm
mb = m.^(-1);                            

V_anf = 16;               %in km/h                      
g=9.81;                   % [m/s^2] Erdbeschleunigung 
rho_L=1.2;                % [kg/m^3] Dichte Luft; beachte Fehler über Temperatureinfluss
ms = 1000;                %Umrechnungsfaktor ms in s

% %Daten Kupplung Verbrenner
% J1v = 1; %Massenträgheitsmoment Antriebsseite in kgm^2
% J2v = 5; %Massenträgheitsmoment der Abtriebseite in kgm^2
% myv = 0.3; %Gleitreibkoeefizient 
% rav = 15; %Außenradius Kupplungsbelag in cm
% riv = 10; %Innenradius Kupplungsbelag in cm
% FKv = 5000; %Anpresskraft in N
% %tkv == tm_V; %Zeit Anpresskraft
% 
% 
% %Daten Kupplung E_Maschine
% J1e = 1; %Massenträgheitsmoment Antriebsseite in kgm^2
% J2e = 5; %Massenträgheitsmoment der Abtriebseite in kgm^2
% mye = 0.3; %Gleitreibkoeefizient 
% rae = 15; %Außenradius Kupplungsbelag in cm
% rie = 10; %Innenradius Kupplungsbelag in cm
% FKe = 5000; %Anpresskraft in N
% %tke == tm_E; 




% Bereich Totzone Haftreibung
epsilon_D=0.03; % m/s muss überschritten werden, um Rollreibung zu genenrieren
t = interpft(t_G, tv_m); 
T_sim=max(t(:,1)); % Simulationsdauer von Messdaten abhängig

disp('Die Ausführung in Simulink kann gestartet werden!')