Теорема Котельникова — Найквиста

Теорема Котельникова — Найквиста (также известная как теорема Найквиста — Шеннона, или теорема отсчётов) — фундаментальное положение теории связи и цифровой обработки сигналов, устанавливающее необходимое и достаточное условие для точного восстановления непрерывного (аналогового) сигнала по его дискретным отсчётам (выборкам). Теорема утверждает, что сигнал с ограниченным спектром (финитной полосой частот) может быть однозначно восстановлен без потерь информации, если частота его дискретизации (частота взятия отсчётов) не менее чем в два раза превышает максимальную частоту, содержащуюся в спектре сигнала.

История

Теорема имеет длительную историю развития, в которой участвовали несколько учёных, работавших независимо друг от друга в разных областях.

Ранние работы

В 1915 году американский математик Эдмунд Тейлор Уиттекер в своей работе «О функциях, представляемых разложениями по интерполяционным формулам» впервые математически сформулировал принцип, согласно которому любая функция с ограниченным спектром может быть восстановлена по своим значениям в равноотстоящих точках. Он вывел интерполяционную формулу, которая позже стала известна как ряд Уиттекера — Шеннона.

Вклад Гарри Найквиста

В 1928 году американский инженер шведского происхождения Гарри Найквист, работавший в Bell Labs, опубликовал статью «Certain Topics in Telegraph Transmission Theory». В ней он исследовал проблему передачи телеграфных сигналов и ввёл понятие «скорость Найквиста» — максимальную скорость передачи символов по каналу без межсимвольной интерференции. Найквист показал, что для передачи сигнала без искажений частота следования импульсов не должна превышать удвоенной полосы пропускания канала. Это был первый шаг к формулировке теоремы в контексте связи.

Формулировка Владимира Котельникова

В 1933 году советский учёный Владимир Александрович Котельников, работавший в Московском энергетическом институте, опубликовал работу «О пропускной способности эфира и проволоки в электросвязи». В ней он впервые строго математически доказал теорему для случая сигналов с ограниченным спектром и предложил практическую схему восстановления сигнала с помощью идеального фильтра нижних частот. В русскоязычной научной и инженерной традиции теорема носит его имя. Котельников также обобщил теорему на случай сигналов с полосой частот, не начинающейся с нуля.

Работы Клода Шеннона

В 1948 году американский математик и инженер Клод Шеннон, также сотрудник Bell Labs, в своей знаменитой статье «Математическая теория связи» (A Mathematical Theory of Communication) заново открыл и популяризировал теорему в контексте цифровой связи. Шеннон привёл её в виде, который стал стандартным в западной литературе, и связал с понятием энтропии и пропускной способности канала. В англоязычном мире теорема чаще всего называется теоремой Найквиста — Шеннона или теоремой отсчётов Шеннона.

Математическая формулировка

Пусть имеется непрерывный сигнал \( x(t) \), спектр которого (преобразование Фурье) равен нулю для всех частот \( |f| \ge F_{\text{max}} \). Такой сигнал называется сигналом с ограниченным спектром (финитной полосой частот). Теорема утверждает, что сигнал \( x(t) \) может быть точно восстановлен по своим дискретным отсчётам \( x(nT_s) \), взятым с частотой \( f_s = 1/T_s \), если выполняется условие:

\[ f_s \ge 2F_{\text{max}} \]

Величина \( 2F_{\text{max}} \) называется частотой Найквиста (или скоростью Найквиста). Частота \( f_s \) называется частотой дискретизации.

Формула восстановления (ряд Котельникова)

Если условие выполняется, исходный сигнал восстанавливается с помощью интерполяционной формулы:

\[ x(t) = \sum_{n=-\infty}^{\infty} x(nT_s) \cdot \text{sinc}\left( \frac{\pi}{T_s}(t - nT_s) \right) \]

где \( \text{sinc}(u) = \frac{\sin(u)}{u} \) — функция sinc (синус-кардинал). Эта сумма представляет собой свёртку дискретных отсчётов с идеальным фильтром нижних частот, имеющим прямоугольную частотную характеристику.

Ключевые понятия

Частота Найквиста

Частота Найквиста \( f_N = f_s/2 \) — это максимальная частота, которая может быть однозначно представлена в дискретном сигнале при данной частоте дискретизации. Если в сигнале присутствуют частоты выше \( f_N \), они не могут быть корректно восстановлены.

Алиасинг (наложение спектров)

Если частота дискретизации \( f_s \) меньше удвоенной максимальной частоты сигнала (\( f_s < 2F_{\text{max}} \)), возникает явление алиасинга (от англ. aliasing — наложение). Высокочастотные компоненты спектра «маскируются» под низкочастотные, что приводит к необратимым искажениям при восстановлении. Например, при дискретизации вращающегося колеса с недостаточной частотой оно может казаться вращающимся в обратную сторону или стоящим на месте.

