Обобщенный цифровой проект Фильтра Баттерворта
[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 является нормированной частотой среза, в которой ответ значения фильтра равен (приблизительно-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' и для текстового отображения и для графиков
Разработайте обобщенный Фильтр Баттерворта с нормированной частотой среза 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.
