Теорема отсчётов: что считает калькулятор
Оцифровка любого сигнала — от звука до радиочастотного тракта — начинается с дискретизации: непрерывную волну заменяют последовательностью отсчётов, снятых с постоянной частотой fs. Теорема отсчётов Котельникова–Найквиста отвечает на главный вопрос: достаточно ли часто мы берём эти отсчёты, чтобы потом восстановить сигнал без потерь. Калькулятор принимает частоту дискретизации fs и частоту полезного сигнала f, вычисляет частоту Найквиста fs/2, проверяет условие корректной оцифровки и — если оно нарушено — показывает, на какую ложную «частоту-двойник» свернётся сигнал. Тема тесно связана с работой АЦП и с построением цифрового антиалиасингового фильтра перед преобразователем.
fНайквиста = fs / 2 условие: f < fs/2 fмин = 2·fmax
При наложении: fнабл = |f − fs·round(f/fs)| (свёртка спектра)
Смысл прост: чтобы однозначно описать колебание частотой f, за один его период нужно снять минимум два отсчёта. Отсюда и знаменитая «двойка» — минимальная частота дискретизации fмин = 2·fmax, где fmax — самая высокая частота в спектре сигнала. Всё, что лежит ниже границы Найквиста fs/2, оцифровывается честно и восстанавливается однозначно. Всё, что выше, — уже нельзя отличить от низкочастотного двойника: информация теряется безвозвратно, и никакая последующая цифровая обработка её не вернёт.
Важно понимать, что теорема гарантирует не только достаточность отсчётов, но и точный способ восстановления: идеальный непрерывный сигнал получается сверткой отсчётов с функцией sinc(x) = sin(πx)/(πx) — то есть интерполяцией идеальным ФНЧ с частотой среза fs/2. На практике роль такого восстановителя играет реконструирующий фильтр на выходе ЦАП. Поэтому две крайние точки — вход (антиалиасинг перед АЦП) и выход (сглаживание после ЦАП) — всегда работают в паре относительно одной и той же границы Найквиста. Если хотя бы одна из них отсутствует, в спектре появляются паразитные образы, кратные fs.
Наложение спектра: откуда берётся «двойник»
Если частота сигнала превышает fs/2, спектр «складывается» относительно частоты Найквиста, как лист бумаги по линии сгиба. Реальная частота f отражается внутрь рабочего диапазона и наблюдается как ложная fнабл = |f − fs·round(f/fs)|. Это и есть наложение, или алиасинг. Классический бытовой пример — «обратное вращение» колёс автомобиля в кино: частота кадров (та же fs, только для видео) слишком мала для скорости вращения спицы, и глаз видит движение назад. В измерительной технике алиасинг куда коварнее: ложная линия в спектре выглядит абсолютно так же, как настоящая, и по одному лишь результату БПФ-анализа отличить их невозможно. Единственная защита — аналоговый ФНЧ перед АЦП, срезающий всё выше fs/2 ещё до дискретизации; его крутизну подбирают через порядок фильтра.
Отсчётов на период и запас по дискретизации
Отношение fs/f показывает, сколько отсчётов приходится на один период сигнала. Ровно два — это теоретический минимум, работающий лишь для идеального синуса и идеального восстановителя; на практике так не делают. При двух точках на период фаза и амплитуда восстанавливаются на грани, а любой реальный фильтр требует переходной полосы. Поэтому берут запас: аудио оцифровывают с 2,2–2,3-кратным перекрытием (отсюда 44,1 кГц при полосе 20 кГц), а в измерениях и осциллографии применяют оверсэмплинг в 5–10 раз и выше — это упрощает антиалиасинговый фильтр и повышает эффективное разрешение, о чём говорит расчёт SNR и ENOB. Обратная сторона медали — рост объёма данных: удвоение fs удваивает и битрейт потока, который нужно хранить и передавать.
Типовые частоты дискретизации
Стандартные значения fs выбраны так, чтобы с запасом перекрыть полосу сигнала конкретной задачи. Ниже — самые ходовые частоты, их частота Найквиста и характерное применение.
| fs | Найквист fs/2 | Где применяют |
| 8 кГц | 4 кГц | Телефония, голос (полоса до 3,4 кГц) |
| 16 кГц | 8 кГц | Широкополосный голос, VoIP, распознавание речи |
| 44,1 кГц | 22,05 кГц | Audio CD, музыка (слух до 20 кГц) |
| 48 кГц | 24 кГц | Видео, кино, профессиональный звук |
| 96 кГц | 48 кГц | Студийная запись, мастеринг, hi-res audio |
| 192 кГц | 96 кГц | Архивная запись, измерения, ультразвук |
Где применяют — и где не годится
Расчёт по Найквисту нужен везде, где аналоговый мир встречается с цифровым: при выборе АЦП и его тактовой частоты, при проектировании аудиотрактов и цифровых синтезаторов частоты (DDS), в осциллографии, радиоприёме и телеметрии. Он же подсказывает граничную частоту антиалиасингового фильтра и связан с истинным спектральным разрешением при БПФ. Ограничения теоремы важно помнить. Во-первых, она предполагает сигнал с ограниченным спектром — идеально «обрезанных» спектров в природе нет, поэтому запас по fs и реальный фильтр обязательны. Во-вторых, простая формула fмин = 2·fmax относится к сигналу в полосе от нуля; для узкополосного радиосигнала работает бандпасс-сэмплинг (subsampling), где алиасинг используют осознанно, чтобы перенести полосу вниз. В-третьих, для дельта-сигма АЦП частота Найквиста относится к выходному, децимированному потоку, а не к сверхвысокой частоте модулятора.
Пример
На значениях по умолчанию (fs = 1000 Гц, f = 700 Гц) частота Найквиста равна fs/2 = 500 Гц. Сигнал 700 Гц лежит выше неё, значит условие f < fs/2 нарушено и возникает наложение. Ложная наблюдаемая частота: fнабл = |700 − 1000·round(700/1000)| = |700 − 1000| = 300 Гц. На период приходится fs/f = 1000/700 ≈ 1,43 отсчёта — меньше двух, отсюда и потеря. Чтобы оцифровать 700 Гц честно, нужна частота выборки fмин = 2·700 = 1400 Гц и выше.
См. также: Зеркальная частота · АЦП ↔ напряжение · ADC SNR / ENOB · RC-фильтр (антиалиасинг) · Разрешение БПФ · Битрейт потока · все калькуляторы
|