Прототип аналогового lowpass фильтра Бесселя
[z,p,k] = besselap(n)
[z,p,k] = besselap(n) возвращает полюса и усиление порядка-n Прототип аналогового lowpass фильтра Бесселя. n должно быть меньше чем или равно 25. Функция возвращает полюса в длине n вектор-столбец p и усиление в скалярном kZ пустая матрица, потому что нет никаких нулей. Передаточная функция
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.
