Лекция 11

АНАЛИЗ СПЕКТРАЛЬНОЙ ПЛОТНОСТИ
МОЩНОСТИ

Спектральная плотность мощности позволяет судить о частотных свойствах случайного процесса. Она характеризует его интенсивность при различных частотах или, иначе, среднюю мощность, приходящуюся на единицу полосы частот.

Картину распределения средней мощности по частотам называют спектром мощности. Прибор, при помощи которого измеряется спектр мощности, называется анализатором спектра. Найденный в результате измерений спектр называется аппаратным спектром.

Работа анализатора спектра основана на следующих методах измерений:

• 
  методе фильтрации;
  
• 
  методе преобразования по теореме Винера-Хинчена;
  
• 
  методе Фурье-преобразования;
  
• 
  методе с использованием знаковых функций;
  
• 
  методе аппаратного применения ортогональных функций.
  

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

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

Запоминание анализируемых реализаций позволяет ускорить аппаратурный анализ и автоматизировать его.

1. МЕТОД ФИЛЬТРАЦИИ

– средняя мощность стационарного случайного процесса .

– средняя мощность в полосе в окрестности .

– спектральная плотность мощности.

В случае, когда полоса частот конечна, но очень узка, то спектральную плотность мощности можно считать постоянной в этой полосе:

.

Функциональная схема измерителя спектральной плотности мощности методом фильтрации приведена на рис. 1

Рис. 1. Функциональная схема измерителя спектральной плотности мощности

методом фильтрации: h_узк(τ) и h_фнч(t – τ) – импульсная переходная характеристика линейного полосового (узкополосного) фильтра ЛПФ и фильтра низких частот
ФНЧ-интегратора; КД – квадратичный детектор;
РУ – регистрирующее устройство
В этом методе первоначально вырезается узкая полоса спектра исследуемого процесса x(t). Затем выполняются те же операции, что и при измерении средней мощности стационарного случайного процесса: выделенная часть сигнала возводится в квадрат, интегрируется и фиксируется регистрирующим устройством.

Напряжение , где Т – длительность реализации или продолжительность анализа, соответствует оценке спектральной плотности мощности анализируемого случайного процесса.

Установим связь между напряжением и истинной в вероятностном смысле спектральной плотностью мощности .

Величина является случайной от измерения к измерению. Найдем математическое ожидание этой величины при следующих допущениях:

• 
  все звенья не содержат источников шумов;
  
• 
  детектор – безынерционное двухполупериодное квадратичное устройство;
  
• 
  усреднителем служит ФНЧ, устойчивый, коммутируемый в моменты t = 0 и t = T; его полоса пропускания намного уже, чем входного полосового фильтра;
  
• 
  анализируется стационарный в широком смысле процесс со средним значением, равным 0.
  

Выходное напряжение интегратора

,

где – импульсная переходная характеристика фнч,

,

.

Известно, что для узкополосного высокодобротного фильтра с симметричной относительно резонансной частоты амплитудно-частотной характеристикой

,

где – функция времени, медленно изменяющаяся по сравнению с . С учетом сказанного, выходное напряжение ЛПФ

.

Учтем, что можно представить произведением двух интегралов:

.

Отсюда

Введем обозначения , тогда

Усредним обе части уравнения, учтя при этом, что

,

,

:

Введем обозначения

,

а также с учетом того, что

,

окончательно получаем

.

Таким образом, среднее по ансамблю напряжение пропорционально истинной спектральной плотности мощности , причем коэффициент пропорциональности С(t,T) зависит от продолжительности T усреднений и момента t фиксирования отсчетов регистрирующим устройством.

При анализе одной реализации эргодического стационарного процесса значения , отсчитываемые в моменты t=T, флуктуируют около математического ожидания ; причем отклонения уменьшаются с увеличением постоянной времени сглаживающего фильтра.

О рациональном выборе ширины полосы узкополосного
фильтра и продолжительности интегрирования

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

Интервал интегрирования Т должен быть существенно большим интервала корреляции τ_k узкополосного процесса : . При постоянном Т уменьшение может нарушить это условие, что вызовет значительные флуктуации величины , и статистическая надежность результатов будет низкой. Дисперсия с увеличением (при фиксированном Т) уменьшается, а смещение оценки растет.

