Прототип аналогового фильтра нижних частот Бесселя
[z,p,k] = besselap(n)
[z,p,k] = besselap(n) возвращает полюса и коэффициент усиления порядка -n Прототип аналогового фильтра нижних частот Бесселя. n должно быть меньше или равно 25. Функция возвращает полюса в длину n вектор столбца p и коэффициент усиления в скаляре k. z является пустой матрицей, поскольку нет нулей. Передаточная функция:
)) ⋯ (s − p (n))
besselap нормализует полюса и усиление таким образом, что при низкой частоте и высокой частоте прототип Бесселя асимптотически эквивалентен прототипу Баттерворта того же порядка [1]. Величина фильтра составляет менее при единичной частоте отсечки Startc = 1.
Аналоговые фильтры Бесселя характеризуются групповой задержкой, которая является максимально плоской при нулевой частоте и почти постоянной по всей полосе пропускания. Групповая задержка на нулевой частоте равна
1/n
Проектирование аналогового фильтра нижних частот Бесселя 6-го порядка. Отображение его величины и фазовых откликов.
[z,p,k] = besselap(6); % Lowpass filter prototype [num,den] = zp2tf(z,p,k); % Convert to transfer function form freqs(num,den) % Frequency response of analog filter
besselap находит корни фильтра из таблицы подстановки, созданной с помощью программы Symbolic Math Toolbox™.
[1] Rabiner, L. R. и B. Золото. Теория и применение цифровой обработки сигналов. Энглвуд Клиффс, Нью-Джерси: Прентис-Холл, 1975, стр. 228-230.