Преобразуйте БИХ-полосовой фильтр комплекса к БИХ-полосовому фильтру комплекса с различными характеристиками
[Num,Den,AllpassNum,AllpassDen] =
iirbpc2bpc(B,A,Wo,Wt)
[Num,Den,AllpassNum,AllpassDen] =
iirbpc2bpc(B,A,Wo,Wt) возвращает числитель и векторы знаменателя, Num и Den соответственно, целевого фильтра, преобразованного от комплексного полосового прототипа путем применения комплексной полосы пропускания первого порядка, чтобы объединить полосовое преобразование частоты.
Это также возвращает числитель, AllpassNum, и знаменатель, AllpassDen, allpass, сопоставляющего фильтр. Прототип lowpass фильтр дан с числителем, заданным B и знаменателем, заданным A.
Это преобразование эффективно помещает две функции исходного фильтра, расположенного на частотах Wo1 и Wo2, на необходимых целевых местоположениях частоты, Wt1 и Wt2 соответственно. Это принято, что Wt2 больше, чем Wt1. В большинстве случаев функциями, выбранными для преобразования, являются ребра полосы полос пропускания фильтра. В целом возможно выбрать любую функцию; например, ребро полосы задерживания, DC, глубокий минимум в полосе задерживания или другие единицы.
Относительные положения других функций исходного фильтра не изменяются в целевом фильтре. Это означает, что возможно выбрать две функции исходного фильтра, F1 и F2, с F1, предшествующим F2. F1 функции будет все еще предшествовать F2 после преобразования. Однако расстояние между F1 и F2 не будет тем же самым до и после преобразования.
Это преобразование может также использоваться для преобразования других типов фильтров; например, комплексные фильтры метки или резонаторы могут быть изменены местоположение на двух отличных желаемых частотах в любом месте вокруг модульного круга; например, в адаптивной системе.
Разработайте прототип действительный БИХ lowpass эллиптический фильтр с усилением приблизительно-3 дБ в 0.5π рад/выборка.
[b,a] = ellip(3,0.1,30,0.409);
Создайте комплексную полосу пропускания от 0.25π до 0.75π.
[bc,ac] = iirlp2bpc(b,a,0.5,[0.25 0.75]);
Переместите полосу пропускания в между –0.3π и 0.1π.
[num,den] = iirbpc2bpc(bc,ac,[0.25 0.75],[-0.3 0.1]);
Сравните эти три просачивается FVTool.
hvft = fvtool(b,a,bc,ac,num,den); legend(hvft,'Prototype','Positive complex band','Target')
| Переменная | Описание |
|---|---|
B | Числитель прототипа lowpass фильтр |
A | Знаменатель прототипа lowpass фильтр |
Wo | Значения частоты, которые будут преобразованы от прототипного фильтра |
Wt | Желаемые местоположения частоты в преобразованном целевом фильтре |
Num | Числитель целевого фильтра |
Den | Знаменатель целевого фильтра |
AllpassNum | Числитель фильтра отображения |
AllpassDen | Знаменатель фильтра отображения |
Частоты должны быть нормированы, чтобы быть между-1 и 1 с 1 соответствием половине частоты дискретизации.
