Комплекс нулей, полюсов и усиления полосовое преобразование частоты
[Z2,P2,K2,AllpassNum,AllpassDen]
= zpkrateup(Z,P,K,N)
[Z2,P2,K2,AllpassNum,AllpassDen]
= zpkrateup(Z,P,K,N) возвращает нули, Z2, полюса, P2, и фактор усиления, K2, целевого фильтра, преобразовываемого от любого прототипа путем применения Nth-порядок rateup преобразование частоты, где N сверхдискретизировать отношение. Преобразование создает N равные копии прототипа фильтруют частотную характеристику.
Это также возвращает числитель, AllpassNum, и знаменатель, AllpassDen, из allpass, сопоставляющего фильтр. Исходный фильтр lowpass дан с нулями, Z, полюса, P, и фактор усиления, K.
Относительные положения других функций исходного фильтра не изменяются в целевом фильтре. Это означает, что возможно выбрать две функции исходного фильтра, F1 и F2, с F1, предшествующим F2. F1 функции будет все еще предшествовать F2 после преобразования. Однако расстояние между F1 и F2 не будет тем же самым до и после преобразования.
Спроектируйте прототип действительный фильтр полуполосы IIR с помощью стандартного эллиптического подхода:
[b, a] = ellip(3,0.1,30,0.409); z = roots(b); p = roots(a); k = b(1); % Upsample the prototype filter 4 times [z2,p2,k2] = zpkrateup(z, p, k, 4); % Compare prototype filter with target filter fvtool(b, a, k2*poly(z2), poly(p2));
| Переменная | Описание |
|---|---|
Z | Нули прототипа фильтр lowpass |
P | Полюса прототипа фильтр lowpass |
K | Фактор усиления прототипа фильтр lowpass |
N | Целочисленное отношение повышающей дискретизации |
Z2 | Нули целевого фильтра |
P2 | Полюса целевого фильтра |
K2 | Фактор усиления целевого фильтра |
AllpassNum | Числитель фильтра отображения |
AllpassDen | Знаменатель фильтра отображения |
Частоты должны быть нормированы, чтобы быть между-1 и 1 с 1 соответствием половине частоты дискретизации.