Добавил:
Закончил бакалавриат по специальности 11.03.01 Радиотехника в МИЭТе Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Литература / бакулев радиолокация распозн

.pdf
Скачиваний:
39
Добавлен:
10.09.2023
Размер:
15.09 Mб
Скачать

Так как шум и смесь

сигнала с шумом распре­ делены на выходе детек­

тора огибающей (дО) по законам Рэлея и обоб­

щенного Рэлея (Райса),

вероятности

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

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D = F

1+~

---'-

 

 

 

 

 

 

 

 

 

 

 

 

 

= 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 и усредним Л по этим параметрам. При этом будем считать, что они изменяют свои

значения дискретно, т.е. т=т 12,••• ,т,,... ,тп , 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