2D КИХ-фильтр с помощью 1D метода окна
h = fwind1(Hd,win)h на основе желаемой частотной характеристики Hd. fwind1 возвращает h как вычислительная молекула, которая является соответствующей формой, чтобы использовать с filter2. fwind1 использует одномерное окно win сформировать приблизительно циркулярное симметричное двумерное окно с помощью метода Хуана.
fwind1 работает с 1D окнами только; используйте fwind2 работать с двумерными окнами.
В этом примере показано, как спроектировать приблизительно циркулярный симметричный двумерный полосовой фильтр с помощью 1D метода окна.
Создайте векторы частотного диапазона f1 и f2 использование freqspace. Эти векторы имеют длину 21.
[f1,f2] = freqspace(21,'meshgrid');Вычислите расстояние каждого положения от центральной частоты.
r = sqrt(f1.^2 + f2.^2);
Создайте матричный Hd это содержит желаемый полосовой ответ. В этом примере желаемая полоса пропускания между 0,1 и 0.5 (нормированная частота, где 1.0 соответствует половине частоты дискретизации, или радианы).
Hd = ones(21); Hd((r<0.1)|(r>0.5)) = 0;
Отобразите идеальный полосовой ответ.
colormap(parula(64)) mesh(f1,f2,Hd)
Спроектируйте 1D окно. Этот пример использует Окно Хэмминга длины 21.
win = 0.54 - 0.46*cos(2*pi*(0:20)/20);
Постройте 1D окно.
figure plot(linspace(-1,1,21),win);
Используя 1D окно, спроектируйте фильтр, который лучше всего производит эту частотную характеристику
h = fwind1(Hd,win);
Отобразите фактическую частотную характеристику этого фильтра.
freqz2(h)
fwind1 берет одномерную спецификацию окна и формирует приблизительно циркулярное симметричное двумерное окно с помощью метода Хуана,
где w (t) является одномерным окном, и w (n1, n2) является получившимся двумерным окном.
Учитывая два окна, fwind1 формирует отделимое двумерное окно:
fwind1 вызовы fwind2 с Hd и двумерное окно. fwind2 вычисляет h использование обратного преобразования Фурье и умножения двумерным окном:
[1] Лим, Джэ С., двумерная обработка сигналов и обработка изображений, Englewood Cliffs, NJ, Prentice Hall, 1990.