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. По иным вопросам, например если надо исправить заблокированное для перевода слово, обратитесь к редакторам через форму технической поддержки.

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