prony

Метод Prony для создания фильтра

Синтаксис

[Num,Den] = prony(impulse_resp,num_ord,denom_ord)

Описание

[Num,Den] = prony(impulse_resp,num_ord,denom_ord) возвращает числитель Num и знаменатель Den коэффициенты для причинной рациональной системной функции с импульсной характеристикой impulse_resp. Системная функция имеет порядок числителя num_ord и знаменатель заказывает denom_ord. Длины Num и Den num_ord+1 и denom_ord+1. Если длина impulse_resp меньше самого большого порядка (num_ord или denom_ord), impulse_resp дополнен нулями. Войдите 0 в num_ord для системной функции все-полюса. Для все-нулевой системной функции введите 0 для denom_ord.

Примеры

свернуть все

Отфильтруйте Ответы с помощью Метода Прони

Скрипт Open Live Script

Подбирайте 4-й порядок модель IIR к импульсной характеристике фильтра lowpass. Постройте исходные и Prony-спроектированные импульсные характеристики.

d = designfilt('lowpassiir','NumeratorOrder',4,'DenominatorOrder',4, ...
    'HalfPowerFrequency',0.2,'DesignMethod','butter');

impulse_resp = filter(d,[1 zeros(1,31)]);
denom_order = 4;
num_order = 4;
[Num,Den] = prony(impulse_resp,num_order,denom_order);

subplot(2,1,1) 
stem(impz(Num,Den,length(impulse_resp)))
title 'Impulse Response with Prony Design'

subplot(2,1,2)
stem(impulse_resp)
title 'Input Impulse Response'

Подбирайте 10-й порядок модель FIR к импульсной характеристике фильтра highpass. Постройте исходные и Prony-спроектированные частотные характеристики.

d = designfilt('highpassfir', 'FilterOrder', 10, 'CutoffFrequency', .8);

impulse_resp = filter(d,[1 zeros(1,31)]);
num_order = 10;
denom_order = 0;
[Num,Den] = prony(impulse_resp,num_order,denom_order);

fvt = fvtool(Num,Den,d);
legend(fvt,'Prony','Original')

Больше о

свернуть все

Системная функция

system function является z-преобразованием импульсной характеристики h [n]:

$H (z) = \sum_{n = - \infty}^{\infty} h [n] z^{- n}$

rational system function является отношением полиномов в z ^–1. Условно полиномом числителя является B (z), и знаменателем является A (z). Следующее уравнение описывает причинную рациональную системную функцию порядка числителя, который q и знаменатель заказывают p:

$H (z) = \frac{\sum_{k = 0}^{q} b [k] z^{- k}}{1 + \sum_{l = 1}^{p} a [l] z^{- l}}$

где a [0] = 1.

Ссылки

Парки, Томас В. и К. Сидни Беррус. Создание цифровых фильтров. Нью-Йорк: John Wiley & Sons, 1987, стр 226–228.

Темы

Параметрическое моделирование

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

Документация Signal Processing Toolbox

Поддержка

Памятка переводчика

1. Если смысл перевода понятен, то лучше оставьте как есть и не придирайтесь к словам, синонимам и тому подобному. О вкусах не спорим.

2. Не дополняйте перевод комментариями “от себя”. В исправлении не должно появляться дополнительных смыслов и комментариев, отсутствующих в оригинале. Такие правки не получится интегрировать в алгоритме автоматического перевода.

3. Сохраняйте структуру оригинального текста - например, не разбивайте одно предложение на два.

4. Не имеет смысла однотипное исправление перевода какого-то термина во всех предложениях. Исправляйте только в одном месте. Когда Вашу правку одобрят, это исправление будет алгоритмически распространено и на другие части документации.

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

Документация