2-D конечная импульсная характеристика с помощью частотной выборки
h = fsamp2(Hd)Hd, и возвращает коэффициенты фильтра в матрице h. Фильтрация h имеет частотную характеристику, которая проходит через точки в Hd. fsamp2 проектирует двумерные конечные импульсные характеристики фильтры на основе желаемой двумерной частотной характеристики, дискретизированной в точках на Декартовой плоскости.
h = fsamp2(f1,f2,Hd,[m n])m-by- n Конечная импульсная характеристика путем согласования характеристики фильтра в точках в векторах 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, pp. 213-217.