
Литература / М.И.Финкельштейн Основы радиолокации
.pdf(3.4.13) изменяются (рис. 3.15). Изменение фазы свидетель ствует о том, что эквифазная поверхность (т. е. фронт волны) искривлена по сравнению со сферическим фронтом одно точечного источника, расположенного в центре цели (точ ка 0). Вместе с тем при изменении направления на цель оп-
|
ределяется |
направление |
||||
|
нормали |
к фронту |
волны. |
|||
|
Следовательно, |
в |
данном |
|||
|
случае возникают |
ошибки |
||||
|
(так |
называемые |
угловые |
|||
|
шумы). |
|
|
|
центр |
|
|
3. |
Кажущийся |
||||
|
отражения. Фаза отражен |
|||||
|
ного сигнала от двухточеч |
|||||
|
ной цели равна 4л£)А+ф, |
|||||
|
гдеф определяется по фор |
|||||
|
муле (3.4.13). В то же вре |
|||||
Рис. 3.15. Зависимость амплиту |
мя для одноточечной цели, |
|||||
находящейся |
на |
расстоя |
||||
ды и фазы сигнала, отраженного |
нии |
D, |
эта |
фаза |
равна |
|
от двухточечной цели, от направ |
4л£)/Х. Двухточечную цель |
|||||
ления при D~ const |
можно поэтому рассматри вать как эквивалентную одноточечную, расположенную на расстоянии D 4- А£), где AD определяется по формуле
АГ» |
А. |
% |
, |
^2 |
, / 2л |
, . л' |
AZ) |
= — ф — — arctg |
е2+ег |
tg I---- |
L sin 0 |
||
|
4л |
4л |
|
\ |
|
Такая эквивалентная одноточечная цель определяет ка жущийся центр отражения. При изменении соотношений ме жду Ег и Е2 величина AZ? соответственно изменяется. В част ности, при изменении от Et ~ 0 до £2 — 0 величина А£) принимает значения от L sin 0/2 до — L sin 0/2, т. е. кажу щийся центр отражения перемещается вдоль линии Цъ Ц2 (рис. 3.13, а).
На рис. 3.16 показана двухточечная цель при малом угле 0, когда < ОЦхЦ ~ 90°. Пунктиром показана линия рав-
Рис. 3.16. К выводу по ложения кажущегося центра отражения
150
ных фаз одноточечной цели, соответствующей середине двух точечной, а сплошной линией двухточечной цели. При этом кажущийся центр отражения Ц лежит на нормали к линии равных фаз и смещен относительно середины двухточечной цели на ДА. Для определения ДА воспользуемся уравне
нием |
нормали |
в |
полярной |
системе |
координат tg Д0 = |
|||||
= (dD/dQ)/D (здесь D играет роль |
радиуса-вектора), |
от |
||||||||
сюда смещение центра |
ДА — D tg ДО = dDidft. Дифферен |
|||||||||
цируя |
(3.4.14), |
находим |
|
|
|
|
|
|
||
|
|
|
|
ДА = б/П/б/0 = |
|
|
|
|||
|
£ |
|
|
(1— /п)/(1Ч-/п) |
■ X |
|
||||
|
2 |
|
, |
|
|
|
|
’ |
|
|
|
|
|
|
|
/2л |
\ |
|
|||
|
|
1 + 1(1-/«)/(!+/n)]2tg2 |
— A Sin е) |
|
||||||
|
|
|
|
|
cos 0 |
\ Л |
/ |
|
||
|
|
|
X |
~ |
|
• |
|
|
||
|
|
|
|
|
\ |
|
|
|||
|
|
|
|
/2л |
|
|
|
|||
|
|
|
|
cos2 |
——Lsin 0 |
|
|
|
||
|
|
|
|
\ |
X |
|
/ |
|
|
|
При sin 0 = 0; |
7./2А; Х/А, ... |
(максимум ЭОП двухточечной |
||||||||
цели) |
|
|
|
___ L |
|
__ |
|
|
|
|
|
|
АГ |
|
|
УОщ |
|
Оц2 |
|
||
|
|
|
|
|
|
УаЦ1 + У°Ц2 |
|
|||
Максимальное смещение соответствует нулям ДОР двух |
||||||||||
точечной цели |
[см. |
(3.4.7)], когда |
sin 0 = Х/4А; ЗХ/4А, |
..., |
||||||
и равно |
У |
|
У°Ц2 |
п П |
|
|
+~У°Ц2 |
|
||
АГ |
___ _£. |
|
|
|
|
|||||
|
тах~ 2 |
У^-У^ |
|
~ 2 |
У^~У^2 ’ |
|
т. е. ДАтах > А/2, что означает смещение центра отраже ния за геометрические размеры двухточечной цели (случай ац1 = сц2 надо исключить, так как при этом ДА == 0).
Легко видеть, что и для направлений, расположенных между максимумами и нулями ДОР, возможно соотношение ДАтах > А/2. Например, при т = Уац1/ац2 = 2 и sin 0 = = А/5А имеем ДА = ЗА/2.
Флуктуация кажущегося центра отражения приводит к погрешностям измерения дальности и угловых координат. При этом можно считать, что погрешность измерения угло вых координат вызывается флуктуациями фронта отражен ной волны.
В заключение отметим, что доплеровский сдвиг частоты, который пропорционален производной фазе, для случая сложной цели так же, как угловые координаты и дальность, является флуктуирующей величиной.
151
4. Скорость флуктуаций ЭОП. Для решения практиче ских задач требуется знание основных параметров флукту аций, возникающих при движении сложных целей. Измене ние взаимного расположения блестящих точек целей и РЛС приводит к сложной интерференционной картине, которая может быть описана посредством многолепестковой переме щающейся ДОР. С другой стороны, из-за изменения расстоя ния между блестящими точками и РЛС возникают вторичные доплеровские биения, которые также объясняют наличие флуктуаций. В зависимости от решаемой задачи использует ся тот или иной метод описания данного явления. Проиллю стрируем сказанное на примере определения частоты флук туаций отраженного сигнала двухточечной цели.
Разность фаз колебаний отдельных точечных отражате лей согласно (3.4.5) <p1>2 ~ 4nL sin 0А. В процессе движе ния двухточечной цели происходят случайные изменения угла 0 (/) , что приводит к флуктуациям, т. е. случайным изменениям результирующей амплитуды £р [формула (3.4.4)] от (Ег + Е2) до |£j — £21Для определения мгно венной частоты флуктуаций достаточно определить абсо лютное значение производной фазы:
р |
_ |
1 I |
^ф1,2 |
I ___ |
1 |
I ^Ф1.2 |
I I ^0 I ___ |
фл ~~ |
I |
dt |
I “ ~2л" I |
I I ~dT I ~ |
|||
|
|
|
= -TLl^r|lcosel- |
(3.4.15) |
|||
|
|
|
|
I dt |
I |
|
|
Наиболее |
интересен |
случай |
0 « 0, |
когда |
<ЗЛ1б>
Полученный результат можно легко найти на основе по нятия ДОР. На рис. 3.17, а показана точечная цель, дви жущаяся прямолинейно. Из-за смещения цели относитель но РЛС меняется ЭОП, что и приводит к флуктуациям отра женного сигнала. Пусть, например, отражающие свойства самолета соответствуют ДОР некоторой эквивалентной пла стины длиной а — L. Тогда согласно (3.2.26) ширина боко вых лепестков A0fe = X/2L. Так как движение цели отно сительно РЛС можно характеризовать ее мгновенным вра щением с угловой скоростью dftldt, то частота флуктуаций
^НЧг1/Л0‘=1Чг12£Л’
I dt 11 I dt I что соответствует формуле (3.4.16).
152
Такой же результат получается и вследствие вторично го эффекта Доплера. На рис. 3.17, б показан общий случай движения двухточечной цели относительно РЛС. Радиаль ные скорости целей обозначены ир1 и ир2. Движение отдель ных целей можно заменить поступательным движением цент ра системы из двух целей со средней скоростью (ир1 4- vp2)/2
Vpf - vp2
2
VP1+VPL \u.
2
Vpl-VpZ _
z
в)
Рис. 3.17. К объяснению флуктуаций отраженного сигнала при дви жении цели
и дополнительным движением целей со скоростями (ир1 —
— ир2)/2 в противоположных направлениях (рис. 3.17, в), т. е. мгновенным вращением системы вокруг ее центра с уг ловой скоростью
I |
dft |
I _ (оР1—^рг)/2 |
Opt — РРа |
' |
(3 4 |
) |
|
I |
dt |
I |
L cos 0/2 ~ |
L |
\ |
||
Различие радиальных скоростей ир1 и vp2 приводит к |
|||||||
тому, что доплеровские сдвиги FД1 |
и Гд2 |
также будут раз |
ными. Результирующие колебания, возникающие в антенне РЛС при смешении сигналов, отраженных от каждой из Целей, представляют собой вторичные доплеровские биения с частотой флуктуаций, равной [см. (3.4.17)1
Гфл = ДГдв=Гд1-/?дг«^1-^-|,
Л, I al I
т. е. снова получена формула (3.4.16).
Проведем теперь анализ полученного результата. Вер немся к основной формуле (3.4.15). Обычно dftldt и 0, а зна чит, и ^фл — весьма медленно меняющиеся случайные функции времени, так что можно принять d&ldt ~ const и cos 0 const. Отсюда dqrJdt ~ 2лГфЛ (0 ~ const и <р12 т ъ 2лГфЛ (0 t.
153
При условии Ег ф Е2 достаточно, например, выполне ния неравенства Е2< ®,7Elt чтобы в формуле (3.4.4) мож но было пренебречь величиной по сравнению с Е1 и да
лее |
(на основании разложения в |
ряд ]/1 + х « 1 4- х/2 |
|
при |
^«1) получить |
|
|
|
Ер(0«£1[1 + Acos2n^(()t 1. |
(3.4.18) |
|
|
L |
J |
|
т. е. отраженный сигнал модулирован медленно меняющей ся случайной функцией времени.
Заметим, что поверхность такой сложной цели, как, на пример, самолет, можно рассматривать как совокупность пар различно расположенных блестящих точек. Так как рас стояние между ними L может меняться от 0 до Lmax, то от раженный сигнал модулирован не одной частотой-, а спектром частот, лежащих в интервале от 0 до Ефлтах. В большинстве случаев можно принять 9» 0 и считать
^Фдтах = ^14г|- |
(3-4-19* |
л I at I |
|
Для ориентировочной оценки этой величины рассмотрим случай, когда самолет, летящий со скоростью v = 300 м/с, совершает вираж с ускорением 19,6 м/с2 (2g). При этом уг ловая скорость вращения dQIdt = 2g/v = 6,5 X 10-2 рад/с (« 3,7%). Если Lmax = 30 м и Л = 10 см, то в соответст вии с (3.4.19) Ефлтах — 39 Гц. В этом случае минимальный период флуктуации ЕфЛт1п = 1/ЕфЛтах ~ 25 мс. Если, на
пример, отраженный сигнал состоит |
из N = 10 импульсов |
с периодом повторения Тп —1 мс, |
то общая длительность |
пачки NTn < ^.4mln и можно считать, что в пределах длительности пачки случайных изменений амплитуды им пульсов не происходит. Однако они будут происходить от пачки к пачке («дружно» флуктуирующая пачка).
Таким образом, максимального своего значения, состав ляющего несколько десятков герц, частота флуктуаций до стигает только при маневрах самолета. Если самолет не ма
неврирует, т. |
е. когда |
колебания |
самолета составляют |
± (1 ... 2)° со |
скоростью |
dQldt = 1 |
... 2%, частота флук |
туации заметна меньше.
Важным обстоятельством является зависимость частоты флуктуаций от длины волны. Как видно из (3.4.16), частота флуктуаций становится ниже при увеличении длины волны.
Из рис. 3.17, а легко видеть, что 0 « vtID, откуда по формуле (3.4.16) получаем Ёфл « 2Lv!\D, т. е. частота
флуктуаций растет с уменьшением расстояния.
154
3.5.ВЕРОЯТНОСТНАЯ ОЦЕНКА ЭОП
1.Плотность распределения амплитуды и фазы отражен ного сигнала. Задача определения характеристик отражения реальной цели, являющейся совокупностью многих элемен тарных отражателей, связана с большими трудностями. Следует иметь в виду, что при всяких изменениях относитель
ного положения цели и РЛС будут меняться расстояния до
элементарных отражателей и и* |
ЭОП Поэтому здесь необ |
|||
ходимо воспользоваться ста |
|
|||
тистической |
моделью |
слож |
|
|
ной цели и определить ее ве |
|
|||
роятностные характеристики. |
|
|||
Представим модель |
слож |
|
||
ной цели в виде совокупности |
|
|||
большого числа случайно рас |
|
|||
положенных |
независимых и |
|
||
равноценных |
отражателей, |
|
||
амплитуда |
и |
фаза колебаний |
|
|
которых являются случайны |
|
|||
ми величинами. Однако в со- |
|
|||
ставе отражателей может на |
Рис. 3.18. Векторное представ |
|||
ходиться |
один, доминирую |
ление отраженных сигналов от |
||
щий над остальными |
и даю |
сложной цели |
щий сильный стабильный не случайный отраженный сигнал (блестящая точка). Это
случай слабо флуктуирующей цели. Если же стабильный сигнал блестящей точки отсутствует, то цель является сильно флуктуирующей.
Как следует из (3.4.1), вторичное поле /г-го элемен тарного отражателя в комплексном виде характеризуется выражением
|
ёй = £„ е>“' = *Е е“ |
е>®', |
(3.5.1) |
причем фаза |
изменяется в пределах от — л до л, а сиг |
||
нал блестящей |
точки примем в виде |
|
|
|
ёо = ЕоеК |
|
(3.5.2) |
На комплексной плоскости (рис. 3.18) вектор Ео отложен |
|||
вдоль вещественной оси, а векторы |
имеют |
случайные |
значения амплитуды и фазы. Амплитуды их можно разло жить на ортогональные вещественные и мнимые элементар ные составляющие Eh cos и jEft sin После суммиро вания п таких элементарных составляющих получим ам
155
плитуды ортогональных составляющих суммарного случай ного процесса (разложение на квадратурные составляющие):
£si = 2 cos (pfe |
= £2 cos ф2, |
(3.5.3) |
A=! |
|
|
£22 = 2 £&sin(Pfc |
=^ssincp2. |
|
R=1 |
|
|
При этом случайная комплексная амплитуда |
|
|
£2 = £2е^2 |
(3.5.4) |
характеризует вектор суммарного поля п элементарных от ражателей.
Составляющая находится в фазе с вещественной ча стью Ео и поэтому именуется фазной. Составляющая Е2 =
— Ex?, сдвинутая по фазе относительно Ео на 90°, называ ется внефазной. Кроме того, на рис. 3.18 показана резуль тирующая вещественная составляющая £\ = EQ + £2i. Результирующий же вектор
ер = Ер ei“f = Ер е1фр е1й<,
где Ер и фр — случайные величины, характеризующие дли ну и фазу результирующего вектора поля вторичного излу чения.
Действительное мгновенное значение
ер = Ер cos (cot + ФР) — Ео cos at + E% cos ((at + фД
(3.5.5)
где e0 = Eq cos — стабильный сигнал или когерентная составляющая, а е2 = f2 cos (at + ф2) — случайный сиг нал или некогерентная составляющая.
Заметим, что составляющие E^i и £22 случайной ком плексной амплитуды £2 могут принимать как положитель ные, так и отрицательные значения. Так как нет основания для преобладания сигналов того или иного знака, то сред
ние значения |
соответственно |
равны |
|
£si=£22 = 0; |
E^Eq+E^Eq', £2 = £22 = О. |
(3.5.6) |
|
Что касается дисперсии, то здесь также нет оснований |
|||
ожидать различий, откуда |
|
|
|
(£i-£i)2 - |
=• = о1. |
(3.5.7) |
Составляющие Ег и £2 совершают случайные колебания относительно своих средних значений, причем эти колеба-
156
ния распределены по гауссовскому закону. Это обстоятель ство связано с тем, что каждая составляющая обязана сво им происхождением большому числу элементарных отража телей. Вместе с тем, согласно центральной предельной тео реме, сумма большого числа независимых случайных вели чин, среди которых нет явно преобладающих над совокуп
ностью остальных, распределена по гауссовскому |
закону. |
||
Таким образом, имеем |
|
|
|
= |
— ехр |
(Ех ~ Ео)2 ' |
|
|
2ла2 |
2сг| |
(3.5.8) |
w (£2)—1— ехр |
/72 |
||
д а |
|
*]/2лсг2 *2а
Найдем теперь закон распределения параметров резуль тирующего вектора, а именно: его амплитуды Ер и фазы <рр (см. рис. 3.18). Будем исходить из того, что ортогональные
составляющие |
и Е2 являются независимыми случайными |
величинами. Как известно из теории вероятностей, любая нормально распределенная двумерная случайная величина путем поворота осей координат может быть приведена к си стеме двух нормально распределенных независимых слу чайных величин. При этом двумерная плотность распределе
ния вероятностей равна |
произведению |
|
w (Elt Е2) = w (Ег) w (Е2) = |
|
|
1 |
(Ех - Ео)2 + Ej |
(3.5.9) |
-------ехр |
|
2а|
Так как нас интересуют законы распределения вероят ностей результирующей амплитуды £р и фазы <рр, то целе сообразно перейти от прямоугольной системы координат, -в которой представлены величины Ег и Е2» к полярной си стеме координат и величинам Ер и <рр.
Связь между этими системами координат следующая:
Ер = VEl + Е%; |
фр = arctg (Е^/Е^ |
Ех = Ер cos фр; |
Е2 = Ер sin <рр. |
Для нахождения двумерной плотности распределения ве роятностей в полярной системе координат воспользуемся известным из теории вероятностей соотношением
(£„, <рр) = w (Elt |
, |
(3.5.10) |
О (Ер, фр)
157
где второй |
множитель — функциональный |
преобразова |
||||||
тель (так называемый якобиан), равный |
|
|
||||||
|
|
dEi |
dEt |
cos фр — Ep sin фр |
||||
d(Elt E2) |
_ dEp |
дфр |
||||||
|
|
|
= BD. |
|||||
d(Ep, <pp) |
dE2 |
dE2 |
sin фр |
|
p |
|||
|
|
dEp |
dq)p |
Bp cos фр |
||||
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
(3.5.11) |
|
Подставляя (3.5.11) в (3.5.10), |
получаем после элемен |
|||||||
тарных |
преобразований |
|
|
|
|
|||
& (Вр, |
|
Р |
/ |
— Еа \ |
|
F |
||
Фр) = —ехр (------ — ) ехр [ |
\ |
cos фр |
||||||
|
|
2лог?, |
\ |
2ст| |
/ |
ог| |
||
|
|
|
|
|
|
|
(3.5.12) |
Теперь остается найти одномерные плотности распреде ления величин Ер и фр. При этом надо иметь в виду, что значения фр могут лежать в интервале от — л до л, а значе
ния Ер соответствуют |
Ер > 0. |
В результате |
получим |
|||
я |
w(Ep, ФрМфр- —%- [ехр—(В^+ В§)/2о^] х |
|||||
w(Ep)= I |
||||||
J |
|
Оу |
|
|
|
|
|
|
X /0 (Ер Е0/а1), |
|
(3.5.13) |
||
где |
|
|
|
r>*£ Eq |
|
|
|
v |
|
л |
|
|
|
/ |
|
_Р—L соз фп |
dq> |
|
||
/-EpE»\=J_ |
С |
е “s |
(3.5.14) |
|||
|
I oj ) |
2л |
J |
|
|
|
|
|
|
|
|||
|
|
|
— Л |
|
|
|
/0 — функция Бесселя нулевого порядка мнимого аргумен та (модифицированная функция Бесселя).
Функция 1Q (х) |
характеризуется тем, что при х = О |
она равна единице, |
а при х » 1 справедливо асимптотиче |
ское выражение /0 |
(х) ~ ех/]/г2лх. Функция распределения |
вероятностей (3.5.13) именуется обобщенным законом рас пределения Рэлея или распределением Райса.
Часто удобно пользоваться безразмерными величинами v = Ер/в?, а = Eq/o^- Для перехода от плотности распре деления случайной величины Ер к плотности распределения функционально связанной с ней величины v — Ер/в% вос пользуемся известным правилом:
w(v)= w(E„)|-^P-|-ue-^‘+a’)/2/0(w). (3.5.15) I dv I
158
Аналогично для двумерной плотности (3.5.12) получаем
выражение |
|
|
' |
|
0,(0, фр) = J_(,e -(О’+а«)/2е«’«*Фр |
(3.5.16) |
|||
|
2л |
|
|
|
Закон распределения фазы определяется путем интегри |
||||
рования по v в интервале (0, |
оо) . |
|
|
|
|
со |
|
|
|
|
(фр)=J w (v> Фр) |
= |
|
|
|
о |
|
|
|
== f J_re*)/2-+a(« |
eaDCOS<pPdo. |
|
||
J |
2л |
|
|
|
о |
|
|
|
|
После несложных |
преобразований имеем |
|
||
до (фр) = ~~ е“а’/2 + —[1 |
+ Ф (a cos фр)1 |
х |
||
2Л |
2 Д/2^ |
|
|
|
|
Xe-(a’s,n’p)* /2, |
(3.5.17) |
2 р2
где Ф (г) =—— I e~xi/2dx—интеграл вероятности.
Т/2л .) * о
Важным для практики случаем является отсутствие ста бильной составляющей, когда Ео = 0, а — 0. При этом на основании (3.5.15) закон распределения относительной ам плитуды будет иметь вид
w (с») = не-0’/2. |
(3.5. 18) |
Это простое распределение Рэлея, которое характерно для амплитуды (огибающей) случайного процесса, имеюще
го нормальный закон распределения. Среднее |
значение от |
носительной амплитуды |
|
и = Ju2 e~t’*/2tfo = |
(3.5.19) |
На рис. 3.19, а изображены кривые распределения (3.5.15) и (3.5.17), построенные для различных значений со ставляющей сигнала. По мере роста сигнала (т. е. вели чины а = £0/(Tz) закон распределения нормализуется. Это легко уяснить из векторной диаграммы рис. 3.18. Действи тельно, по мере роста вектора Ео результирующий вектор
159