firlpnorm

Наименьшее количество P-нормы оптимальный КИХ-фильтр

Синтаксис

b = firlpnorm(n,f,edges,a)
b = firlpnorm(n,f,edges,a,w)
b = firlpnorm(n,f,edges,a,w,p)
b = firlpnorm(n,f,edges,a,w,p,dens)
b = firlpnorm(n,f,edges,a,w,p,dens,initnum)
b = firlpnorm(...,'minphase')
[b,err] = firlpnorm(...)

Описание

b = firlpnorm(n,f,edges,a) возвращается фильтр числителя заказывают n, который представляет наилучшее приближение частотной характеристике, описанной f и a в наименьшем-количестве-Pth смысла нормы. P установлен в 128 по умолчанию, который чрезвычайно эквивалентный норме бесконечности. Векторный edges задает частоты ребра полосы для многополосных проектов. firlpnorm использует неограниченный алгоритм квазиньютона, чтобы разработать заданный фильтр.

f и a должны иметь то же число элементов, которое может превысить число элементов в edges. Это позволяет вам задать фильтры с любым контуром усиления в каждой полосе. Однако частоты в edges должны также быть в векторном f. Всегда используйте freqz, чтобы проверять получившийся фильтр.

Примечание

firlpnorm использует нелинейную стандартную программу оптимизации, которая не может сходиться в некоторых случаях проекта фильтра. Кроме того, алгоритм не является подходящим для определенного крупного заказа (порядок> 100) проекты фильтра.

b = firlpnorm(n,f,edges,a,w) использует веса в w, чтобы взвесить ошибку. w имеет одну запись на точку частоты (та же длина как f и a), который говорит firlpnorm, сколько акцента, чтобы поставить минимизацию ошибки около каждой частоты указывают относительно других точек. Например,

b = firlpnorm(20,[0 .15 .4 .5 1],[0 .4 .5 1],... 
[1 1.6 1 0 0],[1 1 1 10 10])

разрабатывает фильтр lowpass с пиком 1,6 в полосе пропускания, и с акцентом, сделанным минимизации ошибки в полосе задерживания.

b = firlpnorm(n,f,edges,a,w,p) то, где p является двухэлементным вектором [pmin pmax], позволяет вам задать минимальные и максимальные значения p, используемого в наименьшем-количестве-pth алгоритма. Значением по умолчанию является [2 128], который по существу приводит к L-бесконечности, или Чебышеву, норме. pmin и pmax должны быть четными числами. Алгоритм проекта начинает оптимизировать фильтр с pmin и перемещается к оптимальному фильтру в смысл pmax. Когда p установлен в 'inspect', firlpnorm не оптимизирует получившийся фильтр. Вы можете использовать эту функцию, чтобы осмотреть начальное нулевое размещение.

b = firlpnorm(n,f,edges,a,w,p,dens) задает плотность сетки dens, используемый в оптимизации. Количество узлов решетки [dens*(n+1)]. Значение по умолчанию равняется 20. Можно задать dens как одноэлементный массив ячеек. Сетка равномерно распределена.

b = firlpnorm(n,f,edges,a,w,p,dens,initnum) позволяет вам определить первоначальную оценку коэффициентов числителя фильтра в векторном initnum. Это может оказаться полезным для трудных задач оптимизации. Нулевой полюсом редактор в программном обеспечении Signal Processing Toolbox™ может использоваться для генерации initnum.

b = firlpnorm(...,'minphase') то, где 'minphase' является последним аргументом в списке аргументов, генерирует КИХ-фильтр минимальной фазы. По умолчанию, фильтры смешанной фазы проекта firlpnorm. Определение входной опции 'minphase' заставляет firlpnorm использовать различный метод оптимизации, чтобы разработать фильтр минимальной фазы. В результате различной используемой оптимизации фильтр минимальной фазы может привести к немного отличающимся результатам.

[b,err] = firlpnorm(...) возвращает наименьшее-количество-pth ошибки приближения err.

Примеры

свернуть все

Lowpass фильтруют с пиком 1,4 в полосе пропускания.

b = firlpnorm(22,[0 .15 .4 .5 1],[0 .4 .5 1],[1 1.4 1 0 0],...
[1 1 1 2 2]);
fvtool(b)

Получившийся фильтр является lowpass, с желаемыми 1,4 пиками в полосе пропускания (заметьте эти 1.4, заданные в векторе a).

Фильтр минимальной фазы Highpass оптимизирован для с 4 нормами.

b = firlpnorm(44,[0 .4 .45 1],[0 .4 .45 1],[0 0 1 1],[5 1 1 1],...
[2 4],'minphase');
fvtool(b)

Это - минимальная фаза, highpass фильтр.

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

fvtool(b,'polezero')

Ссылки

Saramaki, T, конечный импульсный проект фильтра ответа, руководство для цифрового сигнала ProcessingMitra, S.K. и J.F. Wiley-межнаука редакторов кайзера, Нью-Йорк, 1993, глава 4.

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

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

| | | | | | |

Введенный в R2011a

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