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