Вычислительные анализаторы спектра

Бурное развитие цифровой вычислительной техники дало толчок для создания цифровых анализаторов спектра, которые по совокупно сти дискретных выборок вычисляют (5.19) при замене интеграла на сумму из N выборок.

где a (s) — отсчеты весовой функции a (t), выделяющей участок реализации функции f{t) с длительностью Г_р = (JV- 1) Аt.

Получение спектра процесса численными методами в виде конечных сумм дискретных решетчатых функций (5.19) называют дискретным преобразованием Фурье (ДПФ). Оно отличается от спектра, полученного интегральными преобразованиями. Спектр ДПФ периодически размножен: повторяется по оси частот с периодом, равным частоте отсчетов /₀тс = 1/Д?. Это устраняется соответствующим выбором дискретных частот со_/?, при этом числовое преобразование Фурье называют конечным преобразованием Фурье (К11Ф).

Алгоритм быстрого преобразования Фурье (БПФ) построен так, что вычисление частотных компонентов выполняют не делением больших последовательностей на меньшие, а в обратном порядке, начиная с N исходных последовательностей, содержащих по одному члену в каждой, увеличивая затем вдвое число членов в последовательности и уменьшая вдвое число последовательностей при каждой операции умножения.

Цифровые анализаторы спектра (ЦАС) могут работать в режиме определения спектра по единственной группе учитываемых выборочных значений N_n |/| е {/(s)}, где s = 0,.

1,…, N — 1, и в режиме периодической обработки чередующихся групп N [f] е [f{nQ + 5)}, где п определяет номер обрабатываемой группы, a Q — смещение соседних обрабатываемых групп. Режим периодической обработки развивающейся во времени последовательности на соприкасающихся (Q = N) или пересекающихся (1 < Q < N) группах выборочных значений принято называть режимом анализа спектра в реальном времени. Отличительной особенностью анализа в реальном времени является то, что он проводится без потери информации. При этом различают циклический (на соприкасающихся группах) и скользящий (на пересекающихся группах) спектральные анализы.

На практике исследуемые группы выборочных значений вводятся в ЦАС непрерывным потоком с интервалом дискретности At, значение длительности которого определяется в соответствии с теоремой отсчетов верхней граничной частотой сигнала Q_B: At = я/П_н, так что количество выборочных значений сигнала, поступающего на вход ЦАС в единицу времени, равно AN = Q_B /п.

Во избежание переполнения блоков оперативной памяти вычислительного анализатора при сколь угодно длительном его функционировании в режиме реального времени скорость ввода информации f (s) не должна в среднем превышать скорости ее обработки, т. е. скорости формирования отсчетов спектральной функции 5(со). При этом, поскольку в ЦАС осуществляется групповая обработка информации и в формировании каждого отсчета спектра участвуют все N выборочных значений сигнала данной группы, всегда будет существовать задержка готовности результатов не менее чем на N тактов ввода. Фактическое время задержки полной готовности результатов всегда несколько больше указанного за счет дополнительных затрат на обработку учитываемых выборочных отсчетов в каждой группе. Конкретные значения, устанавливающие связь между допустимой скоростью ввода информации и требуемым быстродействием операционных устройств ЦАС, работающих в реальном времени, зависят от режима анализа спектра (циклический или скользящий), режима функционирования, структурной организации и от реализуемых алгоритмов обработки. Скользящий анализ спектра используется редко, поэтому будем рассматривать в дальнейшем только циклический.

Структурная схема ЦАС для циклического анализа спектра представлена на рис. 5.17. В режиме реального времени с разделением ввода и обработки сначала все N отсчетов обрабатываемой группы N_n f в течение времени Т_в через устройство ввода-вывода информации (УВВ) но мере их поступления пересылаются в оперативное запоминающее устройство (ОЗУ). Последним тактом ввода к ОЗУ подключается арифметическое устройство (АУ), которое осуществляет обработку зафиксированной группы отсчетов в соответствии с реализуемым алгоритмом ДПФ или БПФ. Служебная информация о весовых коэффициентах exp (~j2nps)/N вводится в АУ из постоянного запоминающего устройства (ПЗУ). Вывод результатов обработки может осуществляться непосредственно из АУ либо через ОЗУ. Синхронизация работы всех операционных блоков ЦАС выполняется командами, вырабатываемыми в устройстве управления (УУ).

Рис. 5−17. Структурная схема ЦАС для циклического анализа спектра В режиме реального времени с разделением ввода и обработки информации полная обработка группы выборочных значений N_n [f] должна быть выполнена до поступления первого отсчета следующей группы N ,[/], т. е. в течение одного интервала дискретности At. Таким образом, полное время обработки всех N отсчетов должно удовлетворять условию Г₍₎ < At (рис. 5.18).

Поскольку обработка информации связана с выполнением большого количества элементарных операций обработки, каждая из которых состоит из умножения и суммирования, АУ анализатора должно обладать весьма высоким быстродействием. Так, если обработка информации в реальном времени осуществляется по алгоритму ДПФ и требует, как известно, выполнения N² элементарных операций, то время выполнения одной элементарной операции определяется соотношением.

