%% Hauptskript
%Berechnung der optimalen Dämpfung bei Kompfort- und Sicherheitsauslegung
%bei einem 1/4- Fahrzeug
global M m ca cr ba vf Amp Psi OM
Vertikaldynamik_OptimaleDaempfung_Daten;

%% Erstellen der Straße
OM0=1;      % Wieso ist Omega0=1???
DeltaOM=(OMmax-OMmin)/(n-1);
OM=OMmin:DeltaOM:OMmax;
Phi=Phi0*OM0.*(OM./OM0).^(-2);
Amp=(2*Phi*DeltaOM).^(1/2);
%RandStream.setDefaultStream(RandStream('mt19937ar','seed',5489)); %Gleiche Zufallszahlen benutzen
Psi=2*pi*rand(size(OM));
u0=sum(Amp.*sin(Psi));

%% Definition der statischen Lasten 
Fa0=M*g; %statische Last der Aufbaufeder
Fr0=(M+m)*g; %statische Radlast

%% Anfangsbedingungen
y0=[u0;0;u0;0];
t=[0 te];

%% Berechnung der Varianzen für unterschiedliche Dämpferwerte
sigma_s=zeros(size(Iba));   % Variablen vorreservieren um Rechnung effizienter zu machen
sigma_k=zeros(size(Iba));   % Variablen vorreservieren um Rechnung effizienter zu machen

for h=1:length(Iba)
    ba=Iba(h);
    disp(['rechne für ba=',num2str(ba)])
    [t,y]=ode45(@VertikaldynamikOptimaleDaempfungFCN,t,y0);
    st=(y(:,1)-y(:,3))/sR;  % Berechnung des normierten Federweges
    
end