z = [0.01];     % Zähler der ÜF
n = [1 0.02];   % Nenner der ÜF
Gs = tf(z,n);   % G = z/n wird als Übertragungsfunktion definiert
%
sim('Sprung',1000); % Aufruf und Ausführung der Simulation "Sprung.mdl"
% mit Simulationszeit t = 1000s
%
figure(1) 
    hold on;
    plot(step(Gs),'r')               % Sprungantwort aus step()-Befehl
    plot(Sprung(:,1),Sprung(:,2),... % Einheitssprung und Sprungantwort aus
        Sprung(:,1),Sprung(:,3))     % der Simulation mit "Sprung.mdl"
    grid on;
    title('Sprungantwort: Regelstrecke ohne PID-Regler')
    xlabel('Zeit in s')
    ylabel('Amplitude')
    legend('Sprungantwort aus step()',...
        'Einheitssprung',...
        'Sprungantwort aus Simulation',...
        'location','NorthEast')
    axis([0, 100, 0, 1.5])