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

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

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

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

Для Р-представимых операторов плотности задача статистического описания поля сводится к отысканию классического функционала распределения вероятностей для амплитуд монохроматических плоских волн а(р, ш). С этим функционалом взаимооднозначно связан харак­

теристический

функционал *>, определяемый выражением

Ф 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, 2) =

= J ü i(P . cùi—cù2)g-*j(p,

14, —ші)еГз(р, гг , 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

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