Теорема Котельникова

Содержание

DSPL-2.0 — свободная библиотека алгоритмов цифровой обработки сигналов

Распространяется под лицензией LGPL v3

Страница проекта на SourceForge libdspl-2.0 on SourceForge

Обнаружили ошибку? Выделите ее мышью и нажмите ctrl+enter
Вводные замечания

В предыдущем разделе мы рассмотрели модель дискретного сигнала s_д(t) как результат произведения исходного аналогового сигнала s(t) и решетчатой функции ш_T(t). Мы также проанализировали спектральные свойства решетчатой функции и выяснили, что преобразование Фурье \mathcal{F}\big[ ш_T(t) \big] = \Omega ш_{\Omega}(\omega), где \Omega = 2\pi/T — циклическая частота дискретизации.

Согласно результатам предыдущего параграфа, спектр S_д(\omega)

equation 1
(1)
дискретного сигнала пропорционален сумме копий спектральных плотностей S(\omega) исходного аналогового сигнала s(t), отстоящих друг от друга на частоту дискретизации \Omega рад/c.

Мы уже отмечали, что (1) справедливо вне зависимости от величины \Omega и формы спектральной плотности исходного сигнала.

В данном разделе мы рассмотрим фундаментальную теорему, которая формулирует условия однозначного представления непрерывного сигнала s(t) набором равноотстоящих дискретных отсчетов, взятых с интервалом T с.

Дискретизация сигнала с ограниченной полосой

Пусть аналоговый видеосигнал[1] s(t) имеет спектральную плотность S(\omega), ограниченную полосу B рад/c, как это показано на рисунке 1а. Спектральная плотность S(\omega) сигнала с ограниченной полосой равна нулю, если |\omega|\geq B/2. Для определенности будем полагать, что S(\pm B/2) = 0, т.е. на границах полосы BS(\omega) строго равна нулю.

Тогда спектр дискретного сигнала S_д(\omega) при различном соотношении полосы B исходного аналогового сигнала и частоты дискретизации \Omega, согласно

(1), имеет вид как это показано на рисунке 1.

При \Omega < B (рисунок 1б) копии спектральной плотности S(\omega) частично перекрываются по частоте, как это показано штриховкой, и смешиваются в результате суммирования копий S(\omega) при дискретизации сигнала.

Спектральная плотность дискретного сигнала  при различном соотношении полосы сигнала  и частоты дискретизации : 
	
  
	а — спектральная плотность исходного сигнала;  б — частота дискретизации  меньше полосы ;  	в — частота дискретизации  равна полосе ;  г — частота дискретизации  больше полосы 
Рисунок 1. Спектральная плотность дискретного сигнала S_д(\omega) при различном соотношении полосы сигнала B и частоты дискретизации \Omega:
а — спектральная плотность исходного сигнала;
б — частота дискретизации \Omega меньше полосы B;
в — частота дискретизации \Omega равна полосе B;
г — частота дискретизации \Omega больше полосы B

В результате, мы не можем выделить спектральную плотность исходного сигнала из S_д(\omega) из-за эффекта наложения смещенных копий S(\omega), который носит название алиасинга. Мы еще будем детально анализировать эффект алиасинга в следующих разделах.

При \Omega = B (рисунок 1в) копии спектральной плотности перестают перекрываться по частоте, потому что мы потребовали S(\pm B/2) = 0. Тогда появляется возможность выделить исходную S(\omega) из спектра дискретного сигнала. Для этого достаточно умножить S_д(\omega) на H(\omega) в виде частотной характеристики идеального фильтра нижних частот:

equation 2
(2)
показанной на рисунке 1в пунктирной линией.

При увеличении частоты дискретизации \Omega > B (рисунок 1г), копии S(\omega), входящие в спектр дискретного сигнала, разносятся по частоте, и мы также можем выделить исходную S(\omega) путем умножения на H(\omega) вида (2).

Теорема Котельникова

Таким образом, рассмотрев различные соотношения полосы исходного сигнала B и частоты дискретизации \Omega, мы можем окончательно сформулировать следующую теорему:

Теорема (Котельникова). Если спектральная плотность видеосигнала s(t) ограничена полосой B рад/c, т.е. S(\omega) = 0 при |\omega| \geq B/2, тогда он может быть представлен своими равноотстоящими дискретными отсчетами s(nT), взятыми с периодом Т = 2 \pi/B c, как:

