Откуда на входе нелинейного сигнала взялась гармоника

Гармоники в нелинейных сигналах возникают из-за искажений, вызванных нелинейностью передаточных характеристик систем. В линейных цепях сигнал на выходе повторяет форму входного, но с изменением амплитуды и фазы. В нелинейных системах – усилителях, диодах, транзисторах – зависимость между входом и выходом описывается полиномом или экспоненциальной функцией. Например, для квадратичной нелинейности вида y = a₀ + a₁x + a₂x² при подаче синусоидального сигнала x(t) = A sin(ωt) на выходе появляются компоненты на частотах 0, ω и 2ω. Вторая гармоника (2ω) – прямое следствие члена a₂x².

В реальных устройствах нелинейность проявляется при превышении пороговых значений входного сигнала. Для кремниевых диодов порог открытия составляет ~0,7 В, а для GaAs-транзисторов – ~1,2 В. При амплитудах выше этих значений начинается генерация гармоник. Коэффициент нелинейных искажений (THD) для усилителей класса AB достигает 0,1–0,5% при номинальной мощности, но возрастает до 5–10% при перегрузке. Для минимизации гармоник рекомендуется использовать отрицательную обратную связь с глубиной не менее 20 дБ или применять предварительное ограничение сигнала.

В цифровых системах гармоники возникают из-за квантования и нелинейности АЦП/ЦАП. Например, 12-битный АЦП с диапазоном ±1 В имеет шаг квантования ~0,5 мВ. При входном сигнале с амплитудой 0,8 В искажения проявляются как гармоники на частотах n·f_вх, где n – целое число. Для снижения эффекта применяют дизеринг – добавление шума с амплитудой ~0,5 LSB, что размывает нелинейность и уменьшает THD на 10–15 дБ.

В оптических системах нелинейность возникает из-за зависимости показателя преломления от интенсивности света (эффект Керра). В одномодовом волокне при мощности сигнала >10 мВт начинается генерация третьей гармоники на длине волны ~500 нм. Для подавления используют волокна с низким коэффициентом нелинейности (n₂ < 2,2·10⁻²⁰ м²/Вт) или распределенное усиление с шагом <100 км.

Как нелинейные элементы искажают форму исходного сигнала

Искажения проявляются в виде ограничения амплитуды, появления новых спектральных составляющих и изменения фазовых соотношений. В полупроводниковых приборах нелинейность обусловлена экспоненциальной вольт-амперной характеристикой p-n-перехода. При подаче сигнала с амплитудой 0,7 В на кремниевый диод форма тока на выходе становится асимметричной: положительная полуволна усиливается, а отрицательная – срезается. Это приводит к генерации чётных гармоник, особенно заметных в аудиотехнике, где даже 0,1% нелинейных искажений воспринимаются как «грязь» в звуке.

В ферромагнитных материалах нелинейность связана с гистерезисом и насыщением. При намагничивании сердечника трансформатора магнитная индукция перестаёт линейно зависеть от напряжённости поля при превышении определённого уровня. Для электротехнической стали марки 3413 насыщение наступает при индукции ~1,8 Тл. Если входной сигнал вызывает колебания индукции в диапазоне 1,5–2,0 Тл, на выходе появляются нечётные гармоники (3-я, 5-я), амплитуда которых может составлять до 5% от основной частоты. Это приводит к искажению формы напряжения в силовых сетях, особенно при работе с несинусоидальными нагрузками.

Для количественной оценки искажений используют коэффициент гармоник (THD), который рассчитывается как отношение среднеквадратичного значения гармоник к основной частоте. В высококачественных аудиосистемах THD не должен превышать 0,01%, в то время как в импульсных источниках питания допустимы значения до 5%. При проектировании схем с нелинейными элементами критически важно учитывать рабочую точку: смещение её в область меньшей нелинейности (например, за счёт выбора напряжения смещения в усилителях) снижает уровень гармоник на 15–20 дБ.

Компенсация искажений возможна с помощью обратной связи или предыскажений. В усилителях мощности применение отрицательной обратной связи глубиной 20 дБ уменьшает THD в 10 раз, но требует запаса по усилению. Альтернативный подход – использование цифровых корректоров, которые в реальном времени анализируют спектр сигнала и подавляют гармоники с помощью адаптивных фильтров. Например, в системах связи стандарта LTE-Advanced применяются алгоритмы предварительного искажения сигнала, компенсирующие нелинейность усилителей мощности с точностью до 0,5%.

В импульсных устройствах нелинейность проявляется в виде интермодуляционных искажений, когда при подаче двух и более частот на входе на выходе возникают комбинационные составляющие вида f₁ ± f₂. В радиопередатчиках это приводит к появлению паразитных каналов излучения, нарушающих электромагнитную совместимость. Для подавления таких эффектов используют фильтры нижних частот с крутизной спада 40 дБ/декаду или квадратурные модуляторы с цифровой коррекцией фазы. В системах с ШИМ-модуляцией нелинейность ключевых элементов (например, MOSFET) вызывает субгармоники, которые устраняются синхронизацией тактовой частоты с несущей.

