maxflat

Обобщенный цифровой проект Фильтра Баттерворта

Синтаксис

[b,a] = maxflat(n,m,Wn)
b = maxflat(n,'sym',Wn)
[b,a,b1,b2] = maxflat(n,m,Wn)
[b,a,b1,b2,sos,g] = maxflat(n,m,Wn)
[...] = maxflat(n,m,Wn,'design_flag')

Описание

[b,a] = maxflat(n,m,Wn) lowpass Фильтр Баттерворта с числителем и коэффициентами знаменателя b и a порядков n и m, соответственно. Wn является нормированной частотой среза, в которой ответ значения фильтра равен 1/2  (приблизительно-3 дБ). Wn должен быть между 0 и 1, где 1 соответствует частоте Найквиста.

b = maxflat(n,'sym',Wn) симметричный КИХ-Фильтр Баттерворта. n должен быть даже, и Wn ограничивается подынтервалом [0,1]. Функция повышает ошибку, если Wn задан за пределами этого подынтервала.

[b,a,b1,b2] = maxflat(n,m,Wn) возвращает два полинома b1 и b2, продукт которого равен полиному числителя b (то есть,   b = conv(b1,b2)). b1 содержит все нули в z = -1, и b2 содержит все другие нули.

[b,a,b1,b2,sos,g] = maxflat(n,m,Wn) возвращает представление разделов второго порядка фильтра как матрица фильтра sos и усиление g.

[...] = maxflat(n,m,Wn,'design_flag') позволяет вам контролировать проект фильтра, где 'design_flag'

  • 'trace' для текстового отображения таблицы проекта используется в проекте

  • 'plots' для графиков значения фильтра, групповой задержки, и нулей и полюсов

  • 'both' и для текстового отображения и для графиков

Примеры

свернуть все

Разработайте обобщенный Фильтр Баттерворта с нормированной частотой среза 0.2π rad/s. Задайте порядок числителя 10 и порядок знаменателя 2. Визуализируйте частотную характеристику фильтра.

n = 10;
m = 2;
Wn = 0.2;

[b,a] = maxflat(n,m,Wn);
fvtool(b,a)

Алгоритмы

Метод состоит из использования формул, полиномиального корневого открытия и преобразования полиномиальных корней.

Ссылки

[1] Selesnick, Иван В. и К. Сидни Беррус. “Обобщенный Цифровой Проект Фильтра Баттерворта”. IEEE® Transactions на Обработке сигналов. Издание 46, Номер 6, 1998, стр 1688–1694.

Расширенные возможности

Смотрите также

| |

Представлено до R2006a

Для просмотра документации необходимо авторизоваться на сайте