CC BY
41
Далее имеем
f(y) Е [PN(y|xm' )/PN(y|xm)]X!
ImVm J (22)
для всех y є Yn , p> 0, X> 0,
поскольку из определения (17) следует, что для
у є Лm правая часть неравенства (22) больше 1, а
при у £ Л m она по крайней мере больше 0.
Подставляя выражение (22) для f(y) в (21), получаем
PEm ^Е [Pn (У xm)]1_XP { Е [Pn (y|xm' W .
У ImVm J (23)
p > 0, X> 0.
Поскольку X и p представляют собой произвольные положительные числа, можно выбрать Х=1/(1+р), что дает неравенство
PEm ^Е [PN(yXm)]1/(1+P) { Е [PN (y|x m' )]1/(1+Р)} ,
у [mVm J
р > 0.
(24)
Ясно, однако, что аддитивная граница представляет собой ее частный случай, соответствующий значению r=1 в неравенстве (24).
Литература: 1. Бурдаков С.Н., Верещак А.П., Гурьев В.Е., Кривенко С.А. Простая верхняя граница вероятности ошибки диагностирования/ / Радиоэлектроника и информатика. 2001. № 1. С. 4-6.
Поступила в редколлегию 27.04.2000
Рецензент: д-р техн. наук, проф. Пресняков И.Н.
Бурдаков Сергей Николаевич, начальник отдела АО НИИРИ. Адрес: Украина, 61054, Харьков, ул. Академика Павлова, 271, тел. 26-52-60.
Верещак Александр Петрович, канд. техн. наук, директор АО НИИРИ. Адрес: Украина, 61054, Харьков, ул. Академика Павлова, 271, тел. 26-52-00.
Гурьев Владимир Ефимович, начальник отдела АО НИИРИ. Адрес: Украина, 61054, Харьков, ул. Академика Павлова, 271, тел. 26-13-10.
Кривенко Станислав Анатольевич, канд. техн. наук, доцент, начальник сектора АО НИИРИ. Научные интересы: радиотехнические системы технической диагностики. Адрес: Украина, 61054, Харьков, ул. Академика Павлова, 271, тел. 26-95-36.
УДК 621.396.96’06
МОДЕЛЬ РАССЕЯНИЯ ВОЛН НА НЕОДНОРОДНОСТЯХ АТМОСФЕРЫ
КАРТАШОВ В.М., САКАЛО С.Н.
Разрабатывается модель рассеивающих объектов систем акустического и радиоакустического зондирования атмосферы в области пространственных спектров. Модель адекватно описывает особенности рассеяния волн на решетках и может использоваться для анализа свойств зондирующих сигналов.
Определение вида и характеристик рассеянного сигнала в задачах зондирования атмосферы акустическими и электромагнитными волнами сопряжено с решением достаточно сложных волновых задач [1,2], что требует специальной подготовки, значительных усилий и не всегда приводит к физически прозрачным результатам. Именно этими обстоятельствами объясняются ошибки и заблуждения, нередко встречающиеся в известной литературе, о которых говорится также в [2].
Для разработчиков систем зондирования атмосфе -ры при решении задач анализа и синтеза зондирующих сигналов, синтеза оптимальных алгоритмов приема, оценки точностных характеристик системы целесообразно иметь более простой и физически наглядный модельный подход, основанный на понятиях и процедурах, используемых в теории систем, который воспроизводит характерные особенности процесса рассеяния и принимаемого сигнала. В данной статье разрабатывается структурнофизическая модель, отвечающая этим условиям и позволяющая с помощью достаточно простых и привычных для инженера преобразований, осно-
ванных на аппарате теории линейной фильтрации, определять вид сигнала, рассеянного различными объектами, и его основные характеристики. Возможность рассмотреть рассеяние волн с позиций линейной фильтрации обусловлена линейностью уравнений Максвелла в среде без потерь, а также линейностью волновых уравнений, описывающих распространение и рассеяние звука в диапазоне небольших значений амплитуды.
С целью выяснить характерные для задач зондирования атмосферы особенности процесса рассеяния, учитываемые далее при создании соответствующих моделей, запишем и проанализируем основные соотношения, которые определяют рассеяние электромагнитной волны на звуке. Выражение [3] для
рассеянного поля е1 (?ь t), полученное в борновс-ком приближении при достаточно общих предположениях, исходя из нестационарного волнового уравнения, имеет вид
?1Д)
1 д 2l(r,.t)
4ле2 а2 ,
ї( ?1.0=Шs
(
r,t -
|?1
V
(
r,t -
|Ї1
V
d3r
?1 - ?
(1)
c
c
где ? — радиус-вектор точки пространства; t —
время; е — вектор напряженности электрического поля; s — изменения диэлектрической проницаемости среды, вызванные звуковой волной; c — скорость распространения света. Формула (1) определяет поле в некоторой точке ?1. Интегрирование здесь производится в пределах области рассеяния,
РИ, 2001, № 2
15
содержащей значимые (отличные от нуля) значения s и e.
Будем считать излучаемые электромагнитную и звуковую волны в пределах рассеивающего объема сферическими волнами, исходящими из точки
q = 0, где расположены центры соответствующих антенн. В этой точке разместим также центр приемной апертуры. Если пренебречь изменением направления поляризации в пределах области рассеяния, то для падающего электромагнитного поля можно записать
e(r, t) = — [г - c(t - t^]exp [- j(rat - keO], (2)
r
где E — комплексная огибающая электрического поля волны; ю — несущая частота; ke =ю /с — волновое число; to — промежуток времени между моментами излучения акустического и радиосигнала (время задержки). В (2) и далее векторную величину e записываем в виде скаляра. Акустические колебания представим в виде
s(r> 0 = -[г - cs^exP[- jM - ksr)] , (3)
зоцдирования атмосферы, в которых полезный сигнал формируется в результате брэгговского рассеяния волн на пространственно-распределенных неоднородностях в виде решетки, дальнейшие рассуждения представим в общем виде. Выражение, определяющее рассеянные акустический и электромагнитный сигналы, запишем в следующей форме:
y(0 = J. (5)
Для акустических систем под e в (5) будем понимать зондирующий сигнал, под s(t, г) — естественную неоднородность показателя преломления для звуковых волн; y(t) — общее обозначение рассеянного сигнала для акустических и радиоакустических систем.
Правая часть (5) в соответствии с равенством Парсеваля может быть переписана в виде
y«=(<1 (<-т
, (6)
где k = 2 л /г — волновое число;
здесь S — комплексная огибающая акустического сигнала; cs — скорость распространения звука; ks = 2 л / X s — волновое число для звука; q — несущая частота звука.
Перепишем выражение (1) с учетом (2) и (3) [3], выполнив в нем интегрирование по угловым переменным, дающее постоянный множитель G:
і “
ei(0,t) = -G J E2r - c(t - <oK
r| 1+-C-1-cst
Se(t,k/2),Ss(t,k),Sy(t,k/2) — пространственные
спектры соответствующих функций, зависящие от времени.
Зависимость пространственных спектров Se(t,k/2),Ss(t,k) отt соответствует изменению начальной фазы пространственного колебания, что приводит к появлению дополнительного фазового множителя в спектре бегущих синусоидальных волн:
:exp(- j(ra-Q)t + j
2ke - ks I 1 + -
r>dr
(4)
s
c
Как видно из (4), первое слагаемое в показателе экспоненты не зависит от дальности r и соответствующий сомножитель может быть вынесен за знак интеграла. Этот сомножитель определяет несущую частоту рассеянного сигнала, которая сдвинута относительно излучаемого сигнала на частоту звука q . Второе слагаемое в показателе экспоненты содержит параметр расстройки условия Брэгга
q = 2ke -ksI1 +— I. При существенном отличииq
от нуля данный сомножитель вследствие осцилляций подынтегрального выражения “зануляет” значение интеграла.
Проанализировав полученное выражение (4), приходим к выводу, что при фиксированном значении t интеграл в (4) представляет собой корреляционный интеграл функций s и e, описывающих зондирующие акустический и электромагнитный сигналы.
Учитывая общность процессов рассеяния, характерных для акустических и радиоакустических систем
Ss(t,k) = Js(t,r)e_jkrdr = e_jkcstSs(k)~ e“^sfc),
Seft,-l = J e(t,2r>_jkrdr = LeJe(r>-** *
V 2 / -X 2 -X
!e - j»t 2
Se
і J ■
(7)
где cs — скорость перемещения неоднородности; Q и ю — временные частоты процессов s и e .
В (7) показатели степени экспонент записаны в предположении k = ks для первого и k/2 = ke — для второго уравнений, где^ = 2л / X s; ke = 2л / X e — пространственные частоты неоднородности и зондирующего сигнала. При этих условиях kscs = Q , kec = ю .
Тогда (6) принимает вид
y(t) = 4-eJ Se(|k(k)dk . (8)
—<X У У
Как видно из (8), для монохроматических процессов e и s частота принимаемого колебания при
выполнении условия Брэгга 2k e = ks сдвинута
16
РИ, 2001, № 2
относительно излучаемого сигнала на частоту q , для радиоакустических систем сдвиг частоты рассеянного радиосигнала равен частоте звука.
Из (6), (8) следует, что пространственный спектр рассеяного сигнала определяется произведением спектров взаимодействующих сигналов:
Sy(k/2) = Se(k/2)sS(k). (9)
Для вещественных сигналов справедливо выражение
Sy (k/2) = Se(k/2)S*(k) = Se (k / 2& (k).
Если принять во внимание, что рассеивающая волну решетка движется, то пространственное резонансное рассеяние имеет место при условии
2ke = ks(і + cs /c), и выражение (9) соответственно должно быть записано в виде
Sy(k/2) = Se(k/2)S*[k(i + cs/c)]. (10)
Поскольку отношение cs /c для радиоволн в атмосфере очень мало (~ І0 6 ), то для многих случаев, в частности при анализе энергетических аспектов рассеяния радиоволн, отличие между выражениями (9), (10) можно считать несущественным.
Sy (t, k / 2) представляет собой взаимный энергетический пространственный спектр сигналов e(t,2r) и s(t,r). Как известно [4], для вещественных
сигналов Re|Sy(t,k/2)J — четная, а Im[Sy (t, k/2)J -нечетная функции частоты. Вклад в интеграл (6) дает только вещественная часть Sy (t, k /2), поэтому
y(0
і “
— І Re
2п о
S
y
dk.
Применяя равенство Парсеваля к строго полученному соотношению (4), имеем выражение, совпадающее с (8) с точностью до амплитудного множителя к, который зависит только от расстояния до области рассеяния, и множителя, характеризующего сдвиг сигналов в пространстве:
е^Од) = Ke-j“-^t J Se
k
2
- jk
х e
c
2
(t _t0) “cst
J dk'
Выражение (9) является физически достаточно содержательным и позволяет наглядно интерпретировать многие особенности брэгговского рассеяния на решетках применительно к акустическому и радиоакустическому зондированию атмосферы.
В соответствии с (9) каждая составляющая пространственного спектра падающего сигнала выбирает из спектра неоднородности компоненту, отвечающую условию Брэгга, и рассеивается на ней.
РИ, 2001, № 2
Для получения заметного по амплитуде рассеяного сигнала в пространственном спектре неоднородности должны присутствовать компоненты с достаточной амплитудой, лежащие в полосе частот сигнала.
Следовательно, рассеивающая решетка представляет собой некоторый фильтр, характеризующийся в пространстве волновых чисел комплексным коэффициентом передачи Ss (k), а в пространстве расстояний — импульсной характеристикой.
Понятно, что алгоритм формирования отраженного сигнала (5) соответствует линейному фильтру с
импульсной характеристикой sh (t, r), которая является зеркальным отображением в пространстве неоднородности sh (t, r) = s(t, r0 - r), где r0 — координата наиболее удаленной точки рассеивающей области.
Интеграл пространственной свертки, описывающий работу фильтра, имеет следующий вид:
yft) = Je(t,ct-2rk(t,r)dr . (11)
Спектры неоднородности и импульсной характеристики фильтра связаны соотношением
Ss(t,k) = S*4(t,k>-jkr° .
Значение параметра r0 фильтра определяет задержку рассеянного сигнала в пространстве (во времени), для радиоакустических систем r0 = cst.
Импульсные характеристики рассматриваемых объектов в атмосфере в общем случае являются случайными характеристиками, независимые переменные которых — пространство и время. Поэтому их можно представить в виде суммы регулярной
sh0(t,r) и флуктуационной shf-(t,r) составляющих:
sh (t, r) = sh0 (t, r) + shf (t, r). Регулярную (когерентную) составляющую импульсной характеристики определим как математическое ожидание
(sh(t,r)) = sho(t,r), она формирует когерентную составляющую в рассеянном сигнале; случайную составляющую shf(t,r) = sh(t,r) - sh0(t,r) будем характеризовать корреляционной функцией
Kti,t2,ri,r^ = (shf(ti,r0shf(t2,r^) .
Эти две статистические характеристики полностью описывают случайный процесс s^ (t, r), определяющий рассеяние волн в атмосфере, значения которого в соответствии с [2] распределены по нормальному закону. Наличие регулярной составляющей
sM(t, r) характерно для радиоакутического зондирования (РАЗ), в акустическом зондировании (АЗ) s h0 (t> r) = 0 (регулярную компоненту естественного поля неоднородностей, обуславливающую рефракцию волн, в данной модели не учитываем).
17
Как и в sh(t,r), в Ss(k), а также в Ssh(k), можно выделить усреднением по ансамблю регулярную Sso(k), а также флуктуационную Ssf(k) составляющие.
Если учитывать только детерминированную составляющую акустического сигнала в РАЗ, то данный фильтр — фильтр с переменными (во времени и в пространстве) параметрами; в АЗ фильтр, описывающий рассеяние, имеет случайные параметры. Импульсные характеристики объектов (каналов), зависящие только от времени, принято считать системными характеристиками феноменологических моделей. Алгоритмы (5), (11) соответствуют реальному физическому механизму формирования отраженного сигнала и, следовательно, они должны быть отнесены к структурно-физическим моделям.
Физическая сторона работы фильтра состоит в следующем. Для используемой модели взаимодействующие зондирующий сигнал и неоднородность можно представить в виде суперпозиции плоских бегущих волн, распространяющихся строго параллельно:
s(t,r) = f Ss(kS)exp|j(k'sr - k'scsO-dkS ,
e(t, r) = JSe(ke)exp[j(ker - kect)]ike .
Каждая гармоническая составляющая спектра сигнала взаимодействует только с той спектральной компонентой неоднородности, для которой выполняется условие Брэгга 2ke = ks(l + cs/c). Через другие составляющие спектра неоднородности гармоника сигнала в приближении однократного рассеяния проходит не рассеиваясь, поскольку за счет интерференции дифрагированные волны гасятся.
Для рассеянного сигнала можно соответственно записать
y(t,r) = JSy(ke)exp[j(ker - 1асФа . (12)
Полагая в (12) r = 0 , получаем выражение, определяющее процесс преобразования рассеянного сигнала, представленного в пространственной форме, во временное колебание:
y(0=J Sy(ke)exp[- jkectJdke. (13)
Как следует из (13), спектр временного рассеянного сигнала, формируемого на приемной антенне, полностью определяется его пространственным спектром.
Дополним разработанную модель деталями, позволяющими уточнить частотный состав принимаемого сигнала с учетом доплеровских сдвигов. Каждая гармоническая составляющая пространственного спектра зондирующего колебания в результате рассеяния на соответствующей гармонике неоднородности, движущейся со скоростью cs, согласно формуле двойного эффекта Доплера [2] приобре-
тает частоту юр = ю [(l-cs /c)/(l + cs /с)] . Полученный при рассеянии сдвиг частоты составляет
ю -Юр = 2ю cs /c[l + cs /c]_l. Принимая во внимание условие Брэгга 2k e = k's(l + cs/c), запишем: ю = cke = cks(l + cs/c)/2 .
Тогда ю' - ю'р = k'scs = Q', т.е. сдвиг частоты каждой составляющей пространственного спектра сигнала равен временной частоте соответствующей гармоники неоднородности.
Соответствие между данными частотами описывается формулой
П' = 2k'ecs/(l + cs/c). (14)
Используя выражения (9) и (14), можно записать формулу для спектра доплеровских частот Sg (Q') рассеянного сигнала:
Sg (Q') = Sy [Q'(l + cs /c)/2cs] =
= Se[Q'(l + cJcV2Cs]S^(Q7cs). (15)
Таким образом, полученные результаты позволяют представить процесс рассеяния зондирующего сигнала на естественных или искусственно созданных неоднородностях атмосферы как прохождение его через некоторый фильтр, характеризующийся в пространстве волновых чисел комплексной частотной характеристикой (или соответствующей ей импульсной характеристикой). Сопутствующие рассеянию изменения вида и параметров сигнала интерпретируются как его частотные и фазовые искажения, определяемые особенностями неоднородности и комплексной частотной (или импульсной) характеристикой.
Литература: 1. Красненко Н.П. Акустическое зондирование атмосферы. Новосибирск: Наука, 1986. 168 с. 2. Каллистратова М.А., Кон А.И. Радиоакустическое зондирование атмосферы. М.: Наука, 1985. 200 с. 3. Кон А.И., Татарский В.И. Частотный спектр сигнала при радиоакустическом зондировании атмосферы // Изв. АН СССР. Физика атмосферы и океана. 1980. Т.16, №3. С. 219-228. 4. Тихонов В.И. Статистическая радиотехника. М.: Радио и связь, 1982. 623 с.
Поступила в редколлегию 02.11.2000
Рецензент: д-р техн. наук, проф. Сухаревский О.И.
Карташов Владимир Михайлович, канд. техн. наук, докторант ХТУРЭ. Научные интересы: методы дистанционного зондирования атмосферы. Увлечения: спорт, автомобиль. Адрес: Украина, 61166, Харьков, пр. Ленина, 14, тел. 40-95-87.
Сакало Сергей Николаевич, канд. техн. наук, декан радиотехнического факультета ХТУРЭ. Научные интересы: методы синтеза и анализа антенн. Увлечения: автомобиль. Адрес: Украина, 61166, Харьков, пр. Ленина, 14, тел. 40-94-88.
18
РИ, 2001, № 2
