Разработайте аналог 5-го порядка lowpass фильтр Бесселя приблизительно с постоянной групповой задержкой до 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)')