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

книги из ГПНТБ / Курикша, А. А. Квантовая оптика и оптическая локация (статистическая теория)

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

для выяснения основных особенностей оптимальной об­

работки при рассматриваемых здесь свойствах

сигналов

и помех.

 

 

 

 

 

Спектр поля

на апертуре будем представлять в виде

суперпозиции

плоских

волн

(см. §. 2.1)

 

у[(г, ш) =

-

ik / 4

g - \

а (р, ш) е У г р г / (р) d9,

(5.3.1)

где размерность у (г, ш) выбрана так, что \y(r, t)\~ сов­ падает с плотностью потока мощности; а(р,ш) опера­ тор уничтожения, соответствующий элементарной пло­

ской волне с направлением прихода

р, f(p) функция,

учитывающая

ослабление влияния

поля

при увеличе­

нии угла Ѳ между нормалью -к апертуре и

р (в § 1.4 бы­

ло выбрано

| f (р) |2 =cos Ѳ). Операторы

а (р, со) и

а+ (р, ео) = а(р,—со) подчиняются перестановочным соот­

ношениям (1.4.6).

Сопоставляя (5.3.1) и (5.2.11), замечаем, что роль суммирования по индексу в (5.3.1) играет интегрирова­ ние по р, а со играет роль индекса, значения которого

одинаковы и в â(p, со), и в у(г, со). Поэтому ядро K(Si, S 2 ) в данном случае имеет вид

 

 

К (Г, , Г,, со, <в') =

Ка

(Г,, га ) 5 (со' -

со),

где

учетом

того, что в гл. 1 выбрано

\f (p)|2 =cos Ѳ)

 

 

 

К (г г\—

ftc°3

А ( * | Г , г2 |)1

,е о г,\

Эта

функция

уже встречалась в гл. 1

(функция корре­

ляции

фона)

и в гл. 2 [см. (2.1.10)]. Было

показано, что

эта функция имеет спектр, постоянный в области про­ странственных частот \p\<k и равный нулю вне этой области. Пренебрегая краевыми эффектами вблизи гра­

ницы апертуры и учитывая, что функции, к

которым

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

Km(vlt

Tz),

имеют спектр, сосредоточенный в области

p<k,

бу­

дем считать ядро (5.3.2) эквивалентным ô-функции

(см.

гл.

2):

 

 

 

КАт1г)^2-фЛ(г,-гг).

(5.3.3)

160

При этом ядро оператора, обратного оператору с ядром (5.3.2), запишем в виде

L (г,. га, Ш „ œs) «

8 (г, - г„) § К

- «>,).

(5.3.4)

Подставляя (5.3.4)

в (5.2.21) и считая

сигналы

£Уѵ,и>)

достаточно узкополосными, чтобы пренебречь изменени­ ем энергии кванта в пределах полосы сигнала, получаем

