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

книги / Теория и расчет импульсных и цифровых оптико-электронных систем

..pdf
Скачиваний:
4
Добавлен:
12.11.2023
Размер:
18.87 Mб
Скачать

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

Флуктуационные шумы приемного тракта включают в себя шумы фотоприемника [дробовые шумы фототока, тепловые шумы,

дополнительные шумы

фотоэлектронных умножителей (ФЭУ)

и лавинных фотодиодов

(ЛФД), связанные со статистическим

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

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

Действие естественных источников излучения (Солнца, пла­ нет, звезд, облаков, атмосферы, поверхности земли), создающих мешающий фон, может проявляться как в увеличении шумов фотоприемника, так и в появлении ложных импульсных сигналов. Последние наблюдаются в случаях достаточно резких изменений интенсивности фоновых помех вследствие различного рода пере­ мещений поля приема прибора относительно источников фоно­ вого излучения (сканирования, поступательного движения при­ бора, угловых вибрационных колебаний и т. д.). Можно отметить, что при достаточно малых длительностях рабочих сигналов (напри­ мер, порядка 10—100 нс) влияние рассматриваемых фоновых помех во многих случаях может быть сведено к минимуму мето­ дами линейной фильтрации (выбором передаточной функции тракта с соответствующей полосой непрозрачности по низким частотам). При этом одновременно будет достигаться уменьшение (либо прак­ тическое устранение) влияния составляющих шума с убывающим по частоте энергетическим спектром (избыточные шумы) и шума вибраций (микрофонного эффекта).

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

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

41

препятствовать обнаружению рабочего сигнала,

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

на относительно малых дистанциях до объекта локации.

При последующем рассмотрении диффузных

помех (а также

влияния турбулентности среды) будем ориентироваться на наибо­ лее распространенный класс оптико-электронных систем, рабо­ тающих в атмосфере. Задачи оптической локации объектов и передачи сигнала в водной среде имеют свои особенности, свя­ занные со спецификой характеристик и параметров среды распро­ странения излучения.

Вопросам рассеяния оптического излучения в атмосфере по­ священо большое количество публикаций, основные из которых приводятся, например, в монографиях [7, 12, 131. Для инженер­ ных задач важной является возможность оперативной прибли­ женной оценки характеристик диффузных помех применительно к типовым схемам локационного наблюдения. В этом плане можно отметить работы 115, 39], в которых для случая нестационарного облучения среды в приближении однократного рассеяния полу­ чены соотношения, связывающие временную функцию диффуз­ ной помехи с оптическими характеристиками атмосферы и пара­ метрами приемно-излучающего устройства.

Проведенные оценки точности приближения однократного рассеяния показали, что при типовых параметрах оптических полей ОЭС влияние рассеяния высших порядков становится суще­ ственным лишь при относительно больших плотностях среды, например при оптических толщинах атмосферы больше 15—20.

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

Формирование диффузной помехи в результате отражения (обратного рассеяния) зондирующего импульса от атмосферы можно рассматривать как операцию линейного преобразования оптического сигнала (подобно случаю отражения сигнала от про­ тяженных объектов). В соответствии с этим атмосферу* как источ­ ник диффузной помехи удобно описывать импульсной харак­ теристикой обратного рассеяния ga (£), понимая под ней облу­ ченность, создаваемую в точке приема (на приемном объективе) при зондирующем импульсе в виде б-функции: Фиб (t) = W„ 6 (/), причем W9 = 1 Дж.

Характеристика gn (t) может быть представлена в обобщенном

виде

 

 

 

W “ Л (О Т Ю/Л

 

(2.71)

где г\ (/) — функция, зависящая

только от

параметров

среды;

Y (2) — функция, зависящая от

геометрии

и параметров

опти­

ческих полей излучения и приема.

При этом Y (t) = 0, если t <3 t0, где t0 — время, соответствую­ щее границе так называемой «теневой зоны», из которой рассеян­

42

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

Например, для однородной. атмосферы в случае коллимиро­ ванных несовмещенных равномерных (конусных) полей излучения и приема с параллельными оптическими осями при пренебреже­ нии размерами апертур приемно-излучающего устройства («то­ чечные» излучатель и приемник) импульсная характеристика обратного рассеяния, согласно работе [15], имеет вид

g* (t) =

[ f

-

arcsin -T ~

' Г V

1 -

(-f

У

 

 

 

 

 

 

 

 

(2.72)

где

а — показатель рассеяния

атмосферы;

<ря — среднее

значе­

ние

индикатрисы

рассеяния

в

направлении

назад;

 

— показа­

тель ослабления

атмосферы.

 

 

 

 

 

 

Если зондирующий импульс характеризуется мощностью излу­

чения Ф„ (0, то

диффузная

помеха (облученность

на

приемном

объективе) может быть представлена в виде интегральной свертки

£д(0==t J (t —т) gR(т) da.

и

Отметим, что теперь исходное пренебрежение размерами апер­ тур приемоизлучателя не препятствует нахождению потока диф­ фузной помехи на фотоприемнике обычным способом, т. е. путем умножения £ д (^) на эффективную площадь приемного объектива.

При непараллельных (пересекающихся) оптических осях полей излучения и приема рассеивающий объем среды является ограни­ ченным и импульсная характеристика обратного рассеяния со­ стоит из нескольких участков, описываемых различными функ­ циями [39], причем значения временных интервалов, соответ­ ствующих этим функциям, зависят от соотношений между углами излучения, приема и схождения оптических осей.

Типичная форма импульсной характеристики ga (t) имеет вид одновершинной симметричной кривой, положение tm максимума которой в общем случае зависит от геометрических параметров схемы локационного наблюдения и показателя ослабления атмо­ сферы 8j. При этом в условиях относительно высокой прозрач­ ности среды временная структура рассеянного назад излучения определяется в основном геометрией наблюдения. В частности,

при параллельных оптических осях и L0et = 0,01, где L0 =

ctj2 —

дистанция, соответствующая

границе

«теневой

зоны»,

можно

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

положить tm~ 1,5/0 [15].

По мере повышения показателя ослабления максимум функ­

ции £д (t) сдвигается в сторону начала

импульса,

т. е. крутизна

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

43

влияние на формирование диффузной помехи оказывают пара­ метры среды.

В случае идеализированной совмещенной схемы локационного наблюдения (совпадающие оптические оси излучения и приема) выражение для импульсной характеристики ga (t) (в приближении однократного рассеяния) может быть представлено в виде

гд№ = - ^ е х р ( - е1сО,

(2.73)

где Кв — коэффициент, зависящий от характеристик полей излу­ чения и приема.

В частности, при равномерных согласованных конусных полях излучения и приема Кн — 2, что следует из сопоставления выра­ жений (2.72) и (2.73) при t0 = 0. Это же значение Кн сохраняется и для других форм диаграмм направленности, например при гаус­ совом поле излучения, если в пределах эффективного угла излу­ чения, в котором сосредоточена основная часть энергии излуче­ ния зондирующего сигнала, поле приема может считаться рав­ номерным.

Как видно из формулы (2.73), кривая gд (/) для совмещенной схемы наблюдения является монотонно убывающей (экспонен­ циально-квадратической) функцией. Характеристика диффузной помехи соответственно будет иметь максимум в момент времени /м, находящийся в пределах длительности зондирующего импульса. Заметим, что формула (2.72) переходит в выражение (2.73) при t > /01т. е‘ при достаточно больших значениях t характеристика диффузной помехи не зависит от расстояния между параллель­ ными , оптическими осями полей излучения и приема (базы прибора).

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

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

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

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

Применительно к работе ОЭС в условиях сильно замутненной атмосферы, когда возникает необходимость учета рассеяния выс­ ших порядков, для определения характеристик диффузных по­ мех обычно прибегают к моделированию на ЭВМ процесса рассея­ ния излучения (движения фотона в рассеивающей среде) на основе метода Монте-Карло [28].

Экспериментальные исследования не обнаружили скольконибудь существенного различия энергетики рассеяния лазерного (когерентного) и некогерентного излучения.

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

Ослабление сигнала в атмосфере характеризуется (в преде­ лах применимости закона Бугера) величиной прозрачности среды

где 2L — общая длина пути, проходимого излучением в атмо­ сфере.

При работе в однородной атмосфере тв = ехр (—-2^1). Пока­ затель ослабления гх является суммой показателей рассеяния о и поглощения Кп> причем в большинстве практических случаев основной вклад в ослабление сигнала вносится рассеянием (при работе в видимой и ближней ИК-областях спектра).

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

Когда длина трассы намного превышает размеры наибольших неоднородностей атмосферы, плотность вероятности амплитуд

45

лучистого сигнала может быть охарактеризована логарифмиче­ ски — нормальным законом

Wc (W =

ехр

(In / э — (In / э ) ) 2

1

(2.74)

К

J

 

 

 

где ом = <(1п / 3 — <1п /э» 2> — дисперсия флуктуаций

логарифма

силы излучения; ( ) — знак математического ожидания.

Флуктуации амплитуды сигнала, обусловленные

турбулент­

ностью атмосферы, имеют низкочастотный энергетический спектр, максимум которого при горизонтальном наблюдении обычно не превышает 5—15 Гц. Эти флуктуации являются помехой мульти­ пликативного типа и при взаимодействии с другими видами ампли­ тудной нестабильности сигнала приводят к изменению итогового закона распределения амплитуд и определенному ухудшению условий обнаружения сигнала. В этом заключается основное проявление влияния турбулентности среды при энергетическом приеме.

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

фронта оптического

сигнала (адаптивные оптические си­

стемы).

 

Помехи, обусловленные отражением рабочего излучения

от различного рода

ландшафтных образований, могут суще­

ственно затруднить выделение полезных сигналов от наземных объектов.

Статистические свойства подобных помех в настоящее время еще мало изучены. Несколько исследован вопрос о характеристи­ ках помех, связанных с отражением рабочего излучения от взвол­ нованной водной поверхности [5], которая может считаться ста­ тистически однородной. При этом основную опасность здесь обычно представляет зеркальная составляющая отраженного излучения, влияние которой проявляется лишь при достаточно больших углах визирования относительно горизонта.

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

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

Основные методы повышения помехозащищенности оптико­ электронных приборов рассмотрены в монографии [621.

46

Глава 3

МЕТОДЫ ОПТИМАЛЬНОГО ОБНАРУЖЕНИЯ ИМПУЛЬСНЫХ ОПТИЧЕСКИХ СИГНАЛОВ

§ 3.1. ОСНОВНЫЕ ТЕОРЕТИЧЕСКИЕ ПОЛОЖЕНИЯ

При приеме оптических импульсов прежде всего решается задача обнаружения, сущность которой в классической поста­ новке сводится к принятию на основе определенного правила решения о наличии или отсутствии рабочего сигнала в реализации случайного процесса, воздействующего на вход приемного устрой­ ства ОЭС. Полное описание случайного процесса £ (/) достигается заданием совокупности многомерных плотностей вероятностей.

Выберем N произвольных значений аргумента случайной функции £ (/) на интервале наблюдений Т и рассмотрим значения | (t) в эти моменты времени. Любое конкретное значение Jtft слу­

чайной

величины £ (th) называется выборочным,

а совокупность

значений xh — выборкой. Число

N выборочных

значений назы­

вают объемом

выборки.

вероятностей:

 

Зададим следующую систему

 

Wi

(хъ

A) dx — вероятность

того, что хх < £ (А) < хх 4- dx^

W2 (*ь

х2, А» A) dx1dx2 — вероятность того, что хх < £ (А) <3

<4 хх 4- dx-ь х2 <3 £ (А) <3 х2 + dx2\

 

WN

{хъ

х2,

...» xN, А» А» •••»

Av) dxtdx2 ... dxN — вероятность

того, что

 

 

 

 

 

<

I (А)

<

*i 4* dxi;

 

 

х2 <

£ (А)

<

х2 4- dx2\

 

 

Хи

%(Av)

Хи 4* dxjy.

 

 

Функция WN(xi, x2, ..., xN, А» At •••» A4 называется N-мерной плотностью вероятности случайного процесса £ (/).

Плотность вероятности может зависеть от состояния S рас­

сматриваемого случайного процесса.

В этом случае переходят

к условным плотностям вероятности случайного процесса

А»

(А'Ь А 1*5); W 2 (х ъ х2, А» А |*5);

..., W ц {Хъ Х2, ...,

А» .... Av|*5) либо к условным интегральным функциям распре­

деления

Рг (хъ А 1*5); Рг \хъ x2t А. А 1*5); ...; Ря(хъ х2) ..., xNf

At А>

Av |*5).

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

Например, совместная плотность вероятности выборочных

значений с

независимыми элементами хъ х2у .... xN равна

 

N

Г „

(х |S) = WN (xlt хг..........IS) = П Wt(хк15).

 

ft—1

47

Для выборки из дискретного распределения

N

PN(x=sxь . . ., Х = ЛГдг|5)= n ^ fris .

6=1

Функции правдоподобия являются удобными средствами иссле­ дования, особенно при изучении дискретных систем. Однако в ряде задач полезнее непрерывное описание случайных процессов. Для этого переходят от функций правдоподобия к функционалам

правдоподобия [1 ].

 

при

Выберем значения аргумента случайной функции £ (/)

*i <*2 < ••• < * *

таким образом,

чтобы tm tm_i — А,

т =

= 2, 3, ...» N; А =

Tl(N — 1), где

Т — длительность выборки.

Рассмотрим предел, к которому стремится функция правдо­ подобия в случае, когда А -*• О, N -> оо. Характеристика случай­ ного процесса, получаемая путем такого предельного перехода,

называется функционалом правдоподобия, т. е.

 

3T[x(t)\S]=\m WA*u *2» •••»

^1» ^2» •••» IN I

(3*1)

A-*-0

 

 

N-+oo

Следует отметить, что функционалы правдоподобия полностью учитывают все статистические свойства случайного процесса.

При синтезе ОЭС обнаружения и выделения сигналов исполь­ зуются отношения правдоподобия

Л = WN(х| Sk)/WN|Sn) либо A = ST[x(t)\Sh\!T [х (t) |Sn).

При этом правило принятия решения записывается в следу­ ющем виде:

нк

Л С,

(3.2)

Нп

 

где ffh и Нп — соответствующие гипотезы о наличии состояния Sk либо Sn соответственно; С — величина, характеризующая поро­ говый уровень в зависимости от выбранного критерия качества.

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

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

втабл. 3.1.

Втабл. 3.1 приведены следующие обозначения: р и q — апри­ орные вероятности наличия и отсутствия сигнала соответственно;

Яоо и Ли — платы, соответствующие поощрениям за принятие правильных решений; П01 и Л10 — платы, соответствующие штра­ фам за принятие ошибочных решений; Рл. т — условная вероят­

ность ошибки первого

рода (вероятность ложной

тревоги);

Рю (С) — интегральная

функция распределения

отношения

48

 

 

 

 

 

Таблица 3.1

к»

Критерий кячества

 

 

с

п/п

 

 

1

Критерий минимума среднего риска

Я

Лад — Лоо

2

Максимум

апостериорной вероят­

Р

Лго — /711

 

 

q/p

 

ности

 

 

 

3

Максимум

правдоподобия

 

 

1

4

Минимаксный

q*

Л01 — Лоо

Р*

л

1 0 - Л „

5

Неймана-Пирсона

РМ (Л) =

1 — Рл. т

правдоподобия при гипотезе Н0 (отсутствие сигнала); р* и q* — значения априорных вероятностей наличия и отсутствия сигнала, которым соответствует абсолютный максимум байесовского риска.

Отношениеq*/p* определяется, согласно работе [26], решением трансцендентного уравнения вида

I п\\П + с [1 - Р10 (q*/p*)C] = Рц (Cq*!p%

(3.3)

где Ри (Cq*lp*) — интегральная функция распределения отноше­ ния правдоподобия при гипотезе Нх (наличие сигнала).

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

Величина отношения сигнал/шум (по току или напряжению) на выходе приемного тракта линейной системы, характеризуемого

комплексным

коэффициентом передачи

К (/со),

определяется

зависимостью

 

 

 

 

 

 

00

 

 

 

 

 

- j - J

S (/а) к

(/и)

da

 

 

И = г---- Ч

---------------------- тТ/Г»

<3-4)

 

 

0 (“ Ж

(/“ ) I2

J

 

 

 

 

 

 

где 5 (/©) =

F [S (t) ] A J

s (t) er№ dt — спектральная функция

сигнала s (t);

— по

 

 

 

 

G (©) — энергетический спектр входных шумов;

t0 — значение времени, когда сигнал на выходе приемного тракта достигает максимального значения.

49

При использовании оптимального (с позиций обнаружения сигнала) линейного фильтра достигается максимально возможное отношение сигнал/шум

 

 

 

<3-5>

Если на выходе действует нормальный белый шум с

G (со) = <3,

то выражение (3.5) примет вид

 

Г

00

- 1 1 / 2

 

И = | т | - |

|S(/(o)|=dcoJ .

(3.6)

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

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

§ 3.2. ОБНАРУЖЕНИЕ ДЕТЕРМИНИРОВАННОГО СИГНАЛА НА ФОНЕ БЕЛОГО ГАУССОВА ШУМА В УСЛОВИЯХ АДДИТИВНОСТИ СИГНАЛА И ШУМА

Для определения оптимального алгоритма обнаружения де­ терминированного сигнала при постоянной спектральной плот­ ности шума G воспользуемся отношением правдоподобия

 

Ю] ^ 7°

(3.7)

где

[х (t) |s] = £д exp j— g- J[* (t) s (i)fdt , — функционал

правдоподобияпри наличии полезного сигнала s (/), а &~[x(t)\Q] = = &д ехр — ^ -J x2(/)d/J — функционал правдоподобия при от­

сутствии сигнала; kA — некоторый множитель, не зависящий от выборки хъ х2, ..., xN, но зависящий от А.

Для последующих расчетов удобно использовать логарифм отношения правдоподобия

in Л = J *2 (/) dt

 

Нх

(3.8)

 

^ In С.

т

т

но

 

 

 

После несложных преобразований приходим к записи правила принятия решения в следующем виде:

j

s

(/)

Л-

(0dt

 

4 - Сin

J

S2

(0dt.

(3.9)

 

 

 

и

+ - Y

 

 

50