Аналог функции Бесселя lowpass фильтрует прототип
[z,p,k] = besselap(n)
[z,p,k] = besselap(n) возвращает полюса и усиление аналога функции Бесселя порядка-n lowpass прототип фильтра. n должен быть меньше чем или равен 25. Функция возвращает полюса в длине вектор-столбец n p и усиление в скалярном k. z является пустой матрицей, потому что нет никаких нулей. Передаточная функция
besselap нормирует полюса и усиление так, чтобы в низкой частоте и высокой частоте прототип функции Бесселя был асимптотически эквивалентен прототипу Баттерворта того же порядка [1]. Значение фильтра является меньше, чем в частоте среза единицы Ωc = 1.
Аналоговые фильтры Бесселя характеризуются групповой задержкой, которая является максимально плоской на нулевой частоте и почти постоянной в полосе пропускания. Групповая задержка на нулевой частоте
Разработайте аналог функции Бесселя 6-го порядка lowpass фильтр. Отобразите его значение и фазовые отклики.
[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. Золото. Теория и Приложение Цифровой обработки сигналов. Englewood Cliffs, NJ: Prentice Hall, 1975, стр 228–230.
