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

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

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

где

 

 

 

 

«Р (8|2) = )П (Я) I г) = £ * (фс) |<л-|г> |2 ,

(5.1.7)

 

 

X

 

 

Р ( г ) =

Jp« (Я)р(г|Я)гіЯ,

 

 

^(S) = ^

j p n « P № ) r ( S , Х)гіЯ.

(5.1.8)

Минимум среднего риска должен обеспечиваться вы­

бором л- и

к(6\х).

 

 

Легко видеть, что 0<;cp(ô|z)<; 1 и что минимум

риска

получается

тогда, когда при заданном z функция cp (ô | z) =

= 1 для того б, при котором тг(<5) минимальна, и

равна

нулю для остальных о. Это требование выполняется, если

x=z и л(б|г,) удовлетворяет условиям, сформулирован­ ным для ф(б|г) . Следовательно, измеряя z, можно по­ лучить минимум риска, и, следовательно, z является до­ статочной статистикой.

Этот же результат можно получить из следующих простых соображений [88].

Воспользуемся статистическим толкованием операто­ ра плотности (5.1.5), т. е. будем считать формально, что смешанное состояние получено из множества чистых со­

стояний

с

операторами плотности

\z >

(z\

случайным

выбором

с

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

p(z\l)

для

исходов. Прини­

мая это

толкование,

считаем,

что

система

находится

в одном

из чистых

состояний

|z>

< z\

и

лучшее, что

можно

получить, измеряя

параметры

системы, — это

полностью определить ее состояние. Для этого достаточ­ но измерить набор величин z, являющийся, следователь­ но, достаточной статистикой. Эту статистику можно за­ тем использовать в классической процедуре теории решений.

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

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

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

Обсуждаемое требование не выполняется в важной для практики задаче разделения сигналов от близкорас­ положенных целей, решению которой для классических сигналов уделялось и продолжает уделяться большое внимание. В квантовомеханическом случае пока не уда­ лось найти оптимальный приемник для этой задачи. Возможные постановки задачи синтеза и ограничения, связанные с кватовой структурой поля, будут обсужде­ ны в § 5.5.

Итак, описаны общие • подходы к задаче проверки статистических гипотез о квантовомеханической системе. В случае, если такой системой является поле, заполняю­ щее, вообще говоря, все пространство, немаловажно и то, какая часть поля считается доступной для наблюде­ ния. В работе К- В. Хелстрома ,[86] в качестве доступной для наблюдений системы рассматривалось поле в огра­ ниченном объеме — ловушка с затвором со стороны входной апертуры. В этой части подход К. В. Хелстрома близок к использованному П. А. Бакутом в (84, 85].

В работе автора ,[88] доступным для наблюдения счи­ талось изменяющееся поле на плоской апертуре. Эта ситуация ближе, по-видимому, к имеющимся на прак­ тике. Такой же случай был рассмотрен и К. В. Хелстромом в [90].

Результаты синтеза оптимального приемника для различных вариантов наблюдения поля (в объеме, на апертуре и т. д.), разумеется, получаются различными. Однако для поля в ловушке и на апертуре эти результа­ ты, как справедливо отмечено в (91], могут быть преобра­ зованы друг в друга, поскольку эти поля связаны взаим­ нооднозначно.

Процедура синтеза во всех случаях содержит два этапа. Сначала по оператору плотности поля во всем пространстве необходимо построить оператор плотности доступной для наблюдения части поля (подсистемы по отношению к полю), а затем .применить к подсистеме опи-

151

санную процедуру нахождения оптимального способа принятия решения.

Для большей наглядности сначала рассмотрим более простой случай дискретной совокупности мод, а затем задачу оптимальной регистрации поля на апертуре.

5.2. Оптимальный приемник для дискретной совокупности мод

Рассмотрим сначала случай, при котором исследуе­ мое поле представляет собой дискретную совокупность конечного числа (т.) мод. Будем считать, как и раньше, что оператор плотности поля представлен в виде

Р = ] . . ^ Р ( Ы ) П І « І ) Ы * Ч

С5 -2 -1 )

где |ш ) — вектор

состояния у-й моды с

заданной ам­

плитудой о,-; p({aj})

—весовая функция

Р-представле-

ния, совпадающая с классической функцией распреде­ ления для комплексных амплитуд; интегрирование про­ изводится по комплексным плоскостям CZJ.