Рис. 5.18. Временная диаграмма работы ЦАС При заданном времени выполнения одной элементарной операции т_оп можно определить верхнюю частоту для анализа спектра в реальном времени.

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

т.е. сокращается в 2N/ log₂ N раз. Верхняя граница рабочего диапазона частот при этом.

Однако и в этом случае рабочий диапазон частот на практике не превышает сотен герц.

Снизить требования к быстродействию операционных блоков ЦАС (расширить верхнюю границу частотного диапазона) можно также, используя процессоры БПФ с каскадной структурой организации параллельной работы q = log₂A арифметических устройств (АУ_;.) и q — 1 блоков па;

Рис. 5.19. Структурная схема ЦАС с каскадной структурой организации параллельной работы мя ги (БП;) емкостью 2' комплексных слоев каждый. Структурная схема такого устройства приведена на рис. 5.19.

Такая структура организации ЦАС позволяет полностью обработать одну группу выборочных значений за N тактов пересылки. Таким образом, время выполнения элементарной операции определяется соотношением.

где верхняя граничная частота при заданном т.

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

Существенное расширение диапазона рабочих частот в ЦАС можно получить при совмещении ввода информации и ее обработки, достигаемом в результате использования буферного ОЗУ либо активного использования пауз между выборочными значениями исследуемого сигнала.

В ЦАС с буферной памятью независимо от реализуемого алгоритма обработки (БПФ или ДПФ) и структурной организации процессора обрабатывается одна группа выборочных значений N_n f в течение времени ввода следующей группы отсчетов Л^г" _{+ 1}[/] в буферное ОЗУ, т. е. в течение времени T_Q = NAt. Затем АУ переключается на обработку информации, накопленной в буферном ОЗУ, а освободившееся ОЗУ используется в режиме буфера для накопления следующей группы выборочных значений N_n+2[f и т. д. Схема такого устройства отличается от приведенных па рис. 5.17 и 5.19 только наличием буферного оперативного запоминающего устройства (БОЗУ) и двух двухпозиционных коммутаторов (ДК), осуществляющих переключение УВВ и АУ в каналах обмена информацией с ОЗУ и БОЗУ. Применение БОЗУ увеличивает возможное время обработки одной группы выборочных значений в N раз, что в такой же пропорции уменьшает необходимое быстродействие операционных блоков АУ.

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

а верхняя граничная частота возрастает до значения.

и может достигать единиц мегагерц.

Циклический анализ спектра в реальном времени можно осуществить, применив ДПФ с активным использованием пауз между вводимыми выборочными значениями. Это позволяет существенно снизить требования к быстродействию АУ без применения буферной памяти и коммутации потоков выходной информации.

Таким образом, основными характеристиками вычислительных анализаторов спектра являются диапазон анализируемых частот Q_B, число вычисляемых значений спектральной функции N и динамические характеристики, определяемые АЦП и входными цепями.

Цифровые анализаторы спектра кроме спектральных характеристик обычно вычисляют и статистические.

В качестве примера можно привести цифровой анализатор спектра СК4−71, предназначенный для измерения спектральных, корреляционных и статистических характеристик сигналов в реальном масштабе времени.

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

На практике отдельные характеристики процессов обычно измеряются с помощью специализированных приборов: анализаторов спектра, измерителей корреляционных и статистических характеристик, измерителей нелинейных искажений и параметров цепей. Анализатор СК4−71 представляет собой такой тип аппаратуры, в которой специфические функции многочисленных приборов моделируются с помощью программ: для изменения характера функционирования достаточно вызова соответствующей программы без аппаратурного переустройства системы. Комплекс программ анализатора спектра позволяет сочетать в одном приборе практически все функциональные возможности, необходимые для всестороннего анализа различных сигналов (рис. 5.20).

Рис. 5.20. Функциональные возможности цифрового анализатора спектра.

Работа цифрового анализатора спектра основана на вычислительном принципе определения параметров сигналов. Структурная схема анализатора СК4−71 приведена на рис. 5.21.

Рис. 5.21. Структурная схема анализатора спектра СК4−71.

Входные аналоговые сигналы по одному (А) или двум (А, Б) каналам поступают на соответствующие усилители с переменным коэффициентом усиления, которые приводят различные пределы входных сигналов (от 0,125 до 8 В) к постоянному значению, необходимому для нормального функционирования последующих трактов. Далее сигналы поступают на фильтр нижних частот (ФНЧ), где происходит выделение полосы частот, подлежащей анализу. По команде оператора ФНЧ может быть выключен. С выхода фильтров сигналы поступают на АЦП, где преобразуются в параллельный 10-разрядный двоичный код. Возможна работа как одного, так и обоих каналов. В последнем случае выборки мгновенных значений сигнала берутся одновременно в обоих каналах, что позволяет сохранить в цифровом коде информацию о фазовых соотношениях сигналов, необходимую для измерения взаимных характеристик. Частота выборки определяется кварцевым генератором и может изменяться оператором в пределах от 0,2 Гц до 100 кГц. Эта частота определяет отсчетный масштаб прибора во временной и частотной областях.

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

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

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