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) в область
,
будут вносить вклад и другие копии
. Это явление, как и
при обработке одномерных сигналов, носит название эффекта наложения, поскольку
высокочастотные компоненты
будут «маскироваться» под низкочастотные
компоненты
.