CC BY
73
ОБНАРУЖЕНИЕ ОПТИЧЕСКИХ СИГНАЛОВ ПРИ ПРИЁМЕ ПОТОКА ФОТОЭЛЕКТРОНОВ С НЕИЗВЕСТНОЙ ФОРМОЙ ПЛОТНОСТИ
Парфёнов В.И., Кириллов В.С.
Воронежский государственный университет
Аннотация
Предложен алгоритм оптимального обнаружения оптического сигнала с неизвестной формой плотности потока по статистике фотоэлектронов. Рассчитаны характеристики приёма методом статистического моделирования.
Ключевые слова: дробовой шум, пуассоновская аппроксимация, априори неизвестная интенсивность, оценка максимального правдоподобия, критерий идеального наблюдателя, метод статистического моделирования.
Задача приёма слабых оптических сигналов на фоне шумов представляет значительный интерес при разработке и исследовании систем связи, оптической локации (в координаторах цели лазерных головок самонаведения, авиационных, артиллерийских и танковых лазерных прицелах, биноклях-дальномерах) [1, 2 и др.]. Оптический сигнал в виде светового поля распространяется в канале связи, далее на приёмной стороне оптическое поле собирается оптикой приёмника и детектируется фотодетектором. Основным способом приёма оптического сигнала является так называемое прямое детектирование (энергетический приём), при котором оптический сигнал направляется на фотодетектор и на его выходе образуется электрический сигнал (фототок). В отличие от приёма сигналов в радиосвязи и радиолокации, когда шумы и помехи практически всегда можно считать независимыми, детектирование оптических сигналов сопровождается неустранимым дробовым шумом, статистические характеристики которого тесно связаны со значением самого сигнала. Причём дисперсия дробового шума в каждой точке оси времени пропорциональна интенсивности сигнала в той же точке. При этом выходной сигнал фотоприёмника можно представить в виде аддитивной комбинации двух составляющих, одна из которых пропорциональна оптическому сигналу, а другая - случайная - является дробовым шумом, рождённым самим этим сигналом. Однако эту информацию можно учесть и использовать при синтезе обнаружителей оптических сигналов на фоне комплекса шумов (дробовых, тепловых и др.) только в случае, если имеются априорные сведения о параметрах и форме сигналов. В противном случае приходится полагать, что форма оптического сигнала (точнее, его плотность потока как функция времени) априори неизвестна.
Как показано в [1], форма оптических импульсов может оказывать достаточно существенное влияние как на эффективность их обнаружения, так и на точность оценки неизвестных параметров, например, временного положения. Однако в большинстве известных работ (см., например, [3, 4]) подобная задача решалась в предположении о том, что форма плотности потока сигнала и шума является известной. Итак, предположим, что наблюдению на интер-
вале времени [0, Т] доступна некоторая реализация, представляющая собой либо аддитивную смесь оптического сигнала с неизвестной формой интенсивности и фоновым излучением (гипотеза Н1), либо только фоновое излучение (гипотеза Н0). Далее в соответствии с теорией фотоэффекта [5] примем, что поток первичных фотоэлектронов и, соответственно, поток одноэлектронных импульсов на выходе фотоумножителя или лавинного диода подчиняется распределению Пуассона. Что касается статистики первичных электронов, то при приёме когерентного излучения эта модель адекватна реальности. Если же принимается тепловое излучение, то при условии, что длительность времени измерения превышает время когерентности, соответствующее этому случаю биноминальное распределение фотоотсчётов с отрицательным показателем также переходит в пуассоновское [6]. Однако следует отметить, что реальные фотодетекторы из-за сложности явлений размножения фотоэлектронов искажают исходные распределения фотоэлектронов [7]. Следовательно, подобную пуассоновскую модель следует рассматривать как первое приближение. Что касается фонового и теплового излучений, то, учитывая, что изменение во времени их интенсивности происходит существенно медленнее, чем изменение формы используемых сигналов, в дальнейшем будем считать, что значения плотности потока (интенсивности) фона и соответствующего ему дробового шума могут быть заранее оценены. С учётом этого в дальнейшем полагаем, что плотность потока фонового излучения известна и постоянна. Что касается темнового тока фотодетектора, то он также может быть учтён путём соответствующей коррекции значений фонового излучения [8].
Перейдём теперь к нахождению алгоритма обнаружения оптического сигнала с неизвестной формой интенсивности 11 (/) при наличии фонового излучения с известной постоянной интенсивностью 10. В соответствии с принятой нами моделью поток фотоэлектронов подчиняется закону Пуассона. Это значит, что вероятность насчитать N фотоэлектронов на интервале времени Дt вблизи точки tt определяется как
Р(N,^)А] ехр[-1(Ґ,)Аґ]:
(1)
N!
где ) = 1* ^) +10 при выполнении гипотезы Н1
либо ) = 10 при выполнении гипотезы Н0.
Совместная вероятность выпадения N фотоэлектронов на интервалах длительности At в точках t1, t2, ... tM (функция правдоподобия) в силу независимости пуассоновских отсчётов будет определяться произведением вида
Р(N ...,мм,^..хм) = Пр(N.,ґ,).
(2)
При выполнении гипотезы Н1 из (1), (2) получим следующее выражение для логарифма функции правдоподобия:
м
ІПр = ІП[1, (ґ, ) + 1о ] -
І=1
-Іп(Аґ)-[1, (Ґ,) + 1о]Аґ-ІП(N!)}.
Аналогично, при выполнении гипотезы Н0 логарифм функции правдоподобия:
м
Іп Ро = £{^ Іп (1о)-Іп (Аґ )-1о Аґ - Іп (N !)}.
,=1
Из этих выражений находим логарифм отношения правдоподобия
Л = Іп Р1 - Іп Ро =
1, (Ґі )
1 + -
1о
-1, (ґІ )Аґ \.
(3)
Заметим, что это выражение получено в предположении, что интенсивность оптического сигнала 1* (t) известна.
Для нахождения решающей статистики в случае априори неизвестной формы интенсивности сигнала поступим следующим образом. В выражении (3) 1 (t,■), О = 1,2,... М) - неизвестные величины.
Определим оценки этих величин, исходя из следующей системы уравнений максимального правдоподобия:
ЭЛ
Э1, (ґ,)
= о.
Из этой системы уравнений находим
£ N
1 ^)=А-10. At
Учитывая неотрицательность интенсивности, положим
К (Ґі ) =
N7 Аґ -1о, >1о Аґ,
N. >1о Аґ.
Подставляя эту оценку в (3), находим следующее выражение для логарифма отношения правдоподобия:
Л=І і *
х м (N. -1о Аґ)
( N. ^
Іп . 1 1
1 КоАґ; ]
+ 1о Аґ >х
(4)
где и(х) = 0 при х > 0 и и(х) = 1 при х < 0 .
Таким образом, алгоритм обнаружения оптического сигнала с неизвестной формой интенсивности при наличии фонового излучения должен заключаться в формировании величины Л (4) и сравнении её с некоторым порогом к, определяемым из заданного критерия эффективности. В дальнейшем в качестве критерия эффективности будем рассматривать критерий идеального наблюдателя [9], в соответствии с которым необходимо минимизировать полную (среднюю) вероятность ошибки Ре. Средняя вероятность ошибки Ре определялась как (а + Р) / 2, где а - вероятность превышения порога обнаружения решающей статистикой (4) при выполнении гипотезы Н0 (отсутствие полезного сигнала на входе приёмника). Вероятность Р - это вероятность непревышения порога обнаружения решающей статистикой (4) при выполнении гипотезы Н1 (наличие полезного сигнала на входе приёмника). Определение вероятности ошибки Ре осуществлялось с помощью статистического моделирования на ЭВМ по 50000 реализаций сигнала с фоном (гипотеза Н1) и только фона (гипотеза Н0). При моделировании в случае выполнения гипотезы Н0 осуществлялось формирование пуассонов-ского потока с плотностью 10 , рассчитывалось количество фотоэлектронов на интервале наблюдения, формировалась достаточная статистика в соответствии с формулой (4) и осуществлялось сравнение с набором пороговых значений. При превышении порога фиксировалась ошибка. Вероятность а(к) определялась как отношение количества превышений порога к в N испытаниях к количеству этих испытаний. Аналогичным образом вычисляется вероятность Р(к), за исключением того, что при формировании пуассоновского потока его плотность предполагалась равной 10 +1* ^). После получения вероятности Ре (к) = (а(к) + Р(к)) / 2 для набора пороговых значений к , определялся тот порог, при котором эта вероятность минимальна. Эта вероятность и представляет собой искомую среднюю вероятность ошибки.
Рассматривались две задачи: обнаружение сигнала с известным временем прихода (для определённости полагаемым равным нулю) и со случайным (равномерно распределённым на интервале наблюдения) временем прихода. При моделировании полагалось, что форма интенсивности сигнала имеет экспоненциальный вид, а именно:
1 * ^^-% ]и О1 - ^),
І=1
,=1
о
где t0 - момент прихода сигнала (равный либо 0, либо принимающий случайные значения от реализации к реализации на интервале наблюдения [0; Т-Т0]), где Т0 - длительность сигнала по уровню 0,9 от его полной энергии.
Были исследованы зависимости вероятности средней ошибки Ре от параметров: V = уТ,
gn = 10Т и отношения сигнал/шум:
Т1 {t)и
Ч = 1^----^ ------,
0 \Т Vgn
где и = 1Т.
Здесь параметр V характеризует скорость изменения плотности потока во времени; чем больше параметр V , тем быстрее изменяется плотность потока. Параметр gn - это интегральная плотность потока фонового излучения.
На рис. 1 приведены зависимости Ре (М) при априори известном времени прихода = 0, а на рис. 2
- при случайном времени прихода оптического сигнала. Кривые 1 на этих рисунках построены при V = 10, gn = 20, ч = 1, кривые 2 - при V = 10, gn = 20, ч = 0,5, кривые 3 - при V = 5, gn = 20, Ч = 1, кривые 4 - при V = 5, gn = 20, ч = 0,5 и, наконец, кривые 5 - при V = 10, gn = 10 , ч = 1.
Рис. 1. Зависимость вероятности полной ошибки от числа интервалов разбиения М (10 = 0 )
Из анализа рис. 1 и 2 можно сделать следующие выводы. Естественно, средняя вероятность обнаружения при случайном времени прихода больше по сравнению со случаем априори известного и постоянного времени прихода. Кроме того, в обоих случаях (в первом - более явно, во втором - менее) в зависимостях Ре (М) наблюдаются минимумы. Следовательно, в алгоритме (4) величина М не должна выбираться слишком большой, как следует из стандартных умозаключений. Малые значения параметра М — 1 + 3 приводят к завышенным в несколько раз вероятностям ошибки, что достаточно очевидно, так как при этом интервалы разбиения являются
слишком большими и в результате на оценку интенсивности сигнала будет оказывать существенное влияние фоновое излучение. Если же М велико, то увеличивается количество оцениваемых параметров, что также приводит к росту вероятности ошибки. Видно, что для рассматриваемых случаев оптимальное значение М, обеспечивающее минимум средней вероятности ошибки, должно быть порядка 4^7.
Рис. 2. Зависимость вероятности полной ошибки от числа интервалов разбиения М (10 - случайно)
На рис. 3 и 4 изображены зависимости Ре (ч) также при t0 = 0 и при случайном времени прихода, соответственно.
На рис. 5 и 6 приведены зависимости Ре (gn) при аналогичных условиях наблюдения.
Рис. 3. Зависимость вероятности полной ошибки от отношения сигнал/шум ч (10 = 0)
Рис. 4. Зависимость вероятности полной ошибки от отношения сигнал/шум ч (10 - случайно)
Рис. 5. Зависимость вероятности полной ошибки от интегральной интенсивности фона ^ (10 = 0)
Рис. 6. Зависимость вероятности полной ошибки от интегральной интенсивности фона ^ (10- случайно)
Рис. 3 и 4 построены при gn = 2, а рис. 5 и 6 -при ч = 1. На всех этих рисунках кривые 1 построены при М = 1, кривые 2 - при М = 4 и кривые 3 -при М = 30 . Как и ранее, подтверждается наличие оптимального значения параметра М ~ 4 , обеспечивающего минимум средней вероятности ошибки. Кроме того, следует отметить следующий факт. При случайном времени прихода наименьшую среднюю вероятность ошибки обеспечивает выбор М = 1, если мало отношение сигнал/шум или мала интенсивность шума 10. Однако подобные ситуации практического интереса не представляют, так как при этом средняя вероятность ошибки принимает слишком большие значения (больше, чем 0,1).
Отметим также, что, как и всегда [10, 11] в оптических задачах, средняя вероятность ошибки зависит не только от отношения сигнал/шум, но и от интенсивности фонового излучения (при постоянном отношении сигнал/шум).
Таким образом, с помощью компьютерного моделирования получены зависимости средней вероятности ошибки от различных параметров оптических сигнала и шума в случае неизвестной формы плотности потока сигнала. Следует отметить, что в известных работах по данной тематике синтез оптимальных структур обнаружителей при неизвестной форме плотности потока сигнала не был выполнен.
Влияние формы плотности потока сигнала в основном исследовалось либо на качественном уровне [8], либо определялся проигрыш по каким-либо характеристикам (например, ОСШ) при отклонении формы потока от предполагаемой [1]. В данной работе выполнены синтез и анализ оптимального алгоритма обнаружения оптического сигнала при отсутствии информации о форме плотности потока этого сигнала. Анализ приведённых результатов моделирования свидетельствует о следующем. Повышение эффективности обнаружения (наличие минимума в зависимости Pe (M) при определённом значении параметра M) не сильно зависит от отношения сигнал/шум. Подобная зависимость более явно проявляется при больших интенсивностях шума 10 и при достаточно быстром изменении интенсивности сигнала во времени, т.е. при достаточно больших величинах параметра v = gT.
Литература
1. Лебедько, Е.Г. Системы оптической локации / Е.Г. Ле-бедько. - СПб: НИУ ИТМО, 2012. - 129 с.
2. Козинцев, М.П. Основы импульсной лазерной локации / М.П. Козинцев, В.М. Белов [и др.] - М.: Изд. МГТУ им. Н.Э. Баумана, 2006. - 512 с.
3. Березин, В. В. Обнаружение и оценивание координат изображений точечных объектов в задачах астронавигации и адаптивной оптики / В.В. Березин, А.К. Цы-цулин // Вестник ТОГУ. -2008. - № 1(8). - C. 11-20.
4. Куцов, Р. В. Обнаружение и оценка вектора скорости движения неоднородного пространственно-протяжённого объекта по изображению / Р.В. Куцов, А.П. Трифонов // Автоматика и телемеханика. - 2009. - № 8.-C. 96-109.
5. Глаубер, Р. Оптическая когерентность и статистика фотонов / Р. Глаубер. - М.: Мир, 1966. - 451 с.
6. Гудмэн, Дж. Статистическая оптика / Дж. Гудмэн. -М.: Мир, 1988. - 527 с.
7. Рогачёв, В.А. Классификация режимов работы фоторегистраторов в системах локации / В. А. Рогачёв // Инфокоммуникационные технологии. - 2008. - № 3. -C. 42-46.
8. Соболев, В. С. Максимально правдоподобные оценки параметров оптических сигналов с учётом дробового и фонового шумов / В.С. Соболев, Ю.А. Полещук // Компьютерная оптика. - 2006. - Т. 30. - С. 98-106.
9. Тихонов, В.И. Оптимальный приём сигналов / В.И. Тихонов. - М.: Радио и связь, 1983. - 320 с.
10. Гальярди, Р.М. Оптическая связь / Р.М. Гальярди, Ш. Карп; пер. с англ. - М.: Связь, 1978. - 424 с.
11. Федосеев, В.И. Приём пространственно-временных сигналов в оптико-электронных системах / В.И. Федосеев. - М.: У ниверситетская книга, 2011. - 232 с.
References
1. Lebetco, E.G. Optical location systems / E.G. Lebetco. -Saint-Petersburg: “SRU ITMO” Publisher, 2012. - 129 p.
- (In Russian).
2. Kozintsev, M.P. Foundations of impulse laser location / M.P. Kozintsev, V.M. Belov [et al.]. - Moscow: Bauman MSTU Publisher, 2006. - 512 p. - (In Russian).
3. Berezin, V.V. Point objects image coordinates detection and estimation in astronavigation and adaptive optics / V.V. Berezin, A.K. Tsitsulin // TOGU messenger. - 2008.
- N 1(8). - P. 11-20. - (In Russian).
4. Kutsov, R.V. Nonhomogenious space-lengthy object velocity vector detection and estimation according to its image / R.V. Kutsov, A.P. Trifonov // Automatics and telemechanics. - 2009. - N 8. - P. 96-109. - (In Russian).
5. Glauber, R. Optical coherence and photons statistics / R. Glauber. - Moscow: “Mir” Publisher, 1966. - 451 p. -(In Russian).
6. Goodman, D. Statistical optics / D. Goodman. - Moscow: “Mir” Publisher, 1988. - 527 p. - (In Russian).
7. Rogachev, V.A. Operating condition classification of photoregistrator in location systems / V.A. Rogachev // In-focommunication technology. - 2008. - N 3. - P. 42-46. -(In Russian).
8. Sobolev, V.S. Optical signals parameters maximum likelihood estimations in shot and background noises presence / V.S. Sobolev, U.A. Poleshuk // Computer optics. - 2006. -V. 30. - P. 98-106. - (In Russian).
9. Tichonov, V.I. Optimal signals recept / V.I. Tichonov. -Moscow: “Radio i svyaz” Publisher, 1983. - 320 p. - (In Russian).
10. Galiardi, R.M. Optical communication / R.M. Galiardi, Sh. Karp. - Moscow: “Svyaz” Publisher, 1978. - 424 p. -(In Russian).
11. Fedoseev, V.I. Space-time signal reception in optical-electronical systems / V.I. Fedoseev. - “Universitetskaya kniga” Publisher, 2011. - 232 p.- (In Russian).
OPTICAL SIGNAL DETECTION ACCORDING TO THE PHOTOELECTRON FLOW WITH UNKNOWN DENSITY WAVEFORM
V.I. Parfenov, V.S. Kirillov Voronezh State University
Abstract
The optical signal with unknown flow density optimal detection algorithm according to photoelectron counting statistics is offered. The reception characteristics is calculated by the method of mathematical modeling.
Key words: shot noise, poisson approximation, a priori unknown intensity waveform, maximum likelihood estimation, an ideal observer criterion, statistical modeling method.
Сведения об авторах
Парфёнов Владимир Иванович, 1961 года рождения. В 1983 году с отличием окончил Воронежский государственный университет (ВГУ) по специальности «Радиофизика и электроника». Доктор физико-математических наук (2002 год), профессор, профессор кафедры радиофизики ВГУ. Парфёнов В.И. - специалист в области приема и обработки информации. В списке научных работ В.И. Парфёнова порядка 140 статей, 1 монография, 6 авторских свидетельств и патентов.
E-mail: vip@phys. vsu. ru .
Vladimir Ivanovich Parfenov (b. 1961) graduated with honours from Voronezh State University (VSU), majoring in Radiophysics and electronics. He received his Doctor in Physics & Maths (2002) degrees, professor of radiophysics department. Parfenov V.I. is a specialist in information reception and processing. He is co-author of 140 scientific papers, 1 monograph and 6 inventions and patents.
Кириллов Владислав Сергеевич, 1986 года рождения, в 2009 году окончил Воронежский государственный университет (ВГУ) по магистерской программе «Статистическая радиофизика». Кириллов В.С. - аспирант кафедры радиофизики ВГУ. Область научных интересов: разработка и исследование оптических систем передачи информации. В списке научных работ В.С. Кириллова 5 статей.
E-mail: vkirillow@mail.ru .
Vladislav Sergeevich Kirillov (b. 1986) graduated from “Statistical radiophysics” master program of Voronezh State University (VSU) in 2009. Kirillov V.S. is a postgraduate student of VSU radiophysics department. His research interests is transmission information optical systems development and analysis. He is co-author of 5 scientific papers.
Поступила в редакцию 2 марта 2012 г.
