zpkbpc2bpc
Комплекс нулей, полюсов и усиления полосовое преобразование частоты
Синтаксис
[Z2,P2,K2,AllpassNum,AllpassDen]
= zpkbpc2bpc(Z,P,K,Wo,Wt)
Описание
[Z2,P2,K2,AllpassNum,AllpassDen]
= zpkbpc2bpc(Z,P,K,Wo,Wt) возвращает нули, Z2, полюса, P2, и фактор усиления, K2, целевого фильтра, преобразованного от комплексного полосового прототипа путем применения комплексной полосы пропускания первого порядка, чтобы объединить полосовое преобразование частоты.
Это также возвращает числитель, AllpassNum, и знаменатель, AllpassDen, из allpass, сопоставляющего фильтр. Исходный фильтр lowpass дан с нулями, Z, полюса, P, и фактор усиления, K.
Это преобразование эффективно помещает две функции исходного фильтра, расположенного на частотах Wo1 и Wo2, на необходимых целевых местоположениях частоты, Wt1 и Wt2 соответственно. Это принято, что Wt2 больше Wt1. В большинстве случаев функциями, выбранными для преобразования, являются ребра полосы полос пропускания фильтра. В целом возможно выбрать любую функцию; например, ребро полосы задерживания, DC, глубокий минимум в полосе задерживания или другие единицы.
Относительные положения других функций исходного фильтра не изменяются в целевом фильтре. Это означает, что возможно выбрать две функции исходного фильтра, F1 и F2, с F1, предшествующим F2. F1 функции будет все еще предшествовать F2 после преобразования. Однако расстояние между F1 и F2 не будет тем же самым до и после преобразования.
Это преобразование может также использоваться в преобразовании других типов фильтров; например, комплексные фильтры метки или резонаторы могут быть изменены местоположение на двух отличных желаемых частотах в любом месте вокруг модульного круга; например, в адаптивной системе.
Примеры
Спроектируйте прототип действительный БИХ-полуленточный фильтр с помощью стандартного эллиптического подхода:
[b, a] = ellip(3,0.1,30,0.409);
Создайте комплексную полосу пропускания от 0,25 до 0,75:
[b, a] = iirlp2bpc(b,a,0.5,[0.25,0.75]); z = roots(b); p = roots(a); k = b(1); [z2,p2,k2] = zpkbpc2bpc(z,p,k,[0.25, 0.75],[-0.75, -0.25]);
Проверьте результат путем сравнения прототипного фильтра с целевым фильтром:
fvtool(b, a, k2*poly(z2), poly(p2));
Сравнение фильтров в FVTool показывает результаты в качестве примера. Используйте функции в FVTool, чтобы проверять коэффициенты фильтра или другие исследования фильтра.
Аргументы
| Переменная | Описание |
|---|---|
Z | Нули прототипа фильтр lowpass |
P | Полюса прототипа фильтр lowpass |
K | Фактор усиления прототипа фильтр lowpass |
Wo | Значение частоты, которое будет преобразовано от прототипного фильтра |
Wt | Желаемое местоположение частоты в преобразованном целевом фильтре |
Z2 | Нули целевого фильтра |
P2 | Полюса целевого фильтра |
K2 | Фактор усиления целевого фильтра |
AllpassNum | Числитель фильтра отображения |
AllpassDen | Знаменатель фильтра отображения |
Частоты должны быть нормированы, чтобы быть между-1 и 1 с 1 соответствием половине частоты дискретизации.