Хотя по форме представление (5.2.1) диагонально,

совокупность амплитуд {а3} не

является искомой

доста­

точной статистикой. Операторы

a,j не эрмитовы, a

p({a.j})

нельзя рассматривать как вероятностную меру на мно­ жестве состояний |{«j}> из-за их неортогональности. Нужно искать диагонализирующее представление, свя­ занное с собственными функциями каких-либо эрмито­ вых операторов.

Рассмотрим два случая: гауссова сигнала и регуляр­ ного сигнала со случайной фазой при наличии аддитив­ ной гауссовой помехи с некоррелированными значениями

амплитуд

ctj.

 

 

 

 

 

 

 

В первом

случае

 

 

 

 

 

Р({аМ

=

» » d e i n e n е х Р

J -

S

« W * J .

(5-2.2)

где

— матрица,

обратная корреляционной

матрице

 

 

Rjk

=

(аV , ) =

Nfa

+

Rjhc,

(5.2.3)

Nj — среднее число квантов помехи в у'-й моде, ||i?jftcll— корреляционная мятрица полезного сигнала.

152

Перейдем от комплексных амплитуд d,- к новому на­ бору амплитуд ß„ с помощью унитарного линейного пре­ образования, матрица которого ЦУ^ІІ составлена из соб­

ственных векторов

матрицы \\Rji,\\. В § 1.2 было показа­

но, что величины

>

' і

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

Операторы bv, b*. соответствующие амплитудам ß*v, ßv, обладают всеми свойствами операторов уничтожения

и

рождения.

Собственные

векторы

операторов

(а,}

являются

одновременно

собственными

векторами

операторов ѵ} и могут быть представлены в виде про­ изведений собственных векторов |ßv ):

|{«*}>=П|«*)=1Ш>=П|р,>.

При

такой замене переменных

оператор

плотности

(5.2.1)

с весовой функцией Р({щ})

вида (5.2.2)

преобра­

зуется в произведение операторов плотности для типов колебаний с амплитудами ßv

Р =

І Г Р , = П Г - ^ -

exp

v )<ß„Kßv . (5.2.4)

где nv

— собственные

значения

матрицы Ц ^ Ц . Диагона-

лизация р получается

диагонализацией сомножителей"р„.

Подставляя в (5.2.4) разложение (1.2.15) |ßv ) по соб­ ственным векторам оператора числа квантов

=ѴѢ = E vhv*kffîu,

(5.2.5)

iі..'.ft

 

получаем

11—220

153

 

К" К"

"С + "ѵ )

!J]Pv(")l'0(«J- (5.2.6)

Согласно (5.2.6) оператор рѵ днагонален в представ­ лении чисел заполнения. Таким образом, для рассматри­ ваемого случая достаточной статистикой является сово­ купность чисел заполнения /?.ѵ для комплексных типов

колебаний с амплитудами ßv , связанными с амплитудами üj колебаний в исходных модах унитарным линейным преобразованием, диагоиализирующим корреляционную матрицу.

Совокупность {/гѵ} является, очевидно, достаточной

статистикой для принятия решения и тогда, когда допу­ скаемым гипотезам соответствуют гауссовы поля с ком­

мутирующими корреляционными

матрицами ||/?jh(X,)ll

(À —индекс гипотезы), поскольку

эти матрицы имеют

общую систему собственных векторов. Логарифм функ­ ции правдоподобия для совокупности измеренных вели­

чин {/гѵ|

есть

 

 

 

 

 

 

 

 

 

 

 

v *

 

 

 

 

 

 

— In [1 - I - Яѵ (A)] j — m In it.

(5.2.7)

Формула

(5.2.7)

определяет

способ объединения из­

меренных

величин

ѵ. Этот способ

отличается от

клас­

сического,

при

котором

в L входит

сумма /гѵ ѵ .

Разли­

чие, очевидно,

исчезает

при

v>l.

 

 

Рассмотрим теперь случай, при котором каждая мо­ да представляет собой суперпозицию регулярного сиг­ нала и гауссова фона, причем сигнал во всех модах имеет общую случайную фазу с равномерным распреде­ лением в интервале (0; 2л). Для упрощения предполо­ жим, что интенсивность фона для всех рассматриваемых мод одинакова. Реально это предположение обычно вы-

154

полняется из-за относительной узкополосности регуляр­ ных сигналов. В этом случае

/ > ( М ) = / .

 

 

1

+

\ ІІ

 

 

 

N

 

 

 

 

 

 

(5.2.8)

где N, как и в

(5.2.3), средняя интенсивность

фона (в чи­

слах квантов).

 

 

 

 

 

Перейдем к новым переменным

 

 

ßv =

2

1

 

(5-2-9)

считая SJl ) :=tj,

а

остальные ^.ѵ )

