% Eingabe der Systemparameter
omega_n = input('Geben Sie die Kreisfrequenz omega_n ein: ');
K_P = input('Geben Sie den Verstärkungsfaktor K_P ein: ');
N = 500;

% Dämpfungsgrade
D = [0.1, 0.3, 0.7, 1.0, 2.0, 3.0];

% Zeitachse definieren
t = linspace(0, 4, N); % Anpassen für angemessenen Zeitbereich

% Erstellen der Sprungantworten für jeden Dämpfungsgrad
figure; % Erstellt ein neues Figure-Fenster
hold on; % Hält das aktuelle Figure-Fenster für mehrere Plots offen
for i = 1:length(D)
    % Berechnung der Systemantwort
    sys = tf(K_P * [omega_n^2], [1, 2*D(i)*omega_n, omega_n^2]);
    [y, t] = step(sys, t);
    
    % Zeichnen der Sprungantwort
    plot(t, y, 'DisplayName', sprintf('D = %.1f', D(i)));
end

% Anpassungen des Diagramms
xlabel('Zeit (s)');
ylabel('Ausgang');
title('Sprungantworten eines PT2-Systems für verschiedene Dämpfungsgrade');
legend('show'); % Zeigt die Legende an
grid on; % Fügt ein Gitter hinzu
hold off; % Beendet den "hold on"-Status