книги из ГПНТБ / Курикша, А. А. Квантовая оптика и оптическая локация (статистическая теория)
.pdfплотности поля, являющегося суперпозицией регулярно го сигнала и гауссовой помехи, также Р-представим.
Для Р-представимых операторов плотности задача статистического описания поля сводится к отысканию классического функционала распределения вероятностей для амплитуд монохроматических плоских волн а(р, ш). С этим функционалом взаимооднозначно связан харак
теристический |
функционал *>, определяемый выражением |
Ф h (s)] = |
(exp {i j [ч * (s) a (s) + 7) (s) a* (s)} ds), (1.4.12) |
где через s для краткости обозначена совокупность пе
ременных |
р, ш, а |
скобки |
( ) означают усреднение |
по |
|||||
всем a (s). |
|
|
|
|
|
|
|
||
Характеристический |
функционал удобно |
представить |
|||||||
как |
среднее |
значение |
соответствующего |
оператора. |
|||||
Естественно определить этот оператор, заменив a(s) |
и |
||||||||
a*(s) |
на |
â+(s) |
и â(s) |
соответственно. При |
этом |
чтобы |
|||
для |
Р-представимых операторов плотности результат |
||||||||
был |
тождественен |
(1.4.12), нужно, как нетрудно |
убе |
диться, взять нормальную форму получающегося опе ратора:
Ф [ Ч ( в ) ] - < е ^ ( ' і Г ; ! + ( ' ) ' ' Ѵ ^ ( ' ) Л > . |
(1.4.13) |
Сформулируем теперь более частные |
допущения |
о свойствах сигнала в рамках его классического описа
ния |
с помощью |
распределения р{а). В качестве доступ |
||
ного |
для наблюдения |
сигнала будем |
рассматривать |
|
поле |
на плоской апертуре регистрирующего прибора, |
|||
представляющее |
собой |
аддитивную |
смесь фонового |
излучения и полезного сигнала.
Фоновое излучение будем представлять в виде супер позиции некоррелированных плоских волн с равномер ным угловым распределением- интенсивности в пределах рассматриваемого сектора. При рассмотрении поля на апертуре естественно учитывать только составляющие электромагнитного поля, лежащие в плоскости аперту ры, поскольку только они определяют поток мощности через апертуру. Корреляционные функции для этих со
ставляющих |
электрического |
и магнитного полей |
имеют |
|
*' Это совместный |
'-характеристический функционал для |
Rea(s) |
||
[с аргументом |
Rer|(s)] |
и lmï)(s) |
[с аргументом Inrri(s)]. |
|
50
весьма сложный вид *>. Чтобы не усложнять задачу, будем использовать скалярную модель фонового излу чения с функцией корреляции вида
(У (г,, ш.) у* (га . ш2)) = 2иУѴф К ) |
| е ' : / г , ( |
Г і ^ Г 2 ) р cos Ыр X |
|
X s к - « g = 2 ^ ф К) 8 к - |
U . J ^ ^ Ц - р ^ , |
||
|
|
|
(1.4.14) |
где M |
)—спектральная плотность |
потока мощности, |
J,(CÛ
приходящегося на единичный телесный угол. Множи тель cos Ѳ в этой формуле введен для учета ослабления вклада в поток мощности за апертурой 'плоских волн, приходящих под большими углами к нормали. Физи чески такое ослабление должно происходить, однако его учет в модели фонового излучения для дальнейшего рассмотрения не является принципиально необходимым.
Заметим, что величина Л^ф(со)А,2 (\к=2л/к— длина волны) может быть истолкована как средняя энергия элементарной плоской волны, т. е. энергия, приходя щаяся на одну степень свободы электромагнитного по ля. В этом легко убедиться, рассматривая сначала поле в конечном объеме с какими-либо граничными условия ми, а затем устремляя размеры объема к бесконечности. Отсюда следует, в частности, что для равновесного излучения
ЛГф (ш): |
Яш |
1 |
|
|
где %(Ù — энергия кванта, Ѳ—температура (в энергети ческих единицах).
Для простоты и единообразия рассмотрения весьма заманчиво считать полезный сигнал также суперпози цией независимых плоских волн. Это допущение спра ведливо в очень многих важных для практики случаях. Преждевсего, это относится к классическому случаю точечного источника излучения и апертуры, размеры которой малы_по сравнению с размерами -первой зоны Френеля ]/*Ä,r (Л- — длина волны, г — расстояние от источника до центра апертуры). Условие это обычно выполнено в астрономии, однако при наблюдении за
*' Эти корреляционные функции имеют такой же вид, как и для теплового излучения полупространства [17].
4* |
51 |
объектами и приеме сообщений на Земле или в около земном космосе оно может нарушаться. Учет сферич ности волны, излученной точечным источником, незна чительно усложняет теорию и не является поэтому обременительным.
Если источник можно представить в виде совокупно сти точечных излучателей, расположенных близко друг •к другу, то для суммарного поля на апертуре (расстоя ние от центра апертуры до источника велико) справед ливо следующее приближение:
у ( Г , ю ) = = 2 а ^ ^ А Ä |
£ Ѵ г - ( Р о г ) ^ ° s a f i l k p i r |
|
|
i |
i |
|
UJ — комплексная |
(1.4.15) |
где |
амплитуда волны у-го излучателя |
|
(с |
учетом ослабления |
при распространении); г,-— ра |
диус-вектор, проведенный из центра апертуры к у'-му
излучателю; р^=—rj/г,-; |
р0 — орт, характеризующий на |
||
правление на условный |
центр источника; |
г0 — расстоя |
|
ние до него. |
|
|
|
Интерференция |
волн |
от элементарных |
излучателей |
здесь происходит |
так же, как и на очень |
больших уда- |
-лениях, когда волны плоские, а сферичность волны 'про является так, как если бы источник был точечным. При переходе к непрерывному распределению излучателей сумму в (1.4.15) следует заменить интегралом по на правлениям прихода волн.
Втипичном случае некогерентного излучения или рассеяния (например, рассеяние сильно шероховатой поверхностью) волны элементарных излучателей имеют независимые случайные фазы, принимающие любые зна
чения из интервала 0—2л с одинаковой вероятностью. Суперпозицию большого числа таких волн можно счи тать подчиняющейся гауссову распределению с нулевым средним значением и функцией корреляции
R (г,, г„, m) = (yjr,, ш) у* (r2 > со)) =
= S i * j I2 exp {lk [\ г, - г,I - I r 2 - г,|]} « i
~ e x P { - 1 - £ r [ r ' - r 2 - ( p ° r ' ) 2 + (p °r *)2 l} X
X S l ^ l ' e x p I t A p i f r - r J ] . |
(1.4.16) |
/ |
|
52
Моменты поля с неодинаковым числом перемножаемых значений у я у* будут, очевидно, равны нулю согласно предположению о свойствах фаз элементарных излу чателей.
Переходя к непрерывному распределению интенсив ности сигнала по направлениям, получаем
где I (р) — плотность потока энергии, приходящей из единичного телесного угла.
Из полученных формул видно, что при рассмотрении статистических характеристик поля, полностью опреде ляемых в гауссовом приближении функцией корреляции, условия представимости поля в виде суперпозиции ста тистически независимых плоских волн становятся менее жесткими, чем при рассмотрении реализаций. Множи тель перед интегралом в (1.4.17) можно отбросить, если
показатель экспоненты |
мал всюду |
на |
апертуре |
для |
|||
—г2 | порядка радиуса -корреляции поля. Этот радиус |
|||||||
определяется |
быстротой спадания |
интеграла |
с ростом |
||||
|гі—Г21 и, поскольку |
интеграл является |
преобразова |
|||||
нием |
Фурье |
углового |
распределения |
интенсивности, |
|||
имеет величину порядка |
ѴѲо, где Ѳо — характерный |
угло |
|||||
вой |
размер |
источника. |
С учетом |
сказанного |
сферич |
||
ностью волн можно пренебречь, если |
|
|
|
|
|||
|
|
|
|
|
|
(1.4.18) |
где D — диаметр апертуры, L — линейный размер источ ника.
Это условие обычно выполняется в оптической лока ции. Если размер апертуры имеет величину порядка А./Ѳо, которая требуется для определения очертания и разме ров источника, то условие (1.4.18) можно заменить следующим:-
Ѳ0 > ѴЩг0. |
(1.4.19) |
Для устройств связи при использовании в излучаю щей системе коллимации светового пучка или ОКГ представление поля в виде суперпозиции независимых плоских волн оправданно в еще большей степени, чем для совокупности точечных излучателей.
53
На разрешающую способность и точностные харак теристики оптических систем существенно влияют флюк туации, связанные с рассеянием света на турбулентных неоднородностях среды. Эти флюктуации 'представляют интерес для ряда приложений оптики (астрономия, гео дезия и т. д.). Не ставя задачу полного обзора работ в этой области, ограничимся описанием свойств и .коли чественных характеристик флюктуации.
Флюктуации фазы ср(г) поля на апертуре, располо женной перпендикулярно направлению распространения волны, можно рассматривать как гауссово случайное поле со стационарными приращениями. Если турбулент ность подчиняется закону Колмогорова — Обухова, то структурную функцию фазы при плавном изменении параметров турбулентности можно записать в -виде [18]
A, (I |
г, - г , |
I) |
<? (УГ> |
« |
ЗА» |
I |
|
|
І X |
||
|
|
|
|
|
Г, - |
Г , |
5/3 |
||||
|
|
|
|
X\C\(z)dz, |
|
|
|
|
|
(1.4.20) |
|
|
|
|
|
о |
|
|
|
|
|
|
|
где k—2n/% |
(% — длина |
волны); L — длина |
пути, про |
||||||||
ходимого плоской волной |
в турбулентной |
атмосфере; |
|||||||||
Сп — параметр, |
характеризующий |
глубину |
флюктуации |
||||||||
показателя |
преломления |
(вблизи |
земной |
поверхности |
|||||||
С п ~ 2 - 1 0 - 8 |
см1 '»). |
|
|
|
|
|
|
|
|||
Для приземной трассы формулу (1.4.19) можно упро |
|||||||||||
стить, считая |
|
Cn (2)=const. Для |
наклонной |
трассы |
|||||||
с углом места # = 7 н-10° и более |
|
|
|
|
|
|
|||||
|
|
г і - |
|
r21) ~ 10-°/г2| г, — г2 15 / 2 |
sec 9-, |
|
(1.4.21) |
АДІ
где X и /' должны быть измерены в сантиметрах. Фазовые флюктуации приводят к ошибкам измере
ния угла прихода волны. Если определить измеряемое значение угла как средний по апертуре наклон фронта волны по отношению к плоскости апертуры, то средне квадратичное значение ошибки, обусловленной фазовы
ми флюктуациям'И, есть [18] |
|
|
|
ol~D9(d)/kacr, |
(1.4.22) |
где d— размер |
апертуры в направлении |
отсчета угла. |
На приземной трассе длиной в 1 км среднеквадра |
||
тичная ошибка |
может достигать 2 + 7". |
На наклонной |
54
трассе [20] днем при & = |
90°, |
ы, |
сФ « |
1-4-2". |
|
При фазовых флюктуацнях существенно снижается разрешающая способность приемных систем с большой апертурой. Приближенно можно считать, что при доста точно больших размерах апертуры разрешающая спо собность определяется радиусом г„ когерентности поля, определяемым из условия
по формуле iAlcpi~<À/2r„. Экспериментальные измерения разрешающей способности [21] в видимом диапазоне
дают значения Лісрі в интервале |
от долей угловой секунды |
до десяти секунд в зависимости |
от состояния атмосферы. |
Например, для Паломарской обсерватории наиболее ве
роятное значение Дср=2", что при |
Х = 0,5 миом соответст |
|
вует |
/'к— 2,5 см. |
|
Интерференция волн, рассеянных на иеоднородно- |
||
стях, |
приводит к амплитудным |
флюктуациям поля. |
Радиус корреляции этих флюктуации равен примерно размеру золы Френеля. Закон 'распределения іна доста точно коротких трассах близок к логарифмически нор мальному і[22], а дисперсия оыі2 логарифма интенсив ности растет пропорционально расстоянию, как это следует из теории, основанной на методе плавных воз мущений. На больших расстояниях от источника экспе
риментальная |
зависимость аыі2Щ |
насыщается |
и даль |
|
ше начинает |
медленно спадать. Максимум |
a i n / 2 |
s £ l , что |
|
соответствует |
относительной дисперсии |
интенсивности |
^1,7. Для наклонной трассы [19, 23] получены значения
оъ.г~ 0,4 |
(cTj/I Ä= 0,4). |
Предпринят ряд Попыток теорети |
|
ческого описания этих |
закономерностей. |
|
|
При |
больших размерах апертуры |
амплитудные |
флюктуации принимаемого сигнала ослабляются благо даря усреднению по апертуре. Поэтому в дальнейшем при рассмотрении приемников, собирающих энергию со
всей |
апертуры, эти |
флюктуации не будем |
учитывать. |
В |
заключение |
кратко остановимся на |
характере |
флюктуации поля в фокусе объектива, связанных с фа зовыми искажениями на апертуре ![24]. Кроме атмосфер ного рассеяния причиной этих искажений могут быть неточности изготовления оптической системы. Считая
55
флюктуации фазы, подчиняющиеся гауссову распреде лению, легко получить выражение для функции корре ляции тюля в фокусе объектива
<£*Ф ( Р і ) Еф (р,)) ^ |
j j |
Е* (г,) Е (г2 ) ехр |
о ф (|г , |
- г2 J) - ik (р, - |
р |
2 ) г, - ikp2 (г, - r2 ) I rfr.rfr,, (1.4.23) |
|
где £ ( г ) —поле на |
апертуре, р = Г ф / / , |
г ф — радиус-вектор |
|
в фокальной плоскости, /—фокусное |
расстояние. |
Рассматривая (1.4.22) при условии, что радиус коге рентности искаженного поля на апертуре мал по срав нению с размерами апертуры, легко показать, что ра диус фокального пятна от точечной цели связан с ра
диусом когерентности |
соотношением |
ро=''ф//~Я/гк, |
|||
а |
радиус корреляции поля в фокальной |
плоскости |
равен |
||
pK |
= 2yd (d— диаметр апертуры). |
|
|
|
|
|
Представляя аналогично (1.4.22) функцию взаимной |
||||
корреляции полей в фокусе объектива |
от двух удален |
||||
ных точечных целей с угловым расстоянием |
іДѲ, нетруд |
||||
но |
показать, что при |
указанных условиях |
эти |
поля |
в каждой точке фокальной плоскости некоррелированы, если Следовательно, сигналы от целей, разре шаемых идеальной оптической системой, несмотря на наличие искажений, суммируются в фокальной плоско сти по интенсивности. Указанные свойства флюктуации будут использованы в дальнейшем.
ОХАРАКТЕРИСТИКИ ПРИЕМНИКОВ, ФОРМИРУЮЩИХ КВАДРАТИЧНЫЕ ФУНКЦИОНАЛЫ ПОЛЯ
Втипичном оптическом приемнике принимаемый сиг нал— поле, попадающее на апертуру оптической систе мы, подвергается следующим преобразованиям. Снача ла системой линз, зеркал, призм, светофильтров произ водится пространственная и спектральная фильтрация сигнала, описываемая линейным интегральным 'преобра зованием входного поля по координатам и по .времени. Затем сигнал поступает на фотодетектор, преобразую щий мощность светового потока, в . электрическийток.
56
Упрощенно преобразование поля фотодетектором можно
представить (более |
детальный анализ будет |
проведен |
|
в гл. 3) как формирование квадрата |
модуля |
комплекс |
|
ной амплитуды поля |
с 'последующей |
линейной |
фильтра |
цией, учитьшающей инерционность фотодетѳктора. Очень часто после детектирования сигнал подвергается допол
нительной линейной |
фильтрации для накопления |
полез |
||||
ной составляющей, |
выделения модуляции светового по |
|||||
тока и т. д. |
|
|
|
|
|
|
Описанная |
последовательность преобразований мо |
|||||
жет быть записана |
в форме квадратичного функционала |
|||||
(интегрального |
аналога |
'квадратичной формы) |
от |
поля |
||
на 'входе приемника. Представление |
выходного |
сигнала |
||||
в виде квадратичного |
функционала |
справедливо для |
весьма широкого класса 'способов обработки. В данной главе будут рассмотрены характеристики надежности обнаружения и точности измерений параметров сигнала для таких опособов обработки. Это позволит установить общий для идеальных и реальных приемников характер всех основных зависимостей, определяющих качество
приема, |
так |
что в дальнейшем можно |
будет, |
избегая |
|
повторения, |
ссылаться на результаты данной главы. |
||||
|
2.1. Статистические |
характеристики |
|
||
|
квадратичных функционалов |
поля *) |
|
||
Пусть |
{Fj} — совокупность |
квадратичных |
функцио |
налов поля, измеряемых в приемнике. Каждый такой функционал можно записать в виде
S —со
координатам; Vj(u, г2 , h) обладает следующим свой ством:
|
Vj(r2, |
4 , Г І , |
ti)=V*j(rit |
U, |
r2,t2). |
|
Для |
дальнейшего |
рассмотрения |
удобноперейти |
|||
к спектральной записи |
(2.1.1): |
|
|
|||
|
оо |
|
|
|
|
|
F j = |
J J J J ^ Г " |
Ш " Г " |
У |
' |
У * (Г " Ш ^ |
d r i R F R 2 D ( B . D ( D 2 . |
|
5 О |
|
|
|
|
(2-1-2) |
*> Общин подход, •иопользовашый здесь, предложен в [25].
57
где
|
оо |
|
Vi (г„ |
га , fflj = - ^ Ï - j j " v i (r .. |
r=- '2) exp (Kf, - |
|
— 0 0 |
|
|
— mjjdttdt,. |
(2.1.3) |
Формулы (2.1.1) — (2.1.3) охватывают, в частности, случай произвольной пространственно-временной линей ной фильтрации с последующим детектированием (пере ходом к интенсивности поля) при наличии и при отсут ствии последетекторного сглаживания (фильтрации ин тенсивности) .
Если последетекторное сглаживание отсутствует, то
|
|
Vj(ru |
|
tu Га, |
іг)=ё*і(гі, |
ti)gj(r2, |
іг), |
|
|
|
V j O V œ , , r s , œ s ) = - 1 i - f i r * j ( r t , a i 1 ) 5 r J ( r s i m 3 ) , |
(2.1.4) |
|||||
где |
gj(r, |
t)—весовая |
функция |
фильтра, |
gj(r, |
ш) — ее |
||
спектр. |
|
|
|
|
|
|
|
|
|
При |
наличии |
последетекторного сглаживания |
|||||
|
|
|
|
0 0 |
|
|
|
|
Vi |
(г„ t1,ra,tj |
= |
^vj(p,t)g*i(p,r1,t,t1)gi{p,rat, |
|
Q dçdt |
|||
|
|
|
|
— 0 0 |
|
|
|
|
|
|
|
|
|
0 0 |
|
|
|
|
Vi |
( Г „ «D„ Г 2 , (Da) = |
Ç j * О, ( p , 0 g*,- ( p , |
Г, f, U>,) X |
||||
|
|
|
|
X f t f P . r ^ f . œ J r f p A . |
|
(2.1.5) |
где Uj(p, 0 —весовая функция; p—параметр, по кото рому производится последетекторное сглаживание (им •может быть направление при дифракции на апертуре, номер канала, координата в фокусе линзы и т. п.); gj(P> г, t, f) —импульсная переходная функция додетекторного фильтра.
Если параметры этого фильтра постоянны и
gj(P, г, t, fi)=gj(p, г, i—h),
то |
|
Vj (ri, |
Г2, Cû2) = |
= J ü i(P . cùi—cù2)g-*j(p, |
14, —ші)еГз(р, гг , — cû2)rfp, |
|
(2.1.6) |
58
где üj(p, |
<о), g}{p, |
r, tu) — спектры |
u(p, |
f) и |
gj{p, |
r, t) |
||||||||
соответственно. |
|
|
|
|
|
|
|
|
|
|
|
|||
Подставляя |
в |
(2.1.2) вместо у*(г,ш) |
оператор |
у~(г, ш) |
||||||||||
в виде |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
' |
у (г, о) = |
- |
/А ] / Ж |
J S (р, «о) е' 7 г р г f (р) |
ф , |
|
||||||||
где /(р) —-множитель, |
учитывающий ослабление |
волн |
||||||||||||
в зависимости |
от угла |
прихода, и вводя для сокращения |
||||||||||||
записи четырехмерный вектор s=.(p, ш), получаем |
|
|||||||||||||
|
Т3 |
= |
J Jf/j (s„ s j a + (s,) a (s,) |
rfs.dSj, |
|
(2.1.7) |
||||||||
где^ (s,. s |
2 |
) |
= |
- |
n{zf2 |
|
Г |
(?,) |
f (PJJI^ (r-• ».• |
|||||
|
|
|
X e x p ( — ik1p1r1 |
-f- iK?2Tù |
s |
rfr,c?ra. |
|
|
(2.1.8) |
|||||
|
|
|
|
|
|
|||||||||
Чтобы |
функционалы |
были |
все одновременно |
изме |
||||||||||
римы, должны быть выполнены |
условия |
|
|
|
|
|||||||||
?f t ] = |
J I Й> (s,)a(s2 ) cfsjrfs, [ [Uj (s,, s) £ 4 |
(s, |
s2)— |
|||||||||||
|
|
|
|
|
- L / , l ( S 1 , S ) t / j ( S , S , ) ] c?s = |
0 |
|
|
|
|||||
|
|
|
|
|
|
|
при |
}фе. |
|
|
|
|
(2.1.9) |
|
Отсюда |
следует, |
что |
интегральные |
операторы |
с |
ядра |
||||||||
ми Uj, Uu должны |
коммутировать. |
|
|
|
|
|
||||||||
В том |
случае, |
когда область S — плоская |
апертура |
регистрирующего прибора и ее размеры велики по сравнению с длиной волны для всего рассматриваемого спектра частот, условие (2.1.9) можно упростить. При интегрировании по р произведения UjUk получим инте грал
î e ^ - ^ | / ( p ) | « r f p .
Если считать, что волны,' приходящие на апертуру под углом Ѳ к нормали ослабляются пропорционально cos Ѳ, то этот интеграл, уже встречавшийся при выводе выражения (1.4.14) для функции корреляции фонового излучения, равен
Функция К(гі—г2) |
имеет пространственный спектр, |
равномерный в области пространственных частот р<к и обращающийся в иѵль при p>k. Эта функция быстро
59