книги из ГПНТБ / Курикша, А. А. Квантовая оптика и оптическая локация (статистическая теория)
.pdfдля выяснения основных особенностей оптимальной об
работки при рассматриваемых здесь свойствах |
сигналов |
||||
и помех. |
|
|
|
|
|
Спектр поля |
на апертуре будем представлять в виде |
||||
суперпозиции |
плоских |
волн |
(см. §. 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 |
= |
Sù |
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
со |
|
|
|
|
|
где |
|
|
|
|
|
|
|
|
|
(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 |
