2D КИХ-фильтр с помощью выборки частоты
h = fsamp2(Hd)h = fsamp2(f1,f2,Hd,[m n])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] Лим, Джэ С., Двумерная Обработка сигналов и Обработка изображений, Englewood Cliffs, NJ, Prentice Hall, 1990, стр 213-217.