выбирая

из условий

взаимной ортогональности. При этом (5.2.7) преобразует­ ся в произведение распределений для различных ßv, при­ чем каждый сомножитель будет зависеть только от |ßj . Применяя в каждом сомножителе разложение | ßv> по |/гѵ> , легко убедиться, что если p(ßv ) зависит только от Iß,|, то достаточной статистикой является совокупность величин «v = |ßJ 2 . Этот результат верен и для совокуп­ ности m взаимно ортоганальных сигналов, что следует из вида использованного преобразования (5.2.9).

Распределение pv(/z) приобретает вид

о

у e-iei'rfШ |2

L ( 0 )

/

 

) е х

р / _ _

^ _ \ t

Л е

а\Щ

( 1 + л О " - и Ln

\

/Ѵ(/Ѵ+1)У

Р \

N+lj

 

 

 

 

 

 

 

(5.2.10)

где

L^0 ) (х) — полином Лагерра;

Му =

|

|2.

 

Это распределение уже было рассмотрено

в гл. 2

[см.

(2.2.1)].

При М^>1

и

М » 1

оно

превращается

в обобщенное релеевское, а при ІѴ«СІ в

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

Произведение распределений рч ѵ ) определяет опти­ мальный способ объединения результатов измерений «ѵ

в тех случаях, если принимаемое решение касается всех сигналов в совокупности.

I I *

155

В рассмотренном случае доступными для наблюдения являлись амплитуды {щ}. Большой интерес 'представляет отыскание оптимальной обработки, когда наблюдается суперпозиция мод, представляющая поле в определенной области пространства например:

y(s) = IlCj(s)aj,

s^S0.

(5.2.11)

i

Поле в области 5 о не определяет {а;} однозначно и ха­ рактеризует некоторую подсистему поля в целом. Пер­ вым шагом должно быть получение оператора плотности для этой подсистемы.

Пусть У Ѵ (s)—ортогональная система функций в S0. Раз­ ложим у (s) по этой системе (считая, что это возможно):

 

 

= £ Т А (*)•

( 5 - 2 - 1 2 )

 

 

V

 

ТѴ =

^

[ y(s)v\{s)ds

=

 

 

•So

 

/

S0

j

 

где Яѵ — нормировочный коэффициент.

Потребуем, чтобы операторы*'уѵ , у* подчинялись та­ ким же коммутационным соотношениям, как и Oj, ay:

ІТѴ .

YjT] =

S C*Vj ^^ftôjft =

 

 

 

/.ft

 

= " W И K { s >

' S s ) U ï ( S l ) v * » ( S : ) d S l d s ° = S v | i '

( 5 ' 2 - 1 4 )

V ^ S,,

 

 

 

где

 

С*з (s.) C3 (s2).

 

К (s„ sa) = £

(5.2.15)

*' Напомним, что соглаоно введенным в гл. 1, обозначениям оператору âj соответствуют собственные значения a*j, поэтому в опе­

раторном аналоге соотношений (5.2.11), (5.2.12) коэффициенты за­ менены на комплексно-сопряженные.

156

Условие (5.2.14) выполняется, если vv(s)

определяют­

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

 

 

(5.2.16)

т. е. являются собственными функциями

интегрального

оператора с ядром K(si, sz). Разложение

(5.2.12) по

этим функциям возможно в силу теоремы

Гильберта —

Шмидта [23]. Заметим, что при этих условиях амплитуды Yv оказываются некоррелированными, если были не коррелированы амплитуды а,-.

Дополним систему функций vv (s) до полной ортого­ нальной системы функций во всем пространстве (дока­ зательство возможности дополнения см. в (3]). Такое построение можно осуществить, введя ортогональную си­ стему функций, полную по отношению к функциям, представимым истокообразно с ядром K(si, S2) в подпро­ странстве So, дополняющем 50 до всего пространства, и полагая эти функции равными нулю на So, а функции

У ѵ (s) считая равными нулю на So. Если подпространст­ ва S0 и So имеют разную размерность, то vv (s) следует домножить на о-функцию с аргументом, обращающимся в нуль на So. Полученную таким образом дополнитель­ ную систему функций обозначим {г^0> (s)}, а коэффициен­ ты в разложении поля по этим функциям — 7'°'.

Оператор плотности

полл в Sg получаем из

(5.2.1),

переходя к новой

системе мод и беря по всем уѵ ( 0 )

I - J p .

({ъ})Ц\ъ)(ъ\а2ъ,

(5.2.17)

где

 

 

 

p.({r.})=J

- 5 Р ( { « Л { Т Л Г } ] } ) П ^ -

