Преобразуйте БИХ фильтр lowpass к БИХ-M-ленточному-фильтру
[Num,Den,AllpassNum,AllpassDen] =
iirlp2mb(B,A,Wo,Wt)
[Num,Den,AllpassNum,AllpassDen]=iirlp2mb(B,A,Wo,Wt,Pass)
[Num,Den,AllpassNum,AllpassDen] =
iirlp2mb(B,A,Wo,Wt) возвращает числитель и векторы знаменателя, Num и Den соответственно, целевого фильтра, преобразованного от действительного lowpass, моделируют путем применения Mth-порядок действительный lowpass к действительному несколько полосовое отображение частоты. По умолчанию функция DC сохранена в ее исходном местоположении.
[Num,Den,AllpassNum,AllpassDen]=iirlp2mb(B,A,Wo,Wt,Pass) позволяет вам задавать дополнительный параметр, Pass, который выбирает между использованием “Мобильности DC” и “Мобильностью Найквиста”. В первом случае функция Найквиста остается в своем исходном местоположении, и функция DC свободна перемещаться. Во втором случае функция DC сохранена на исходной частоте, и функция Найквиста подвижна.
Это также возвращает числитель, 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);
Создайте действительный многополосный фильтр с двумя полосами пропускания.
[num1,den1] = iirlp2mb(b,a,0.5,[2 4 6 8]/10);
Создайте действительный многополосный фильтр с двумя полосами задерживания.
[num2,den2] = iirlp2mb(b,a,0.5,[2 4 6 8]/10, 'stop');Сравните ответы величины фильтров с помощью FVTool.
hvft = fvtool(b,a,num1,den1,num2,den2); legend(hvft,'Prototype','Two passbands','Two stopbands')
| Переменная | Описание |
|---|---|
B | Числитель прототипа фильтр lowpass |
A | Знаменатель прототипа фильтр lowpass |
Wo | Значение частоты, которое будет преобразовано от прототипного фильтра |
Wt | Желаемые местоположения частоты в преобразованном целевом фильтре |
Pass | Выбор ( |
Num | Числитель целевого фильтра |
Den | Знаменатель целевого фильтра |
AllpassNum | Числитель фильтра отображения |
AllpassDen | Знаменатель фильтра отображения |
Частоты должны быть нормированы, чтобы быть между 0 и 1 с 1 соответствием половине частоты дискретизации.
[1] Франкитти, J.C., “Все-передайте интерполяция фильтра и проблемы преобразования частоты”, Тезис магистра наук, Отдел Электротехники и Вычислительной техники, Университета Колорадо, 1985.
[2] Feyh, G., Й.К. Франкитти и К.Т. Муллис, “Все-передают интерполяция фильтра и проблема преобразования частоты”, Продолжения 20-я Конференция Asilomar по Сигналам, Системам и Компьютерам, Пасифик-Гроуву, Калифорния, стр 164-168, ноябрь 1986.
[3] Муллис, К.Т. и Р. А. Робертс, Цифровая обработка сигналов, разделяет 6.7, Чтение, Масса., Аддисон-Уэсли, 1987.
[4] Feyh, G., В.Б. Джонс и К.Т. Муллис, “Расширение Алгоритма Шура для преобразований частоты”, Линейные схемы, Системы и Обработка сигналов: Теория и Приложение, К. Дж. Бирнс и др. Редакторы, Амстердам: Elsevier, 1988.