
% Tworzenie modelu
s = tf('s');

%% Czlon inercyjny I rzedu 1/bs+c
A=0;
B=1;
C=1;
sys = (1/(A*s^2+B*s+C))
figure(1)
subplot(2,1,1)
bode(sys,s);
hold on;
subplot(2,1,2)
nyquist(sys,s)
%% Czlon inercyjny II rzedu 1/(a^2*s+bs+c)
A=1;
B=2;
C=1;
sys = (1/(A*s^2+B*s+C))
figure(2)
subplot(2,1,1)
bode(sys,s);
hold on;
subplot(2,1,2)
nyquist(sys,s)
%% Czlon oscylacyjny
A=1;
B=1;
C=1;
sys = (1/(A*s^2+B*s+C))
figure(3)
subplot(2,1,1) 
bode(sys,s);
hold on;
subplot(2,1,2)
nyquist(sys,s)
%% Cz³on calkujacy
%% Cz³on ró¿niczkuj¹cy
%% Cz³on calkujacy z inercj¹
A=1;
B=1;
C=0;
sys = (1/(A*s^2+B*s+C))
figure(6)
subplot(2,1,1)
bode(sys,s);
hold on;
subplot(2,1,2)
nyquist(sys,s)
%% Cz³on ró¿niczkuj¹cy z inercj¹
A=0;
B=1;
C=1;
sys = (s/(A*s^2+B*s+C))
figure(7)
subplot(2,1,1)
bode(sys,s);
hold on;
subplot(2,1,2)
nyquist(sys,s)
%% Cz³on opóŸniaj¹cy
A=0;
B=1;
C=0;
sys = exp(-s*B);
figure(8)
subplot(2,1,1)
bode(sys,s);
hold on;
subplot(2,1,2)
nyquist(sys,s)