2-D конечная импульсная характеристика с помощью частотного преобразования
h = ftrans2(b,t)h который соответствует одномерному конечная импульсная характеристика b использование преобразования t. b должен быть одномерным фильтром типа I (четная симметричная нечетная длина), который можно вернуть fir1, fir2, или firpm в программном обеспечении Signal Processing Toolbox. Матрица преобразования t содержит коэффициенты, которые определяют используемое частотное преобразование.
Использование ftrans2 спроектировать приблизительно циркулярно симметричный двумерный полосно-пропускающий фильтр с полосой пропускания от 0,1 до 0,6 (нормированная частота, где 1,0 соответствует половине частоты дискретизации, или, по радианам). Начиная с ftrans2 преобразует одномерную конечную импульсную характеристику фильтр, чтобы создать двумерный фильтр, сначала спроектирует одномерную конечную импульсную характеристику полосно-пропускающий фильтр с помощью функции Signal Processing Toolbox firpm.
colormap(jet(64))
b = firpm(10,[0 0.05 0.15 0.55 0.65 1],[0 0 1 1 0 0]);
[H,w] = freqz(b,1,128,'whole');
plot(w/pi-1,fftshift(abs(H)))
Использование ftrans2 с преобразованием МакКлеллана по умолчанию для создания требуемого приблизительно кругового симметричного фильтра.
h = ftrans2(b); freqz2(h)
Преобразование ниже задает частотную характеристику двумерного фильтра, возвращенного ftrans2.
где B (в) - преобразование Фурье одномерного фильтра b:
и T ('1,' 2) является преобразованием Фурье матрицы преобразования t:
Возвращенный фильтр h - обратное преобразование Фурье H (
[1] Lim, Jae S., Двумерная обработка сигналов и изображений, Englewood Cliffs, NJ, Prentice Hall, 1990, pp. 218-237.