2-D КИХ-фильтр с использованием частотной выборки
h = fsamp2(Hd)Hdи возвращает коэффициенты фильтра в матрице h. Фильтр h имеет частотную характеристику, которая проходит через точки в Hd. fsamp2 конструирует двумерные КИХ-фильтры на основе желаемой двумерной частотной характеристики, дискретизированной в точках на декартовой плоскости.
h = fsamp2(f1,f2,Hd,[m n])mоколо-n Фильтр FIR путем сопоставления отклика фильтра в точках векторов f1 и f2. Частотные векторы f1 и f2 имеют нормированную частоту, где 1,0 соответствует половине частоты дискретизации, или δ радиан. Результирующий фильтр как можно ближе подходит к требуемому отклику в смысле наименьших квадратов. Для достижения наилучших результатов должно быть не менее m*n требуемые частотные точки. fsamp2 выдает предупреждение, если указано менее m*n точки.
Использовать fsamp2 для конструирования приблизительно симметричного двумерного полосного фильтра с полосой пропускания от 0,1 до 0,5 (нормированная частота, где 1,0 соответствует половине частоты дискретизации, или радиан).
Создание матрицы Hd который содержит требуемый полосовой отклик. Использовать freqspace для создания частотных векторов f1 и f2.
[f1,f2] = freqspace(21,'meshgrid');
Hd = ones(21);
r = sqrt(f1.^2 + f2.^2);
Hd((r<0.1)|(r>0.5)) = 0;
colormap(jet(64))
mesh(f1,f2,Hd)
Создайте фильтр, проходящий через этот ответ.
h = fsamp2(Hd); freqz2(h)
fsamp2 вычисляет фильтр h принимая обратное дискретное преобразование Фурье желаемой частотной характеристики. Если желаемая частотная характеристика является действительной и симметричной (нулевая фаза), результирующий фильтр также является нулевой фазой.
[1] Lim, Jae S., двумерная обработка сигналов и изображений, Englewood Cliffs, NJ, Prentice Hall, 1990, стр. 213-217.