%System o transmitancji L/M
%sys=tf(L,M)
%T=2,4,6

%inercja I rzedu k/Ts+1
%k=1;

%odpowiedz skokowa:
L=[1];
M1=[2 1];
M2=[4 1];
M3=[6 1];
t=0.01:0.1:50;

sys1=tf(L,M1);
sys2=tf(L,M2);
sys3=tf(L,M3);
subplot(3,1,1);
step(sys1,t);
hold on;
grid on;
step(sys2,t);
step(sys3,t);


%odpowiedź impulsowa:
L=[1];
M1=[2 1];
M2=[4 1];
M3=[6 1];
t=0.01:0.1:50;
sys1=tf(L,M1);
sys2=tf(L,M2);
sys3=tf(L,M3);
subplot(3,1,2);
impulse(sys1,t);
hold on;
grid on;
impulse(sys2,t);
impulse(sys3,t);

%charkaterystyka amplitudowo fazowa
w = 0.01:0.01:100;

subplot(3,1,3);

nyquist(L,M1,w);


grid on;
hold on;