Преобразуйте БИХ lowpass, фильтр к БИХ объединяют M-ленточный-фильтр
[Num,Den,AllpassNum,AllpassDen] =
iirlp2mbc(B,A,Wo,Wc)
[Num,Den,AllpassNum,AllpassDen] =
iirlp2mbc(B,A,Wo,Wc) возвращает числитель и векторы знаменателя, Num и Den соответственно, целевого фильтра, преобразованного от действительного прототипа lowpass путем применения th-порядка M действительный lowpass, чтобы объединить мультиполосовое преобразование частоты.
Это также возвращает числитель, AllpassNum, и знаменатель, AllpassDen, allpass, сопоставляющего фильтр. Прототип lowpass фильтр дан с числителем, заданным B и знаменателем, заданным A.
Это преобразование эффективно помещает одну функцию исходного фильтра, расположенного на частоте Во, на необходимых целевых местоположениях частоты, Wt1..., WtM.
Относительные положения других функций исходного фильтра не изменяются в целевом фильтре. Это означает, что возможно выбрать две функции исходного фильтра, F1 и F2, с F1, предшествующим F2. F1 функции будет все еще предшествовать F2 после преобразования. Однако расстояние между F1 и F2 не будет тем же самым до и после преобразования.
Выбор функции, подвергающейся этому преобразованию, не ограничивается частотой среза исходного фильтра lowpass. В целом возможно выбрать любую функцию; например, ребро полосы задерживания, DC, глубокий минимум в полосе задерживания или другие единицы.
Это преобразование может также использоваться для преобразования других типов фильтров; например, фильтры метки или резонаторы могут быть легко реплицированы во многие необходимые местоположения частоты. Хорошее приложение было бы адаптивной тональной схемой отмены, реагирующей на изменяющийся номер и местоположение тонов.
Разработайте прототип действительный БИХ lowpass эллиптический фильтр с усилением приблизительно-3 дБ в 0.5π рад/выборка.
[b,a] = ellip(3,0.1,30,0.409);
Создайте комплексный многополосный фильтр с двумя полосами пропускания.
[num,den] = iirlp2mbc(b,a,0.5,[-7 -5 6 8]/10);
Сравните ответы значения фильтров с помощью FVTool. iirlp2mbc реплицирует желаемую функцию в 0,5 в фильтре lowpass в четырех местах в многополосном фильтре.
hvft = fvtool(b,a,num,den); legend(hvft,'Prototype','Target')
| Переменная | Описание |
|---|---|
B | Числитель прототипа lowpass фильтр. |
A | Знаменатель прототипа lowpass фильтр. |
Wo | Значение частоты, которое будет преобразовано от прототипного фильтра. Это должно быть нормировано, чтобы быть между 0 и 1 с 1 соответствием половине частоты дискретизации. |
Wc | Желаемые местоположения частоты в преобразованном целевом фильтре. Они должны быть нормированы, чтобы быть между-1 и 1 с 1 соответствием половине частоты дискретизации. |
Num | Числитель целевого фильтра. |
Den | Знаменатель целевого фильтра. |
AllpassNum | Числитель фильтра отображения. |
AllpassDen | Знаменатель фильтра отображения. |
