2-D конечная импульсная характеристика с 1-D оконного метода
fwind1 функция проектирует 2-D конечная импульсная характеристика с помощью оконного метода. fwind1 использует спецификацию 1-D окна для разработки 2-D конечная импульсная характеристика на основе желаемой частотной характеристики. fwind1 работает только с 1-D окнами. Использовать fwind2 для работы с 2-D окнами.
Можно применить 2-D конечная импульсная характеристика к изображениям при помощи filter2 функция.
Этот пример показывает, как спроектировать приблизительно циркулярно симметричный двумерный полосно-пропускающий фильтр с помощью 1-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)
Спроектируйте 1-D окно. Этот пример использует окно Хэмминга длины 21.
win = 0.54 - 0.46*cos(2*pi*(0:20)/20);
Постройте график 1-D окна.
figure plot(linspace(-1,1,21),win);
Используя 1-D окно, спроектируйте фильтр, который лучше всего производит эту частотную характеристику
h = fwind1(Hd,win);
Отображение фактической частотной характеристики этого фильтра.
freqz2(h)
fwind1 функция принимает одномерную спецификацию окна и формирует приблизительно круговое симметричное двумерное окно с помощью метода Хуана,
где w (t) - одномерное окно, а w (n1, n2) - полученное двумерное окно.
Учитывая два окна, fwind1 функция образует разделяемое двумерное окно:
fwind1 функция вызывает fwind2 с желаемой частотной характеристикой Hd и двумерное окно. fwind2 функция вычисляет h использование обратного преобразования Фурье и умножения на двумерное окно:
[1] Lim, Jae S., Двумерная обработка сигналов и изображений, Englewood Cliffs, NJ, Prentice Hall, 1990.