2D КИХ-фильтр с помощью 2D метода окна
fwind2 функционируйте проектирует 2D КИХ-фильтры с помощью метода окна. fwind2 использует 2D спецификацию окна, чтобы спроектировать 2D КИХ-фильтр на основе желаемой частотной характеристики. fwind2 работает с 2D окнами только. Использование fwind1 создать 2D КИХ-фильтр из 1D окна.
Можно применить 2D КИХ-фильтр к изображениям при помощи filter2 функция.
В этом примере показано, как спроектировать приблизительно циркулярный симметричный двумерный полосовой фильтр с помощью 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.