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