1.4. Дискретизация непрерывных двумерных сигналов

Почти все дискретные последовательности получаются как результат представления некоторых реальных непрерывных сигналов. Существует много различных способов представления непрерывных сигналов - разложение в ряды Фурье и Тейлора, разложение по нетригонометрическим ортогональным функциям и т. д., однако значительно чаще других способов используется периодическая дискретизация, что частично связано с простотой ее реализации. В настоящем разделе будет рассмотрена взаимосвязь характеристик непрерывных сигналов и получаемых из них путем периодической дискретизации дискретных последовательностей. Мы сделаем это дважды - сначала для частного случая периодической дискретизации по прямоугольному растру, а затем в более общем случае использования других растров дискретизации.

1.4.1. Периодическая дискретизация по прямоугольному растру

Из нескольких способов обобщения одномерной периодической дискретизации на двумерный случай наиболее простым является периодическая дискретизация в прямоугольных координатах, которую мы будем для простоты называть прямоугольной дискретизацией. Если  - двумерный непрерывный сигнал, то дискретный сигнал , полученный из него путем прямоугольной дискретизации, имеет вид

,                (1.115)

где  и  - положительные вещественные константы, известные как горизонтальный и вертикальный интервалы дискретизации. Расположение отсчетов на плоскости  показано на рис. 1.23. Исходя из вида сформированной таким образом последовательности, необходимо ответить на два вопроса: можно ли по значениям  восстановить сигнал  и как связаны Фурье-преобразования сигналов  и ?

Рис. 1.23. Расположение отсчетов на плоскости  для случая прямоугольной дискретизации.

Прежде всего определим двумерное преобразование Фурье для непрерывных сигналов:

,                        (1.116)

.             (1.117)

Поскольку , то, используя (1.117), можно записать

.                 (1.118)

Далее преобразуем это выражение так, чтобы получить обратное преобразование Фурье для дискретных сигналов. Начнем с подстановки  и ; это даст нам правильную форму экспоненциальных множителей. Получим

.                 (1.119)

Двойной интеграл по всей плоскости  можно разложить на бесконечную последовательность интегралов, каждый из которых имеет квадратную область интегрирования площадью . Пусть  представляет собой квадратную область

, .

Тогда (1.119) можно записать следующим образом:

.

Заменяя  на  и  на , можно устранить зависимость пределов интегрирования от  и  и получить

              (1.120)

Второй экспоненциальный множитель в (1.120) равен 1 для всех значений целочисленных переменных , ,  и . Теперь выражение (1.120) записано в форме, совпадающей с обратным преобразованием Фурье, откуда можно заключить, что

                                              (1.121)

иначе

.                                    (1.122)

Выражение (1.122) и дает нам искомую взаимосвязь между преобразованиями Фурье непрерывного и дискретного сигналов. Правую часть этого выражения можно рассматривать как периодическое продолжение, или дополнение функции , дающее периодическую функцию .

Если непрерывный сигнал  является сигналом с ограниченным спектром, выражение (1.122) допускает дальнейшее упрощение. Преобразование Фурье  сигнала с ограниченным спектром имеет нулевые значения вне некоторой конечной области на плоскости . Для простоты предположим, что периоды дискретизации  и  выбраны достаточно малыми, так что

 для , .                                                     (1.123)

Тогда выражение (1.122) упрощается:

 при  и .                     (1.124)

Значения  вне этой области определяются периодичностью .

На рис. 1.24,а представлено изображение преобразования Фурье сигнала с ограниченным частотным спектром. Периодическое повторение преобразования Фурье дает периодическую функцию, показанную на рис. 1.24,б.

Рис. 1.24. а - преобразование Фурье непрерывного сигнала с ограниченным частотным спектром; б - периодическое повторение этого преобразования.

Пока  удовлетворяет уравнению (1.123), его можно восстановить по , обратив выражение (1.124):

.                  (1.125)

 

Следовательно, в данном случае непрерывный сигнал  можно восстановить по дискретному сигналу. Для иллюстрации этого положения выразим  через его преобразование Фурье

                      (1.126)

Здесь для удобства использованы обозначения  и . Выразим теперь  через :

      (1.127)

Вместе взятые выражения (1.115), (1.125) и (1.127) образуют основу двумерной теоремы отсчетов. Эта теорема утверждает, что непрерывный сигнал с ограниченным спектром может быть восстановлен по значениям его отсчетов. Для обеспечения условия (1.123) интервалы дискретизации  и  должны быть достаточно малыми, либо, что то же самое, частоты дискретизации  и  должны быть достаточно велики.

Непрерывный сигнал с неограниченным частотным спектром также можно подвергнуть дискретизации, однако в этом случае выражения (1.124) и (1.125) несправедливы, так как при периодическом повторении (1.122) в область ,  будут вносить вклад и другие копии . Это явление, как и при обработке одномерных сигналов, носит название эффекта наложения, поскольку высокочастотные компоненты  будут «маскироваться» под низкочастотные компоненты .