Разработайте аналоговый lowpass фильтр Бесселя пятого порядка с приблизительно постоянной групповой задержкой до рад/сек. Постройте график величины и фазовых откликов фильтра с помощью 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)')