30 lines
983 B
Matlab
30 lines
983 B
Matlab
% 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 |