Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Штагер, Е. А. Рассеяние волн на телах сложной формы

.pdf
Скачиваний:
58
Добавлен:
19.10.2023
Размер:
7.41 Mб
Скачать

Здесь

_i

 

/ = / ( ? ) = 4 -

cos2<p, (5.57)

ij.

 

 

р, cp — полярные координаты радиус-вектора, проведен­ ного из начала координат до границы области т (кх, X),

полученной проектированием моделирующего объема V на плоскость переменных Xv, X .

Иногда представляет интерес статистический анализ не всех измерений, а только части их, отобранной из всего ансамбля по определенному правилу. В данном случае важно выяснить, как изменяются флуктуации уг­ ловых координат в наборе измерений, полученном из полного ансамбля в результате отбора измерений с нор­ мированной интенсивностью, большей q. Определим вероятность правильной пеленгации в этом подансамбле.

Для решения поставленной задачи сформулируем не­ сколько иначе уже полученный результат для вероятно­ сти правильной пеленгации в полном ансамбле измере­ ний. Согласно теореме Байеса, вероятность правильной пеленгации в полном ансамбле измерений равна произ­ ведению двух вероятностей, одна из которых является вероятностью попадания случайного измерения в рас­ сматриваемый подансамбль, а вторая является вероят­ ностью правильной пеленгации в выделенном подансамб­ ле измерений. Интегрируя W (/0) по /о в пределах от q до бесконечности, находим, что вероятность попадания случайного измерения в подансамбль І0> q равна е_(С Следовательно, полная вероятность правильной пеленга­ ции равна

Q - e-'/Qj,

где Qi — вероятность правильной пеленгации в подан­ самбле реализаций, удовлетворяющем условию f 0^>qr Сравнивая этот результат с (5.56), приходим к выводу,, что

e-9/Р’

(5.58)

Qi = 1

d<o.

/ ( 1

+ /Р=)

 

Из (5.56) и (5.58) непосредственно видно,,что функ­ ции Q и Qi не зависят от ■О*, Фу. Это обстоятельство по­

169

зволяет доказать следующие общие свойства функций

Q и Qi.

1. Вероятность правильной пеленгации тела сложной формы инвариантна относительно линейных преобразо­

ваний

его формы вдоль осей X, У, Z (в соответствии

с § 5.3

предполагается, что ось Z проходит через центр

тяжести моделирующего объема V, а ориентация оси X выбрана из условия (5.12)).

Действительно, величины Q и Qі зависят не от вида области V в системе координат XYZ, а от вида области V в координатах Kx,Kv,kz. Но, согласно (5.46), при изме­ нении масштабов вдоль осей X, Y, Z координаты произ­ вольной точки пространства в переменных МДуДг не из­ меняются.

Как частный вывод из доказанного утверждения следует независимость функций Q и Qi от линейных размеров тела сложной формы.

Отмеченные свойства функций Q, Qі существенно упрощают анализ конкретных задач. Например, если рассматривается вероятность правильной пеленгации те­ ла сложной формы, причем моделирующий объем имеет форму эллипсоида, и наблюдатель расположен на одной из осей симметрии, не делая никаких выкладок, можно заключить, что размеры эллипсоида на величину Q не влияют. Действительно, линейным преобразованием двух координатных осей эллипсоид преобразуется в шар, а радиус шара на величину Q не влияет, так как Q не зависит от линейных размеров тела сложной формы.

Из (5.56), (5.58) также вытекает следующее свой­ ство.

2. Чем больше интенсивность плотности потока энер­ гии, тем больше вероятность правильной пеленгации.

Действительно, если принят сигнал интенсивности q, то можно утверждать, что имеется реализация из подан­ самбля, удовлетворяющего условию h ^ q . Но, как видно из (5.58), вероятность правильной пеленгации в таком подансамбле стремится к единице с ростом q. Прибли­ жение к предельному значению монотонное.

3. Чем выше пороговый уровень q, начиная с которо­ то приемная система фиксирует сигналы, тем меньше ве­ роятность правильной пеленгации. Максимальная веро­

170

ятность правильной леленгции достигается при q — О, когда измеряются все потоки энергии.

Для доказательства продифференцируем (5. 56) по q. Тогда получим

'1 e-Q.fp’

^ = е-<

— L

С

/

- d v — 1

(5.59)

дд

2л;хд.;^ѵ

J

 

 

Выполняя интегрирование, можно показать, что

2я

ипоскольку функции / и р положительно определены, первый член в квадратных скобках формулы (5.49) меньше или равен единице. Таким образом,

что и требовалось доказать.

е. при ма­

При 9 = 0,

согласно

(5.59), dQ/dq = 0, т.

лых пороговых

уровнях

q отклонение Q от

максималь­

ного значения — величина

второго порядка малости.

Уменьшение величины Q с ростом q не противоречит отмечавшемуся ранее увеличению вероятности правиль­ ной пеленгации для потоков большой интенсивности. Действительно, хотя с увеличением уровня ограничения принимаемые потоки энергии реже отклоняются от на­ правления на тело сложной формы, одновременно уве­ личивается и число случаев, когда антенна вообще «не видит» объекта. Монотонное уменьшение Q означает, что второй эффект всегда сильнее первого. При доста­ точно больших значениях q имеем:

Qj-M, Q-. е-?,

(5.60)

т. е. вероятность правильном пеленгации стремится к вероятности обнаружения сигнала.

В заключение приведем примеры конкретных расче­ тов на основе общих соотношений, полученных в гл. 5.

171

5.6.Примеры расчетов

Пример 1. Пеленгация объектов с простой формой моделирующего объема

Проведем расчеты для случаев, когда моделирующий объем V имеет форму эллипсоида, прямоугольника или цилиндра. Будем предполагать, что наблюдатель находится на одной из осей сим­ метрии пеленгуемого тела, так что условие (5.12) выполняется при произвольной ориентации осей X, Y. Однако, чтобы упростить ра­ счеты, ориентацию этих осей также будем выбирать вдоль осей симметрии моделирующего объема V.

Наиболее просто вычисляется вероятность правильной пеленга­

ции Q (q)

в случае,

когда

моделирующий объем — тело вращения

(ось симметрии совпадает с осью Z). Из

(5.57) имеем

 

 

 

 

 

Р= =

1.

Щу = ,'4 =

и2.

/ (?)

= І/и-2-

 

 

 

Подынтегральная

функция

в (5.56)

от

переменной

интегрирования

не зависит, и для

Q (q)

имеем

 

 

 

 

 

 

 

 

 

Q (И =

 

 

 

 

 

 

 

(5.61)

В качестве иллюстрации

проведем

вычисление

;і~2 для

шара.

В системе

координат

ХЛ., Ху,

\ z его

радиус равен единице, и

 

 

^

 

 

 

^ rfSt ^ (г sin !) cos ?)V2 sin SW/-

"ПР

 

Объем

шара в системе

координат Хл., Xv, X,

равняется

4г./3,

и, со­

гласно

(5.47).

 

 

 

 

Их

 

 

 

 

 

 

 

 

 

 

 

 

 

=

5.

 

 

 

(5.62)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Hydv}.

 

 

 

 

 