equation 3
(3)
Для доказательства данной теоремы перейдем во временну́ю область, взяв обратное преобразование Фурье для правой и левой частей выражения (2):
equation 4
(4)
Учтем, что произведение S_д(\omega) H(\omega) в частотной области эквивалентно свертке сигналов во времени. Тогда (4) принимает вид:
equation 5
(5)
где  h(t) =  \mathcal{F}^{-1}\big[H(\omega)\big] — импульсная характеристика фильтра нижних частот H(\omega).

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

equation 6
(6)
Тогда при \Omega = 2\alpha получаем:
equation 7
(7)
Подставим (7) в (5), а также учтем представление дискретного сигнала через решетчатую функцию , рассмотренное в предыдущем разделе. Тогда получим:

equation 8
(8)
Поменяем операторы суммирования и интегрирования, а также используем фильтрующее свойство дельта-функции:
equation 9
(9)
что и требовалось доказать.

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

Отсчеты исходного дискретного сигнала показаны на рисунке 2а. Интерполирующая функция \operatorname{sinc}(\pi t/T) показана на рисунке 2б. Данная функция равна единице при t = 0, и равна нулю во все остальные моменты дискретизации t_n = nT.

Восстановление непрерывного сигнала : 
		
    
		а —исходный дискретный сигнал; б — функция ;   	в — интерполяция смещенными функциями ;    г — восстановленный аналоговый сигнал 
Рисунок 2. Восстановление непрерывного сигнала s(t):
а —исходный дискретный сигнал; б — функция \operatorname{sinc}(\pi t/T);
в — интерполяция смещенными функциями  \operatorname{sinc}(\pi (t-nT)/T);
г — восстановленный аналоговый сигнал s(t)

В соответствии с (9), интерполирующая функция смещается на момент взятия каждого отсчета, и масштабируется по амплитуде как это показано на рисунке 2в. В результате суммирования всех смещенных во времени интерполирующих функций происходит восстановление исходного непрерывного сигнала s(t) (2г). При этом значения сигнала s(t) в моменты дискретизации t = nTn = 0, \, \pm 1, \, \pm 2 , \ldots не меняются.

Краткая историческая справка

В русскоязычной литературе рассмотренная теорема носит имя В.А. Котельникова [1], однако в мировом сообществе она именуется в честь К. Шеннона [2, 3], а также Г. Найквиста [4] , и Э. Уиттакера [5]. Часто теорему называют просто теоремой дискретизации или теоремой отсчетов[2], или используют именование содержащее несколько фамилий [6].

Прошло без малого сто лет с опубликованной В.А. Котельниковым работы, а споры вокруг первенства доказательства не утихают. Приведем историческую справку.

В 1915 году Э. Уиттакер опубликовал работу [5], в которой вводил интерполяционный ряд вида (3), названный им кардинальным рядом. Однако он не ставил цели найти однозначного представления функций при помощи дискретных отсчетов, а стремился заменить «плохие», с точки зрения анализа функции (имеющие бесконечные разрывы или быстрые осцилляции), рядом (3).

В своей работе Уиттакер пишет: «Заметим, что мы можем построить бесконечное множество интерполяционных рядов проходящих через заданные точки f(x) . . . Но данные функции не удовлетворяют свойству кардинальной функции, в частности полного подавления периодических компонент с периодом меньшим 2w. Между тем все эти функции являются решением задачи "Найти аналитическое выражение для функции, имеющей равные значения для аргумента  aa + wa -wa + 2w\ldots": которая является фундаментальной задачей теории интерполяции».[3]

Таким образом, Э. Уиттакер вводит интерполяционный ряд, но не приводит теоремы об однозначном представлении функции рядом (3), а напротив, говорит, что интерполяционных рядов может быть бесконечно много.

Другой работой, которую можно считать предвестником теоремы Котельникова является статья Г. Найквиста [4], опубликованная в 1928 году. В своей работе, Найквист рассматривает вопрос безошибочной передачи телеграфного сигнала по узкополосному каналу связи с возможностью исключения межсимвольной интерференции. Главный вывод, который делает Найквист: «Частотный диапазон, который должен быть передан для определения одной полосы численно равен сигнальной скорости передачи»[4]. Найквист также не формулирует никаких теорем о представлении сигналов дискретными отсчетами.

