fsamp2
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 точки.
Примеры
Создайте 2-D конечная импульсная характеристика с помощью частотной выборки
Использование 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.