Калькулятор алиасинга сигнала

Частота Найквиста	24 кГц
Исходный тон	Выше границы

Почему высокий тон появляется на низкой частоте

Инженер направляет стробоскоп на настольный вентилятор с цветной меткой на лопасти

Ультразвуковая помеха в записи. Тон 30 кГц при частоте дискретизации 48 кГц получает алиас 18 кГц. На спектре записи может появиться линия внутри слышимого диапазона, хотя исходное колебание было выше него. Так произойдёт, если аналоговый вход пропустит этот тон: реальный фильтр перед преобразованием способен его ослабить. Число 18 кГц не доказывает, что именно такую частоту излучал источник.

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

Калькулятор рассматривает один действительный тон и равномерные идеальные отсчёты. Введите частоту тона f и частоту дискретизации fₛ в одинаковых единицах, здесь в кГц. Главный ответ есть неотрицательная эквивалентная частота в диапазоне от нуля до fₛ/2. Дополнительное число показывает частоту Найквиста, то есть верхнюю границу этого диапазона.

\begin{aligned}r&=f\bmod f_s\\f_a&=\min(r,f_s-r)\\f_N&=\frac{f_s}{2}\end{aligned}

Сначала частота сдвигается на целое число частот дискретизации, затем остаток выше половины диапазона отражается вниз. Например, остаток 30 при fₛ = 48 заменяется на 48 − 30 = 18. Эквивалентный способ получить ответ: найти расстояние от f до ближайшего целого кратного fₛ. Множитель 2 в границе Найквиста относится к отсчётам во времени, а не к двум каналам или двум битам.

Когда спектральная линия сохраняет своё место

Измерить обычный тон. Для 10 кГц при дискретизации 48 кГц ответ остаётся 10 кГц, а граница равна 24 кГц. Именно этот тон не переносится в другую частоту. Но такой расчёт ничего не говорит о других помехах на входе, точности амплитуды или шуме устройства. Отсутствие переноса одного компонента не является проверкой всей системы записи.

Разобраться с датчиком вибрации. Пусть входное колебание имеет частоту 1300 Гц, а данные собираются со скоростью 1000 отсчётов в секунду. Для формы это 1,3 и 1 кГц. Получаем алиас 0,3 кГц, то есть 300 Гц. Низкая линия в файле может быть следствием дискретизации, а не отдельной механической вибрацией. Чтобы различить причины, нужно проверить полосу датчика и аналогового тракта, а также повторить измерение с подходящими условиями.

Частота дискретизации описывает, как часто берётся отсчёт, а частота тона описывает, как часто повторяется само колебание. Запись «48 кГц» означает 48 000 отсчётов в секунду, а не 48 000 бит или байт. У многоканального прибора полезно уточнить частоту для одного канала: суммарная скорость передачи данных не заменяет период его измерений.

При частоте ниже fₛ/2 алиас равен исходному тону. Выше границы отображение складывает последовательные участки спектра в первую полосу. В отражённых участках меняется порядок частот: повышение исходного тона может уменьшать его алиас. Поэтому нельзя считать, что повышение частоты дискретизации всегда плавно уменьшает показанное число; при смене полосы оно ведёт себя иначе.

Нулевой алиас и граница Найквиста

Наблюдать синхронный генератор. Для 96 кГц и 48 кГц получаем нулевой алиас. Между соседними измерениями проходят ровно два периода, поэтому отсчёты попадают в одну и ту же фазу. Последовательность постоянна: при одной фазе она нулевая, при другой имеет ненулевой уровень. Ответ 0 кГц не означает, что генератор остановился или что его амплитуда обязательно исчезла.

Проверить два отсчёта на период. При тоне 24 кГц и дискретизации 48 кГц алиас равен 24 кГц, ровно границе Найквиста. Для синусоиды с начальной фазой ноль все отсчёты будут нулевыми. При фазе четверть периода они чередуются между +A и −A. Одна и та же частота дискретизации даёт совершенно разную видимую амплитуду, поэтому равенство не гарантирует восстановление неизвестного сигнала.

\sin(\pi n+\varphi)=(-1)^n\sin\varphi

