2-D КИХ-фильтр с помощью преобразования частоты
h = ftrans2(b,t)h соответствует одномерному FIR-фильтру b с помощью преобразования t. b должен быть одномерным фильтром типа I (четный симметричный фильтр с нечетной длиной), который может быть возвращен fir1, fir2, или firpm на панели инструментов обработки сигналов. Матрица преобразования t содержит коэффициенты, определяющие используемое преобразование частоты.
Использовать ftrans2 создание приблизительно кругового симметричного двухмерного полосного фильтра с полосой пропускания от 0,1 до 0,6 (нормированная частота, где 1,0 соответствует половине частоты дискретизации, или δ радиан). С тех пор ftrans2 преобразует одномерный фильтр FIR для создания двухмерного фильтра, сначала проектируя одномерный полосовой фильтр FIR с помощью функции панели инструментов обработки сигналов 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.
T (start1, start2),
где B (λ) - преобразование Фурье одномерного фильтраb:
− jstartn
и T (start1, start2) - преобразование Фурье матрицы преобразованияt:
jü 1n1e − jstart2n2.
Возвращенный фильтр h - обратное преобразование Фурье для H (start1, start2):
ej
[1] Lim, Jae S., двумерная обработка сигналов и изображений, Englewood Cliffs, NJ, Prentice Hall, 1990, стр. 218-237.