Преобразование частоты с нулевым полюсным усилением
[Z2,P2,K2] = zpkftransf(Z,P,K,AllpassNum,AllpassDen)
[Z2,P2,K2] = zpkftransf(Z,P,K,AllpassNum,AllpassDen) возвращает нули, Z2, полюса, P2 и коэффициент усиления, K2, преобразованного цифрового фильтра нижних частот. Прототип фильтра нижних частот имеет нули, Z, полюса, Pи коэффициент усиления, K. Если AllpassDen не указано, по умолчанию будет установлено значение 1. Если ни то ни другое AllpassNum ни AllpassDen , то функция возвращает входной фильтр.
Спроектируйте прототип реального полуполосного фильтра БИХ с использованием стандартного эллиптического подхода:
[b, a] = ellip(3,0.1,30,0.409); [AlpNum, AlpDen] = allpasslp2lp(0.5, 0.25); [z2, p2, k2] = zpkftransf(roots(b),roots(a),b(1),AlpNum,AlpDen);
Проверьте результат, сравнив фильтр прототипа с целевым фильтром:
fvtool(b, a, k2*poly(z2), poly(p2));
После преобразования фильтра вы получите ответ, показанный на рисунке, где полоса пропускания была сдвинута к нулю.
| Переменная | Описание |
|---|---|
Z | Нули прототипа фильтра нижних частот |
P | Полюса прототипа фильтра нижних частот |
K | Коэффициент усиления фильтра нижних частот прототипа |
FTFNum | Числитель фильтра отображения |
FTFDen | Знаменатель фильтра отображения |
Z2 | Нули целевого фильтра |
P2 | Полюса целевого фильтра |
K2 | Коэффициент усиления целевого фильтра |