Allpass фильтруют для комплексного преобразования сдвига
[AllpassNum,AllpassDen] = allpassshiftc(Wo,Wt)
[AllpassNum,AllpassDen] = allpassshiftc(0,0.5)
[AllpassNum,AllpassDen] = allpassshiftc(0,-0.5)
[AllpassNum,AllpassDen] = allpassshiftc(Wo,Wt) возвращает числитель, AllpassNum, и знаменатель, AllpassDen, векторы allpass, сопоставляющего фильтр для выполнения комплексного сдвига частоты частотной характеристики цифрового фильтра произвольной суммой.
[AllpassNum,AllpassDen] = allpassshiftc(0,0.5) вычисляет фильтр allpass для того, чтобы сделать Гильбертово преобразование, 90 степеней против часовой стрелки вращение исходного фильтра в частотном диапазоне.
[AllpassNum,AllpassDen] = allpassshiftc(0,-0.5) вычисляет фильтр allpass для того, чтобы сделать обратное Гильбертово преобразование, т.е. 90 степеней по часовой стрелке вращение исходного фильтра в частотном диапазоне.
Спроектируйте фильтр allpass, точно вращающий целый фильтр по сумме, заданной местоположением выбранной функции от исходного фильтра, Wo=0.5, и его необходимое положение в целевом фильтре, Wt=0.25:
Wo = 0.5; Wt = 0.25; [AllpassNum, AllpassDen] = allpassshiftc(Wo, Wt);
Вычислите частотная характеристика отображения просачиваются полный спектр:
[h, f] = freqz(AllpassNum, AllpassDen, 'whole');| Переменная | Описание |
|---|---|
Wo | Значение частоты, которое будет преобразовано от прототипного фильтра |
Wt | Желаемое местоположение частоты в преобразованном целевом фильтре |
AllpassNum | Числитель фильтра отображения |
AllpassDen | Знаменатель фильтра отображения |
Частоты должны быть нормированы, чтобы быть между-1 и 1 с 1 соответствием половине частоты дискретизации.
Оппенхейм, A.V., Р.В. Шафер и Дж.Р. Бак, обработка сигналов дискретного времени, Prentice-Hall International Inc., 1989.
Датта-Рой, Южная Каролина и Б. Кумар, “На Цифровых Дифференциаторах, Гильбертовых Преобразователях и Полуполосе Фильтры Lowpass”, IEEE® Transactions на Образовании, издании 32, стр 314-318, август 1989.