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