filternorm

2-норма или норма по бесконечности цифрового фильтра

Синтаксис

L = filternorm(b,a)
L = filternorm(b,a,pnorm)
L = filternorm(b,a,2,tol)

Описание

Типичное использование для норм фильтра заключается в цифровом фильтре масштабировании для уменьшения эффектов квантования. Масштабирование часто улучшает отношение сигнал/шум фильтра, не приводя к переполнению данных. Можно также использовать 2-норму, чтобы вычислить энергию импульсной характеристики фильтра.

L = filternorm(b,a) вычисляет 2-норму цифрового фильтра, заданную коэффициентами числителя в b и коэффициенты знаменателя в a.

L = filternorm(b,a,pnorm) вычисляет 2- или норму по бесконечности (inf-норму) цифрового фильтра, где pnorm является либо 2, либо inf.

L = filternorm(b,a,2,tol) вычисляет 2-норму БИХ с заданным допуском, tol. Допуск может быть задан только для БИХ 2-норма расчетов. pnorm в этом случае должно быть 2. Если tol не задан, значение по умолчанию равно 10–8.

Примеры

свернуть все

Вычислите 2-норму БИХ Butterworth с допуском 10-10. Задайте нормированную частоту среза 0.5π рад/с и порядок фильтра 5.

[b,a] = butter(5,0.5);
L2 = filternorm(b,a,2,1e-10)
L2 = 0.7071

Вычислите норму бесконечности КИХ-трансформатора Гильберта порядка 30 и нормированную ширину перехода 0.2π рад/с.

b = firpm(30,[.1 .9],[1 1],'Hilbert');
Linf = filternorm(b,1,inf)
Linf = 1.0028

Алгоритмы

Задан фильтр с H частотной характеристики (e), L p -norm для 1 p < ∞ задается как

H(ejω)p=(12πππ|H(ejω)|pdω)1/p.

Для случая p → ∞, L ∞ -norm является

H(ejω)=maxπωπ|H(ejω)|.

Для случая p = 2 теорема Парсеваля утверждает, что

H(ejω)2=(12πππ|H(ejω)|2dω)1/2=(n|h(n)|2)1/2,

где h (n) - импульсная характеристика фильтра. Энергия импульсной характеристики является квадратной L 2-нормой.

Ссылки

[1] Jackson, L.B. Цифровые фильтры and Signal Processing: with MATLAB Exercises. 3rd Ed. Hingham, MA: Kluwer Academic Publishers, 1996, Chapter 11.

См. также

|

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