Этот пример показывает, как спроектировать приблизительно циркулярно симметричный двумерный полосно-пропускающий фильтр с помощью 2-D оконного метода.

Создайте векторы частотной области значений f1 и f2 использование freqspace. Эти векторы имеют длину 21.

[f1,f2] = freqspace(21,'meshgrid');

Вычислите расстояние между положениями и центральной частотой.

r = sqrt(f1.^2 + f2.^2);

Создайте матрицу Hd который содержит требуемую полосу пропускания. В этом примере желаемая ширина полосы пропускания находится между 0,1 и 0,5 (нормированная частота, где 1,0 соответствует половине частоты дискретизации, или $$\pi$$ радианы).

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)