Частотная характеристика фильтра Бесселя Lowpass

Разработайте аналог 5-го порядка lowpass фильтр Бесселя приблизительно с постоянной групповой задержкой до 104 rad/s. Постройте значение и фазовые отклики фильтра с помощью freqs.

[b,a] = besself(5,10000);
freqs(b,a)

Вычислите ответ групповой задержки фильтра как производная развернутого фазового отклика. Постройте групповую задержку, чтобы проверить, что это приблизительно постоянно до частоты среза.

[h,w] = freqs(b,a,1000);
grpdel = diff(unwrap(angle(h)))./diff(w);

clf
semilogx(w(2:end),grpdel)
xlabel('Frequency (rad/s)')
ylabel('Group delay (s)')

Для просмотра документации необходимо авторизоваться на сайте