В 1933 году В.А. Котельников опубликовал свою статью [1], в которой впервые появляется теорема:

Любую функцию F(t), состоящую из частот от 0 до f_1 периодов в секунду, можно представить рядом

equation 10
(10)
где k — целое число; \omega_1 = 2\pi f_1D_k — постоянные, зависящие от F(t).

И наоборот, любая функция F(t), представленная рядом (10), состоит лишь из частот от 0 до f_1 периодов в секунду.

К сожалению, работа В.А. Котельникова не была переведена на английский язык, и широкая научная дискуссия началась только после публикации К. Шеннона [2], в результате чего, в англоязычной литературе принято называть данную теорему именем К. Шеннона.

Мы же называем данную теорему именем В.А. Котельникова по праву первенства формулировки и доказательства[5], что, впрочем, не умоляет заслуг К. Шеннона внесшего фундаментальный вклад в развитие теории передачи информации.

Выводы

В данном разделе мы ввели рассмотрели фундаментальную теорему теории цифровой обработки сигналов: теорему Котельникова.

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

В следующем разделе мы рассмотрим некоторые эффекты практической дискретизации сигналов и некоторые обобщения теоремы Котельникова.

странице обсуждения статьи

Смотри также
Аналоговые, дискретные и цифровые сигналы

Алиасинг при дискретизации сигналов

Примечания

[1] Напомним, что видеосигналом называется вещественный или комплексный сигнал s(t), чья спектральная плотность S(\omega) сосредоточена в области нижних частот в окрестности \omega = 0.

[2] Sampling theorem

[3] «We may remark in passing that it is possible to construct an infinite number of functions cotabular with f(x) . . . But this function does not possess the property characteristic of the cardinal function, namely, that periodic constituents of period less than 2w 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 aa + wa - wa + 2w\ldots 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], профессор Ганс Дитер Люке пишет: «Вероятно, В.А. Котельников был первым ученым, давшим точную формулировку и доказательство теоремы дискретизации, применительно к теории связи.»

Список литературы
[1] Котельников В.А. О пропускной способности эфира и проволоки в электросвязи. Материалы к 1 Всесоюзному съезду по вопросам технической реконструкции дела связи и развития слаботочной промышленности,1933, Всесоюзный энергетический комитет.

[2] Shannon C. E. A mathematical theory of communication The Bell System Technical Journal, Oct. 1948, Vol. 27, Num. 4, pp. 623–656. doi 10.1002/j.1538-7305.1948.tb00917.x

[3] Shannon C. Communication in the presence of noise Proceedings of the Institute of Radio Engineers, Jan. 1949, Vol. 37, Num. 1, pp. 10–21.

[4] Nyquist H. Certain topics in telegraph transmission theory Transactions of the American Institute of Electrical Engineers, Apr. 1928, Vol. 47, Num. 2, 617–644, doi 10.1109/T-AIEE.1928.5055024

[5] Whittaker E. XVIII.— On the Functions which are represented by the Expansions of the Interpolation Theory. Proceedings of the Royal Society of Edinburgh, 1915, Vol. 35, pp. 181 – 194. doi 10.1017/S0370164600017806

[6] Джерри А. Дж. Теорема отсчётов Шеннона, её различные обобщения и приложения. Обзор. ТИИЭР, 1977, № 11, т. 65, стр. 53--89.

[7] Luke H. D. The origins of the sampling theorem. IEEE Communications Magazine, Apr. 1999, Vol. 37, Num. 4, pp. 106-108, doi 10.1109/35.755459

[8] Баскаков, С.И. Радиотехнические цепи и сигналы. Москва, ЛЕНАНД, 2016, 528 c. ISBN 978-5-9710-2464-4

[9] Гоноровский И.С. Радиотехнические цепи и сигналы Москва, Советское радио, 1977, 608 c.

[10] Bracewell R. The Fourier Transform and Its Applications McGraw-Hills, 1986, 474 c. ISBN 0-07-007-015-6

Последнее изменение страницы: 12.05.2022 (19:41:42)
Страница создана Latex to HTML translator ver. 5.20.11.14