При выборе нелинейных элементов для конкретных задач необходимо учитывать их динамические характеристики. Например, в высокочастотных цепях важен параметр IP3 (точка пересечения третьего порядка), определяющий уровень интермодуляционных искажений. Для малошумящих усилителей (LNA) в приёмниках СВЧ-диапазона IP3 должна быть не менее +10 дБм, чтобы обеспечить линейность при входных сигналах до -20 дБм. В аналоговых вычислительных устройствах нелинейность операционных усилителей компенсируют схемотехническими методами, такими как каскодное включение или использование дифференциальных пар с динамической нагрузкой.

Методы математического описания гармонических составляющих

Гармонические составляющие в нелинейных сигналах возникают из-за искажений, вызванных нелинейностью систем. Основным инструментом их анализа служит разложение в ряд Фурье, позволяющее представить периодический сигнал в виде суммы синусоидальных и косинусоидальных компонент. Для дискретных сигналов применяется дискретное преобразование Фурье (ДПФ), реализуемое через алгоритм быстрого преобразования Фурье (БПФ). Коэффициенты ряда Фурье вычисляются по формулам:

Амплитуда k-й гармоники: $ a_k = \frac{2}{T} \int_{0}^{T} f(t) \cos(k \omega_0 t) dt $,
Фаза: $ b_k = \frac{2}{T} \int_{0}^{T} f(t) \sin(k \omega_0 t) dt $,
Спектральная плотность: $ c_k = \sqrt{a_k^2 + b_k^2} $.

При анализе непериодических сигналов используется преобразование Фурье, переводящее сигнал из временной области в частотную. Для сигналов с конечной энергией применяется интеграл Фурье, где спектральная плотность определяется как $ F(\omega) = \int_{-\infty}^{\infty} f(t) e^{-i \omega t} dt $. Этот метод критичен для систем с непрерывным спектром, например, в радиотехнике и акустике.

Нелинейные искажения часто моделируются с помощью полиномиальных аппроксимаций. Если входной сигнал $ x(t) $ проходит через нелинейную систему с характеристикой $ y = f(x) $, то выходной сигнал можно разложить в ряд Тейлора вокруг рабочей точки. Для слабых нелинейностей ограничиваются первыми тремя членами:

$ y(t) \approx a_0 + a_1 x(t) + a_2 x^2(t) + a_3 x^3(t) $.

При подаче на вход синусоидального сигнала $ x(t) = A \cos(\omega t) $ выходной сигнал будет содержать гармоники на частотах $ \omega, 2\omega, 3\omega $. Амплитуды гармоник зависят от коэффициентов $ a_2 $ и $ a_3 $: вторая гармоника пропорциональна $ a_2 A^2 $, третья – $ a_3 A^3 $. Этот подход применяется для оценки коэффициента нелинейных искажений (КНИ) в аудиосистемах.

$При подаче на вход синусоидального сигнала $ x(t) = A \cos(\omega t) $ выходной сигнал будет содержать гармоники на частотах $ \omega, 2\omega, 3\omega $. Амплитуды гармоник зависят от коэффициентов $ a_2 $ и $ a_3 $: вторая гармоника пропорциональна $ a_2 A^2 $, третья – $ a_3 A^3 $. Этот подход применяется для оценки коэффициента нелинейных искажений (КНИ) в аудиосистемах.$

Для анализа случайных сигналов с гармоническими составляющими используется корреляционный анализ. Автокорреляционная функция $ R(\tau) $ сигнала $ f(t) $ позволяет выявить периодические компоненты даже при наличии шума. Преобразование Фурье от $ R(\tau) $ дает спектральную плотность мощности (СПМ), где пики соответствуют частотам гармоник. Метод эффективен для сигналов с низким отношением сигнал/шум, например, в сейсмических данных или биомедицинских записях.

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

Окно Ханна: минимальные боковые лепестки (-31 дБ), умеренное разрешение.
Окно Блэкмана: сверхнизкие боковые лепестки (-58 дБ), но широкий главный лепесток.

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

Практическая реализация методов требует учета ограничений вычислительных ресурсов. Для сигналов с высокой частотой дискретизации (например, >1 МГц) БПФ может быть заменено на полифазные фильтры или алгоритмы субполосного анализа. В системах реального времени применяются рекурсивные методы, такие как рекурсивное ДПФ или адаптивные фильтры Калмана, обновляющие спектральные оценки по мере поступления новых данных. Для минимизации вычислительных затрат рекомендуется использовать библиотеки с оптимизированными реализациями БПФ (FFTW, Intel MKL).

