Обработка огибающих реализаций сигналов в реальном времени в автономных информационных системах Текст научной статьи по специальности «Физика»

электронное научно-техническое издание

НАУКА и ОБРАЗОВАНИЕ

Эя № ФС 77 - 30569. Государствен идя регистрация №0421100025. ISSN 1994- 0408_

Обработка огибающих реализаций сигналов в реальном времени в автономных информационных системах.

77-30569/342346

# 03, март 2012 Хохлов В. К.

УДК 621.396.96

МГТУ им. Н.Э. Баумана [email protected] [email protected]

В автономных информационных системах (АИС), например акустических и сейсмических, входные сигналы часто являются случайными процессами с неизвестной и изменяющейся средней частотой энергетического спектра и информация о свойствах объектов заключена в статистических характеристиках огибающих входных реализаций [3]. Поэтому в АИС необходимо выделять информативные параметры огибающих и вычислять статистики, инвариантные к средней частоте энергетического спектра.

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

Решить отмеченную проблему можно, если отсчеты огибающей сигнала брать при помощи сопряженного сигналу по Гильберту процесса.

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

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

x(t) = E(t) cos [Q0t + ф(0],

где E(t) - огибающая; wo - средняя частота; 9(t) - случайная фаза.

На основании [1] автокорреляционная функция (АКФ) огибающей Е(^) для узкополосного стационарного нормального случайного процесса может быть представлена как

KБ{т) = П°2

2

1+^М+уК^-3)!!]2

1 + 4 +у 22л(п)2 Г Т

(1)

= а2'

Б (т)]-[1 - г2 (ТЖо М®

где а 2 - дисперсия исходного узкополосного случайного процесса; г0(т) - огибающая нормированной АКФ; Щг0(т)], Е[г0(т)] - полные эллиптические интегралы первого и второго рода соответственно.

Первый член в (1) равен квадрату постоянной составляющей огибающей Е(^)

.2

па

= А (2)

2

Из (1) следует, что среднее значение квадрата огибающей Е(^) равно

^Е = Ке(0) = 2а2,

а дисперсия огибающей Е(^) будет

В = ^Е -А = ( 2 -П ]а2.

Используя (1), (2), можно найти нормированную АКФ огибающей узкополосного нормального случайного процесса

Г

кбет)-А

б

i б

П

4(4 -п)

В

(3)

г (т)+У гЛт)

Как показали расчеты АКФ по формуле (3) с учетом (1), для большинства практических случаев можно ограничиться первым членом в квадратных скобках (3), т. е.

^Ьф^4 (4)

Чтобы устранить влияние нестационарности входного процесса {х(0}, обусловленное изменением средней частоты ш0 на отсчеты АКФ, будем производить отсчеты с помощью сопряженного ему по Гильберту процесса {у(0} (рис. 1).

У) = - — |¿г= ¿{УЬт^о г- <?(г)]. Пт- г

хГО-УЮ

хГ^.,] /1X / / Л УЙ 1

V V' / \ |\ / \ 1 \ в дя

V 1 и ХЩ]

Рис. 1 - Получение отсчетов огибающей реализации сигнала

Как видно из рис. 1 при ф(:) = ®0 г- (р(:^)= пТ1, у(Ь) = 0 и х(^) = Е(^) = (-1)пЕ, .

Таким образом, для получения отсчетов огибающей Е, необходимо производить отсчеты реализации х(^) в моменты перехода через ноль сопряженного по Гильберту процесса у(^). Реализации сопряженных процессов можно получить на выходах широкополосного фазорасщепителя со сдвигом фаз на 90°.

Моменты времени, в которые берутся отсчеты, разделены интервалами Ат, к, где

Ат, к = г, - г к =

а.

На основании [2] дисперсия интервалов Ат, к равна

я(ДгЛ )=-0- [б»-Б, а

- к — а

(6)

где Бп

|/ - к— а

- значение корреляционной функции фазы при сдвиге , - к| л/ю0;

Б9 (о) = —— величина дисперсии фазы [1].

Из уравнения (6) следует, что для нахождения о(Ат,к) необходимо определить корреляционную функцию Вд>(т) при фиксированных значениях аргумента

т = р—, р = 0,1.....(к-2), (к- 1),к.

а

(7)

Выражение для В(т) может быть получено в виде степенного ряда по г0 (г) [1]

ЬИр:/ЛесЬпотац.еёи.ги/ёос/342346.Ь1:т1

3

Г I п +

т

вТ)=1 уу ПпЧ с2п(г)+ПгГ Т)+(т)+(8)