ßv = р = - j * ( " ^ (г. /) t/*y (г, t) I - I /

Г J|t/„(r, 0 | W r

—оо S

(5.3.5)

Величины | ß j 2 образуют достаточную статистику.

Этот результат совпадает с получающимся при клас­ сическом рассмотрении (умножение на ожидаемый сиг­ нал и интегрирование). Подчеркнем, что в рассматривае­ мом высокочастотном приближении (когда длина волны мала по сравнению с размерами апертуры) одновремен­ но измеримыми оказываются величины ßv , соответствую­ щие просто ортогональным (без всякого веса) сигналам

(г> 'Ч и> в

частности, сигналам, не перекрывающимся

по времени

или имеющим неперекрывающиеся спектры,

а также сигналам на неперекрывающихся частях апер­ туры.

Рассмотрим теперь случай гауссова поля, состоящего из излучения фона с равномерным распределением по углам прихода и полезного сигнала в виде пакета пло­ ских волн со средним направлением прихода р0 . Закон изменения сигнала во времени будем считать стационар­ ным случайным процессом со спектральной плотностью среднего числа квантов Л^с(со), промодулированным функцией u(t). При этом функцию корреляции поля у(г, со) записываем в виде

• R (г,. », га, ш) =N (ш) КШ (г,, г2) 8 (ш'—ш)

+

 

+ Rc К . ш 2 ) ехр [іф (rlt ш,) гф 2, со,)],

 

(5.3.6)

где N((Ù) ='іѴф(со)^2 ш среднее число квантов

излуче­

ния фона на одну степень свободы (в элементарной пло­ ской волне); -ф(г, со) характеризует пространственную

161

фазовую модуляцию поля на апертуре;

со

R*c (со,, ш2 ) = - i - j " W c (eu) U (со, - со) (7* (со, - со) й Ц (5 . 3 . 7)

и

с7(со)спектр нормированного закона модуляции «СО- Задача отыскания достаточной статистики сводится в этом случае к'решению уравнения ( 5 . 2 . 2 0 ) , записывае­

мого в виде

0

0

 

R* ( Г , CD, Г ' , CD') V ( Г ' , СО') dco'dr'

 

оJ

Sf

=

 

 

=

ÄJtfe (r,r')V(r',<o)dr'.

(5 . 3 . 8 )

Для отыскания решения прибегнем к ряду упрощений.

Будем считать, что

условия,

при которых

справедливо

равенство

(5 . 3 . 3) ,

выполнены. Тогда

Ѵ(г,

ш)

удастся

представить в виде

 

 

 

 

 

 

Ѵ ( г ,со) = Ѵ ( с о ) е ' ф ( г ' ш ) ,

 

 

причем для Ѵ(со) из (5 . 3 . 8)

с учетом

(5 . 3 . 6)

получим

следующее

уравнение:

 

 

 

 

оо

 

 

 

 

 

 

5 j R*c

(со, со') V (<•>') сісо' =2ф>

(» - N (со)) V (со).

( 5 . 3 . 9 )

о

 

 

 

 

 

 

Это уравнение удается решить для двух крайних случа­ ев: когда спектр /Ѵс(си) значительно уже или значительно шире спектра модуляции іУ(со). Результат для первого случая получается таким же, как и для регулярного сиг­ нала со случайной фазой. Соответствующая функция правдоподобия получается из функции правдоподобия

регулярного

сигнала со

случайной фазой

усреднением

по амплитуде.

 

 

 

 

 

 

В случае

узкого спектра

модуляции

 

 

 

 

+

0 О

 

 

 

R * с (со,, « g «

- ^ - ]/"У ѴС (СО ,) І ѴС (СО 2 )

j0 0

U (со) U* 2 -со,+со)й?со=

= K t f o M t f c K )

J

I и (t) I2

e - ' ( - ' - ^ > ' Ä = •

 

 

—оо

 

 

 

 

 

= ѴЩ^ЖК)

 

U0 (ш, -

ш2).

(5.3.10)

162

Выделяя из Ѵ ( ш ) множитель

j.''/Vc (to)

(V ( ш ) = ѵ / W 0

(u>) Ѵ ( ш ) )

и вынося

с (аз)

из под

 

интеграла,

(5.3.9) перепишем в

виде

 

yVc(co)5j't7o(cû—со') У (со') dec/=

 

 

 

 

 

 

 

 

= 2іфи(й—W (ш)] Ѵ'((о).

 

 

(5.3.11)

 

Аппроксимируем

закон

модуляции

интеисивностей

| « ( / ) | 2 ступенчатой

функцией

 

 

 

 

 

 

 

 

 

 

 

 

со

 

 

 

 

 

 

 

 

 

 

 

 

Ѵ = — 0 0

 

 

 

 

 

где

J ] ( J C ) = 1

при

I X | < 1/2

