Mein MATLAB Forum

function Zug21n
global L m s h g
global fre
fre = 0.04*2*pi;
m = 84000;
s=19.5;L=2;h=4.4;
g = 9.81;
%-- Anfangswerte
x=[0]; z=[0];
dx = 0; dz = 0;
lam=[0];
y0 = [x dx z dz lam];
%--
M = eye(5);
M(5,5)=0;
dy = dgl(0,y0)
options=odeset('Mass',M,'RelTol',1.0e-6,'AbsTol',1.0e-6);
[t,y] = ode15s(@dgl, [0:0.025:200], y0, options);
plot(t,y(2),t,y(4))

function dy = dgl(t,y)
global m
dy = zeros(5,1);
x=y(1);dx=y(2);ddx=dy(2);
z=y(3);dz=y(4);ddz=dy(4);
lam=y(5);
%-- Rad-Schiene
Fx1 = lam*dGleis1(x); Fz1 = -lam;
%-- m*a = F;
dy(1) = dx;   % XPUNKT
dy(2) = (Antriebskraft(t)*cos(winkel(x))-Bremskraft(t)*cos(winkel(x))+Fx1)/m;   %XPUNKT PUNKT
dy(3) = dz;   %Y PUNKT
dy(4) = (Antriebskraft(t)*sin(winkel(x))-Bremskraft(t)*sin(winkel(x))+Fz1)/m;   %YPUNKT PUNKT
dy(5) = ddGleis1(x)*dx^2+dGleis1(x)*ddx-ddz;                                   % RANDBEDINGUNG 2 mal nach t abgeleitet


function z = ddGleis1(x)
global fre
z = -sin(fre*x)*fre*fre;
%z = 0*x;

function z = dGleis1(x)
global fre
z = cos(fre*x)*fre;
%z = 0*x;

function z = Gleis1(x)
global fre
z = sin(fre*x);
%z = 0*x;

function alpha = winkel(x)
alpha = atan(dGleis1(x));

function F = Antriebskraft(t)
%--
  F =  200000;
  if t>= 80
      F = 0;
  end
 
function F = Bremskraft(t)
%--
  F=0*t;
  if t>= 100
      F = 500000;
  end
  %F = F.*signum(-v);