CC BY
26
Труды БГТУ, 2018, серия 3, № 2, с. 71-75
71
УДК 621.383.29:621.514.015
В. И. Никитенок1, С. С. Ветохин2, А. М. Бахарь1, Е. В. Терешко2
1 Белорусский государственный университет 2Белорусский государственный технологический университет
МЕТОДИКА РАСЧЕТА РАБОЧИХ ХАРАКТЕРИСТИК АДАПТИВНОГО ОБНАРУЖИТЕЛЯ СЛАБЫХ ОПТИЧЕСКИХ СИГНАЛОВ
Рассматривается обнаружение слабых оптических сигналов, когда на выходе фотоприемника наблюдается последовательность импульсов, распределенных по закону Пуассона, другими словами, стационарный или простейший пуассоновский поток, представляющий собой либо чистый шум, либо смесь полезного сигнала с шумом. С опорой на элементы теории статистических решений разработана методика получения рабочих характеристик адаптивного обнаружителя. В модели для стабилизации вероятности ложной тревоги используют подстройку порога обнаружения с учетом характеристик классифицированной обучающей выборки. В качестве последней принят «шумовой» пуассоновский поток.
Ключевые слова: слабый оптический сигнал, фотоприемник, адаптивный обнаружитель, вероятность ложной тревоги, вероятность правильного обнаружения.
V. I. Nikitsionak1, S. S. Vetokhin2, A. M. Bakhar1, E. V. Tereshko2
'Belarusian State University 2Belarusian State Technological University
THE METHOD OF CULCULATION OF THE WORKING CHARACTERISTICS OF A WEAK OPTICAL SIGNALS ADAPTIVE DETECTOR
This study examines the case of weak light signals detection when the Poisson series of the detector output pulses is occurred (stationary process), which represents a pure noise or mixture of noise and signal. The method of defining the working characteristics of an adaptive detector, which is elaborated under the theory of statistical decisions, is offered. The method uses the adaptation of the detection threshold to stabilize the value of alarm probability. It applies the characteristics of a background series that is a simple "noise" Poisson process.
Key words: weak optical signal, photo detector, adaptive detector, false alarm probability, true detection probability.
Введение. При обнаружении слабых оптических сигналов на выходе фотоприемника имеется простейший пуассоновский поток (ППП) с интенсивностью (прием чистого шума) или = + (прием смеси сигнала с шумом) [1, 2, 3]. Известно, что в этом случае оптимальное по критерию Неймана - Пирсона принятие решений на обнаружение состоит в сравнении числа принимаемых импульсов за фиксированное время с порогом обнаружения, зависящим от задаваемой вероятности ложной тревоги и интенсивности потока шума [1, 2, 3].
Если в качестве входных воздействий рассматривать последовательность временных интервалов между соседними импульсами ППП, которые распределены экспоненциально [1-4], то оптимальное обнаружение включает их суммирование при фиксированном числе импульсов и сравнение с порогом решения, зависящим как от задаваемой вероятности ложной тревоги, так и от величины [4]. Очевидно, что оба алгоритма обнаружения статистически эквивалентны. Оперирование с непрерывными
случайными величинами облегчает последующий анализ, не сужая общности решаемой проблемы, которая состоит в том, что рабочие характеристики оптимального обнаружителя слабых оптических сигналов получены при условии точно известного параметра экспоненциального распределения. Однако на практике такой случай является исключением из правила. В этих условиях, как показано в публикации [5], реальные характеристики обнаружения могут существенно отличаться от рассчитанных, так как параметр экспоненциального распределения при наличии чистого шума может уйти от нашего контроля, и порог обнаружения определяется неверно. Как отмечалось, в данном случае имеет место априорная неопределенность относительно параметра экспоненциального распределения.
Обратимся к двум известным принципам преодоления априорной неопределенности. Один и них состоит в разработке адаптивных алгоритмов, другой - в применении непараметрических, в том числе ранговых, алгоритмов. В настоящей
статье рассматривается пример применения адаптивного алгоритма обнаружения слабых оптических сигналов с применением методов обучения. В ней дополняются и уточняются результаты, представленные в работе [6]. В случае параметрической априорной неопределенности эти методы базируются на теории статистических решений с использованием классифицированной обучающей выборки [4]. Вид функции распределения известен (экспоненциальный) и в течение периода обучения, который должен предшествовать моменту принятия решения на обнаружение, надо на основании наблюдений сформировать оценку неизвестного параметра V Эта оценка используется затем для вычисления порога обнаружения вместо неизвестного истинного значения ^о интенсивности ППП при приеме чистого шума. Подчеркнем, должна быть твердая уверенность именно в этом: оценка действительно вырабатывается по чистому шуму. Имеет место так называемое обучение с учителем. В ходе обучения в данном случае происходит изменение параметра алгоритма обучения с целью приближения его к оптимальному алгоритму с точки зрения задаваемого критерия качества. Подобные устройства принято называть адаптивными (или самоорганизующимися) [4]. Обратим внимание, что применение методов обучения приводит к системам с характеристиками, близкими к оптимальным с известными параметрами.
Основная часть. Как отмечалось, в адаптивном обнаружителе в целях определения оптимального значения порога обнаружения для оценки интенсивности шума используется классифицированный обучающий ППП импульсов (обучающая выборка), соответствующий приему чистого шума.
Для удобства расчетов перейдем к величине А0, обратной А,0:
Ао = 1/ V (1)
Тогда экспоненциальное распределение последовательности временных интервалов между соседними импульсами ППП при приеме чистого шума имеет вид (плотность вероятности случайной величины) [4]:
/(Г) = 1/ Аое/А°, Г>0, Ао >0. (2)
Среднее значение такой случайной величины равно А0, а ее дисперсия А2. Эффективной несмещенной оценкой параметра А0 экспоненциального распределения (1) является выборочное среднее:
А0 = 1/ «01 , (3)
где т0 - объем классифицированной обучающей выборки, по которой проводится оценка параметра -90.
Как известно [4], дисперсия выборочного среднего всегда (при независимых элементах выборки размера т0) в т0 раз меньше дисперсии исходного распределения:
V{А0> = А2 /т0, (4)
а математическое ожидание
Е{А 0> = А0. (5)
При т0 >> 1 распределение оценки (3) приближается к гауссову с параметрами (4) и (5):
/ (А0) = 1/д/ 2^ {А0> X Xехр{-(А0 -Е{А0>)2/^{А0>)> = = у]0,5 т0 / п / А0 X Xехр{-(А0 - А0)2 / (2А2 /т0)>. (6)
Поскольку точность оценки параметра -90 зависит от объема обучающей последовательности т0, то и показатели качества адаптивного обнаружителя являются функциями этой величины. Поэтому вероятности правильного обнаружения и ложной тревоги адаптивного обнаружителя (Оа, Ка) вычисляют усреднением этих показателей для оптимального обнаружителя (О, К) по неизвестному параметру, получаемому путем оценивания.
Таким образом, уравнение рабочей характеристики адаптивного обнаружителя с учетом выражений (4)-(6) и (8) из [5] можно записать в виде
Ва = £О(т,К;А0 / АсШ^МА = = Г~ф(а0/ А^^т - (1+ А0/ Ас)Ф-1(1 - К))х Ху/0,5т0 / п / А0 х Xехр{-(А0 - А0)2 /(2А2 /т,)^ (7)
где Ф(х) = 1 / (2п)| е- /2dt - интеграл вероятности, Ф-1(х) - функция, обратная интегралу вероятности, Ас =1/ Хс.
Вычислить уравнение (7) аналитически в простой форме не удается. Поэтому воспользуемся известной аппроксимацией интеграла вероятности (ошибка аппроксимации составляет 5%) [7]:
Ф(х) = 1 - аехр{-Ь(й + х)2}, (8)
при х = 0-5,5 и 1 - Ф(х) = 0,5 10-8, где а = 0,65, Ь = 0,443, й = 0,75.
Остановимся более подробно на использовании приближения интеграла вероятности (8).
Из уравнения (8) следует, что его можно использовать для нахождения вероятности пра-
вильного обнаружения, значение которой превышает 0,5:
В = Ф( ) =
= 1 - аехр{-Ъ(а + хВ)2}, В > 0,5, (9)
где нижний индекс у х «привязывает» его к искомому показателю качества обнаружения. В данном случае это вероятность правильного обнаружения.
Очевидно, что если В < 0,5, то вероятность правильного обнаружения надо вычислять через вероятность пропуска 1 - В, которая больше 0,5, и при выполнении условий в формуле (8):
1 - В = Ф(хх_в ) =
= 1 - аехр{-Ъ(а + Х!-в)2}, (10)
откуда
В = аехр{-Ъ(а + х1-В)2}, при В < 0,5. (11)
Вероятность ложной тревоги / может находиться в пределах 10-2-10-8. При этом условия в формуле (8) не выполняются. Поэтому / найдем через вероятность правильного необнаружения 1 - /:
1 - / = Ф( х^) = = 1 - аехр{-Ъ(а + хь/)2}, (12)
откуда
/ = 1 -Ф( хг-р) = а ехр{-Ъ(а + хг-р )2}, (13)
х,-/ = ф-1 (Ф( хь/)) = ф-1 (1 - /) =
= у1 1/ Ъ 1п(а / /') - а. (14)
Вероятность правильного обнаружения В в уравнении (7) с учетом формул (9) и (14) равна
В(т,/;/дс) = Ф(£0 /Ъс4т -- (1 + £0/ ^) Ф-1(1 - /)) = = 1 - а ехр{-Ъ(2а1/Ъ 1п(а / /) + + (00/ -V1/Ъ 1п(а / /) + а))2}. (15)
Подставим выражение (15) в (7): Ва = {1 - а ехр{ -Ъ(2 а1/Ъ 1п( а / /) + + (V )(4т1/Ъ 1п(а / /) + а ))2}} х х^/0,5т0 / п / Ф0 х
хехр{-(Ф0 -^0)2 /(2^0 /т0)}^ = = 1 - а^/0,5т0 / п / Ф0 х
х ехр{-Ъ(2а 1/Ъ 1п(а / /) +
+ (О0/ )(Vm-у/1/b ln(a / F) + ¿))2} х Xexp{-(^o -^о)2/(2^0/mo)}d^. (16)
Выражение (16) можно привести к виду Ва = 1 - aj0,5m0 / п / tf 0 exp{-bL2 -Mtf 2} X
X £ exp{-^ + q^Ä (17)
где
l = 2d-JrnuOJF), m = 0,5m0 /
q = 2(M- bLK), ^2 = bK2 + M; K = (1/){4m-yj 1/ b ln( а / F) + d).
Для вычисления формулы (17) воспользуемся табличным интегралом [8]:
2 2
exp{-^ х ± qx}dx = ^л/п /^ exp{q2 / (4/)}. (18)
С учетом (18) выражение (17) приводится к виду
Ва = 1 - а^/0,5m0 / п / Ф0 х
х exp{-bL2 - M£0Ь/л /VbK2 + M х
X exp{(2(Md0 - bLK))2 /(4(bK2 + M))}. (19)
После преобразований формулы (19) имеем
Ва = 1 -а(1 + 2bg2(d + 4m -1/ b ln(a / F))2 / m0)-0'5 x X exp{-b(d(2 + g) + g4m -- (1 + g)^1/ b ln(а / F) )2 / (1 + 2bg2 X X (d + Vm 1/bln(a /F))2 /m0)}, (20)
при Ва > 0,5.
Проведем проверку правильности полученного выражения. Положим, m0 ^ да, и из уравнения (20) получим
lim m0 ^ Ва = 1 - exp{-b(d (2 + g) +
+ g4m - (1 + g )^1/ b ln( а / F ))2} = В. (21)
Замечаем, что формула (21) представляет собой рабочую характеристику оптимального обнаружителя с точно известными параметрами [5]. Это подтверждает правильность выполненных преобразований.
Обратим внимание на то, что в известных авторам работах условную вероятность ложной тревоги адаптивного обнаружителя не рассматривают. Однако ее анализ приводит к интересным для практики результатам. Как отмечалось,
вероятность ложной тревоги адаптивного обнаружителя можно получить, усредняя по неизвестному параметру выражение для вероятности ложной тревоги оптимального обнаружителя:
^(т,К; А0/ А0) = Ф(А0/ А{)(4т --Ф-1(1 - К)) -4т).
Далее для использования аппроксимации интеграла вероятности перейдем к вероятности правильного необнаружения:
Опр.н (т, К; А 0/ Ао) = 1 - Ка (т, К; А0/ А0) =
= 1 -Ф(4т -А 0/ А0(4т -Ф-1(1 - К))) = = Ф(А0 / А0 (л/т -Ф-1 (1 - К)) - 4т) = = Ф(л/т +А0/ А0 X X ^1/ Ь 1п(а / К) - й -4т)). (22)
Вероятность правильного необнаружения адаптивного обнаружителя получаем, усредняя выражение (22) по неизвестному параметру
Одр.н.а = Одр.н (т, К; А0 / А0)/(А0^А0 =
= Ф(^т +А0/ А0^1/Ь 1п(а / К) -й --4т Хц/0,5то / п / А0 X X ехр{-(А0 - А0)2 / (2А2 / mо)>dАо. (23)
Проводя преобразования уравнения (23), аналогичные получению вероятности правильного обнаружения, получаем
Опр.н а = 1 - а(1 + 2Ь^1/ Ь 1п(а / К) -- й -4т)2/ т0)-0,5 X
X ехр{- 1п(а / К )/(1 + + 2Ь^1/ Ь 1п(а / К) - й -4т)2 / т0)>, (24) откуда
К = 1 - Опр.н.а = а(1 + 2Ь^1/ Ь 1п(а / К) -
- й-4т)2 /т0)-0,5 Xехр{- 1п(а/К)/(1 + + 2Ь^1/Ь 1п(а /К) -й-4т)2/т0)>. (25)
В формуле (25) значения т должны выбираться не произвольно, а, как и в публикации [5], с учетом обеспечения заданных О и К при располагаемом £:
т = ((Ф-1 (О) + (1 + £)Ф-1 (1 - К))/ £)2. (26)
При т0 ^ да из выражения (27) получаем
11Шт0 ^ Ка = а ехр{- 1п(а / К )> = К. (27)
Логарифмируя выражение (27), имеем 1п(Ка / а) ~ ~ -1п(а / К) и Ка ~ К, что подтверждает правильность выполненных преобразований.
Заключение. Итак, в адаптивном обнаружителе для стабилизации условной вероятности ложной тревоги осуществляется подстройка порога обнаружения по классифицированной обучающей выборке, которой является «шумовой» ППП.
Представленная методика расчета рабочих характеристик адаптивного обнаружителя слабых оптических сигналов позволила получить сравнительно простые выражения для условных вероятностей правильного обнаружения и ложной тревоги. Они в дальнейшем могут быть применены различными пользователями для конкретного анализа работы адаптивных обнаружителей слабых оптических сигналов, сравнения их с оптимальными.
Литература
1. Шереметьев А. Г. Статистическая теория лазерной связи. М.: Связь, 1971. 264 с.
2. Макарова С. Н. Обнаружение слабых сигналов в оптическом диапазоне волн // Теория и техника радиолокации: Тр. МАИ. 1972. Вып. 207: Машиностроение. C. 198-208.
3. Сигналы и помехи в лазерной локации / под ред. В. Е. Зуева. М.: Радио и связь, 1985. 264 с.
4. Левин Б. Р. Теоретические основы статистической радиотехники. М.: Радио и связь, 1989. 752 с.
5. Никитенок В. И., Ветохин С. С. Об оптимальном обнаружении слабых оптических сигналов // Труды БГТУ. 2015. № 6: Физ.-мат. науки и информатика. С. 88-91.
6. Никитенок В. И. Быстрые непараметрические алгоритмы обнаружения сигналов. Минск: БГУ, 2010. 131 с.
7. Тамм Ю. А., Гомозова Т. М. К аппроксимации интеграла вероятности // Электросвязь. 1970. № 9. C. 77-78.
8. Градштейн И. С., Рыжик И. М. Таблицы интегралов, сумм, рядов и произведений. 7-e изд. СПб.: БХВ-Петербург, 2011. 1237 с.
References
1. Sheremet'ev A. G. Statisticheskaya teoriya lazernoy svyazi [Statistical theory of laser communication]. Moscow, Svyaz' Publ., 1971. 264 p.
2. Makarova S. N. Detection of weak signals in the optical band. Trudy MAI "Teoriya i praktika radiolokatsii" [MAI Proceedings "Theory and practice of radiolocation"], 1972, issue 207: Mechanical Engineering, pp. 198-208 (In Russian).
3. Signaly i pomekhi v lazernoy lokatsii [Signals and hindrance in laser location]. Moscow, Radio i svyaz' Publ., 1985. 264 p.
4. Levin B. R. Teoreticheskie osnovy statisticheskoy radiotekhniki [Theoretical basic of statistical radio technique]. Moscow, Radio i svyaz' Publ., 1989. 752 p.
5. Nikitionak V. I., Vetokhin S. S. On the optimal detection of weak light signals. Trudy BGTU [Proceedings of BSTU], 2015, no. 6: Physics and Mathematics. Informatics, pp. 88-91 (In Russian).
6. Nikitionak V. I. Bystrye neparametricheskie algoritmy obnaruzheniya signalov [Fast nonparametric algorithms of signal detection]. Minsk, BGU Publ., 2010. 131 p.
7. Tamm Yu. A., Gomozova T. M. On approximation of the probability integral. Elektrosvyaz' [Telecommunications], 1970, no. 9, pp. 77-78 (In Russian).
8. Gradshteyn I. S., Ryzhik I. M. Tablitsy integralov, summ, ryadov iproizvedeniy [The tables of integrals, sums, series, and multiplications]. Moscow, Nauka Publ., 2011. 1237 p.
Информация об авторах
Никитенок Виктор Иванович — кандидат технических наук, доцент кафедры физики и аэрокосмических технологий. Белорусский государственный университет (220006, г. Минск, пр-т Независимости, 4, Республика Беларусь). E-mail: [email protected]
Ветохин Сергей Сергеевич — кандидат физико-математических наук, доцент, заведующий кафедрой физико-химических методов сертификации продукции. Белорусский государственный технологический университет (220006, г. Минск, ул. Свердлова, 13а, Республика Беларусь). E-mail: [email protected]
Бахарь Андрей Михайлович - первый заместитель начальника военного факультета. Белорусский государственный университет (220006, г. Минск, пр-т Независимости, 4, Республика Беларусь). E-mail: [email protected]
Терешко Елена Владимировна - ассистент кафедры высшей математики. Белорусский государственный технологический университет (220006, г. Минск, ул. Свердлова, 13а, Республика Беларусь). E-mail: [email protected]
Information about the authors
Nikitsionak Viktar Ivanavich — PhD (Engineering), Associate Professor, the Department of Physics and Aerospace Technologies. Belarusian State University (4, Nezavisimosti Ave., 220006, Minsk, Republic of Belarus). E-mail: [email protected]
Vetokhin Siarhei Siarheevich — PhD (Physics and Mathematics), Associate Professor, Head of the Department for Physical and Chemical Methods of Testing. Belarusian State Technological University (13a, Sverdlova str., 220006, Minsk, Republic of Belarus). E-mail: [email protected]
Bakhar Andrey Mikhaylavich - Deputy Head of the Military Department. Belarusian State University (4, Nezavisimosti Ave., 220006, Minsk, Republic of Belarus). E-mail: [email protected]
Tereshko Elena Vladimirovna - assistant lecturer, the Department of Higher Mathematics. Belarusian State Technological University (13a, Sverdlova str., 220006, Minsk, Republic of Belarus). E-mail: [email protected]
Поступила 10.05.2018