и [ ] ( л ' ) =

0

при | J C | > 1/2.

 

 

 

 

 

 

 

 

 

д

 

 

 

 

Тогда

 

 

 

 

sin со - g -

оо

ц е _ ' ѵ ш Д .

 

 

 

 

[/„ (ш) =

 

^

V

 

 

Используя

 

условие

 

с(со + 2л/Д) «JVC (to),

простой

подстановкой

нетрудно показать, что в этом случае урав­

нению (5.3.11)

удовлетворяют функции

 

 

 

 

 

 

 

 

д

f

 

2*

\

 

 

 

 

 

 

 

 

sin

I Ш — (X—г- I

 

 

v

' 4

»

=

^

 

 

 

 

 

 

е-'"'ш Л .

(5.3.12)

соответствующие

собственным значениям

 

 

 

% = «ѵ

 

J

+NU^-).

(5.3.13)

Постоянную А

 

 

 

ftfx

из условия

нормировки

 

определяем

Окончательное

 

выражение

для V

(г, <о) имеет вид

 

 

 

 

 

 

 

д

 

 

 

 

 

 

 

 

 

Н-

sin-n-(cû — ш )

 

.

 

V . (г. ») =

к

 

(ш — Ш | А )

е ' * ^ » 1 - ^ .

(5.3.14)

где

 

 

 

 

 

ш =

(х2и/Д.

 

 

 

 

 

 

 

 

 

163

Достаточной статистикой является совокупность ве­ личин /гѵи. = 2*ѵ [ і гѵ и ., причем

 

 

со

 

 

 

 

 

 

 

 

 

z

= [ dm { drV*

, (Г, to) у (Г,

с о ) =

 

X

 

 

O

S

'

 

 

 

 

 

 

 

 

 

 

 

à

 

 

 

 

 

 

 

CO

 

ein (Cû — Cü.J -ô-

.,,

S | .

 

 

 

X ± . \ d * \ —

Ц _ і е - « + ( г . » ) + - д х

 

 

 

 

 

oo

 

 

 

 

 

X{y(r,cu)dr

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

со

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

(5.3.15)

 

 

 

 

 

 

 

 

 

 

 

 

y(/) =

-ir j'ei t t 'dcoji ,(r,ai)e_ '"'i , ( r , w ) rfr.

(5.3.16)

 

 

 

 

0

 

 

 

 

 

 

Если для всех

г на апертуре ехр (г, а>

,)Ф(г, со^)]}:^

« s i ,

то

У(/;) в

(5.3.15) можно заліенить на Y (t):

 

 

M ' ) = f^(r.0e~""*( r 't O , l )

rfr.

(5.3.17)

Формирование

достаточной

статистики

 

распадается

на следующие операции: пространственную

корреляцион­

ную

обработку

с

опорными сигналами

ехр[іг|з(г, со )],

пропускание каждого ц-го сигнала через фильтр с ча­ стотной характеристикой вида (sin я)/я и отсчет ннтенсивностей на выходах фильтров в дискретные моменты времени \ =ѵД. Эту достаточную статистику можно, очевидно, использовать как в задачах проверки гипотез и оценки параметров, так и в задачах фильтрации мо­ дулирующего сигнала.

Если рассматривается совокупность сигналов с орто­ гональными законами пространственной модуляции

- j - \ е х Р *— (г -ш ) - Ф* ( г . Ш Л d T = hb>

s

164

то достаточная статистика для всей совокупности сигналов может быть составлена из величин |zv [ i |2 , определяемых формулами (5.3.15), (5.3.16) с различ­ ными Фэ- (г, ш).

5.4. Основные особенности структуры

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

сучетом квантовых эффектов

Вэтом параграфе выделим и обсудим основные осо­ бенности оптимальной обработки поля с учетом кванто­

вых эффектов и сравним эту обработку с получающейся

вклассическом случае.

Вслучае регулярного сигнала со случайной фазой обработка не отличается от классической. Если рассма­ тривается совокупность взаимно ортогональных сигна­ лов, то способ объединения величин лѵ , полученных для

каждого из них (см. (5.3.5)], определяется функцией правдоподобия, записываемой в виде произведения рас­ пределений (5.2.10). При произвольных значениях пара­ метров это выражение получается, очень сложным. Для сильных сигнала и шума ( М » 1 и N^>1) оно переходит в классическое (с обобщенными релеевскими распреде­

лениями). Для У

 

 

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

(5.2.10)

переходит

в пуассоновское

и отношение

правдоподобия

приобрета­

 

Ѵ < С І

 

 

 

 

ет вид

 

 

 

 

 

 

 

 

Л (я,

пт) = П

 

 

 

exp j J

« In

- л Л .

(5.4.1)

Для совокупности

сигналов

достаточной статистикой

является

в этом

случае величина

 

 

В случае гауссова

сигнала с широким

спектром

флюктуации элементы

достаточной

статистики

опреде­

ляются формулой (5.3.15), а способ

их объединения —

165

формулой (5.2.7). Для отношения правдоподобия имеем

In Л [у (г, т =

-±-^У*Цх)Г(Ц

• X

 

—со

 

г*

 

 

 

 

 

Р

Jn

X

 

 

 

 

X dt4L-

^ in

+

 

1 +

/V ( %

 

 

 

 

где все обозначения те же, что и в § 5.3. Выражение для А(у(г, t)] можно существенно упростить, если заменить суммирование по интегрированием по частоте. В ре­ зультате, используя предположение о широкополосное™ флюктуации, получаем

 

 

0 0

/

 

 

 

 

In Л [у (г,*)]*»

J Л

J w(t,t

t')Y(f)dt'

 

 

 

 

 

(0 I 2

M (со)

(5.4.3)

 

 

 

 

1 + N (со)

 

 

 

 

 

где

M((Ù) =Nc(cû)S/%m

— спектральная

плотность сред­

него

числа квантов

сигнала за

секунду; w(t, т ) — и м ­

пульсная переходная функция фильтра, квадрат модуля частотной характеристики которого медленно меняется

со временем и определяется

следующей формулой:

 

 

(5.4.4)

1

+

\аМ*Т+Щ»)

Множитель l/fco обеспечивает измерение выходного сиг­ нала в числах квантов независимо от их энергии.

Если # ( W ) < 1 ,

M(<ä)\u(t) | 2 < 1 , то

 

|Ktf,co)|2

Sftco In

l+\u(t)

, M (со)

(5.4.5)

 

N (со)

Такой вид весовой функции типичен

для случая

весьма

слабых (по абсолютной

величине)

сигналов, когда кван­

товые эффекты играют существенную роль.

 

Если, наоборот,

(со) ^> 1, то

 

 

 

 

1

 

'I и (if)|g

(со)

(5.4.6)

Sftw

N (со) [N (со) +

M (со)]

 

166

Это выражение совпадает с получающимся при класси­

ческом рассмотрении [31]. Множители й ш

могут быть

добавлены к N (со) и M (со), в результате

чего перейдем

к измерению мощности сигнала и шума в абсолютных единицах. Квантовые эффекты в этом случае совершен­ но несущественны.

Аналогично ведет себя второе слагаемое

в (5.4.3).

При

А/(ш)-СІ и

\u(t)2\M(со)

<С 1 оно имеет вид

 

 

со

со

 

 

 

j Ä j | « ( / ) | a A f ( c u ) û f a > ,

 

 

 

—со

О

 

 

так

что отношение правдоподобия полностью

совпадает

с получающимся

для пуассоновского потока

(см. гл. 3).

При А/(со)3>1 это слагаемое совпадает с получающимся при классическом рассмотрении.

Таким образом, отличие временной и спектральной обработок в квантовом случае по сравнению с классиче-

-скими сводится к видоизменению характеристики филь­ тра. Пространственная обработка для одиночного точеч­ ного источника не отличается от классической.

При рассмотрении совокупности точечных источников или протяженного источника первостепенными становят­ ся требования одновременной измеримости сигналов от отдельных элементов источника. Эта проблема не воз­ никает, если сигналы от отдельных элементов ортого­ нальны. Если элементарные сигналы элемента источника не разрешены по модуляции, то условие ортогонально­ сти можно записать в виде

Ц?І ~

Р*.«°) = 4 - j Y / i

( P ^ e * ) r d r œ ЬІК,

(5.4.7)

 

 

s

 

 

где Pj орты направлений

на отдельные элементы ис­

точника (здесь считаем,

что

можно либо

пренебречь

сферичностью

волн, либо

учесть ее так, как это сделано

в (1.4.15)).

 

 

 

 

Без наложения условий ортогональности задача син­ теза сильно усложняется. В дальнейшем будем рассма­ тривать приемник с набором расстроенных угловых ка­

налов, удовлетворяющих условиям (5.4.7), в

качестве

квазиоптимального во всех случаях, включая

случай

протяженной цели.

 

167

Заметим, что для прямоугольной апертуры совокуп­ ность непродетектированных выходных сигналов таких приемных каналов при максимально плотном заполне­ нии ими пространства р образует достаточную статисти­ ку. Это следует непосредственно из теоремы Котельникова. Максимально плотное заполнение получается при расположении каналов в узлах прямоугольной «решет­ ки» с расстоянием (в угловой мере) между линиями X/d, где d— размер соответствующей стороны апертуры.

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

положении каналов

в центрах

и вершинах

правильных

шестиугольников, заполняющих

картинную

плоскость,

со сторонами \,22\/d

(d — диаметр апертуры).

Проанализируем оптимальный способ объединения сигналов в угловых каналах при наблюдении заданной совокупности точечных источников, сигналы которых флюктуируют по закону Гаусса. Будем считать сначала сигнал медленно флюктуирующим и достаточно узкопо­ лосным, чтобы пренебречь изменением спектральной плотности фона в полосе сигнала, и полагать ехрг/фг*» »=ехр icûopr/c для всех рассматриваемых р и г. Как сле­ дует из (5.2.7), (5.3.5), отношение правдоподобия в этом случае можно записать в виде

А і У № е х Р { ^ п < ; + : т л і ; ' ' -

- і » ( і + г т Ѵ ) } '

(5 -4 -8 )

где (без учета сферичности волн)

 

J «К j y{x.t) е ~ "

p3-r dr

1

 

—со S

 

со

Sj|-u(f)|*di

0 0

Mj — среднее число квантов

полезного

сигнала от /-го

элемента; N—N(®o)—среднее

число

квантов фона на

одну степень свободы. Для случаев малых и больших чисел квантов из (5.4.8) получаем предельные соотно-

168

шеиия того же типа, что и для рассмотренного

случая

быстрых

флюктуации: при М<^1, N <g;l

способ

объеди­

нения tij

приближается

к типичному для

пуассоновского,

а

при N^>1—-для

классического случая.

 

 

 

При

протяженной

цели Mj = \x{ßj)

=M(pj)X\/S,

где

M

(p)/S — средняя

плотность

потока

квантов, приходя­

щихся на единицу

телесного

угла.

Вводя формально

угловую

плотность

числа

зарегистрированных

квантов

 

 

 

V(P) =

£ S ( P - K < ) ,

 

 

(5.4.9)

где ні — угловые (направления, на которых были

зареги­

стрированы кванты, можно записать (5.4.8) для протя­ женной цели в виде

Л (г,

(5.4.10)

При такой форме записи результаты для случаев ма­ лого и большого числа квантов (Л/•< 1, u.<С 1 и Л / > 1 ) пол­ ностью совпадают с результатами для пуассоновского и для классического случаев соответственно.

 

Совершенно

аналогичный

результат

получаем для

протяженного

источника

прибыстрых

флюктуациях:

в

(5.4.3)

добавляется

суммирование

по /, a M (со)

в

(5.4.3)

и (5.4.4) заменяется

на M (р, со).

5.5. Оптимальная обработка поля при наличии фазовых искажений

Определим, как.меняется оптимальная обработка по­ ля при наличии фазовых искажений поля на апертуре. Эти искажения могут быть обусловлены рассеянием на турбулентности в среде (атмосфере, например, см. гл. 1) или погрешностями оптической системы. Известно, что такие искажения существенно снижают разрешающую способность оптических приборов при обычно используе­ мых в этих приборах способах обработки поля.. Синтези­ рованный в {92] для классических полей способ обработ­ ки позволяет, в принципе, обеопечить разрешающую спо-

12—220

169

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