2т=1П=0 п!(п + т)! 2 4 12

где Г^п + - у- функция; г0 (т) - огибающая нормированного коэффициента корреляции случайного процесса {^)}; г02 (г) = Г (г) + г2 (г).

Причем

Г (г) = | Бх(<а )со,(с - с0 )тёс I | Бх(со )dto

г, (т) = | (<) ^п (с - с 0 )тёс I | £(с }ёс

]-1 у ]-1

)

У

где £х(с) - энергетический спектр случайного процесса {х(/)}. При прямоугольном энергетическом спектре с шириной полосы Ас

Ас

^ (<) =

^ при |с-с0| < 0 при |с-с0| >

2 ' Ас

Г0 (т) определяется выражением

Г (т)= Вт! Т

Ас У тАс

Если спектральная плотность случайного процесса аппроксимируется

гауссоидой

¿Х(с)= ^0ехР

(с-с0 ) в2

где в =

Ас

то выражение для г0 (г) имеет вид

г0 (т)= ехР

С П2 2\

вТ

V 4 у

Обозначим относительную ширину энергетического спектра Ас /с0 = а , тогда интересующие значения г0 (г) получим в виде:

гл

сл

30

3

для прямоугольного энергетического спектра

Г (т) = Г

V ®су

. I па )( па

Ч рт II Рт

для гауссова спектра

,(т) =

(

п

V ®0;

= ехр[- 4 (ра)1.

(10)

Коэффициентр определяется как

т= p—, p= 0,1,...,(к- 2), (к-l),k.

К-1 ад л

|\ к=£

уИ.1

1

«.к

К'В.Л К. Чу ш

о о.г о/, -эб о, в * о

Рис. 2 - Зависимость корреляционной функции фазы от относительной ширины полосы гауссова энергетического спектра при различных сдвигах между сечениями

Из выражений (6) и (8) видно, что корреляционная функция случайной фазы, а, следовательно, и дисперсия длительностей интервалов между нулями при значениях аргумента, определяемых уравнениями (7), для каждого из рассмотренных видов энергетических спектров являются функциями только относительной ширины полосы энергетического спектра и не зависят от средней частоты а>0.

Для удобства обозначим г - к = п, тогда математическое ожидание интервалов между отсчетами равно

Тгк )= п — .

®0

Нормированный коэффициент корреляции отсчетов огибающей, разделенных интервалами Дт/ к, получим в виде

?Е(^т1к)= ¡Ц^тГАт^Гк (11)

где Ж(Атгк) - плотность распределения вероятностей случайной величины Атг к.

Разложим функцию гЕ(Атг к) в ряд Тейлора в районе точки |г - к| л/ю0, ограничиваясь первыми тремя членами разложения

ГЕ(АТк) = Е1-кПа0 )+ ГЕ (У- кПа0 )АТк +

(12)

+2ге (у-кпа0 +...,

где АТгк = АТгк -м{АТгк) .

Подставив выражение (12) в (11), получим:

ГЕ(АТ1к)= Ф-кП00 )+ 2ГЕ кП®0 )^(АТк). (13)

На основании (4) гЕ (/- кж/®о) будет иметь вид

