Теорема Котельникова
![]() DSPL-2.0 — свободная библиотека алгоритмов цифровой обработки сигналов Распространяется под лицензией LGPL v3
Страница проекта на SourceForge
|
В предыдущем разделе мы рассмотрели
модель дискретного сигнала
как результат произведения исходного аналогового сигнала
и решетчатой функции
. Мы также проанализировали спектральные свойства решетчатой функции и выяснили, что преобразование Фурье
, где
— циклическая частота дискретизации.
Согласно
результатам предыдущего параграфа,
спектр




Мы уже отмечали, что (1)
справедливо вне зависимости от величины и
формы спектральной плотности исходного сигнала.
В данном разделе мы рассмотрим фундаментальную теорему, которая формулирует условия однозначного представления непрерывного сигнала набором равноотстоящих дискретных отсчетов, взятых с интервалом
с.
Пусть аналоговый видеосигнал[1] имеет спектральную плотность
, ограниченную полосу
рад/c, как это показано на рисунке 1а. Спектральная плотность
сигнала с ограниченной полосой равна нулю, если
. Для определенности будем полагать, что
, т.е. на границах полосы
,
строго равна нулю.
Тогда спектр дискретного сигнала при различном соотношении полосы
исходного аналогового сигнала и частоты дискретизации
, согласно
(1), имеет вид как это показано на рисунке 1.
При (рисунок 1б) копии спектральной плотности
частично перекрываются по частоте, как это показано штриховкой, и смешиваются в результате суммирования копий
при дискретизации сигнала.




а — спектральная плотность исходного сигнала;
б — частота дискретизации


в — частота дискретизации


г — частота дискретизации


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

При увеличении частоты дискретизации (рисунок 1г), копии
, входящие в спектр дискретного сигнала, разносятся по частоте, и мы также можем выделить исходную
путем умножения на
вида (2).
Таким образом, рассмотрев различные соотношения полосы исходного сигнала и частоты дискретизации
, мы можем окончательно сформулировать следующую теорему:
Теорема (Котельникова).
Если спектральная плотность видеосигнала ограничена полосой
рад/c, т.е.
при
, тогда он может быть представлен своими равноотстоящими дискретными отсчетами
, взятыми с периодом
c, как:




![h(t) = \mathcal{F}^{-1}\big[H(\omega)\big]](img/eqlin-23.png)

Рассматривая спектральную плотность -функции
в
мы говорили, что





Выражение представляет собой интерполяционную формулу, позволяющую восстановить аналоговый сигнал по его дискретным отсчетам, как это показано на рисунке 2.
Отсчеты исходного дискретного сигнала показаны на рисунке 2а. Интерполирующая функция показана на рисунке 2б. Данная функция равна единице при
, и равна нулю во все остальные моменты дискретизации
.


а —исходный дискретный сигнал; б — функция

в — интерполяция смещенными функциями

г — восстановленный аналоговый сигнал

В соответствии с (9), интерполирующая функция смещается на момент взятия каждого отсчета, и масштабируется по амплитуде как это показано на рисунке 2в. В результате суммирования всех смещенных во времени интерполирующих функций происходит восстановление исходного непрерывного сигнала (2г). При этом значения сигнала
в моменты дискретизации
,
не меняются.
В русскоязычной литературе рассмотренная теорема носит имя В.А. Котельникова [1], однако в мировом сообществе она именуется в честь К. Шеннона [2, 3], а также Г. Найквиста [4] , и Э. Уиттакера [5]. Часто теорему называют просто теоремой дискретизации или теоремой отсчетов[2], или используют именование содержащее несколько фамилий [6].
Прошло без малого сто лет с опубликованной В.А. Котельниковым работы, а споры вокруг первенства доказательства не утихают. Приведем историческую справку.
В 1915 году Э. Уиттакер опубликовал работу [5], в которой вводил интерполяционный ряд вида (3), названный им кардинальным рядом. Однако он не ставил цели найти однозначного представления функций при помощи дискретных отсчетов, а стремился заменить «плохие», с точки зрения анализа функции (имеющие бесконечные разрывы или быстрые осцилляции), рядом (3).
В своей работе Уиттакер пишет: «Заметим, что мы можем построить бесконечное множество интерполяционных рядов проходящих через заданные точки . . . Но данные функции не удовлетворяют свойству кардинальной функции, в частности полного подавления периодических компонент с периодом меньшим
. Между тем все эти функции являются решением задачи "Найти аналитическое выражение для функции, имеющей равные значения для аргумента
,
,
,
": которая является фундаментальной задачей теории интерполяции».[3]
Таким образом, Э. Уиттакер вводит интерполяционный ряд, но не приводит теоремы об однозначном представлении функции рядом (3), а напротив, говорит, что интерполяционных рядов может быть бесконечно много.
Другой работой, которую можно считать предвестником теоремы Котельникова является статья Г. Найквиста [4], опубликованная в 1928 году. В своей работе, Найквист рассматривает вопрос безошибочной передачи телеграфного сигнала по узкополосному каналу связи с возможностью исключения межсимвольной интерференции. Главный вывод, который делает Найквист: «Частотный диапазон, который должен быть передан для определения одной полосы численно равен сигнальной скорости передачи»[4]. Найквист также не формулирует никаких теорем о представлении сигналов дискретными отсчетами.
В 1933 году В.А. Котельников опубликовал свою статью [1], в которой впервые появляется теорема:
Любую функцию , состоящую из частот от 0 до
периодов в секунду, можно представить рядом





И наоборот, любая функция , представленная рядом (10), состоит лишь из частот от 0 до
периодов в секунду.
К сожалению, работа В.А. Котельникова не была переведена на английский язык, и широкая научная дискуссия началась только после публикации К. Шеннона [2], в результате чего, в англоязычной литературе принято называть данную теорему именем К. Шеннона.
Мы же называем данную теорему именем В.А. Котельникова по праву первенства формулировки и доказательства[5], что, впрочем, не умоляет заслуг К. Шеннона внесшего фундаментальный вклад в развитие теории передачи информации.
В данном разделе мы ввели рассмотрели фундаментальную теорему теории цифровой обработки сигналов: теорему Котельникова.
Мы рассмотрели спектральную плотность дискретного сигнала и условия, при которых возможно восстановление аналогового сигнала по имеющимся дискретным отсчетам.
В следующем разделе мы рассмотрим некоторые эффекты практической дискретизации сигналов и некоторые обобщения теоремы Котельникова.
Алиасинг при дискретизации сигналов
[1] Напомним, что видеосигналом называется вещественный или комплексный сигнал , чья спектральная плотность
сосредоточена в области нижних частот в окрестности
.
[2] Sampling theorem
[3] «We may remark in passing that it is possible to construct an infinite
number of functions cotabular with . . . But this function
does not possess the property characteristic of the cardinal function,
namely, that periodic constituents of period less than
are absent.
Such functions are, however, all of thein solutions of the problem
"To find an analytical expression for a function when we know the
values which it has for the values
,
,
,
of its
argument": which is essentially the fundamental problem of the theory
of interpolation»
[4] «The frequency range which must be transmitted to specify one band is numerically equal to the speed of signaling»
[5] В своей статье [7], профессор Ганс Дитер Люке пишет: «Вероятно, В.А. Котельников был первым ученым, давшим точную формулировку и доказательство теоремы дискретизации, применительно к теории связи.»