clear all       % alle Variablen löschen

% Modellparameter
m = 1;              % kg Schlitten und Spulenmasse
L = 2e-3;           % H Induktivität der Spule
R = 2;              % Ohm Widerstand der Spule
FR = 5;             % N Reibingskraft de Führung
KF = 12;            % N/A Kraftkonstante des Motors
KV = 12;            % Vs/m Spannungskonstante des Motors
Tel = L/R;          % s elektrische Zeitkonstante der Wicklung
U0 = 12;            % V Speisespannung
C = 6.3;            % F Kondensatorkapazität
k = 100000;          % N/m Federkonstante


% Simulationsparameter definieren
simulation_opt=simset('Solver', 'ode5', 'FixedStep', 1e-4);
Tsim = 0.60;        % Simulationszeit in s

% Simulation starten
[x, y, t]=sim('tauchspulmotor_mod',Tsim, simulation_opt);

% Simulationsergebnisse darstellen
figure(1); hold off         % erste Grafik aufrufen, alte Darstellung löschen
plot (ts, xs, 'r--');       % Zeitverlauf der Position darstellen
hold on;                    
plot (ts, vs, 'k'); grid;        % Zeitverlauf der Geschwindigkeit darstellen, Gitterlinien


xlabel('t in s');
ylabel('Position x in mm, Geschwindigkeit v in m/s');
title('Anlaufsimulation des Tauchspulmotors');

figure(2); hold off         % zweite Grafik aufrufen
plot (ts, Is, 'b'); grid;        % Zeitverlauf des Stroms darstellen, Gitterlinien

xlabel('t in s');
ylabel('Strom in A');
title('Anlaufsimulation des Tauchspulmotors');