Allpass фильтруют для lowpass к bandstop преобразованию
[AllpassNum,AllpassDen] = allpasslp2bs(Wo,Wt)
[AllpassNum,AllpassDen] = allpasslp2bs(Wo,Wt)
возвращает числитель, AllpassNum
, и знаменатель, AllpassDen
, allpass второго порядка, сопоставляющего фильтр для выполнения действительного lowpass к действительному bandstop преобразованию частоты. Это преобразование эффективно помещает одну функцию исходного фильтра, расположенного на частоте-Wo
, на необходимом целевом местоположении частоты, Wt1 и второй функции, первоначально в +W
o, в новом местоположении, W
t2. Это принято, что W
t2 больше, чем W
t1. Это преобразование реализует "Мобильность Найквиста", что означает, что функция DC остается в DC, но функция Найквиста перемещается в местоположение, зависящее от выбора W
o и W
t.
Относительные положения других функций исходного фильтра изменяются в целевом фильтре. Это означает, что возможно выбрать две функции исходного фильтра, F
1 и F
2, с F
1 предыдущий F
2. После функции преобразования F
2 будет предшествовать F
1 в целевом фильтре. Однако расстояние между F
1 и F
2 не будет тем же самым до и после преобразования.
Выбор функции, подвергающейся lowpass к bandstop преобразованию, не ограничивается только частотой среза исходного фильтра lowpass. В целом возможно выбрать любую функцию; например, ребро полосы задерживания, DC, глубокий минимум в полосе задерживания или другие единицы.
Разработайте фильтр allpass, изменяющий фильтр lowpass с частотой среза в W
o=0.5
к действительному заграждающему фильтру с частотами среза в W
t1=0.25
и W
t2=0.375
:
Wo = 0.5; Wt = [0.25, 0.375]; [AllpassNum, AllpassDen] = allpasslp2bs(Wo, Wt); [h, f] = freqz(AllpassNum, AllpassDen, 'whole'); plot(f/pi, abs(angle(h))/pi, Wt, Wo, 'ro'); title('Mapping Function Wo(Wt)'); xlabel('New Frequency, Wt'); ylabel('Old Frequency, Wo');
В фигуре вы находите функцию filter отображения, как определено примером. Обратите внимание, что ответ нормирован к π:
Переменная | Описание |
---|---|
Wo | Значение частоты, которое будет преобразовано от прототипного фильтра |
Wt | Желаемые местоположения частоты в преобразованном целевом фильтре |
AllpassNum | Числитель фильтра отображения |
AllpassDen | Знаменатель фильтра отображения |
Частоты должны быть нормированы, чтобы быть между 0 и 1 с 1 соответствием половине частоты дискретизации.
Constantinides, A.G., “Спектральные преобразования для цифровых фильтров”, IEEE® Proceedings, издание 117, № 8, стр 1585-1590, август 1970.
Nowrouzian, B. и А.Г. Констэнтинайдс, “Параметры передаточной функции ссылки прототипа в преобразованиях частоты дискретного времени”, Продолжения 33-й Среднезападный Симпозиум по Схемам и Системам, Калгари, Канада, изданию 2, стр 1078-1082, август 1990.
Nowrouzian, B. и Л.Т. Брутон, “Решения Закрытой формы в течение дискретного времени эллиптические передаточные функции”, Продолжения 35-го Среднезападного Симпозиума по Схемам и Системам, изданию 2, стр 784-787, 1992.
Constantinides, A.G., “Проект полосовых цифровых фильтров”, Продолжения IEEE, издание 1, стр 1129-1231, июнь 1969.