Если в шаре

имеется

шаровая

полость

радиуса г0 (/•„ <

1), то

вместо (5.62) получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(5.63)

Когда

источники

распределены

в узком

слое вблизи

поверхно­

сти шара, то го-*-1, р - 2-*-3. Для моделирующего объема, имеющего іформу цилиндра (образующая вдоль оси Z), аналогичные расчеты дают значение р—2 = 6.

Согласно (5.63), полость внутри моделирующего объема при­ водит к уменьшению вероятности правильной пеленгации.

■Формулы (5.61), (5.63), как отмечалось в предыдущем парагра­ фе, одновременно описывают вероятность правильной пеленгации

.для наблюдателя, расположенного на оси симметрии моделирующего

172

объема, который имеет форму эллипсоида с полостью, подобной его поверхности и расположенной симметрично относительно глав­ ных осей.

Если область т(л.ѵ, Л„) имеет вид квадрата со сторонами, парал­ лельными координатным осям, интеграл в (5.56) выражается через

элементарные функции лишь при q = 0.

В этом случае

 

Q (0) — —

arc lg

V-x

 

VyV\ V-x

 

/ і

V ,

 

1

 

(5.64)

. arc lg

 

V \ + M-y

Klv-x- УV 11 4+- Kyі4 J.

 

Формула (5.64) позволяет рассчитывать максимальные вероят­ ности правильной пеленгации тел сложной формы, которые модели­ руются совокупностью точек локального отражения, распределен­ ных в объеме, имеющем форму цилиндра или параллелепипеда. Предполагая, что ость X ориентирована параллельно образующей цилиндра, согласно (5.47) имеем

Р-72 = 3,

