2D КИХ-фильтр с помощью 2D метода окна
h = fwind2(Hd,win)h = fwind2(f1,f2,Hd,win)h = fwind2(Hd,win)h с помощью обратного преобразования Фурье желаемой частотной характеристики Hd и умножение окном win. Hd является матрицей, содержащей желаемую частотную характеристику в равномерно распределенных точках в Декартовой плоскости. fwind2 возвращает h как вычислительную молекулу, которая является соответствующей формой, чтобы использовать с filter2. h одного размера как win.
Используйте fwind2, чтобы разработать двумерные КИХ-фильтры с помощью метода окна. fwind2 использует двумерную спецификацию окна, чтобы разработать двумерный КИХ-фильтр на основе желаемой частотной характеристики Hd. fwind2 работает с двумерными окнами; используйте fwind1, чтобы работать с одномерными окнами.
Для точных результатов используйте точки частоты, возвращенные freqspace, чтобы создать Hd.
Этот пример показывает, как разработать приблизительно циркулярный симметричный двумерный полосовой фильтр с помощью 2D метода окна.
Создайте векторы частотного диапазона 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)
Создайте 2D окно Gaussian с помощью fspecial. Нормируйте окно.
win = fspecial('gaussian',21,2);
win = win ./ max(win(:));Отобразите окно.
mesh(win)
Используя 2D окно, разработайте фильтр, который лучше всего производит желаемую частотную характеристику
h = fwind2(Hd,win);
Отобразите фактическую частотную характеристику этого фильтра.
freqz2(h)
fwind2 вычисляет h с помощью обратного преобразования Фурье и умножения двумерным окном win.
[1] Лим, Джэ С., Двумерная Обработка сигналов и Обработка изображений, Englewood Cliffs, NJ, Prentice Hall, 1990, стр 202-213.