Пример расчета частотного отклика системы
Определить частотную характеристику системы с импульсным откликом:
h(0,0) = 0.25, h(0, 1) = 0.125, h( 1,0) = 0.125, h( 1, 1) = 0.0625.
Частотный отклик:
H(wx,wy) = h(n,m)exp(-jnwx-jmwy) = 0.25+0.125[exp(-jwx)+exp(jwx)+exp(-wy)+exp(jwy)]+
+0.0625[exp(-jwx-jwy)+exp(-jwx+jwy)+exp(jwx-jwy)+exp(jwx+jwy)] = 0.25(1+cos wx)(1+cos wy).
Система является примером двумерного фильтра нижних частот. Частотный отклик системы на плоскости (wx,wy), приведенный на рис. 18.3.1, имеет осевую симметрию с коэффициентом передачи 1 в центре (wx=0, wy=0) со спадом до нуля при wx= p и wy= p.
Рис. 18.3.1. Частотная характеристика ФНЧ.
При разделимости импульсного отклика частотный отклик многомерных систем также является разделимой функцией:
h(k,l)= q(k)g(l) Û Q(wx)G(wy)= H(wx,wy)
Q(wx) = Sk q(k) exp(-jk×wx). G(wy) = Sl g(l) exp(-jl×wy).
Импульсный отклик системы. Выражение (18.3.1), по существу, описывает разложение функции Н(wx,wy) в двумерный рад Фурье с коэффициентами разложения в виде отсчетов импульсного отклика h(k,l), т.е. прямое преобразование Фурье. Очевидно, что обратным преобразованием Фурье с интегрированием в пределах одного периода из частотного отклика H(wx,wy) можно получить импульсный отклик системы:
h(k,l) = H(wx,wy) exp(jkwx+jlwy) dwxdwy. (18.3.2)