Здесь n обозначает номер отсчёта, φ есть начальная фаза, A в примере есть амплитуда. Из последовательности на границе нельзя независимо определить неизвестные амплитуду и фазу. Для обычного сигнала с ограниченной полосой условие восстановления требует частоты дискретизации строго выше удвоенной верхней частоты и выполнения остальных предпосылок. Практический фильтр имеет переходную полосу, а не идеальную вертикальную стенку.

Положение «выше границы» относится к исходному тону. Сам вычисленный алиас всегда находится в первой полосе. Например, частота в три раза выше fₛ/2 также может попасть ровно на её край, однако исходный тон от этого не становится граничным. Разделение этих двух утверждений помогает правильно прочитать ответ и не перепутать частоту источника с частотой линии в данных.

Что необходимо проверить в реальном входном тракте

Перенести известную высокую частоту. При тоне 125 кГц и дискретизации 48 кГц ближайшее кратное равно 144 кГц. Разность 125 − 144 = −19 даёт неотрицательный алиас 19 кГц. Арифметика допускает такой перенос, но обычный аудиовход не обязан пропускать 125 кГц. Для намеренного измерения высокой полосы важны аналоговая полоса АЦП и фильтр, выделяющий нужный диапазон.

Антиалиасинговый фильтр ставят перед той дискретизацией, на которой могут наложиться частоты. Расчёт RC-фильтра помогает оценить ослабление на выбранной частоте. Однако один полюс и одна частота среза сами по себе не гарантируют достаточного подавления: нужен допустимый уровень помех и характеристики всего тракта.

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

Частоты хранятся с конечной численной точностью, а ответ округляется до шести значащих цифр. Это не точность генератора или часов прибора. Если большая частота введена с настолько малой поправкой, что она теряется в численном представлении, восстановить поправку из результата нельзя. Неравномерные отсчёты и изменение периода во время записи также выходят за пределы формулы.

Спектр, фильтрация и условия измерения

Увеличение длины записи, разрядности и входной полосы решает разные задачи. При наложении частот важно понимать, какая информация уже потеряна и на каком этапе её ещё можно сохранить.

Что произойдёт с несколькими тонами одновременно?

Каждый тон отображается отдельно. При дискретизации 48 кГц частоты 10 и 38 кГц обе дают линию 10 кГц. Их вклады могут усиливать или ослаблять друг друга в зависимости от амплитуд и фаз. Складывать исходные частоты и искать один алиас суммы нельзя: реальный спектр содержит отдельные компоненты.

Удалит ли цифровой фильтр алиас?

Если 8 и 40 кГц уже наложились на одну линию 8 кГц при дискретизации 48 кГц, фильтр после этого этапа не узнает, какая часть пришла от каждого источника. Подавлять мешающий входной тон нужно до наложения. Внутренняя передискретизация реального АЦП может менять устройство тракта, но её нельзя описать одним внешним числом fₛ.

Поможет ли больше точек FFT?

При частоте дискретизации 48 кГц преобразование 48 000 отсчётов имеет шаг частотной сетки 1 Гц. Оно позволяет подробнее изучить запись, однако граница Найквиста остаётся 24 кГц. Добавление реальных отсчётов продлевает наблюдение при прежнем шаге времени; дополнение нулями лишь уплотняет сетку отображения. Ни то ни другое не разделяет тоны, давшие одинаковые отсчёты.

Всегда ли нужно брать отсчёты чаще удвоенной частоты тона?

Для известной изолированной высокой полосы применяют полосовую дискретизацию: полезный диапазон намеренно переносится вниз. Нужно исключить перекрытие полос, выбрать аналоговый фильтр и обеспечить входную полосу АЦП. Высшая частота сама по себе не описывает эти условия. Расчёт одного алиаса не проверяет пригодность такой архитектуры для целого сигнала.

Подходит ли этот ответ для комплексных I/Q-отсчётов?

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

Можно ли использовать среднюю частоту при пропущенных отсчётах?

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

Устранят ли дополнительные биты АЦП наложение частот?

Нет. При фиксированном диапазоне больше бит означает меньший идеальный шаг по напряжению. Частота дискретизации при этом может остаться прежней, поэтому алиас и граница Найквиста не изменятся. Разрешение АЦП помогает оценить этот отдельный шаг; его нельзя подменять частотой измерений.

Есть что добавить?

Напишите своё мнение, комментарий или предложение.