
книги из ГПНТБ / Штагер, Е. А. Рассеяние волн на телах сложной формы
.pdfТогда
|
|
|
|
С |
|
|
|
|
|
U7 (Д„) = |
|
|
|
||
|
|
|
|
М+Т |
|
|
|
|
С — постоянная, |
|
|
О + Ао) |
|
|
|
где |
от Д0 не зависящая, |
которая однозначно опре |
|||||
деляется условием нормировки. |
|
с 2М степенями сво |
|||||
|
Учитывая, что распределение Стыодента |
||||||
боды имеет вид |
|
|
|
|
|
|
|
|
1Рст(л-, 2М) |
г ( Л<- '- т ) |
|
, |
|
||
|
Г (уИ) / 2 - м |
|
|
|
|||
|
|
|
|
|
|
||
получим |
|
|
|
С 1 |
+ 2м ) |
|
|
|
|
|
|
|
|
||
|
W (Д0) = (2М) ' W CT(Д0 у 2 М ,2 М ) . |
(5.73) |
|||||
|
Для длинных выборок, когда 2М > |
1, согласно [25], |
распреде |
||||
ление Стыодента стремится |
к нормальному |
|
|
||||
|
W CT(x, |
2Л Г )- |
У 2те |
|
0 ^2/и)]’ |
|
|
т. |
е. |
|
|
|
|
|
|
|
|
W |
|
|
|
|
(5.74) |
Таким образом, отклонения До распределены по нормальному закону с дисперсией 1/2М.
Практическое значение имеют и неоптнмальные методы обработки выборки, так как они могут быть более удобными при технической реализации. В качестве направления на тело сложной формы можно брать среднее арифметическое угловых координат, отрабатываемых
пеленгатором |
|
|
уа |
У 'ЫІМ |
(5.75) |
где М — число измерений, или |
угловую |
координату ѵр, соответ |
ствующую элементу выборки с максимальной величиной интенсив ности поля (согласно § 5.5 при этом вероятность правильной пелен гации также повышается). По аналогии с теорией связи будем называть этот метод обработки наблюдений методом переключений. Формулы (5.50), (5.52) позволяют провести количественное иссле
дование этих |
методов. |
Следуя |
[1], |
приведем основные |
результаты. |
|||
Если |
|
|
'‘а |
|
|
|
р |
|
|
Д д |
= |
‘4 |
Д р |
= |
|
||
и W (Да), W (Др) — соответствующие |
плотности вероятности откло |
|||||||
нений чп н чр |
от средних |
значений, |
то |
эффективность |
различных |
|||
12* |
|
|
|
|
|
|
|
179 |
методов обработки наблюдений можно оценивать но быстроте стремления интегральных распределений
^ \Ѵ ( Л - ) d x
-д
к единице. Результаты расчетов интегральных распределений для всех методой приведены на рис. 46 (случай М = 2) и на рис. 47 (случай Af = 5). Пунктиром обозначено распределение отклонений для одноэлементной выборки (в этом случае все методы обработки тождественно совпадают). В соответствии с теорией, наилучшие
Рис. 46. Интегральное распре деление угловых отклонений при М = 2:
1— м е т о д |
с р е д н е в ы б о р о ч н о г о ; 2 — |
|
м е т о д п е р е к л ю ч е н и й ; |
3—о п т и м а л ь |
н ы й м е т о д .
Рис 47. Интегральное распре деление угловых отклонении при М = 5:
I — м е т о д с р е д н е в ы б о р о ч н о г о ; |
2- м е - |
|
т о д п е р е к л ю ч е н и й ; |
3 -- о п т и м а л ь н ы й |
м е т о д .
результаты дает оптимальный метод обработки (5.71). При М ^ 5 метод переключений предпочтительней средневыборочного. Уже при М = 2 имеет место значительное увеличение точности пеленгации. Так, вероятность попаданий измеренных значений Ѵо, ѵР, ѵ0 в об
ласть, занимаемую телом сложной формы (Д ^ У "3 ), при М = 2 соответственно равна 91%, 94%, 98%, в то время как при единич
ном измерении эта вероятность равна 86%. |
W(Аа) и W'(Ap) при |
|
Асимптотическое |
исследование функций |
|
М - °о показывает, |
что они также стремятся |
к нормальным распре |
180
делениям с дисперсиями ІпЛ'І/А-І и 1/(2 ln Af) соответственно. Таким образом, для достаточно длинных выборок из нсоптнмалыіых мето дов обработки более предпочтительным становится метод средневыборочиого. Однако сходимость к нормальным распределениям
уфункций 117(Д0), ЩДр) очень медленная (как 1/ІпМ).
5.7.Дополнительные замечания
Вгл. 5 упрощена .математическая модель отражений от объектов сложной геометрии. В результате упрости лось описатгие отражающих объектов: их можно опреде лять либо геометрией тел, в которых статистически рав номерно распределены участки локального отражения, либо задавая распределение участков локального отра жения в пространстве. Появились новые возможности и для более полного описания отраженного поля.
Для упрощенной модели отражений (модели Дела но) вопрос о 'распределении амплитуд поля решается
элементарно. Амплитуды распределены по закону Рэлея (т. е. интенсивности — экспоненциально) независимо от геометрии тела сложной формы. Поэтому основным объ ектом исследования в гл. 5 является плотность потока энергии. Модуль этого вектора пропорционален интенсив ности отраженного поля (ЭПР тела сложной формы), а его угловые координаты несут информацию об угло вых координатах тела сложной формы.
В гл. 5 получено выражение для эффективного рас пределения плотностей потоков энергии и проанализиро ваны основные статистические закономерности флуктуа ций этой векторной величины. Главный позитивный вы вод— высокая вероятность правильной пеленгации тела сложной формы (пеленгация считается правильной, когда прямая, проведенная вдоль вектора плотности потока энергии, пересекается с отражающим объемом). Соглас но проведенным расчетам при единичном измерении эта величина колеблется в пределах 70—85%, возрастая по мере сглаживания конфигурации моделирующего объ ема. Наибольшая вероятность правильной пеленгации (из рассмотренных примеров) достигается для тела сложной формы, которое -можно представить совокупно стью участков локального отражения, распределенных в объеме эллипсоида, 'наименьшая — для совокупности
181
участков локального отражения, распределенных по по верхности параллелепипеда.
Вероятностное распределение потоков энергии сим метрично по углам для произвольных тел сложной фор мы. Размах флуктуаций угловых координат, отрабатывае мых пеленгатором, определяется угловым размером тела сложной формы в той плоскости, в 'которой проводится измерение угловой координаты.
Флуктуации угловых координат, отрабатываемых пе ленгатором, и флуктуации интенсивности поля (ЭПР тела сложной формы) статистически зависимы. Если все измерения угловых координат разбить на группы по ин тенсивностям принимаемых сигналов, то наибольшие уг ловые флуктуации будут характерны для группы изме рений с минимальной интенсивностью. Это обстоятель ство можно использовать для увеличения вероятности правильной пеленгации (в примере 2 решена задача об оптимальной обработке выборки независимых измерений плотности потока энергии).
В гл. 5 рассмотрен вопрос о влиянии порогового уров ня приемной аппаратуры на точность пеленгования. По казано, что уровни ограничения, малые по сравнению со средней интенсивностью поля, практически не изменя ют вероятности правильной пеленгации по сравнению с нулевым уровнем, а для больших уровней ограничения вероятность правильной пеленгации стремится к веро ятности обнаружения сигнала.
Все формулы, полученные в гл. 5, допускают двоякое толкование. Их можно рассматривать как статистиче ские характеристики плотности потока энергии, получен ные для фиксированного момента времени. Тогда в фор мулы следует подставить распределение участков ло кального отражения для рассматриваемого момента вре мени. Таким образом можно получить статистические характеристики тела сложной формы для определенного сектора углов наблюдения. Если же предполагать, что проводится усреднение по всем углам наблюдений, в ка честве распределения участков локального отражения следует брать соответствующее усредненное распреде ление.
Заметим, что согласованный характер поворотов участков локального отражения в отдельных реализа
182
циях не повлиял на распределение потоков энергии. Все полученные формулы в равной мере применимы и к со вокупности несвязанных отражающих элементов. Напри мер, 'неоднородности тропосферы, брызги вблизи поверх
ности МОрЯ И Т. Д . При ЭТОМ, ОДНа'КО, ДОЛЖНЫ (ВЫПОЛНЯТЬ
С Я другие ограничения, принятые в теории. В частности, угловые размеры отражающей области должны быть малыми.
6
КОРРЕЛЯЦИОННЫЙ АНАЛИЗ ФЛУКТУАЦИЙ ЭФФЕКТИВНОЙ ПЛОЩАДИ РАССЕЯНИЯ
ИУГЛОВЫХ КООРДИНАТ ТЕЛА СЛОЖНОЙ ФОРМЫ
ВПРИБЛИЖЕНИИ ДЕЛАНО
6.1.Влияние жесткой конструкции тела сложной формы на статистические свойства отраженного поля
Впредыдущей главе изучались вероятностные харак теристики флуктуаций интенсивности отраженного поля (ЭПР) и измеряемых угловых координат тела сложной формы. Полученные результаты позволяют анализиро вать независимые измерения этих параметров. Однако радиолокационные системы осуществляют непрерывные измерения параметров поля, так что отбор независимых измерений нельзя провести без корреляционного анализа процесса. Корреляционные характеристики ЭПР и изме ряемых угловых координат тела сложной формы пред ставляют и самостоятельный интерес: они содержат ин формацию о теле сложной формы, их необходимо учиты вать при конструировании приемной аппаратуры.
Основную теоретическую трудность представляет
анализ измеряемых угловых координат. Впервые расчет корреляционных функций для этих параметров был вы полнен Мачмором [76]. В этой работе модель Делано рассматривалась во времени, и для упрощения расчетов были использованы независимые усреднения по началь ным фазам поля, обусловленным расстояниями до участ ков локального отражения, и по траекториям движения
183
этих участков. Однако количественные соотношения, по лученные в [76], весьма приближенны, так как при вы полнении расчетов было сделано необоснованное предпо ложение о статистической независимости сигнала ошиб ки и амплитуды поля.
Ошибка, допущенная Мачмором, была устранена в работах [6, 16, 37, 77]. Однако в процессе'расчетов было использовано предположение о нормальных изменениях поля во времени*. Это предположение оправдано только для элементарных отражателей, совершающих независи мые хаотические движения (например, брызги вблизи поверхности моря). Поля, отраженные от жестко связан ных конструкций, какими являются тела сложной фор мы, таким свойством не обладают.
Корректный расчет корреляционных функций для ЭПР и измеряемых угловых координат тел сложной фор мы в приближении Делано был выполнен лишь недавно [23]. В основе этого расчета лежит двойственный под ход к временным флуктуациям поля, отраженного от тела сложной формы. С одной стороны, при малых разно сах во времени (или пространстве) взаимное расположе ние участков локального отражения предполагается не изменным. Это предположение обусловлено тем обстоя тельством, что линейные размеры участков локального отражения всегда существенно меньше линейных разме ров тела сложной формы. Но, с другой стороны, в соот ветствии с модельным приближением Делано при боль ших разносах во времени (или в пространстве) положе ния участков локального отражения существенно изме няются. Таким образом, можно считать, что имеется ансамбль тел сложной формы, каждое из которых полу чается случайным перемещением участков локального отражения, но внутри каждой из этих реализаций изме нения во времени происходят только в результате пово ротов всех участков локального отражения. Корреляци онные функции ЭПР и измеряемых угловых координат в [23] получаются как результат усреднения по описанно му выше ансамблю реализаций. Таким же методом про водится корреляционный анализ в шестой главе.
* Нормальность процесса не вытекает из нормального распре деления измерений через большие интервалы времени.
184
6.2. Корреляционные функции ЭПР и угловых координат тела сложной формы
Для определения корреляционных функций ЭПР и измеряемых угловых координат тела сложной формы
будет проведен |
анализ |
плотностей |
потоков |
энергии 5 |
|
в общих чертах |
так же, |
как |
в § 5.3 рассматривались |
||
вероятностные характеристики |
этого |
вектора. |
Согласно |
§ 5.3, вектор 5 удобно рассматривать в ортогональной системе координат X Y Z с началом в центре тяжести тела сложной формы, считая, что распределение масс в объеме пропорционально распределению участков локального отражения. Для определенности будем считать, что
излучаются флуктуации вектора S в плоскости X Z ,
определяемой единичными ортами ех, ег. Тогда, если
Si — плотность потока энергии в первый момент вре
мени в точке |
1, 5-2 — плотность |
потока |
энергии во вто |
|||||
рой |
момент |
времени в |
точке |
2, |
то |
задача |
состоит |
|
в корреляционном анализе двух |
двумерных |
векторов |
||||||
(5 .,, |
5 і2), |
(Soд>S 2z). |
что в |
каждой |
из реализации |
|||
Будем |
предполагать, |
случайного процесса движения наблюдателя и тела сложной формы происходят одинаково, т. е. статисти ческие усреднения производятся только по начальным фазам ф„, амплитудам ап полей, отраженных отдельными
участками тела сложной формы, и положениям |
участков |
||||||
локального отражения. |
Тогда, в полной аналогии с пре |
||||||
дыдущим, векторы |
^ |
и So выражаются через распре |
|||||
деленные по нормальному закону векторы Ап, |
АТ, |
А І , |
|||||
АТ |
и |
статистические |
свойства |
плотностей |
потоков |
||
-» |
—> |
полностью |
определяются |
вторыми |
моментами |
||
Si и S 2 |
|||||||
восьми случайных величин А?х, ATz, |
Д Д . 4 , |
Au-, |
АіД |
||||
Аох, |
A2z. |
|
|
|
|
|
Явные выражения А +, А~ через известные величины
1 8 5
приведены в предыдущей главе |
(формулы (5.6), |
(5.7)), |
вычисления моментов Л ц - Atz, Л а |
Atx и т . д . |
прово |
дятся аналогично расчетам § 5.3. Окончательные резуль
таты приведены в табл. 7. |
Момент Л а |
ЛоІ расположен |
|||||
на пересечении |
строки Л а |
и |
столбца |
А%х и т. |
д. |
Для |
|
краткости все моменты нормированы к величине |
Л а |
Ла |
|||||
и введены следующие сокращенные обозначения: |
|
|
|||||
ЧІ = |
5 |
(х ’ У’ |
z ) c o s 2 k ( % 2 |
— R\)dv, |
|
(6.1) |
|
Чп = |
5 |
3». |
z) sin 2k (R2— /?,) dv, |
|
(6.2) |
||
|
|
(0) = |
4-2 |
, |
|
|
|
где W(x, |
y, z) — вероятностная плотность распределения |
||||||||||
каждого |
из |
участков |
локального отражения в простран- |
||||||||
|
|
|
|
|
|
|
|
|
Т а б л и ц а |
7 |
|
|
Atz |
Atz |
А + |
At,- |
Atz |
Atz |
|
|
|
||
|
A |
\ X |
|
A |
t x |
||||||
Atz |
|
1 |
0 |
|
0 |
0 |
4о |
4 t |
4 t |
4 t |
|
Atz |
|
0 |
1 |
1 |
0 |
0 |
- 4 t |
4о |
- 4 t |
|
-1- |
|
4 |
1' |
|||||||||
|
|
|
|
1 |
|
|
|
|
|
|
|
At- |
|
0 |
0 |
4 t (° ) 0 |
4 t |
4 t |
4 І |
4t |
|||
A u |
|
0 |
0 |
|
0 |
4 t О) - 4 t |
41 - 4 t |
4 t |
|||
Atz |
4(t |
- 4 t |
4 t |
| - ? Г |
1 |
0 |
0 |
|
0 |
||
Z |
4 t |
4a |
4 t |
4І |
0 |
1 |
0 |
|
0 |
||
Atx |
4t |
- 4 t |
4 t |
—4 t |
0 |
0 |
4 t ( 0 |
) |
0 |
||
A u |
4 t |
4i |
4 t |
4t |
0 |
0 |
0 |
4-t(°) |
186
стве (соответственно, W {х, у, z)dv — с точностью до по стоянного множителя средняя энергия, рассеянная уча стками локального отражения в точку наблюдения из объема dv), R0— расстояние между точкой наблюдения и телом сложной формы, R і — расстояние между элемен том объема dv и первой точкой наблюдения в первый момент, R2— расстояние между элементами объема dv
и второй |
точкой наблюдения во |
второй момент времени. |
Функции |
R \{х, у, z ), R2(x, у, z ) |
определяются траекто |
рией движения наблюдателя и вращением многоэлемеигного источника.
Будем определять ориентацию вектора S xex + S zez
полярным углом 7), отсчитываемым от оси Z. Факти чески, в пределах точности расчетов параметр т| сов
падает с |
проекцией |
единичного |
вектора |
t, введенного |
||
в предыдущей главе, |
на ось X, |
т. е. равен измеряемой |
||||
угловой координате тела сложной формы в плоскости X Z . |
||||||
Однако |
здесь целесообразно упростить обозначения, |
|||||
так как |
рассматриваются |
флуктуации плотности |
потока |
|||
энергии |
только в одной |
|
|
|
5 |
|
плоскости. Так как т\ — -ф — |
||||||
|
|
|
|
|
|
Ьг |
малая величина, то, выражая 5 |
через А+ и А~, |
соглас |
||||
но (5.15), |
для корреляционной функции |
флуктуаций из |
меряемых угловых координат тела сложной формы имеем
= |
I |
S.n |
1 |
/ (Л ѵ И г1! + |
А ѵ И гі) (Лг2^г2 + |
А х 2A z 2 ) \ |
||||
1 |
Sz lS„ |
I |
Ч ^ г И г 'і |
+ |
Azl 4 і ) |
(A+2Az2 + |
Az2Az2)l |
|||
|
|
+ J+ |
4 |
|
|
|
Az2Az2) |
|
(6.3) |
|
|
|
( А 2ІЛ2І ' АГіАА) (Аг%Аz+2 |
|
|
||||||
где точками обозначены другие слагаемые. |
сводится |
|||||||||
Вычисление |
отдельных |
слагаемых |
в (6.3) |
|||||||
к восьмикратному интегрированию в бесконечных |
преде |
|||||||||
лах соответствующих функций, |
умноженных |
на |
плот |
ность вероятности W (Лц А2)- Однако нормальное рас
пределение W ( 4 , А 2) очень сложно выражается через матрицу моментов, поэтому предварительно выражают
187
W (Д , Д ) через характеристическую функцию нормаль ного распределения, а затем уже выполняют интегриро вание. Эти расчеты были проведены в [75]. В резуль
тате для |
получена формула |
|
|
||||
----- |
/('Ф 'Д |
+ |
ЧоЧі'У — {</Uo |
— <h Чб У |
|
||
|
----------- |
- |
Ы |
2]2 |
|
||
|
|
|
2 [ № |
|
|||
<1о |
1Н |
"Чог ЧчW |
l n |
D |
- |
(<7о ) 2 - ( < |
|
2 [ № Я+ М * ] |
|
|
|
|
|||
I_ |
_ |
|
(_v U t _+ |
ЧдЧѴ)'2 |
|
||
' |
[ |
' |
“ ( ^ r n u r f ) 3+ (■/«)’]' |
l6'4) |
Корреляция интенсивностей поля (величина, пропор циональная корреляции ЭПР) вычисляется значительно
проще. Согласно § 5.3, модуль вектора 5 приближенно
равен его г -- компоненте (угловые флуктуации S по рядка угловых размеров тела сложной формы), поэтому,
если I ѵ / 2 — интенсивности поля в точках 1 и 2 в пер вый и второй моменты времени соответственно, то кор реляционная функция интенсивности поля в точке наблю дения равна
j 1^2 = S z \ S ] (Д 1Д Г (Ді)2] [(До)2-j- (Дг)2]-
Но для нормально распределенных величин Д |
чет |
|
вертые моменты определяются вторыми: |
|
|
Д Д Д Д (ДД) (ДД) -)- |
|
|
+ (ДД) (ДД) + (ДД) (ДД), |
(6.5) |
|
а вторые моменты, необходимые для вычисления |
/ , / 2. |
|
приведены в табл. 7. В |
результате для корреляции ин |
|
тенсивностей получаем |
|
|
Л Л = [1 + |
(до г + ( Д ) 2] (If2. |
(6.6) |
Совершенно аналогично можно рассчитать и другие моменты величин S xl, S zl, S д.2, S z2. Все они достаточно
просто выражаются через дб, д,7- Например,
Д і Д 2 = ](<?Г)'2 + (дТ)2 + д і ?о+ + й2 gö] {^ ~ ■
188