wp = 2*pi*1e5;
wz = wp;
Qp = 2;
f  = logspace(4, 7, 100);
w  = 2*pi*f;

s = tf('s');

h1 = (1 + (s/wz)^2) / (1 + s/(wp*Qp) + (s/wp)^2);
h2 = 1              / (1 + s/(wp*Qp) + (s/wp)^2);

[m, p] = bode(h1, w);
m = m(1, :);
m(m < 1e-12) = 1e-12;
semilogx(f, 20*log(m), 'b');

axis([ min(f) max(f) -50 20 ]);
xlabel('Frequency  [Hz]');
ylabel('Magnitude  [dB]');
fixfig;

figure;
pzmap(h1);
hold on;
t = linspace(0, 2*pi, 100);
x = cos(t);
y = sin(t);
plot(x*wp, y*wp, 'g');
wp = 1.1*wz;
axis([ -wp wp -wp wp ]);
fixfig;