[I"2 - 4 (1 + /■ * ),

(5.65)

а для параллелепипеда с гранями, параллельными координатным плоскостям,

Ѵ-х

= 3

 

6 •

( 5 . 66)

 

 

1—го

параметр г0 показы­

В формулах (5.65), (5.66),

как и в

(5.63),

вает, во сколько раз внутренняя полость меньше всего моделирую­ щего объема тела сложной формы по линейным размерам.

Значения Q(0), которые получаются из формул (5.61) н- (5.66), приведены в табл. 6 в порядке возрастания. Для четырех тел слож­ ной формы, у которых область т(А.*, Ху) имеет вид квадрата, функции Qi (г/) рассчитывались на ЦВМ по формулам (5.56), (5.57), (5.66). Соответствующие кривые приведены на рис. 41. Номера около кривых обозначают номера в первом столбце табл. 6. Как видно из рисунка, кривые не пересекаются, т. е. расположение пеленгуе­ мых объектов в порядке возрастания Q(0), приведенное в таблице, сохраняется и при q ф 0.

Для иллюстрации общих свойств функции Q(q), по графикам функций Q[(q), приведенным на рис. 41, были рассчитаны две кри­ вые Q(q). Результаты расчетов приведены на рис. 42. Обозначения кривых те же, пунктиром проведена кривая ехр(—q). Как видно из рис. 42, в полном соответствии с выводами § 5.5, функции Q{q) (кривые I и 6) монотонно убывают, стремясь к ехр(—q), производ­ ная от Q(q) обращается в нуль при q = 0.

Когда наблюдатель не расположен на оси симметрии моделирую­ щего объема, вычисления значительно усложняются. Ограничимся перечислением основных результатов, полученных в [44].

Тело сложной формы, для которого моделирующий объем имеет форму эллипсоида, оказывается в некотором смысле особым отра­ жающим объектом. Для него вероятность правильной пеленгации

173

Т а б л и ц а 6

В и д м о д е л и р у ю щ е г о

п/п

о б ъ е м а

1Полый параллелепипед

2Полый цилиндр

3Полый эллипсоид

4Сплошной параллелепипед

5Сплошной цилиндр

6Сплошной цилиндр

7Сплошной эллипсоид

не зависит от ориентации относительно определяется формулами (5.61), (5.63).

В и д о б л а с т и

о

т (\ѵ S )

68,98

74,75

О75,00

78,70

о80,00

80,83

о83,33

наблюдателя, т. е. всегда

Рис. 41. Зависимость

вероятности

Рис. 42. Зависимость вероят-

правильной пеленгации

при иитен-

ности правильной пеленгации Q

снвных отражениях Q, от порого-

от порогового уровня д для

вого уровня д для

симметрично

симметрично ориентированных

ориентированных тел

сложной

тел сложной формы,

формы.

 

 

 

Результаты расчетов величины Q(0) для моделирующего объема, имеющего форму параллелепипеда, приведены на рис. 43. Геометрия задачи следующая: наблюдатель находится в плоскости двух осей симметрии; при О- = 0 координатные оси XYZ совпадают с осями

174

Рмс. 43. Зависимость максимальной

вероятности правильной пелен­

гации 0(0) от

ориентации тела

сложной формы, для которого

моделирующий объем имеет форму параллелепипеда:

1-

5 „ = 0 , 2 ; 2- £0 = 0 , 4 ; 3- £ 0 = 0 , 6 ; 4- Е0 = 1.

Рис. 44. Зависимость максимальной вероятности правильной пелен­ гации 0 (0) от ориентации тела сложной формы, для которого моделирующий объем имеет форму цилиндра:

, I - $о = 0,1; 2 —£о = 1; 3 - So = 10.

симметрии параллелепипеда

(наблюдатель на ccnZ); g0— отношение

Lz к Lx; угол

описывает

поворот

параллелепипеда вокруг осп У.

Очевидно, что если ö изменяется и

интервале [0, я/2], то для go

можно ограничиться интервалом [О,

I]. Максимум Q(0) достигается,

когда наблюдатель находится на продолжении диагонали прямо­ угольника, полученного проектированием моделирующего объема на

плоскость XZ.

максимальной

вероятности

правильной

Результаты расчетов

пеленгации Q(0) для цилиндра приведены на рис. 44.

Здесь

ö —

угол между осью Z и осью цилиндра,

go — отношение

высоты

ци­

линдра к диаметру.

 

функций Q\(q) быстро уменьша­

С ростом q пределы изменений

ются. Это обстоятельство

является

следствием отмеченного в § 5.5

Рис. 45. Зависимость вероятности правильной пеленгации при

интенсивных отражениях

от ориентации

тела сложной формы,

для которого моделирующий

объем

имеет

форму

параллелепи­

/ - 5 О = о , '2; 2 - і 0 =

педа:

і0=

 

 

 

0 ,4 ; 3

0 , 6 ;

• / - £ „ =

11.

общего свойства функции

Q\(q): сростом

г/ величина

Qi (q) стремит­

ся к единице, независимо от формы п ориентации тела сложной

формы. Д ія

иллюстрации на рис. 45

приведены результаты

расчетов

вероятности

правильной пеленгации

при интенсивных отражениях

Qi (<7) для

совокупности участков локального отражения,

распреде­

ленных в

моделирующем объеме в виде параллелепипеда

(q = 0,5,

геометрия

та же, что на рпс. 43).

 

 

Проведенный анализ позволяет сделать вывод, что вероятность правильной пеленгации — весьма устойчивая характеристика тела сложной формы. Несмотря на различия рассмотренных форм, вели­ чина Q(0) изменялась в пределах 69 -н 85%, повороты многоэле­ ментного источника не приводят к изменениям Q(0), большим 6%.

Пример 2. Случай нескольких независимых измерений

Согласно (5.49), потоки энергии флуктуируют около направле­ ния на центр тяжести тела сложной формы. Поэтому можно значи­ тельно увеличить точность пеленгации, если использовать набор не­

зависимых измерений плотностей потока энергии S|,. S2, ..., S m . При

176

пеленгации подвижных объектов такой набор можно получить, про­

водя

измерения через достаточно большие интервалы

времени

[66,

76], для неподвижных объектов — модулированием

частоты

[62, 69] или за счет пространственного разноса приемников. Полу­ ченное в § 5.5 вероятностное распределение плотностей потоков энергии позволяет количественно оценить различные методы обра­

ботки выборки 5|, S2, ..., S м.

Ограничимся двумерной задачей. Будем считать, что участки локального отражения статистически равномерно распределены на отрезке, который расположен в плоскости XZ и имеет угловой раз­

мер 2fl-.v. Тогда, согласно (5.48), коэффициент формы ,ч”."2 = 3, и рас­

пределение разности угловых координат і.х Тх имеет вид (5.52). Действительно, при выводе (5.52) предполагалось, что направле­ ние на центр тела сложной формы известно, и относительно этого направления отсчитывались углы tx. Теперь это направление предпо- —>

латается неизвестным, а іх —это полярная координата вектора S при произвольной ориентации полярной оси.

Введем безразмерный полярный угол

ѵ =

) "З:

(5.67)

Тогда, обозначая, как обычно, среднее значение ч через ч, для распределения случайных параметров ч и /„ из (5.50) имеем

_ і_

 

«^ (■'./„) = ( ^ ) 2 exp { _ / 0 [1 -I- (-----v)2j}.

(5.68)

Задача об определении оценки, обеспечивающей максимум функции правдоподобия для выборки измерений ч, / 0, согласно [8], сводится к определению экстремума функции правдоподобия по параметру

ч. В данном случае, так как измерения независимые, функция прав­ доподобия имеет вид

м

_і_

 

(5.69)

П W(ч,„ І оп) = П ( if) 2 exp { - I on [1 + (ч„ - 7 ) 2]}.

Максимум этой

функции

по

параметру ч достигается, когда

сумма

 

 

м

 

 

 

 

2

/ 0« ( ѵ „ - Д 2

(5.70)

 

 

П=1

 

 

минимальна. Но (5.70) представляет собой полином второй степени

по ‘Л следовательно, функция

(5.70)

имеет один

минимум, в кото­

ром производная обращается

в

нуль. Дифференцируя (5.70)

по ѵ

и приравнивая результат нулю,

для

оптимальной

оценки ѵ0

полу­

чаем

 

 

 

 

 

1 2 З а к а з Js2 І б б

177

V о

2 / олѵл

(5.71)

где суммирование по индексу п проводится от 1 до М.

Точность определения ѵ по формуле (5.71) наиболее подробно описывается распределением отклонении ѵ0 от ѵ. Пусть

2 /o«(v« —ѵ)

--V=

(5.72)

2

Тогда вычисление плотности распределения И7 (Д0) удобно прово­

дить с

помощью характеристической

функции.

Согласно (5.68)

и (5.72),

 

 

 

 

 

 

+00

 

 

 

 

 

 

DO

СО

 

 

Ѳ(в) =

е‘“іо = Srf/

 

 

 

d-!M X

 

 

 

 

0

 

 

 

 

 

 

 

 

 

I

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

1/ 01

7оAl

т

 

Y i o n ^ 2 ^оп(ѵл

о*

 

 

b r -

71

exp

 

 

 

 

 

 

 

 

 

 

 

 

 

2 7 О Л

( ѵ я

 

 

 

 

 

 

 

 

Ъ ‘оп

 

 

 

 

Выполняя интегрирование

по -ц,

получаем

 

 

с о

 

СО

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dloM-

Для плотности распределения случайных отклонений Д0

соответ­

ственно

имеем

 

 

 

 

 

 

 

 

 

 

W( Д0)

 

 

'°da =

 

 

 

5 d l 01. .. ^ ( 2

70л і) 2 ехР £ — О +

до)

2 ,оп

dlом-

оо

Перейдем к переменной интегрирования

х п = (і + ^o) 7on-

178

Соседние файлы в папке книги из ГПНТБ