
книги из ГПНТБ / Штагер, Е. А. Рассеяние волн на телах сложной формы
.pdfсобой. Значения моментов распределения ЭПР тел, представленные в табл. 1, могут быть использованы для определения коэффициентов вариации, асимметрии и эксцесса распределения ЭПР цилиндра, прямоуголь ной пластины и двугранного уголкового отражателя.. Указанные коэффициенты вариации а, асимметрии аі и эксцесса И2 по определению равны
Мо |
М. |
М, о |
(1.71) |
|
М%2 |
а2 — ~ГУ>— 3. |
|
|
м 2 |
|
В формулу (1.71) входят центральные моменты M s, ко торые легко выражаются через начальные моменты т/.
М 2= т2 — т2; М ъ.= /п3— 3тхт2 -J- 2т\\
М 4= т4 — 4т:]т1+ бт.2т\ — 3т\. |
(1.72) |
Значения коэффициентов вариации, асимметрии и экс цесса, рассчитанные по формуле (1.71) с учетом (1.72), представлены в табл. 2. Характер изменения коэффи-
|
|
|
Таблица |
2 |
|
О траж аю щ ее тело |
Коэффициент |
Коэф ф ициент |
Коэффициент |
||
вариации а |
асимметрии |
эксцесса |
а 5 |
||
|
Цилиндр (Л, |
R) |
Л / |
|
і / |
kfib0 |
3 |
khü0 |
У |
37t |
У |
2 |
-J- |
|||
|
|
|
|
||||
Пластина (а х |
h) |
1 / 4йЛ90 |
1 / |
khb0 |
3 |
kh\ |
|
У |
37t |
У |
2 |
— |
|||
Двугранный |
|
1 / 4khb„ |
1 / |
kh\ |
3 |
|
|
уголок (я X Л) |
У |
37t |
V |
2 |
- 4 - |
Ы Я , |
циентов таков, что ни одно из стандартных распреде лений, приведенных в Приложении, не может быть использовано для аппроксимации распределения ЭПР цилиндра, пластины и двугранного уголкового отража теля. Это распределение будет найдено в следующем параграфе.
39
1.8.Распределение вероятности эффективной площади
рассеяния тел простой формы
Наиболее исчерпывающей характеристикой отража тельной способности тела в секторе углов его наблю дения является интегральная F (а) или дифференциаль
ная |
функция |
распределения вероятности |
ЭПР этого |
тела. |
|
|
вероятности |
Интегральная функция распределения |
|||
ЗП Р |
тела по |
своему физическому смыслу |
должна об |
ладать свойствами непрерывности и дифференцируемо
сти. Производная этой функции =№ '(а) являет
ся дифференциальной функцией распределения вероят ности значений ЭПР тела пли плотностью распределе ния вероятности ЭПР. Задача аналитического опреде ления W (а) аналогична известной задаче теории веро ятности, когда ищется плотность распределения неслу чайной функции от случайных аргументов, плотность распределения которых задана [25]. Когда ЭПР тела о (А) является функцией одного угла наблюдения ö, в соответствии с общими правилами
W (а) = IV (9) |
1 |
(1.73) |
|
d- |
|||
|
|
Ж |
|
где 117(A)— плотность |
распределения вероятности угла |
||
наблюдения. |
справедливо, |
когда зависимость |
|
Выражение (1.73) |
ö от а является однозначной. Тогда, выполняя операцию дифференцирования и подставляя в правую часть фор мулы (1.73) значение fl, выраженное через а, получаем искомое выражение для W'(a). Нели зависимость fl от су неоднозначна, то в ней необходимо выделить ветви однозначности Ад и операцию дифференцирования про изводить по каждой ветви, а результаты просуммиро
вать: |
|
|
W(c) = 2 w (b „ ) _1_ |
(1.74) |
|
ft |
I d°k |
|
|
и»* |
|
Когда ЭПР тела зависит от двух углов и зависимости А и у от а неоднозначны, определение величины И7(ц)
40
связано с большими трудностями, которые удается пре одолеть с помощью расчетов на ЦВМ [83].
Рассмотрим распределение вероятности ЭПР цилинд ра, прямоугольной пластины и двугранного уголкового отражателя. В § 1.7 было показано, что зависимости
-_г
ІІУУ
. |
|
ш (УМ№ |
и д . |
ЯгlZ I |
|
||||
а п П IJ111 1 I II I I I 1 |
I II Т И У/ 1Л 1/ |
|
---- — |
|
|||||
О |
10° |
20° |
30° |
40° |
50° |
60° |
70° |
30° |
& |
Рис. 17. К определению распределения |
вероятности Э П Р тел про |
||||||||
|
|
|
стои |
формы. |
|
|
|
|
ЭПР этих тел от угла наблюдения приближенно опре деляются одним и тем же выражением
sin2 ( р !))
(1.75)
(/>»)* ’
где p = kli, причем h является длиной цилиндра, сторо ной пластины или ребром уголка в зависимости от того, какое из этих тел рассматривается. Зависимость у — = ср(€), представленная на рис. 17, показывает, что од ному значению у соответствует несколько углов По этой причине определение \Ѵ(а) следует производить по формуле (1.74).
К сожалению, зависимость Ф = <рі (у) не может быть записана в явном виде, что делает невозможным полу чение аналитического выражения для величины Щ а). Тем не менее можно выявить некоторые свойства этого распределения.
1. Как следует из графика функции у = о (4) (рис. 17); производная dy/db обращается в нуль всякий раз, когда khbnt = т (тг/2), где т = 1, 2, 3. . . (2т - 1). ' Следоваг тельно, в этих точках один или несколько членов сум мы (1.74) обращаются в бесконечность, если одновре менно с этим W (4) не обращается в нуль. В секторе, углов наблюдения функция ' W (9), как правило, всегда
т
больше |
нуля. Таким образом, плотность распределения |
||||
W (у) в |
точках |
уш= а /а макс = Sln2(/«n:/2)/(/«ir/2)2 (/»== |
|||
= 1, 2, |
3...) обращается в бесконечность. |
Нетрудно |
|||
показать, |
что эта особенность W (у) является |
интегри |
|||
руемой и в указанных точках интегральная функция |
|||||
вероятности ЭПР этих |
тел имеет конечное приращение. |
||||
2. |
Рассмотрим соотношение между плотностью рас |
||||
пределения |
вероятности и характеристической |
функцией |
|||
|
|
|
|
со |
|
|
|
W Ы |
= І |
$ ѳ fo) ехР (*Vo) du, |
(1.76) |
|
|
|
|
— со |
|
где 0о — величина ЭПР тела, нормированная к 1 м2. Поскольку функции W (оо) и Ѳ (г)) связаны между со бой Фурье-преобразованнем, значения аргументов (Доо и Дг)), при которых ординаты функций становятся пренебрежимо малыми, удовлетворяют известному соотношению ДооДц-^І [34]. Рассматриваемые тела (цилиндр, плоскость и двугранный уголок) имеют раз меры, значительно превосходящие длину волны поля, поэтому Д0О^>1, а следовательно, Дті<§;1 *. В этом слу чае представляется возможным разложить функцию Ѳ(г|) в ряд Тейлора около точки т]= 0. Выполнив ука занное разложение в формуле (1.76), получаем
со со
|
w |
К ) = і |
2 |
$ ( Ң У ^ -а ц . |
(1.77) |
|
|
|
rt =1 |
— сю |
|
ИзвестноT A |
, что |
ä "^~° (п’ l l |
при ц = |
ГU\ равны соответствующим |
|
моментам п-то порядка, поэтому |
|
||||
|
|
|
со |
со |
|
|
W (°о) = |
|
\ |
(1.78) |
|
|
|
|
П—1 — оо |
|
Величины моментов тп ЭПР рассматриваемых тел бы ли получены в предыдущем параграфе (см. формулу <1-68)). 3 аменяя в формулах (1.68) сомножитель, за-
* Неравенство справедливо при длине волны падающего поля свыше 1 мм, т. е. практически во. всем диапазоне используемых на
Практике радиоволн.
42
висящий от порядка момента, через ап, получаем общее выражение
тП |
а п (°о)макс |
(1.79) |
|
|
’ |
||
где (о„)макс — максимальное значение ЭГ1Р тела. |
|
||
Подставляя выражение |
(1.79) |
в формулу (1.78), имеем |
|
|
( \П |
“ |
|
U7 (а0) Ш 0= |
Т |
$ {i-nY^d-n. |
(1.80) |
П=1 |
-СЖ5 |
|
|
В правой части равенства |
(1.80) |
отсутствуют члены, за |
|
висящие от величины сектора усреднения т0о, |
поэтому |
||
произведение W (оо)'X.kh'&o (левая часть этого |
равенст |
ва) также не зависит от ФоЭто может быть лишь в том случае, если функция плотности вероятности 117(сго) за висит от величины, обратной сектору усреднения. Таким
же |
свойством будет |
обладать и интегральная функция |
F (у) распределения |
вероятности ЭПР рассматриваемых |
|
Т е л |
П р и </ = сг/ом акс- |
|
Определим пределы изменения у, в которых указан ное свойство выполняется. Проще всего определить до пустимые пределы изменения у путем расчета функции F{y), принимая величину k№0 за параметр. Результаты таких расчетов, заимствованные из работы [3], пред ставлены на рис. 18, откуда следует, что в области из-
г,%
Рис. 18. Интегральная функция распределения ЭПР тел простой формы (цилиндра, прямоугольной пластины или двугранного уголка) при разных электрических размерах характерной длины этих тел:
1 - Ш 0 = 1 0 ; 2 - *Л9о = 20; 3 - *Л90 = 40.
43
менепия у от 1 до 0,05 ординаты функции F (у) обратно пропорциональны величине сектора усреднения. Ука занная область приближенно соответствует области из менения ЭПР тел в пределах своего главного лепестка.
Таким образом, если сектор усреднения велик по сравнению с шириной главного лепестка диаграммы от ражения от рассматриваемых тел (kh'&0^> 1), то в пре делах изменения ЭПР этих тел от максимального зна чения до пяти сотых этой величины, ординаты распре деления вероятности значений ЭПР обратно пропорцио нальны величине сектора углов усреднения.
1.9.Корреляционная функция флуктуаций эффективной
площади рассеяния тела
В предыдущих параграфах рассматривались исклю чительно энергетические вероятностные характеристики величины ЭПР тел разной формы. Не менее важными являются характеристики, учитывающие взаимную связь между значениями ЭПР тела при разных углах наблю дения. Такая зависимость характеризуется корреляци онной функцией
Ф (»,!),) = |
^ ^ (а, - Ö ) |
(о, - о) W (в,в2) flf8,rf»2, |
(1.81) |
|
где оі |
и 02 — ЭПР тела соответственно при угле наблю |
|||
дения |
üi и |
Ö2; о — среднее значение ЭПР |
тела; |
|
W(ÖT, |
Ö2) ■— двумерная |
плотность вероятности |
углов |
наблюдения.
Полагая, как и ранее, что флуктуации ЭПР тела образуют стационарный и эргодичный процесс, перечис лим ряд свойств корреляционной функции Ф(Ф,, fl2), которые будут использованы в последующем.
1. Корреляционная функция не зависит от выбора начала и направления отсчета углов ■От и öo, а зависит лишь от их разности •От — fl-2 = ДФ.
2.Если Д0-»-0, то Ф(Д'б) стремится к дисперсии ЭПР тела.
3.Если Д0—>-оо, то Ф (ДЭ1) = 0 , что указывает на от сутствие статистической взаимосвязи между двумя вели чинами ЭПР тела, соответствующими достаточно удален ным друг от друга углам наблюдения.
44
Наряду с прямым методом вычисления корреляци онной функции по формуле (1.81) применяется также метод вычисления с помощью характеристической функ
ции [17, 25]. |
что функция |
у = а ( ‘&)/ат»,с |
может |
|
Предположим, |
||||
быть представлена контурным интегралом: |
|
|
||
|
|
|
|
(1.82) |
|
С |
|
|
|
Тогда, подставляя (1.82) в (1.81), получаем |
|
|
||
Ф (ДЯ) = ^ ^ |
Е (іи) Е (is) Ѳ (и, s, ДЯ) duds, |
(1.83) |
||
где Ѳ ( ') — двумерная характеристическая |
функция. |
|||
Практическое |
использование |
формулы |
(1.83) |
зави |
сит от возможностей определения двумерной характе ристической функции процесса изменения углов наблю дения. Наиболее просто функция Ѳ(и, s) выражается в случае нормальных флуктуаций утла гЯ с нулевым сред ним значением и дисперсией D2:
|
Ѳ(и, s) = exp |
— |
(и2+ s2+ 2Kbus) , (1.84) |
где |
— коэффициент |
корреляции углов наблюдения Я; |
и, s — аргументы двумерной характеристической функции. В этом случае расчеты корреляционной функции ЭПР тела по методу характеристической функции оказывают ся сравнительно несложными.
Рассмотренные выше методы определения корреля ционной функции по формулам (1.81) и (1.83) явля ются наиболее распространенными, но не единственны ми [17]. Выбор того или иного способа расчета корре ляционной функции обусловливается простотой дости жения результата. Аналитическое выражение корреля ционной функции ЭПР цилиндра и прямоугольной пла стины прямым методом получено в работе [51]. Ниже приведем решение этой задачи методом характеристи ческой функции.
Как уже |
отмечалось, амплитуда отраженного поля |
от цилиндра, |
прямоугольной пластины и двугранного |
уголкового отражателя в приближении физической оп тики определяется выражением
Уо ($) = V У(Я) = sin р§!рЪ. |
(1.85) |
45
Следуя методу характеристической функции, функцио нальную зависимость (1.85) необходимо представить интегралом вида (1.82). По аналогии с определением спектра прямоугольного импульса, запишем
Е(и) = к/р при —
Е (и) = 0 |
ц р и |и [> /? . |
(1.86) |
С учетом выражения (1.86) интеграл (1.82) примет вид
р |
|
>’o(ö) = ^ r ^ - y exP (^ « )rf«- |
С1-87) |
- р |
|
Полагаем, что углы наблюдения тела распределены по нормальному закону. В этом случае двумерная ха рактеристическая функция определяется выражением (1.84). Подставляя выражения (1.84) и (1.87) в формулу (1.83), получаем
р р
ф ( К ®2) = і |
ехр [— 1 Г (й2 + ^ + |
duds. |
|
||
|
- р - р |
( 1. 88) |
|
|
Интеграл (1.88) достаточно просто вычисляется, если воспользоваться разложением
ехр ( - Щ-u sK ^) = 2 |
a"s"- |
-89) |
Тогда
ф(&і>82) = 4^ 2 ( - 1)" w i [SK"exp ( - ^ r )rfM]-
п = 0 |
—р |
(1.90)
Рассматриваемые в задаче тела по предположению об ладают характерными размерами, значительно превос ходящими длину волны поля, поэтому р^>1. В этих условиях пределы интегрирования в (1.90) могут быть заменены на бесконечные и выражение, заключенное в квадратные скобки, оказывается пропорциональным мо менту л-го порядка нормального распределения.
46
Произведя эти замены, окончательно получаем
KlmО'“ л |
(1.91) |
|
2n+'n\D2n |
||
|
п =О
Для нормального распределения первые четыре момента (см. Приложение) соответственно равны:
т0= 1; т х= 0; т2= D--, тл = 8£И.
Подставляя значения этих моментов в формулу (1.91), получаем приближенное выражение для искомой корре ляционной функции в виде
ф (аі> — и | 2£>Ѵ2 + 16DV ^ |
(1-92) |
В заключение отметим, что для рассмотренной за дачи метод характеристической функции оказался бо лее удобным по сравнению с прямым методом.
1.10. Критерий дальней зоны при измерении средней эффективной площади рассеяния линейного отражателя
При определении ЭПР тела исходят из «бесконечно большого» расстояния R0 между приемо-передающей антенной и отражающим телом. Именно поэтому вели чина ЭПР тела оказывается не зависящей от расстоя ния Ro. Однако измерения ЭПР производят, как пра вило, на конечном удалении от приемо-передающей ан тенны, т. е. вместо ЭПР тела измеряют величину
0R |
(/отр/^пад)» |
(1.93) |
а вместо средней ЭПР тела — величину
^ ( ® ) ° * (*)<**• |
(1.94) |
Обозначения, принятые в формулах (1.93) и (1.94), сов падают с обозначениями в формулах (1.6) и- (1.16) соответственно.
Различия между значениями ия и а оказываются не значительными, если расстояние между телом и приемо передающей антенной выбрано согласно неравенству
ЯМНН> А 2А, |
(1.95) |
47
где h — наибольший поперечный (к направлению рас пространения поля) размер тела, X— длина волны.
Неравенство (1.95) часто называют критерием даль ней зоны * для линейного отражателя (цилиндра дли ной !г, прямоугольной пластины со стороной h и др.)> эффективная площадь рассеяния которого определяется выражением (1.75).
При измерении средней ЭПР тела на конечном уда лении от приемо-передающей антенны также целесооб разно каким-либо образом определить минимальное рас стояние, при котором измеряемая величина будет мало отличаться от средней ЭПР тела, определенной при бес конечном расстоянии. Задачу рассмотрим применительно к линейному отражателю. Для этой цели сначала опре делим ЭПР линейного отражателя при расположении точечной приемо-передающей антенны на расстоянии R0 от него. Процедура указанного расчета полностью эк вивалентна расчету диаграммы направленности прямо угольной апертуры с квадратичным распределением фа зы по ее стороне, изложенному А. 3. Фрадиным в книге [57]. Повторяя этот расчет с небольшими видоизмене ниями, находим
o* = A - ^ { [ C ( n + ІО+ С (ц — 0] +
+ [5 (ң. -М) + |
5 (ц - £)]}, |
(1.96) |
где А — некоторая размерная |
постоянная, |
которая, на |
пример, для кругового цилиндра равна kah2 (а — радиус цилиндра, h — его длина, k — волновое число); ц =
= h V k!Rbcos Я; \ = \r kR ütg 9; С(-) и 5 (■) — соответ ственно косинус-интеграл и синус-интеграл Френеля.
Подставив выражение (1.96) в формулу (1.94), полу чим исходное выражение для определения средней ЭПР линейного отражателя в сферическом падающем поле. Вычисления по этой формуле для цилиндра выполнены с помощью UBM. При этом предполагалось, что углы наблюдения распределены равномерно, т. е. W (й) = = 1/2’é’o, а сектор усреднения 2б0 для цилиндров разной электрической длины охватывал помимо главного ле
* Вывод этого критерия будет дан в § 2.8.
48