Литература / бакулев радиолокация распозн
.pdf
Так как шум и смесь
сигнала с шумом распре делены на выходе детек
тора огибающей (дО) по законам Рэлея и обоб
щенного Рэлея (Райса),
вероятности
F = J-exp |
--- |
|
dZ = |
|
|||||||||
r,) |
z |
{ |
|
|
|
z" } |
|
|
|||||
|
h |
(]"2 |
|
|
"а-2 |
|
|
|
|||||
|
|
z |
|
|
|
- |
|
z |
|
|
|
||
=exp{-_!l__}, |
|
|
|
|
|
|
|
6) |
|||||
|
|
|
20-: |
|
|
|
|
|
|
|
|
|
|
(fJ |
|
|
|
|
|
|
|
|
|
|
Р11с. 3.10. Коррсляциош1ый (а) 11 фильтровой (6) об |
||
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
-:, |
фюой |
|||
|
-pp --,, 0 |
(ZEJdZ, |
|||||||||||
D= |
|
z |
{ |
|
2 |
|
|
|
|
|
|||
|
Z 1 |
|
наружитсли радиоим11ульса с нсювестной 11ачаныюii |
||||||||||
|
|
|
|
|
|
|
} |
|
|
|
|
|
|
/1 |
О": |
|
'2d; |
|
|
|
|
а-; |
|
||||
где о-~ |
NE |
|
|
|
h |
|
=R |
~ |
|||||
= - |
|
|
|
|
|
|
нормированныи порог. |
||||||
|
|
|
0 |
|
|
- |
|
|
|
|
21n- |
|
|
|
|
|
- ; |
|
|
|
|
|
|||||
|
|
|
2 |
|
|
|
а-= |
|
|
|
F |
|
|
Обнаружение флуктуирующего радиоимпульса. «Белый» шум складывается с радиоимпульсом, имеющим неизвестную начальную фа
зу и флуктуирующую амплитуду: u(t,(j)}=aИт(t)cos[ro0t+\j/(l)-<j>]. Плотно
сти распределения вероятностей амплитуды и фазы задаются соотноше-
ниями w |
(<j>)=l/2л и 1,v(a)=а-, ехр{--а2 } |
- |
закон распределения Рэлея· |
. |
|
0 |
а-; |
2 |
|
. |
|
|
(]"а |
|
|
||
Усредняя Л по а и (j), получаем
щ2,r
Л(у)=Л(у,<j>,а)q,,а = JJЛ(y/a,<j>}1v0(a)1v0(<j>)d<j>da=
= jл(y/a),vo(a)da=I·+о о ~~н~~Jwo(a)da,
где Е энергия сигнала при а =1.
Средняя энергия сигнала е = М{а2Е}=ЕМ{а2}=2а/Е. С учетом
этого вычисляем отношение правдоподобия:
Л(у)= ~ахр |
|
2 |
|
ехр |
|
2 |
|
|
(2aZ) |
(ia. |
|||
{ |
---ат |
} |
{ |
--а,Е-}/0 |
|||||||||
|
|
"J' |
|
|
lv 0 |
|
|
-- |
|
||||
|
|
0 au |
|
2а" |
|
|
|
|
|
N 0 |
|
||
Воспользуемся интегралом |
|
|
|
|
|
|
|
||||||
щ |
x |
ехр{ -ах~}/ |
|
|
) |
|
|
{ |
р2 |
|
} |
|
|
|
(Дх)dх = -ехр |
|
-- |
|
|
||||||||
f |
|
0 |
|
|
2а |
|
|
4а |
|
|
|
||
11 |
|
|
|
|
|
|
|
|
|
|
|
|
|
61
и получим окончательное соотношение |
|
|
|
|
|
|
|
|
|||||||
|
Л(у)= |
No - |
ехр{ |
20-~z2 - } . |
|
|
|
|
|
|
|
|
|
|
|
|
|
N0 +E |
|
N0 (N0 +E) |
|
|
|
|
|
|
|
|
|
|
|
Таким образом, решающее правило может строиться на основании |
|||||||||||||||
сравнения Z2 с порогом: |
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(1.12) |
Структура обнаружителя показана на рис. 3.11. |
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
Определим вероятно |
||||
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
сти ошибок F и D: |
|
||||||
|
|
|
|
|
|
|
|
|
|
|
F= J-exp -- dZ= |
||||
|
|
|
|
|
|
|
|
|
|
|
~ z |
{ |
22 |
} |
|
|
|
|
|
|
|
|
|
|
|
|
/1 |
о?, |
|
~ |
|
|
|
|
|
|
|
|
|
|
|
|
=~{-;} |
|
|
||
|
|
|
|
|
|
|
|
а для вычисления D нужно |
|||||||
|
|
|
|
|
|
|
|
знать |
плотность |
вероятно |
|||||
|
|
|
|
|
|
|
|
сти |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2Z |
|
|
|
|
|
|
|
|
|
|
|
|
и{Z/а,0=1)=-х |
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
N0 E |
|
Р11с. 3.11. Корреляционный (а) и фильтровой (б) обна |
хе |
|
|
2 |
2 |
2 |
|||||||||
|
|
|
|
||||||||||||
ружители радиоимпульса с ••еювест•юй начальной фа |
хр |
{ |
Z -a E } 1O(2aZ) |
||||||||||||
зой и флуктуирующей амплитудоii |
|
|
|
|
N.E |
|
N |
' |
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
о |
|
о |
|
тогда |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
w( Z / 0 = 1) = |
2Z |
2 2 ехр{ |
|
|
|
|
|
|
|
|
|
|
||
|
,,,,h |
N0Е+2баЕ |
1+ 2баЕ 1 |
|
|
|
|
||||||||
|
|
|
|
!N Е |
|
|
|
|
|||||||
и поэтому |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
D = f}V ( Z / 0 = 1) dZ =ехр |
2 |
|
|
|
|
|
|
|
|
|
||||
|
( h |
2 ) |
• |
|
|
|
|
|
|||||||
|
|
|
|
|
о |
N |
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
о |
|
|
|
|
|
|
|
|
Сравнение D и F различных моделей сигналов можно осуществ- |
||||||||||||||
лять только при равенстве их энергий, т.е. при |
б"2 = - F, |
= -1 . Поэтому, |
|||||||||||||
исключая h, получаем |
|
|
|
|
|
|
2Е |
2 |
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
D = F |
1+~ |
---'- |
|
|
|
|
|
|
|
|
|
|
|
|
|
Nн = F l+ч , |
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
где ц = РсlРш.
62
Обнаружение флуктуирующего радиоимпульса с неизвестным
(случайным) временем прихода (-r= tн) и с расстройкой по частоте
(о,= Ф;1)- В этом случае сигнал имеет вид
U(t,<p,a,т,ro)=aum(t-т)cos[(ro0+ro)(l-т)+w(t-т)-<p] ,
где а,<р,т,rо - случайные величины с известными априорными распреде лениями вероятностей.
Тогда нужно найти Л(у/0,<р,и,r,ш). Поскольку известно Л(y/0,<p,a)IJ)." =
= Л(у/0), введем в Л(у/0) новые случайные параметры т, ro и усредним Л по этим параметрам. При этом будем считать, что они изменяют свои
значения дискретно, т.е. т=т 1,т2,••• ,т,,... ,тп , ro=ro 1,ro2,...,rok,'···,rom с вероят
ностями ptl=P{т=Тi}, i=l,2,3, .. .,п; pooJ.-=P{ro=rok}, k= 1,2,3, .. .,111.
Определить Л(v/0) можно как и ранее:
Л(у)= I Рr,РшkЛ(у/t'"Шk),
1,k=I
(3.13)
Таким образом, синтезированный обнаружитель состоит из п кана лов по времени запаздывания (дальности) и т каналов по частоте (ско
рости). В действительности и время запаздывания, и смещение частоты - непрерывные величины, поэтому многоканальный обнаружитель яв ляется квазиоптимальным. Только по мере увеличения числа каналов (в пределе до бесконечности) он будет приближаться к оптимальному. Многоканальность по дальности необходима при использовании КО, в
этом случае для каждого канала формируется свой опорный сигнал,
сдвинутый по сравнению с соседними каналами на величину элемента
разрешения по времени (ти)- В ФО из-за инвариантности согласованного фильтра ко времени прихода сигнала многоканальность по дальности не
нужна. а многоканальность по скорости обеспечивается набором (гре
бенкой) фильтров, расстроенных друг относительно друга на величину
элемента разрешения по частоте (скорости) (8rод).
Обнаружение пачки когерентных радиоимпульсов. Реализации аддитивной смеси «белого» шума с когерентными радиоимпульсами
пачки в каждом периоде повторения имеют ту особенность, что шум в
них не коррелирован, так как время корреляции шума ткор меньше пе
риода повторения Тп:
Тк0,1 = 1/Л/~ф~ Ти< Tn.
63
Плотности распределения вероятностей
пп
}V11 (y/0=0)= пW; (у/0 = 0), Wn(}7/0 = )) =пW; (у/0 = 1),
|
|
|
|
|
i=I |
|
i=I |
|
поэтому |
|
|
|
|
|
|
||
|
|
|
|
|
|
,, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Пw,(y/0 = 1) |
,, |
|
|
Л |
п |
( |
у |
) = }Vп(у/0 = 1) = |
i=I |
=ПЛ-( |
) . |
|
|
|
w,,(y/0 =0) |
11 |
'у |
|
||
|
|
|
|
|
|
пW;(y/0 =0) |
,=I |
|
|
|
|
|
|
|
,=1 |
|
|
|
|
|
|
|
|
|
|
|
п
Соответственно ln Лn(У)LlnЛ;(Y).
i=I
В согласии со случаем обнаружения для модели сигнала с полно
стью известными параметрами получаем
Л,(у)~•+i}p{~}•
следовательно,
откуда
I=,11 ~(No /2)lnT +(1/2)Ltl Е, = ипор.
i=I |
,=1 |
Структуры обнаружителей пачек показаны на рис.
rz-
~
~
rиnop
а) б)
в)
3.12.
Здесь ОФ1
- оптимальный
фильтр для оди
ночного импу
льса; ОФ,, -оп
тимальный фи
льтр для пачки
из п импульсов
(состоит из
ОФ 1 и ~); ~ -
1•11с. 3.12. Схемы обнаруж1rrелей пачки когсрентньiх радиоимпульнакопитель им
сов: а - струкrура КО пачки: 6 - ~,руктур ФО: в - структура КФО |
пульсов; ПУ _ |
64
Учитывая, что Z;= z1,cos<p + z2,siшp, получаем
" |
п |
|
LZ,cos(<p- v,) = L[z1,cos<p + z2;sin<p] =Zrcos(<p- vr), |
||
i=I |
i=I |
|
|
|
fz2; |
|
|
|
где Zr= |
|
; ur= arctg ~ . |
|
|
LZi; |
|
||
|
|
i=I |
что приводит к алгоритму обнаружения:
lz~:ипор./
Соответствующие этому алгоритму структуры обнаружителей по
казаны на рис. 3.15.
а)
б)
Сигнал с неизвестной начальной фазой и флук туирующей амплитудой. В соответствии со случаем об наружения сигнала с флук
туирующей ампли,удой, ана
логично предыдущему имеем
Л"{у,<j>,а)=ех+at:•}x
2a;Z;COS ( q> - V;)}
хехр ------ ,
{ No
Р11с. 3.15. Коррсляцио1шый (а) и фильтровой |
п |
(б) обнаружители па•1ек когерентных радио |
Лп(у,q>,а)=ПЛ;(у,q>,а). |
|
|
импульсов с неизвестной начальноn фазой |
1=\ |
|
66
(3.15)
Структура обнаружителя представлена на рис. 3. 16, а прохождение сигнала через нее иллюстрируется эпюрами напряжений (рис. 3.17), где U1(t) - пачка из трех радиоимпульсов на входе обнаружителя; U2(1) -
радиоимпульсы на выходе согласованного фильтра; Uд(t) - видеоим
пульсы на выходе детектора;
|
|
|
|
Uд(t-Тп} - пачка |
импульсов, |
|||
|
|
|
|
задержанная на период и Ид(t- |
||||
|
|
|
|
2Тп) на два периода; Ur,(t) - |
||||
|
|
|
|
результат накопления. |
||||
|
|
|
|
|
Для отождествления опе |
|||
|
|
|
|
ратора ln/0(x) с конкретным |
||||
|
|
|
|
устройством |
рассмотрим по |
|||
|
|
|
|
ведение ln/0(x) при больших и |
||||
|
|
|
|
малых значениях аргументах. |
||||
|
|
|
|
При х>1 ln/0(x)~x |
- линейная |
|||
|
|
|
|
функция, при х<1 lnl0(x}~x2/4 - |
||||
|
|
|
|
квадратичная функция. По |
||||
|
|
|
|
этому можно считать, что не |
||||
|
|
|
|
линейный |
элемент |
ведет себя |
||
|
|
|
|
как обычный амплитудный де- |
||||
Р11с. 3.17. Граф11ки процесса накопления пачки |
тектор и вкл_ючать в схему рис. |
|||||||
некогсре1т1ых радиоимпульсов |
3.18 |
вместо |
блока |
ln/0(x) ам |
||||
|
|
|
|
плитудный детектор или де |
||||
|
|
|
|
тектоf) огибающей. |
|
|||
|
|
|
|
|
Можно показать, что при |
|||
|
|
|
|
нефлуктуирующей |
амплитуде |
|||
Рис. 3.18. Обнаружитель пачки 11екогере11тных ра- |
импульсов |
|
|
|
||||
ДИОll!\Ш}'JIЪСОВ |
|
) |
vJ,z2n-lexp{- z2 Lп, |
|||||
|
|
|
|
F |
||||
|
|
|
|
2п(п-1)! h |
|
2 ГL, |
||
|
|
|
D |
ехр{- nq} .п |
|
{ |
2 } |
|
|
|
|
|
2 fz" |
ехр |
|
--z х |
|
|
|
|
|
(пстi"-1)12 |
• |
|
2 |
|
|
|
|
|
|
11 |
|
|
|
а) |
б) |
х/,,_1( j;qz)dz, |
|
|
|
|||
где q=цi=PclPш - |
|
|
|
|||||
|
|
|
отношение |
|||||
|
|
|
||||||
Р11с. 3.19. Быстрые (а) и мед.ленные (б) флу1<1уа- |
мощности |
сигнала |
к |
мощно- |
||||
цин амптпуды импульсов пачки |
|
|
сти шума в одном импульсе. |
|||||
68
лем, и обратная ей корреляционная матрица R~~ === Qjk. Причем алго
ритм обращения корреляционной матрицы
,,
LR1iQik =б,k,
i=I
где 81k - так называемый символ Кронекера, который равен 1 при j=k
или О приJ::;:k.
Пусть последовательность выборочных значений {y(t1)J'(t2),y(t3), •• • , y{tn)}={y1J,'2J,'з,... J,'n} образует вектор у~ ===У. Совместные плотности рас
пределения вероятностей выборочных значений можно представmъ в виде
Найдем отношение правдоподобия
Л(У) = w(y /0 =1) |
|
|
|
||
|
w(y/0 =О) |
|
|
|
|
= ехр{ff Q iky(t; )u(tk )-_!_ f f Q iku(t i }u(tk )} 11,'(Уо/ 0 = 1) |
|||||
|
,=1 k=I |
· |
2 J=I k=I |
W(Yo /0 = О) |
|
и перейдем к его логарифму для сравнения с порогом решения lnT: |
|||||
|
11 |
11 |
1 п 11 |
|
~ / 0 - 1) |
ln Л(У)= LLQ ,ky(t1)u(tk )--LLQ.iku(ti )u(tk) +ln |
° = ~ln Т. |
||||
|
i=I k=I |
2 i=I k=\ |
|
w(yo/0-O) |
|
Объединяя слагаемые этого выражения, не зависящие от У, в по |
|||||
роговое напряжение Unop , получаем алгоритм |
|
|
|||
,, |
" |
|
|
|
(3.17) |
LLQ1kY(t1)u(tk):unop · |
|
|
|||
.1=1 |
k=I |
|
|
|
|
,,
Если обозначить WJ =L Q1ku(tk), отождествляя W1 с весовыми коэф
k=I
фициентами фильтра, то алгоритм обнаружения становится более понятным:
11
LWiy(t1) ~ Uпор•
/=1
70
