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

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

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

Figure contains 2 axes. Axes 1 contains an object of type line. Axes 2 contains an object of type line.

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

[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)')

Figure contains an axes. The axes contains an object of type line.

См. также

|