Анализ амплитудного спектра и оценка среднеквадратического значения сигнала на основе дискретного преобразования Фурье (ДПФ), страница 4

  1. Исследование отношения сигнал/шум позволяет осуществлять измерение чувствительности радиоприемных устройств, оценку коэффициента шума отдельных модулей и т.д.
  2. Исследование импульсных и частотных характеристик каналов, устройств и модулей, включая (амплитудно-частотных, фазо-частотных, группового времени запаздывания).
  3. Исследование диапазона перестройки уровня и нелинейных искажений сигнала генератора, а также выходного сигнала усилителя с оценкой параметров его амплитудной характеристики.
  4.  Исследование нестабильности частоты и уровня выходного сигнала генераторов, возбудителей передатчиков, а также нестабильности передаточной характеристики усилителей и трактов.
  5. Оценка параметров модуляции, нелинейных искажений огибающей, стабильности характеристик модуляторов и демодуляторов.

3. Анализ  амплитудного спектра

Приведем алгоритм оценки СКЗ, исходя из предположения, что первая гармоника сигнала по уровню превышает все остальные. Это позволит реализовать нахождение основной гармоники путем простого сравнения между собой спектральных линий. При необходимости алгоритм можно усложнить добавлением процедуры оценки других наиболее значимых составляющих, выделяемых по уровню или по частоте. Приведем алгоритм [147] оценки СКЗ по шагам:

  1. Выбирают частоту дискретизации и с помощью АЦП получают массив u[i/fд]  в объеме N дискретных отсчетов сигнала.
  2. Накладывают на массив дискретных отсчетов временное окно, например,  окно Хемминга или Кайзера и получают взвешенный массив.
  3. По взвешенному массиву данных вычисляют прямое преобразования Фурье, получая комплексный спектр  S[ifд/N]=FFT(u[i/fд]).
  4. В спектре S[ifд/N] определяют номер M компоненты, амплитуда которой максимальна, считая ее первой гармоникой сигнала.
  5. Вычисляют СКЗ сигнала. В зависимости от поставленной задачи можно найти СКЗ1 (а) - только первой гармоники сигнала, СКЗ2 (б) - всех гармоник в рабочей полосе частот с учетом постоянной составляющей сигнала, СКЗ3 (в) - в заданной пользователем полосе пропускания:

а) ; б)  ; в)   , где ui амплитуды спектральных компонент; u0 – постоянная составляющая сигнала; N – число дискретных отсчетов в массиве данных; N1=round(Nf/fд)- номер спектральной компоненты, соответствующий положению первой гармоники с разрешением по частоте fд/Nd -  целое число, учитывающее растекание основной гармоники; H=round(NfH/fд) и L=round(NfL/fд) - номера спектральных линий, соответствующих требуемой верхней и нижней границам полосы пропусканияfH и fL; k – весовой коэффициент, учитывающий влияние вида используемого временного окна на  величину СКЗ сигнала.

Для нахождения амплитудного значения гармонического сигнала в приведенной формуле (а) правая часть умножается на . Значение k зависит не только от вида окна, но и от числа учитываемых спектральных компонент. Например, для окна Кайзера при bk=12 и d=5 оптимальное значение k=1.971309650, для  Ханна k=1.632993162, а для Хемминга k=1.586302719.

С целью апробации и детального исследования работоспособности приведенного алгоритма была создана компьютерная программа моделирования, которая позволяет оценить степень влияния объема выборки, уровня нелинейных искажений по второй гармонике (КНИ2), вида окна (Кайзера, Хэмминга, Ханна), частоты дискретизации, разрядности АЦП, уровня приведенных к входу АЦП шумов (бит), числа (d), спектральных компонент, учитывающих растекание гармоники сигнала на соседние линии сетки частот ДПФ. Рабочий интерфейс созданной интерактивной программы «Погрешность анализа амплитудного спектра на основе БПФ» приведен на рис. 1. Здесь в окошечке «Битность» не стоит метка, поэтому приведенный график соответствует границам методической погрешности для выбранных условий работы. Если провести огибающие линии сверху и снизу, то внутренняя  область будет соответствовать значениям методической погрешности для случайной начальной фазы сигнала.

Рис. 1.  График максимальной относительной погрешности для N=2048