clear


m=500         %[kg]           Eigengewicht des Anhängers
g=9,81        %[m/s^2]        Erdbeschleunigung
c=50          %---            Federsteifigkeit
d=103         %---            Dämpfung
F=500         %[N]            Beladung des Anhängers
w_0=sqrt(c/m) %[1/s]          Eigenkreisfrequenz des ungedämpften Systems
delta=d/(2*m) %[rad/s]=[1/s]  Abklingkonstante
D=delta/w_0   %[-]            Dämpfungsgrad


dy_dt=@(t,y)[y(2);-(d/m)*y(2)-(c/m)*y(1)];


[t,y]=ode45(dy_dt,[0 100.0],[0.0 1.0]);



subplot(1,2,1);
plot(t,y(:,1),'r-',t,y(:,2),'b-');

title('Einmassenschwigner');
xlabel('time');
ylabel('Amplitude');
legend('y1(t) Aufbaumasse','y2(t) Fahrbahnanregung');

subplot(1,2,2);
plot(y(:,1),y(:,2));
title('Phasendiagramm');
xlabel('Auslenkung');
ylabel('Impuls');

plot(t,y,t,(F(t)-c*y(:,1))/m); 
legend('Weg','Geschwindigkeit','Beschleunigung');