Практическое измерение амплитуд и фаз гармоник в реальных цепях

Для измерения амплитуд и фаз гармоник в реальных цепях применяют анализаторы спектра с разрешением не хуже 1 Гц и динамическим диапазоном ≥90 дБ. Пример: анализатор Keysight N9040B с полосой 50 ГГц и погрешностью амплитуды ±0,2 дБ. Перед измерениями калибруют прибор по эталонному генератору (например, Rohde & Schwarz SMB100A) на частоте первой гармоники. Входное сопротивление анализатора должно соответствовать импедансу цепи (обычно 50 Ом), иначе возникают отражения, искажающие результаты.

Методика измерений включает:

Подачу тестового сигнала с известными параметрами (например, синусоида 1 кГц, 1 В амплитуды) на вход цепи.
Фиксацию спектра выходного сигнала в полосе до 10-й гармоники (для аудиоцепей) или до 50-й (для ВЧ-устройств).
Коррекцию результатов с учетом АЧХ измерительного тракта, полученной при калибровке.
Расчет фазового сдвига гармоник относительно основной частоты с помощью двухканального осциллографа (например, Tektronix MSO64) или фазометра (погрешность ≤0,5°).

При работе с импульсными сигналами используют оконные функции (Ханна, Блэкмана) для снижения спектральных утечек.

Типичные ошибки при измерениях:

Неучет собственных нелинейностей измерительного оборудования – даже качественные анализаторы вносят искажения ≥–80 дБc на частотах выше 1 МГц.
Игнорирование температурного дрейфа – изменение температуры на 10°C может сместить фазу 3-й гармоники на 2–3° в полупроводниковых цепях.
Неправильный выбор разрешения по частоте (RBW) – слишком широкий фильтр RBW «размазывает» гармоники, узкий увеличивает время измерений.

Для минимизации погрешностей рекомендуется проводить серию измерений (не менее 5) с последующим усреднением.

В промышленных системах (например, в преобразователях частоты) для автоматизации измерений применяют ПО типа MATLAB с пакетом Signal Processing Toolbox или специализированные программы (Keysight VSA, NI LabVIEW). Алгоритмы постобработки включают:

БПФ с перекрытием сегментов (overlap ≥50%) для повышения точности оценки фазы.
Компенсацию групповой задержки тракта с помощью полиномиальной аппроксимации.
Фильтрацию шумов медианным фильтром или вейвлет-преобразованием.

Результаты экспортируют в форматы CSV или HDF5 для дальнейшего анализа в системах мониторинга (например, SCADA).

Влияние нелинейности на спектр сигнала в усилителях и преобразователях

Спектральный анализ показывает, что нелинейность порождает гармоники на частотах, кратных основной: вторая гармоника (2f₀), третья (3f₀) и т. д. В усилителях мощности на биполярных транзисторах доминирует вторая гармоника, уровень которой может превышать −40 дБ относительно основного тона при выходной мощности 1 Вт. В MOSFET-усилителях третья гармоника часто оказывается сильнее из-за квадратичной зависимости тока стока от напряжения затвор-исток. Для снижения уровня гармоник рекомендуется:

Использовать отрицательную обратную связь (ООС) с глубиной не менее 20 дБ – это снижает THD в 10–100 раз, но требует запаса по усилению.
Выбирать транзисторы с линейной передаточной характеристикой (например, GaN HEMT для СВЧ-диапазона).
Оптимизировать рабочую точку: для биполярных транзисторов – ток коллектора 5–10 мА, для MOSFET – напряжение затвор-исток в середине линейного участка.

В импульсных преобразователях (например, DC-DC) нелинейность обусловлена дискретным характером работы ключей и паразитными индуктивностями/емкостями. При частоте коммутации 1 МГц спектр выходного сигнала содержит гармоники на частотах n·fₛ (где fₛ – частота переключения), а также интермодуляционные продукты вида m·fₛ ± k·f₀ (f₀ – частота модуляции). Уровень гармоник зависит от:

Скважности импульсов: при 50% уровень второй гармоники минимален, но растет при отклонении.
Длительности фронтов: при фронтах 10 нс на частоте 1 МГц уровень гармоник выше −30 дБ, при 1 нс – ниже −50 дБ.
Топологии: в резонансных преобразователях (LLC) гармоники подавляются на 15–20 дБ лучше, чем в жестко коммутируемых (buck, boost).

Интермодуляционные искажения (IMD) критичны в многоканальных системах, где два и более сигналов взаимодействуют в нелинейной среде. Например, в усилителях звуковой частоты при подаче двух тонов (1 кГц и 1,1 кГц) возникают продукты третьего порядка на частотах 2·1 кГц − 1,1 кГц = 0,9 кГц и 2·1,1 кГц − 1 кГц = 1,2 кГц. Уровень IMD3 в усилителях класса D может достигать −30 дБ при максимальной мощности, что требует применения фильтров нижних частот с крутизной не менее 40 дБ/декаду.

