CC BY
47
УДК 621.396.96’06
СТАТИСТИЧЕСКИЕ ХАРАКТЕРИСТИКИ РАССЕЯННЫХ СИГНАЛОВ СИСТЕМ АКУСТИЧЕСКОГО И РАДИОАКУСТИЧЕСКОГО ЗОНДИРОВАНИЯ АТМОСФЕРЫ
КАРТАШОВ В.М.
С использованием структурно-физической модели рассеяния волн в турбулентной атмосфере анализируются вероятностные характеристики рассеянных сигналов систем акустического и радиоакустического зондирования атмосферы. Предлагается алгоритм оценивания параметров турбулентности по пространственным характеристикам сигнала радиоакустических систем.
При синтезе алгоритмов обработки локационных сигналов и интерпретации результатов измерений необходимо иметь достоверные априорные сведения о сигнале, базирующиеся на физических представлениях и механизмах, в результате действия которых рассеянный сигнал формируется. Усредненные энергетические характеристики рассеянных сигналов систем акустического и радиоакустического зондирования атмосферы, получаемые усреднением по ансамблю рассеивающих полей, достаточно подробно изучены [1,2]. Однако воссоздать реальную картину рассеяния сигнала с учетом “мгновенной” структуры рассеивающего поля неоднородностей и фазовых соотношений различных спектральных составляющих в рамках таких пред -ставлений нельзя.
В связи с этим разработана структурно-физическая модель (СФМ) рассеяния [3], учитывающая существующую (не осредненную) реализацию случайного поля неоднородностей, его динамику, особенности брэгговского рассеяния. СФМ позволяет анализировать временные и пространственные характеристики полезного сигнала, особенности его спектра. Аналогично в классической радиолокации введение обобщенной энергетической характеристики цели- эффективной площади рассеяния- упрощает некоторые расчеты и качественные представления, но существенно сужает возможности анализа процесса рассеяния и характеристик рассеянного сигнала. Усреднением по ансамблю от СФМ можно переходить к распространенной в литературе энергетической модели, обратный переход в общем случае невозможен. В данной работе анализируются временные и пространственные статистические характеристики рассеянных волновых полей систем акустического и радиоакустического зондирования атмосферы с использованием СФМ.
Поле неоднородностей, рассеивающих в атмосфере акустический сигнал, в соответствии с СФМ [3], представляет собой совокупность квазисинусои-
РИ, 2000, № 4
дальных линейных решеток, оси которых ориентированы вдоль вектора рассеяния b (направлен вдоль оси x), а пространственный период при обратном рассеянии соответствует половине длины волны X зондирующего сигнала. Амплитуды и начальные фазы решеток случайны и зависят только от структуры реальной среды. Каждую из решеток можно рассматривать как выборку узкополосного случайного процесса єДх). Радиус поперечной корреляции решеток в плоскости y, z -р k =%, число их I в рассеивающем объеме с
поперечным размером Lj : I 2L^J .
Продолжительность существования в атмосфере эквивалентных структурных образований — рассеивающих решеток — значительно влияет на характеристики рассеянного сигнала и должна учитываться, хотя бы косвенно, при интерпретации полученных результатов измерений. В связи с этим оценим время существования эквивалентных структурных образований.
Время жизни х неоднородностей с характерным пространственным размером lj можно представить как т и ljU , где u — скорость деформации (диссипации) неоднородностей с масштабом порядка lj.
Величина u2 может быть определена путем интегрирования спектральной плотности энергии E(k) в
диапазоне от kj = 2^ до <» [4]. Как известно, для изотропной турбулентности в инерционной подобласти [1] E(k) = By/3k , где B — безразмерная
постоянная; у — скорость диссипации энергии турбулентности. Тогда для масштабов l < lj
2 “ у _ у у _ 2/
u2 = 2 J By3k /3 = /By3k /3.
П
Время жизни неоднородности, за которое происходит ее существенное изменение (эволюция), определяется приближенной формулой
l / _ У _ У 2/
ТИ lju = i,iB /2 у//lj//. (1)
- V
Отметим, что (1) с точностью до множителя UB /2 совпадает с оценочной формулой, приведенной в [1].
Таким образом, время жизни неоднородности определяется ее размером (чем больше lj, тем больше х) и интенсивностью турбулентности (обратная зависимость). Для значений lj = 0Дм (соответствует длине рассеиваемой волны X = 0,2м), B = j,4 при у = ю _4 м 2 с _/ (очень слабая турбулентность) получим х = 4,2c . При у = Шм2с_/ имеем х = 0,42c.
Условие “замороженности” рассеивающего поля применительно к зондированию атмосферы можно определить следующим образом: т»т s, где xs -
9
длительность зондирующего сигнала. Следовательно, в диапазоне реальных значений параметров турбулентности и зондирующих сигналов рассеивающие решетки можно считать в первом приближении “замороженными”. В то же время, при
сильной турбулентности и большом значении Ts , с учетом того, что решетка состоит из большого числа неоднородностей, появившихся в разные моменты времени, можно ожидать, что эволюция среды способна оказать значимое влияние на характеристики рассеянного сигнала.
Для подтверждения адекватности СФМ и анализа ее особенностей применительно к акустическому и радиоакустическому зондированию рассмотрим далее влияние турбулентности на характеристики звуковой волны, распространяющейся в атмосфере. Анализ функции когерентности звуковой волны [1], распространяющейся в локально-изотропной среде с колмогоровским спектром показателя преломления, показывает, что поперечный радиус когерентности волнового поля определяется формулой
рс = (o,73C2q2ro)'/'5 * , (2)
где C 2 — структурная характеристика поля показателя преломления є для звуковой волны; q = 2^ — волновое число; ro — пройденное волной расстояние. Под поперечным радиусом когерентности здесь понимаем средний размер области волнового фронта, в пределах которой флуктуации фазы волны находятся в допустимых пределах — на расстоянии рс обеспечивается равенство единице структурной функции флуктуаций фазы D(p c) = 1.
Продольная корреляция волнового поля вдоль направления зондирования оказывается существенно более сильной и влияние турбулентности сказывается в первую очередь на флуктуациях амплитуды волны. Продольный радиус когерентности волнового поля оценивается выражением [2]
С использованием (2), (3) производились оценки продольного и поперечного радиусов когерентности звуковой волны на различных дальностях при следующих значениях параметров:
2 -7 - 23
C2 = 10 7 * * 10 м /3 ; f = 1кГц — частота сигнала. Для
сопоставления ниже приведены также соответствующие значения поперечного размера волны (волнового пучка) ро »<Fo, где ф — ширина диаграммы направленности излучателя. При ф = 50 получаем р о и 0,1ro . В качестве характерного продольного размера можно использовать пространственную протяженность акустического импульса
rs = C • xs, где C — скорость звука, т ^ — длитель-
ность импульса. Для xs = 0,1c имеем rs и 34м .
Результаты расчетов сведены в таблицу.
Анализ полученных результатов позволяет сделать
вывод, что звуковая волна при распространении в
10
r0,м о Cl Р с ,м Lc, м
10 1 144 2,7-105
100 10 36 1,7 -104
500 50 13,8 2,5-103
1000 100 9,1 1,1103
турбулентной атмосфере на некотором расстоянии от источника может быть представлена в виде совокупности когерентных волновых столбов (решеток), причем число этих столбов с расстоянием увеличивается, а поперечный размер их уменьшается.
Амплитуда радиосигнала, отраженного от звуковой волны, распространяющейся в турбулентной атмосфере, с расстоянием существенно уменьшается вследствие того, что когерентное рассеяние радиоволны имеет место только в пределах столба. Сигналы, полученные от отдельных столбов, складываются не когерентно (по мощности). По мере увеличения количества таких столбов в пределах волнового пакета результирующий отраженный сигнал уменьшается. Именно этим обстоятельством объясняется зависимость мощности Pe отраженного радиосигнала от дальности r0 в радиоаку-
_26/
стическом зондировании как Pe ~ % /5 .
Сама же акустическая волна с расстоянием способности к дифракции на турбулентных неоднородностях не теряет. Поскольку в соответствии с СФМ рассеивающие периодические неоднородности распределены вдоль направления распространения волны и имеют небольшой поперечный размер, существенно меньший, чем поперечный радиус когерентности волны, то квазикогерентное рассеяние на каждой неоднородности не ослабевает. А сложение сигналов, рассеянных отдельными неоднородностями, как и на небольших расстояниях, происходит не когерентно. Действительно, из основного уравнения акустической локации [5] следует, что мощность Ps рассеянного сигнала зависит от расстояния, пройденного волной до области рассеяния, как Ps r0 , т.е. влияет только расходимость волны.
Если бы рассеяние акустической волны происходило на протяженных в поперечном направлении образованиях, то это привело бы к дополнительному ослаблению рассеянной волны с расстоянием. Справедливость же основного уравнения акустической локации многократно подтверждалась в экспериментах и сомнений не вызывает, а следовательно, СФМ правильно описывает механизм формирования рассеянной волны.
Использование СФМ позволяет анализировать свойства рассеянного сигнала: законы распределения мгновенных значений, флуктуаций огибающей и т.д. С одной стороны, данная модель позволит получать новые сведения о сигнале, с другой, сопоставляя полученные результаты с априорными сведениями о рассеянных полях, известными в
РИ, 2000, № 4
литературе, можно оценивать адекватность модели либо достоверность имеющихся данных.
Закон распределения мгновенных значений акустического сигнала, рассеянного на турбулентных неоднородностях атмосферы, можно считать нормальным, основываясь на том, что результирующее колебание получается как сумма достаточно большого числа сигналов, отраженных отдельными решетками. Проанализируем свойства огибающей и квадратурных составляющих рассеянного сигнала. Зондирующее колебание будем полагать монохроматическим S(t) = Re[ejrot]. Тогда, в соответствии с изложенными представлениями о рассеянии, принимаемый сигнал можно записать в виде
Sr (t) = ]
Rejz К; exp[jra(t -Ti )]j = Re{K exp[jrot]},
где I — число решеток, находящихся в рассеивающем объеме; К; — амплитудный коэффициент і-й решетки; ют i — член, учитывающий фазу і-й рассеивающей решетки; К = Z К; exp(- jrox;) — комп-
i=1
лексная амплитуда рассеянного сигнала, которая в данном случае по определению равна комплексному коэффициенту рассеяния цели.
Вследствие механического перемещения и эволюции решеток значения К; и т; флуктуируют,
следовательно, параметр К также зависит от времени. Учитывая зависимость К от частоты, запишем: К(t,jffl).
Распределения вероятностей комплексного коэффициента рассеяния К(t, jra) и его модуля являются
важными характеристиками цели. Запишем К(t, jra) в следующем виде:
К (t, jra) = К exp(- j0) = At, A + jBt, ra) =
= Z К; COs(ff>T ;) + jZ К; sin(ff>T ;), i=i i=i
где к = |K| и 0 — соответственно модуль и аргумент комплексного коэффициента рассеяния, которые также являются случайными функциями t и ю ; A = К cos 0 , B = К sin 0 - квадратурные составляющие.
P( K) = ~~2 exp|
_ K
,2ст"
2 ] I. Фаза рассеянного сиг-
нала распределена равномерно на интервале [0, 2 л ],
2
средняя мощность сигнала ст равна дисперсии квадратурных составляющих.
При радиоакустическом зондировании атмосферы, особенно на малых дальностях, статистические свойства сигнала будут иными. На небольших дальностях, когда акустический волновой пакет (АВП) представляет собой единую когерентную структуру, рассеянный сигнал можно описать квазидетерминированной моделью с неопределенными параметрами.
Далее по высоте, когда вследствие действия турбулентности происходят рассмотренные выше изменения АВП, флуктуации огибающей подчиняются обобщенному распределению Релея. Характерная особенность такого рассеяния- существование регулярной (зеркальной) составляющей сигнала, физической основой формирования которой является группирование по фазе рассеивающих решеток (столбов) около некоторого доминирующего значения.
Рассмотрим пространственные характеристики рассеянных сигналов. Приведенные далее известные положения и соотношения из [1,2] хорошо объясняются изложенными модельными представлениями.
Угловой радиус корреляции AF акустического волнового поля, рассеянного совокупностью решеток, определяется угловым размером пучка, дифрагированного отдельной решеткой, AF = . Ли-
нейный радиус корреляции в поперечном направмо / „
лении рассеянного поля р rk = . С уменьше-
нием длины решеток l направленность их излучения уменьшается, следовательно, значения AF и Рrk увеличиваются.
Диапазон частот м , в котором сохраняется корреляция рассеянного в обратном направлении сигнала, определяется шириной пространственного спектра
решетки, т. е. значением параметра 2L , М и ^l . Спектральные составляющие зондирующего сигнала, выходящие за полосу м , рассеиваются на структурных образованиях с другим заполнением, следовательно, корреляция нарушается.
A и B образуются в результате сложения большого числа слабо коррелированных составляющих с ограниченными мощностями (дисперсиями) и их можно считать независимыми и нормально распределенными. Поскольку коэффициенты решеток
К; примерно одного порядка, фазовые сдвиги достаточно велики, то дисперсии A и B одинаковы, а их математические ожидания равны нулю. Распределения вероятностей модуля |К| и огибающей результирующего сигнала являются релеевскими:
Электромагнитное поле, рассеянное акустическим волновым пакетом, имеющим когерентную структуру, представляет собой, как известно, сильно коррелированный в поперечном направлении волновой пучок, образующий на поверхности земли “пятно”. По мере распространения АВП в турбулентной атмосфере пятно “расплывается”, а поперечная пространственная корреляция рассеянного поля уменьшается.
Поскольку расширение пятна рассеянного сигнала вызвано влиянием турбулентности, а его размеры
РИ, 2000, № 4
11
для каждой высоты определяются интегральным значением с2 в слое атмосферы, то этот факт можно использовать для оценки с2. Получим соответствующее выражение.
Средняя интенсивность рассеянного радиосигнала
Jr определяется формулой [2] Jr = ^К°Л. 0 2/ 2^,
/ + R0 / Рc )
где J ro — интенсивность отраженного сигнала в бестурбулентной атмосфере; Ro — поперечный размер области взаимодействия электромагнитного
и звукового пучков на расстоянии Го . Представив интенсивности в виде Jro = PS , Jr = PS , где P — мощность рассеянного сигнала; So , S — площади пятна в бестурбулентной и турбулентной атмосфере соответственно, запишем
J
R
S
= 1+
Ro,
Т S / • (4)
JR0 S0 / Pc
Известны формулы [2] для диаметра пятна отраженного радиосигнала в бестурбулентной атмосфе-
ре d0 = ^2a2 + a2)^2 , где as и ae — радиусы акустической и электромагнитной передающих антенн, и
параметра R0 = j, где k = 2Дe ’ “е
волны радиосигнала. Воспользовавшись известным представлением радиуса когерентности волны
e , ^ e - длина
(2), из (4) можно получить соотношение для определения с2 по измеренному значению S :
C2 = 0,227(S - S^/'er”83X13 •
Таким образом, СФМ адекватно отображает физические процессы, происходящие при рассеянии волн в атмосфере, и верно воспроизводит свойства рассеянных волновых полей. Полученные результаты могут использоваться при разработке феноменологических моделей, синтезе алгоритмов обработки полезных сигналов и интерпретации экспериментальных данных.
Литература: 1. Татарский В.И. Распространение волн в турбулентной атмосфере. М.: Наука, 1967. 548 с. 2. Каллистратова М.А., Кон А.И. Радиоакустическое зондирование атмосферы. М.: Наука, 1985. 200 с. 3. Петров В.А., Карташов В.М. Анализ структурно-физической модели рассеяния волн в турбулентной атмосфере // Радиотехника. 2000. Вып. 114. С. 181-184. 4. Турбулентность: принципы и применения/ Под ред. У. Фроста. М.: Мир, 1980. 535 с. 5. Красненко Н.П. Акустическое зондирование атмосферы. Новосибирск: Наука, 1986. 165 с.
Поступила в редколлегию 17.09.2000
Рецензент: д-р техн. наук Величко А. Ф.
Карташов Владимир Михайлович, канд. техн. наук, докторант ХТУРЭ. Научные интересы: методы дистанционного зондирования атмосферы. Увлечения: спорт, автомобиль. Адрес: Украина, 61166, Харьков, пр. Ленина, 14, тел.:40-95-87.
УДК 621.317
МИНИМАЛЬНАЯ ВЕРОЯТНОСТЬ ОШИБКИ ДИАГНОСТИРОВАНИЯ
БУРДАКОВ С.Н., ВЕРЕЩАКА.П., ГУРЬЕВ В.Е, КРИВЕНКО С.А.
Показывается, что тестер максимального правдоподобия зависит только от диагностируемой радиоэлектронной системы.
Рассмотрим системутехнической диагностики (СТД), представленную в общем виде на рис.1. Входные цифровые данные, характеризующие техническое состояние диагностируемого устройства, задаются, как правило, в двоичной форме, однако могут быть представлены и в другом алфавите из qi2 символов. Информация, поступающая со скоростью один символ каждые Ts секунд, содержит К диагностических признаков. Этот блок затем подается в цифровой имитатор в качестве одного из M возможных диагнозов, обозначаемых через H1, H2,..., HM. Здесь M=qK, а q — объем алфавита диагностических сообщений. Цифровой и аналоговый имитаторы вместе осуществляют отображение множества из M диагностических состояний {Hm} на множество из M сигналов {xm(t)} конечной энергии и конечной длительности T=K Ts.
Диагностируемую радиоэлектронную систему (см. рис.1) можно рассматривать как одно случайное
отображение конечного множества имитируемых сигналов {xm(t)} в принимаемый случайный процесс y(t). Прежде чем сигнал появится на выходе системы, на него накладывается целый ряд искажений.
Рис.1. Система технической диагностики
12