(¡/- кПа0 )= 2(4%)к (1-кПа0 )+ (14) кПа0 к(¡,- кПа0 Я

+ Г\\1-

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

На основании (13) с учетом (4)-(6) и (9)-(10) можно сделать вывод, что значения нормированных коэффициентов корреляция отсчетов узкополосного нормального случайного процесса, взятых с помощью сопряженного по Гильберту процесса, зависят от относительной ширины полосы энергетического спектра процесса а и количества интервалов между нулями, разделяющих эти отсчеты п, а также от вида энергетического спектра, и не зависят от средней частоты энергетического спектра входного процесса.

0.90

В.72

»54

И.ЭС

0.19 1.10

Прямоугольный ^^^ спектр

—> —-

спектр

0.Э2

Я.55

а

0.77

1.00

б

Рис. 3. Зависимости коэффициентов корреляции отсчетов огибающей, взятых с помощью сопряженного по Гильберту процесса от относительной ширины полосы энергетического спектра процесса а при сдвигах п между отсчетами а- п=1, б- п=5

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

0.90

0.72

П.54

0.36

Прямоугольный спек,гР

Гауссов ^Ч. спектр

0.10

о.зг

0.55

0.77

1.00

б

Рис. 4. Зависимости АКФ огибающей от относительной ширины полосы энергетического спектра процесса а при сдвигах п между отсчетами а- п=1,б- п=5

а

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

Анализ данных графиков показывает, что до значений а, равных 0,7-0,8, коэффициенты корреляции отсчетов, взятых по предложенной методике, практически совпадают со значениями отсчетов АКФ огибающей. Таким образом, отсчеты огибающей, выделенные с помощью сопряженного по Гильберту процесса, могут быть использованы в качестве информативных параметров, инвариантных к средней частоте сигналов с относительной шириной полосы до значений а, равных 0,7-0,8.

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

обработку сигналов в реальном масштабе времени, что важно для АИС. При цифровой обработке сигналов сопряженный по Гильберту процесс может быть получен при помощи дискретного преобразования Гильберта. Однако, для существенного упрощения цифровой части тракта обработки сигнала, без применения дискретизатора реализаций сигнала на входе, возможно применение в аналоговом тракте широкополосного фазорасщепителя [4] , который преобразует входную реализацию в два процесса, модули спектров которых равны, а фазы на всех частотах сдвинуты на 90 град.. То есть на выходе фазорасщепителя получают два процесса, сопряженых по Гильберту.

На рис. 5 приведена структурная схема тракта обработки сигнала с фазорасщепителем на входе.

Рисунок 5 - Структурная схема тракта обработки огибающей входной реализации с

фазорасщепителем на входе: 1-приемное устройство; 2- фильтр; 3- широкополосный фазорасщепитель; 4- компаратор; 5- ключ; 6- тракт принятия решения.

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

Заключение

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

Литература

1. Левин Б.Р. Теоретические основы статистической радиотехники: в 3-х кн. - М.: Сов. радио, 1974. - Кн. 1. - 552 с.

2. Хохлов В.К. Обнаружение, распознавание и пеленгация объектов в ближней локации.- М.: Изд-во МГТУ им. Н.Э. Баумана, 2005.-336 с.:ил.

3. Хохлов В.К., Гулин Ю.Ю. Выбор информативных признаков в автономных информационных системах с нейросетевыми трактами обработки сигналов Вестник МГТУ им. Н.Э. Баумана. Сер. «Приборостроение». - 2003. - № 3. - с. 70-83.

4. Авраменко А.А., Галямичев Ю.П., Ланнэ А.А. Электрические линии задержки и фазовращатели. - М.:- Связь. 1973.

electronic scientific and technical periodical

SCIENCE and EDUCATION

_EL № KS 77 -3()56'J..VaU421100025. ISSN 1994-jMOg_

Real-time waveform envelope processing in autonomous information systems

77-30569/342346

# 03, March 2012 Hohlov V.K.

Bauman Moscow State Technical University [email protected] [email protected]

The article describes the possibility of obtaining waveform envelope readings of input implementations using the Hilbert conjugated process which can be used as informative parameters of the signals that are invariant to the energy spectrum medium frequency with a relative width of the band to 0.7-0.8. In this case automatic tracking of the sample rate is performed according to the changing medium frequency of the input implementation energy spectrum.

Publications with keywords: waveform envelope, Hilbert transform, average frequency, relative width of the band, energy spectrum

Publications with words: waveform envelope, Hilbert transform, average frequency, relative width of the band, energy spectrum

References

1. Levin B.R. Teoreticheskie osnovy statisticheskoi radiotekhniki [Theoretical foundations of statistical radio engineering]. Moscow, Sov. radio Publ., 1974, book 1. 552 p.

2. Khokhlov V.K. Obnaruzhenie, raspoznavanie ipelengatsiia ob"ektov v blizhnei lokatsii [Detection, identification and direction finding of objects in the nearest location]. Moscow, Bauman MSTU Publ., 2005. 336 p.

3. Khokhlov V.K., Gulin Iu.Iu. Vybor informativnykh priznakov v avtonomnykh informatsionnykh sistemakh s neirosetevymi traktami obrabotki signalov [Selecting informative features in the autonomous information systems with neural network signal processing paths]. VestnikMGTUim. N.E. Baumana. Ser. Priborostroenie [Herald of the Bauman MSTU. Ser. Instrumentation], 2003, no. 3, pp. 70-83.

4. Avramenko A.A., Galiamichev Iu.P., Lanne A.A. Elektricheskie linii zaderzhki i fazovrashchateli [Electrical delay lines and phase shifters]. Moscow, Sviaz' Publ., 1973. 111 p.