- •Сигналы и линейные системы
- •Тема 13. Многомерные сигналы и системы.
- •Содержание
- •Введение.
- •13.1. Двумерные и многомерные сигналы [9].
- •13.2. Двумерные системы.
- •13.3. Частотные характеристики сигналов и систем.
- •13.4. Дискретизация двумерных сигналов [9].
- •13.5. Многомерный спектральный анализ [9].
- •Литература
13.3. Частотные характеристики сигналов и систем.
Частотный отклик системы. Допустим, что двумерная ЛИС-система имеет импульсный отклик h(kx,ly). Подадим на вход системы сигнал вида комплексной синусоиды:
s(n,m) = exp(jnxx+jmyy),
где x и y – значения частоты сигнала соответственно по координатам x и y. Принимая x = 1, y = 1 и выполняя двумерную свертку (13.2.5), получаем:
z(n,m) = h(k,l) exp[jx(n-k)+jy(m-l)] =
= exp(jnx+jmy) h(k,l) exp(-jkx-jly) = H(x,y) exp(jnx+jmy).
H(x,y) = h(k,l) exp(-jkx-jly). (13.3.1)
Таким образом, выходной сигнал представляет собой комплексную синусоиду с теми же значениями частоты, что и у входного сигнала, с изменением амплитуды и фазы за счет комплексного множителя H(x,y), который носит название частотного отклика (частотной характеристики) системы. Для дискретных сигналов частотный отклик периодичен с периодом 2 по обеим частотным переменным:
H(x+2k,y+2l) = H(x,y).
Пример расчета частотного отклика системы.
Определить частотную характеристику системы с импульсным откликом:
h(0,0) = 0.25, h(0,1) = 0.125, h(1,0) = 0.125, h(1,1) = 0.0625.
Частотный отклик:
H(x,y) =h(n,m)exp(-jnx-jmy) = 0.25+0.125[exp(-jx)+exp(jx)+exp(-y)+exp(jy)]+
+0.0625[exp(-jx-jy)+exp(-jx+jy)+exp(jx-jy)+exp(jx+jy)] = 0.25(1+cos x)(1+cos y).
Система является примером двумерного фильтра нижних частот. Частотный отклик системы на плоскости (x,y), приведенный на рис. 13.3.1, имеет осевую симметрию с коэффициентом передачи 1 в центре (x=0, y=0) со спадом до нуля при x= и y=.
Рис.
13.3.1. Частотная характеристика ФНЧ.
h(k,l)= q(k)g(l) Q(x)G(y)= H(x,y)
Q(x) = k q(k) exp(-jkx).
G(y) = l g(l) exp(-jly).
Импульсный отклик системы. Выражение (13.3.1) описывает разложение функции Н(x,y) в двумерный рад Фурье с коэффициентами разложения в виде отсчетов импульсного отклика h(k,l), т.е. прямое преобразование Фурье. Очевидно, что обратным преобразованием Фурье с интегрированием в пределах одного периода из частотного отклика H(x,y) можно получить импульсный отклик системы:
h(k,l) = H(x,y) exp(jkx+jly) dxdy. (13.3.2)
Пример расчета импульсного отклика фильтра.
Определить импульсный отклик идеального фильтра низких частот с прямоугольной частотной характеристикой вида: H(x,y) = 1 при |x|a<, |y|b<; H(x,y) = 0 в остальных случаях.
Импульсный отклик: h(k,l) = exp(jkx+jly) dx dy.
Система разделима: h(k,l)= exp(jkx) dx exp(jly) dy= .
Пример расчета неразделимого импульсного отклика.
Определить импульсный отклик идеального кругового фильтра нижних частот:
H(x,y) = 1 при x2+y2 <R2<2; H(x,y) = 0 в остальных случаях.
Вычисления по круговой области целесообразно выполнять в полярных координатах: =,
arctg(y/x), = arctg(m/n), при этом выражение 13.3.2 перепишется в следующем виде:
h(n,m) = exp[j cos(-)] dd=
= Jo( ) d= (R/2 J1(R ) / ,
где Jo(…), J1(…)- функции Бесселя 1-го рода 0-го и 1-го порядков соответственно.
На рис. 13.3.2 приведена пространственная форма импульсного отклика фильтра, расчет которой проведен при R = 1 с ограничением по N = 10 и M = 10, и сечения отклика по координате m.
Рис. 13.3.2. Круговой низкочастотный фильтр (справа - сечения по координате m).
Свойства двумерного преобразования Фурье. Вышеприведенные преобразования импульсного отклика в частотный отклик и наоборот представляют собой двумерные дискретные преобразования Фурье с прямоугольным растром дискретизации информации, эквивалентные одномерным преобразованиям. На двумерные преобразования с прямоугольным растром переносятся и другие свойства одномерных систем. В частности:
1. Фурье-преобразования сигналов.
S(x,y) = n m s(n,m) exp(-jnx-jmy). (13.3.3)
s(n,m) =S(x,y) exp(jnx+jmy) dxdy. (13.3.4)
2. Теорема о свертке.
z(n,m) = h(n,m) ** s(n,m) H(x,y) S(x,y) = Z(x,y).
z(n,m) = c(n,m) s(n,m) C(x,y) ** S(x,y) = Z(x,y).
3. Основные свойства Фурье-преобразования.
1) Линейность (в том числе для любых комплексных чисел a и b):
аs(n,m)+bz(n,m) aS(x,y)+bZ(x,y).
2) Пространственный сдвиг:
s(n-N,m-M) S(x,y) exp(-jNx-jMy).
3) Дифференцирование:
dS(x,y)/dx -jn s(n,m),
dS(x,y)/dy -jm s(n,m),
d2S(x,y)/(dx dy) -nm s(n,m).
4) Комплексное сопряжение:
х*(n,m) S*(-x,-y).
Вещественная и мнимая части Фурье-образов последовательностей s(n,m):
S(x,y) = S*(-x,-y).
Re [S(x,y)] = Re [S(-x,-y)].
Im [S(x,y)] = -Im [S(-x,-y)].
5) Теорема Парсеваля:
n m s(n,m) s*(n,m) = S(x,y) S*(x,y) dx dy.
В частности, при s(n,m) = s(n,m):
n m |s(n,m)|2 = |S(x,y)|2 dx dy,
где левая часть уравнения представляет собой полную энергию дискретного сигнала s(n,m), a функция |S(n,m)|2 - спектральную плотность энергии сигнала.