automatyka2

k1 = 0.1;
T1 = 2;

k2 = 0.5;
T2 = 1;

% obiekt inercyjny I rzędu

licz_1 = [0, k1];
mian_1 = [T1, 1];

licz_2 = [0, k2];
mian_2 = [T2, 1];

[y_1, x_1, czas_1] = step(licz_1, mian_1)
[y_2, x_2, czas_2] = step(licz_2, mian_2)
plot(czas_1, y_1, czas_2, y_2)
title("Obiekt inercyjny I rzędu, odp. skokowa")
legend(['k=' num2str(k1) ', T=' num2str(T1)], ['k=' num2str(k2) ', T=' num2str(T2)]);


licz_1 = [0, k1];
mian_1 = [T1, 1];

licz_2 = [0, k2];
mian_2 = [T2, 1];

[y_1, x_1, czas_1] = impulse(licz_1, mian_1)
[y_2, x_2, czas_2] = impulse(licz_2, mian_2)
figure()
plot(czas_1, y_1, czas_2, y_2)
title("Obiekt inercyjny I rzędu, odp. impulsowa")
legend(['k=' num2str(k1) ', T=' num2str(T1)], ['k=' num2str(k2) ', T=' num2str(T2)]);


ksi = 2;

licz_1 = [0, 0, k1];
mian_1 = [T1^2, 2*ksi*T1, 1];

licz_2 = [0, k2];
mian_2 = [T1^2, 2*ksi*T1, 1];

[y_1, x_1, czas_1] = step(licz_1, mian_1)
[y_2, x_2, czas_2] = step(licz_2, mian_2)
figure()
plot(czas_1, y_1, czas_2, y_2)
title(strcat("Obiekt inercyjny II rzędu (oscylacyjny), odp. skokowa, ksi =", num2str(ksi)))
legend(['k=' num2str(k1) ', T=' num2str(T1)], ['k=' num2str(k2) ', T=' num2str(T2)]);


licz_1 = [0, 0, k1];
mian_1 = [T1^2, 2*ksi*T1, 1];

licz_2 = [0, k2];
mian_2 = [T1^2, 2*ksi*T1, 1];

[y_1, x_1, czas_1] = impulse(licz_1, mian_1)
[y_2, x_2, czas_2] = impulse(licz_2, mian_2)
figure()
plot(czas_1, y_1, czas_2, y_2)
title(strcat("Obiekt inercyjny II rzędu (oscylacyjny), odp. impulsowa, ksi =", num2str(ksi)))
legend(['k=' num2str(k1) ', T=' num2str(T1)], ['k=' num2str(k2) ', T=' num2str(T2)]);

Ti = 1;
licz_1 = [0, 0, k1];
mian_1 = [Ti*T1, Ti, 0];

licz_2 = [0, k2];
mian_2 = [Ti*T2, Ti, 0];


[y_1, x_1, czas_1] = step(licz_1, mian_1)
[y_2, x_2, czas_2] = step(licz_2, mian_2)
figure()
plot(czas_1, y_1, czas_2, y_2)
title(strcat("Obiekt całkujący rzeczywisty, odp. impulsowa, Ti =", num2str(Ti)))
legend(['k=' num2str(k1) ', T=' num2str(T1)], ['k=' num2str(k2) ', T=' num2str(T2)]);








#####################






clear

k1 = 0.1;
T1 = 2;

k2 = 0.5;
T2 = 1;

% obiekt inercyjny I rzędu

licz_1 = [0, k1];
mian_1 = [T1, 1];

licz_2 = [0, k2];
mian_2 = [T2, 1];

bode(licz_1, mian_1)
hold on
bode(licz_2, mian_2)
legend(['k=' num2str(k1) ', T=' num2str(T1)], ['k=' num2str(k2) ', T=' num2str(T2)]);


figure()
nyquist(licz_1, mian_1);
hold on
nyquist(licz_2, mian_2);
legend(['k=' num2str(k1) ', T=' num2str(T1)], ['k=' num2str(k2) ', T=' num2str(T2)]);