2. ПРИНЦИПЫ ПОСТРОЕНИЯ АНАЛИЗАТОРОВ
СПЕКТРА С ПОЛОСОВЫМИ ФИЛЬТРАМИ

Используются три способа аппаратного анализа спектров:

• 
  одновременный;
  
• 
  последовательный;
  
• 
  комбинированный;
  

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

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

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

Гетеродинный анализатор спектра мощности

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

Рис. 2. Гетеродинный анализатор спектра мощности.
ВУ – входное устройство; К – калибратор; С – смеситель;
ГКЧ – генератор качающейся частоты; ГР – генератор развертки;
ПГ – перестраиваемый гетеродин; ЛПФ – линейный полосовой
фильтр с полосой пропускания ; КД – квадратичный детектор;
И – интегратор; УС – усилитель; Осц – осциллограф
Выделяется тот участок спектра шириной , для составляющих которого выполняется условие

,

в случае использования фильтра с частотой среза F_C

.

3. МЕТОД ОПРЕДЕЛЕНИЯ СПЕКТРАЛЬНОЙ ПЛОТНОСТИ
МОЩНОСТИ ПО КОРРЕЛЯЦИОННОЙ ФУНКЦИИ

Спектральная плотность мощности стационарного случайного процесса и его корреляционная функция связаны согласно теореме Винера – Хинчена парой преобразования Фурье:



.

Для действительных стационарных случайных функций

.

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

,

.

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

Эти искажения можно уменьшить, если использовать следующую оценку

,

где – весовая функция, которую часто называют «окном» (корреляционным окном). Выбор «окна» зависит от характера определяемого спектра мощности.

4. АНАЛИЗ СПЕКТРАЛЬНОЙ ПЛОТНОСТИ
МОЩНОСТИ НА ОСНОВЕ ФУРЬЕ–ПРЕОБРАЗОВАНИЯ
РЕАЛИЗАЦИИ СЛУЧАЙНОГО ПРОЦЕССА

Спектр случайного процесса х(t) определяется

.

Текущий спектр

.

При аппаратном спектральном анализе в качестве оценки спектральной плотности мощности выбирается величина

,

.

Различным реализациям эргодического случайного процесса соответствуют различные спектральные плотности мощности. Истинная в вероятностном смысле спектральная плотность является средним по ансамблю значений, полученных при многократных измерениях, каждое из которых проведено за одинаковый интервал (0,Т).

.

Теоретически рассматриваемому методу присуща бесконечно большая разрешающая способность.

Используя формулу Эйлера, имеем

,

,

.

5. ФУНКЦИОНАЛЬНАЯ СХЕМА ПРИБОРА ДЛЯ АНАЛИЗА СПЕКТРА МОЩНОСТИ ПО ФУРЬЕ–ПРЕОБРАЗОВАНИЮ РЕАЛИЗАЦИИ СЛУЧАЙНОГО ПРОЦЕССА

Рассмотренная выше оценка I(f,T) – это непараметрическая оценка, так как в данном случае оценивается не числовой параметр, характеризующий спектральную плотность мощности G_x(f), а производится оценка всего хода этой функции. Так как число оцениваемых параметров бесконечно, то дисперсия оценки каждой ординаты не убывает с увеличением интервала Т. Это следует также из того, что с увеличением продолжительности интервала Т уменьшается корреляция между значениями спектральной плотности мощности, относящимися к соседним частотам.

Из сказанного следует, что рассмотренная оценка является несостоятельной и не позволяет получить надежное значение спектральной плотности мощности при частоте f, т. е. надежно определить тонкую структуру спектра мощности исследуемого случайного процесса.

В случае, когда исследуемый стационарный случайный процесс представлен записью одной реализации х(t) за весьма большой отрезок времени Т и обладает устойчивостью на этом отрезке, эффективной является следующая методика спектрального анализа.

Запись реализации х(t) разбивают на достаточно большое число N отрезков равной длительности Т₀ (также достаточно большой). На каждом отрезке определяется оценка спектральной плотности мощности

.

Затем находят среднее арифметическое всех измерений:

