
книги из ГПНТБ / Курикша, А. А. Квантовая оптика и оптическая локация (статистическая теория)
.pdfгде |
|
|
|
|
«Р (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 входит |
сумма /гѵ /йѵ . |
Разли |
|||
чие, очевидно, |
исчезает |
при |
rä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