A=[1 88888 3.9506e9]
B=3.9506e9
f=linspace(0,3e4,200);
w=2*pi*f;
H=freqs(B,A,w);
figure(1)
plot(f,abs(H));
grid
xlabel f
ylabel '|H|'
title 'Modulová charakteristika'
figure(2)
plot(f,phase(H));
grid
xlabel f
ylabel 'PH'
title 'Fázová charakteristika'
t=grdelays(B,A,w);
figure(3)
plot(f,t);
grid
xlabel f
ylabel 't'
title 'Skupinové zpoždìní'