Для количественной оценки нелинейности используют параметры:

IP3 (точка пересечения третьего порядка): определяет уровень входного сигнала, при котором мощность интермодуляционных продуктов сравнивается с мощностью основных тонов. В СВЧ-усилителях IP3 должна быть на 10–15 дБ выше максимального уровня сигнала.
P1dB (точка компрессии на 1 дБ): мощность, при которой коэффициент усиления падает на 1 дБ. Для линейных усилителей P1dB должна превышать номинальную выходную мощность на 3–6 дБ.
SFDR (динамический диапазон без паразитных составляющих): разница между уровнем сигнала и максимальной гармоникой. В АЦП SFDR достигает 90–100 дБ при частоте дискретизации 100 МГц.

Подавление гармоник в преобразователях достигается методами:

Широтно-импульсной модуляции (ШИМ) с синхронизацией несущей частоты и модулирующего сигнала – снижает уровень гармоник на 20–30 дБ.
Использования многоуровневых инверторов (например, 3-уровневых) – уменьшает амплитуду гармоник в 2–3 раза по сравнению с 2-уровневыми.
Активной компенсации: введение противофазных гармоник через дополнительные ключи или цифровую обработку.

В аналоговых усилителях ключевую роль играет выбор схемотехники. Усилители с распределенным усилением (например, на лампах бегущей волны) демонстрируют THD ниже 0,01% в полосе до 1 ГГц, но требуют высоковольтного питания. Твердотельные усилители на GaAs или SiGe обеспечивают THD 0,1–0,3% при частотах до 10 ГГц, но критичны к разбросу параметров транзисторов. Для минимизации искажений рекомендуется:

Применять дифференциальные каскады – подавляют четные гармоники на 20–40 дБ.
Использовать токовое зеркало в цепях смещения – стабилизирует рабочую точку при изменении температуры.
Ограничивать размах сигнала на 1–2 дБ ниже P1dB – предотвращает насыщение и рост гармоник.

Цифровые методы коррекции нелинейности включают предискажение сигнала и адаптивную фильтрацию. В системах связи (например, LTE) предискажение на основе обратной модели усилителя снижает уровень гармоник на 10–15 дБ. Адаптивные фильтры (LMS, RLS) отслеживают изменения нелинейности в реальном времени, но требуют вычислительных ресурсов. Для эффективной работы алгоритмов необходима калибровка с использованием тестовых сигналов (например, многотоновых или шумоподобных) и обратной связи по спектру выходного сигнала.

Примеры возникновения гармоник в аудиосистемах и радиотехнике

В аудиосистемах гармоники возникают при перегрузке усилителей или динамиков. Например, при подаче синусоидального сигнала частотой 1 кГц на усилитель класса AB с коэффициентом нелинейных искажений (THD) 0,1% на выходе появятся вторая (2 кГц) и третья (3 кГц) гармоники с амплитудами ~0,07% и ~0,05% от основного тона соответственно. В ламповых усилителях THD может достигать 1–2%, но гармоники распределяются неравномерно: вторая гармоника доминирует (~1,5%), а высшие (4–5 кГц) составляют менее 0,3%. Для минимизации искажений рекомендуется использовать усилители с THD <0,05% (например, класса D на основе ШИМ) и динамики с линейной АЧХ до 20 кГц.

Источник гармоник	Типичные частоты гармоник	Методы снижения
Транзисторные усилители (класс AB)	2f, 3f, 4f (f – основная частота)	Отрицательная обратная связь, выбор транзисторов с низким V_CE(sat)
Радиопередатчики (нелинейность модуляции)	f ± nF (F – частота модуляции, n=1,2,3)	Фильтрация после модулятора, использование линейных усилителей мощности
Аналоговые микшеры (диоды в цепях смесителей)	mf₁ ± nf₂ (m,n – целые числа)	Применение балансных смесителей, снижение уровня входных сигналов

В радиотехнике гармоники критичны для систем связи. Например, в FM-передатчиках при модуляции сигналом 15 кГц на несущей 100 МГц возникают гармоники на 100 МГц ± 15n кГц (n=1,2,3). При мощности передатчика 50 Вт уровень второй гармоники (200 МГц) может достигать –40 дБ относительно несущей, что требует установки фильтров нижних частот с затуханием ≥60 дБ на 200 МГц. В приёмниках гармоники гетеродина (например, 10 МГц) создают ложные каналы приёма на 20 МГц и 30 МГц, устраняемые полосовыми фильтрами с крутизной спада ≥24 дБ/октаву.