.

Эта оценка является состоятельной, так как

.

Блок-схема аналого-цифрового анализатора спектральной плотности мощности приведена на рис. 3. Она включает n генераторов гармонических сигналов ГГС₁ГГС_n , у каждого из которых имеется два выхода для колебаний и , соответственно. Эти колебания перемножаются на коды выборочных значений исследуемого процесса х(t) в умножителях У_1а,У_1б  У_na, У_nб, полученные произведения сглаживаются интеграторами И_1a , И_1б  И_na, И_nб. В результате на выходах последних получаем

,

.

Напряжения u_ia и u_iб через схему электронного переключателя ЭП_i подаются в устройство преобразования напряжение–код, состоящее из устройств сравнения и и временного селектора на триггерах Тг_ia и Тг_iб и клапанах К_ia и К_iб.

Рис. 3. Аналого-цифровой спектроанализатор на основе Фурье-преобразования:
ВУ – входное устройство; АЦП – аналого-цифровой преобразователь;
ГТИ – генератор тактовых импульсов; ДУ – делитель частоты;
– умножители; ГГС₁ – ГГС_n – генераторы гармонических сигналов;
интеграторы; ЭП₁ – ЭП_n – электронные переключатели;
ГПН – генератор пилообразного напряжения; НПН – нелинейный
преобразователь напряжения; – устройства сравнения;
Тг_ia , Тг_iб – триггеры управления временными селекторами;
К_ia , К_iб – клапаны временных селекторов; СИ_i – счетчик импульсов
временных селекторов; Т – таймер продолжительности измерений

Напряжения измеряются путем их сравнения с нелинейно изменяющимся напряжением U_н =, формируемым нелинейным преобразователем напряжения НПН из пилообразного напряжения генератора ГПН. Временной интервал t от начала процесса сравнения до его завершения, т. е. когда u_н станет равным u_ia (либо u_iб), измеряется счетно-импульсным методом: на время t переключается в рабочее состояние триггер Тг_ia (либо триггер Тг_iб), в результате чего открывается клапан К_ia (К_iб) и сигналы опорной частоты таймера Т поступают на соответствующий счетчик СИ_i . В течение эксперимента в счетчике СЧ_i накапливаются данные согласно выражениям:

,

,

.

6. СПЕКТРОАНАЛИЗАТОР НА ОСНОВЕ
ДИСКРЕТНОГО ПРЕОБРАЗОВАНИЯ ФУРЬЕ

Дискретное преобразование Фурье реализуется по формуле:

,

где ; m – дискретный аналог частоты, k – дискретный аналог времени.

Спектр мощности, называемый мгновенным спектром входного сигнала, равный квадрату модуля по Фурье-образу, вычисляется по формуле

.

Блок-схема устройства, реализующая дискретное преобразование

Фурье, приведена на рис. 4.


Рис. 4. Спектроанализатор на дискретном преобразовании Фурье:
НЧФ – низкочастотный фильтр; СЛП – схема линейного пропускания;
АЦП – аналого-цифровой преобразователь; АУ – арифметическое
устройство ДПФ; У – усреднитель; ЦАП – цифро-аналоговый
преобразователь; Осц – осциллограф
Схема линейного пропускания устанавливает продолжительность эксперимента (временное окно), т. е. длительность исследуемой реализации. Дискретизация сигнала осуществляется в схеме АЦП, частота дискретизации f_дискр 2f_макс, где f_макс – максимальная частота в спектре сигнала х(t). Низкочастотный фильтр имеет частоту среза, равную частоте Найквиста (2f_макс). При этом частотная характеристика фильтра должна иметь в области частоты Найквиста крутой излом, однако требования к крутизне уменьшаются по мере увеличения частоты f_дискр по сравнению с частотой Найквиста.

Следует также отметить, что в ряде случаев при обработке сигнала х(t), имеющего периодические компоненты, конечная длительность временного окна приводит к тому, что в спектре сигнала каждая линия сопровождается боковыми лепестками.

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

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

Более подробно алгоритмы БПФ, организация параллельно-поточных вычислений при спектральном анализе данного типа мы будем рассматривать на семинарских занятиях.