Изменение частоты среза ФНЧ

Добрый день.

Построил АЧХ и ФЧХ эллиптического фильтра. Все вроде норм:

[z, p, k] = ellipap(5, 0.5, 20);
w = 0:0.01:2;
[b, a] = zp2tf(z, p, k);
h = freqs(b, a, w);
figure(1)
semilogy (w, abs(h))
figure(2)
plot (w, unwrap(angle(h)))

Но когда попытался масштабировать частотную ось, результат вышел совсем иной:

[z, p, k] = ellipap(5, 0.5, 20);
[b, a] = zp2tf(z, p, k);
f0 = 2e3;
[b, a] = lp2lp(z, p, 2*pi*f0);
f = 0:1:5e3;
h = freqs(b, a, 2*pi*f);
figure(1)
semilogy (f, abs(h))
figure(2)
plot (f, unwrap(angle(h)))

Можете пояснить в чем ошибка? При построении опирался на примеры из Сергиенко, Цифровая обработка сигналов, стр. 118-124.