%% wartosci ogolne
n1=7;
n2=9;
k=n1/n2
m1=n1
m2=4*n2
t=2*n2+4
obiekt=tf([k],[1,m2,m1,1],'IODelay',t)
%% charakterystyki czasowe
figure(1)
impulse(obiekt)
figure(2)
step(obiekt)
%% kupfmuller
s=tf('s')
transkupf=tf([0.778],[10.8 1])
modelkupf=transkupf*exp(-24.5*s)
figure(3)
step(obiekt)
hold on
step(modelkupf)
hold off
figure(4)
impulse(obiekt)
hold on
impulse(modelkupf)
hold off