Антиалиасинговый фильтр

Для предотвращения алиасинга перед дискретизацией аналоговый сигнал пропускают через антиалиасинговый фильтр — аналоговый фильтр нижних частот, который подавляет все частоты выше половины частоты дискретизации. Это обязательный этап в любом аналого-цифровом преобразователе (АЦП).

Практические аспекты

Реальные сигналы

На практике сигналы редко имеют строго ограниченный спектр. Большинство реальных сигналов (речь, музыка, видео) содержат бесконечный спектр, хотя его энергия быстро убывает с ростом частоты. Поэтому теорема применяется с некоторыми допущениями:

• Выбирается частота дискретизации, заведомо превышающая удвоенную максимальную значимую частоту сигнала (например, для аудио CD — 44,1 кГц при максимальной частоте 20 кГц).
• Используются антиалиасинговые фильтры с конечной крутизной спада, что вносит некоторую погрешность.

Идеальный фильтр

Формула восстановления Котельникова требует использования идеального фильтра нижних частот с бесконечной крутизной спада, что физически нереализуемо. На практике применяются фильтры с конечной импульсной характеристикой (КИХ) или бесконечной импульсной характеристикой (БИХ), которые аппроксимируют идеальный фильтр с заданной точностью.

Обобщение на комплексные сигналы

Теорема обобщается на комплексные сигналы (например, квадратурные сигналы в радиосвязи). Для комплексного сигнала с полосой частот \( B \) (от \( f_0 - B/2 \) до \( f_0 + B/2 \)) частота дискретизации должна быть не менее \( B \), а не \( 2f_{\text{max}} \). Это используется в системах с квадратурной модуляцией.

Применение

Теорема Котельникова — Найквиста является основой для всех современных цифровых систем передачи и обработки информации:

• Цифровая запись и воспроизведение звука (CD, DVD-Audio, MP3, потоковые сервисы). Частота дискретизации 44,1 кГц (для CD) выбрана на основе теоремы, чтобы охватить слышимый диапазон 20 Гц — 20 кГц.
• Цифровая фотография и видео (CCD и CMOS матрицы). Пространственная дискретизация изображения также подчиняется теореме: разрешение матрицы должно быть не менее чем вдвое выше максимальной пространственной частоты деталей сцены.
• Цифровая связь (модемы, сотовая связь, Wi-Fi, спутниковая связь). Теорема определяет требования к частоте дискретизации в приёмниках и передатчиках.
• Радиолокация и гидролокация (обработка сигналов отражённых импульсов).
• Медицинская визуализация (МРТ, КТ, УЗИ).
• Цифровые осциллографы и анализаторы спектра.

Интересные факты

• Теорема Котельникова — Найквиста является частным случаем более общей теоремы отсчётов для сигналов с ограниченным спектром в любых ортогональных базисах.
• В 1950-х годах советский математик Александр Харкевич обобщил теорему на случай сигналов, спектр которых не начинается с нуля (полосовые сигналы).
• Нарушение условий теоремы (недостаточная частота дискретизации) является причиной эффекта «муара» на фотографиях и видео при съёмке объектов с мелкой периодической структурой.
• В некоторых системах (например, в радиолокации с синтезированной апертурой) намеренно используют дискретизацию с частотой ниже Найквиста (так называемый «сжатый сенсинг» или компрессивное зондирование), но это требует специальных методов восстановления.

Критика и ограничения

Теорема является идеализацией и не учитывает ряд практических факторов:

• Конечная длительность сигнала. Реальные сигналы имеют конечную длительность, что приводит к бесконечному спектру. Теорема строго применима только к финитным во времени и по спектру сигналам, что противоречит принципу неопределённости Гейзенберга — Габора.
• Шум квантования. При аналого-цифровом преобразовании неизбежно возникает ошибка квантования, которая не учитывается теоремой.
• Неидеальность фильтров. Реальные фильтры не могут полностью подавить частоты выше \( f_N \), что приводит к остаточному алиасингу.
• Джиттер (дрожание) тактового сигнала. Нестабильность момента взятия отсчётов вносит дополнительную погрешность.

Несмотря на эти ограничения, теорема остаётся краеугольным камнем теории цифровой обработки сигналов и является обязательной для изучения в курсах радиотехники, связи и информатики.

Источники

• Котельников В. А. «О пропускной способности эфира и проволоки в электросвязи» (1933).
• Nyquist H. «Certain Topics in Telegraph Transmission Theory» (1928).
• Shannon C. E. «A Mathematical Theory of Communication» (1948).
• Whittaker E. T. «On the Functions Which are Represented by the Expansions of the Interpolation Theory» (1915).
• Марпл-мл. С. Л. «Цифровой спектральный анализ и его приложения» (1987).
• Оппенгейм А. В., Шафер Р. В. «Цифровая обработка сигналов» (1975).