< 5 - 2 - 1 8 )

Таким образом, весовую функцию р . ({YJ) в новом представлении получаем из имеющейся p({<*j}) по обыч­ ным правилам преобразования плоти 'сти вероятностей. Дальнейшая процедура синтеза не отличается от рас­ смотренной в начале данного параграфа.

157

Операции формирования достаточной статистики удобно представить в виде преобразований наблюдаемо­ го поля у (s). Согласно (5.2.5), (5.2.13) элементы доста­ точной статистики для гауссова сигнала есть /?.ѵ = оѵ + Ьѵ - где

К =

Yl V

' =

S

°/» TT \ Î J

( s ) ü

j { s )

d s =

\ V

* M

d s -

 

i

 

 

 

i

 

s0

 

 

 

s0

 

(5.2.19)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для

опорного

сигнала

Vv (s)

получаем

следующее

уравнение:

 

 

 

 

 

 

 

 

 

 

 

 

 

f Я (5, 5') К

ѵ

 

 

 

J

j »,

(5)

ft

 

 

ift

 

 

 

( 5 ' ) els' =/,ft,/ s0

o*

(s') /?*

 

- i -

^ (s') ds' =

 

 

 

 

=

n,

f/C(s,s')Vv (0^',

 

 

 

(5.2.20)

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

(s,

s')

=

(y* (s)

y (s/)) ;

R i h

=

(у*Д/<>;

 

К (s, s')

определено

(5.2.15).

 

 

 

 

 

 

 

Уравнение

(5.2.20) определяет

такое

 

преобразование

п оля, при котором диагонализируются одновременно кор­

реляционная матрица

поля и матрица

коммутаторов

v ,

b*]. Нормировка

Vv (s) задана

условием

 

 

 

 

[6 , Ь*] =

J f К (s,,S = ) V * v

(s,) V v (s3)

ds.ds, =

1.

 

Для

совокупности

взаимно

ортогональных

(с весом

L(sits2))

регулярных сигналов

{tVv (s)} со случайными

фа­

зами достаточной статистикой является совокупность величин /гѵ = pvps,=v, где

j j L * (s,, s2 ) L / * V (S,) (/ (sa) ds,ds2

1/2 (5.2.21)

- S0

158

L (s,, s,) = £ Vj (s,) v*j (s2)/Ä~ — ядро

интегрального опера-

i

Kis^s*).

тора, обратного оператору с ядром

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

пользование приближений

может

привести

к на­

рушению условий одновременной

измеримости

вели­

чин, образующих достаточную

статистику. В связи с этим

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

Опровергнуть это утверждение, конечно, невозможно. Можно лишь противопоставить его категоричности сле­ дующие соображения. Разумеется, не удается построить прибор, измеряющий одновременно неизмеримые величи­ ны. Поэтому фактически прибор, построенный для изме­ рения величин, выбранных по'результатам приближенно­ го решения задачи синтеза, будет измерять не их, а не­ которые другие одновременно измеримые величины. Эти

последние будут в той или иной мере

(в зависимости от

степени приближения)

близки

к тем,

которые

получи­

лись бы при точном решении задачи.

 

 

Это замечание должно в какой-то мере оправдать те

приближения,

которые

будут

использованы в

.§ 5.3.

5.3.

Достаточные статистики

поля

 

на

апертуре

 

 

 

 

Применим полученные результаты к задаче отыска­ ния оптимальной обработки поля на апертуре. В каче­

стве

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

у(г.

ш) для всех со>0 и T Œ S (S— площадь апертуры).

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

159

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