Например, Бобцов

Пуассоновская модель звездного неба и задача обнаружения звезд оптико-электронным прибором

УДК 621.397.274; 621.391.1: 519.27
ПУАССОНОВСКАЯ МОДЕЛЬ ЗВЕЗДНОГО НЕБА И ЗАДАЧА ОБНАРУЖЕНИЯ ЗВЕЗД ОПТИКО-ЭЛЕКТРОННЫМ ПРИБОРОМ

© 2011 г. В. И. Федосеев, доктор техн. наук Научно-производственное предприятие “Геофизика-Космос”, Москва E-mail: coop@geofizika-cosmos.ru

Рассматриваются достоинства и недостатки использования пуассоновской модели звездного неба на примере задачи обнаружения звезд оптико-электронным прибором.

Ключевые слова: звездное небо, пуассоновская модель, обнаружение звезд.

Коды OCIS: 070.4560; 070.6110; 200.3050

Поступила в редакцию 23.06.2010

В современной оптической технике имеется класс аппаратуры, объектом наблюдения которой является звездное небо – отдельные звезды или группы звезд. При проектировании такой аппаратуры и исследовании ее характеристик возникает необходимость в представлении звездного неба в виде математической модели. Довольно часто в качестве такой модели используется пуассоновская модель [1, 2]. Согласно этой модели звездное небо представляется случайным полем точек-звезд, распределение которых по небесной сфере описывается законом Пуассона. Практика показывает, что такая модель может служить простым и удобным инструментом для выявления зависимостей и тенденций на качественном уровне и приближенных оценок; для точных количественных оценок она является слишком грубой. В настоящей работе на примере решения одной частной задачи обнаружения звезд оптико-электронным прибором демонстрируются достоинства и недостатки пуассоновской модели звездного неба.
Рассматривается следующая задача обнаружения. Оптико-электронным прибором проводится наблюдение участка звездного неба. Характеристика обнаружения прибора – зависимость вероятности обнаружения звезды от ее величины Pобн(m) – имеет вид, представленный на рис. 1. Для ярких звезд m ≤ m1 вероятность Pобн(m) близка к 1, для слабых звезд m ≥ m0 вероятность Pобн близка к 0, в зоне же m1 > m > m0 вероятность Pобн(m) существенно отличается как от 1, так и от 0. Эту зону назовем зоной

неуверенного обнаружения. Особенность рассматриваемой задачи в том, что протяженность этой зоны значительна и может составлять 1–2 зв. вел. Такая ситуация имеет место, например, в приборах ориентации космических аппаратов по звездам, где размер этой зоны складывается из ряда составляющих погрешностей энергетических измерений приборов, нестабильности их параметров в условиях эксплуатации, точности знания величин рабочих звезд и др. Прибор должен обнаружить определенное число рабочих звезд, необходимое для решения той или иной задачи (распознавание группы звезд – 3–4 звезды, определение углов трехосной ориентации – не менее 2-х звезд, углов двухосной ориентации – не менее одной звезды и др.). Необходимо определить вероятность обнаружения прибором требуемого количества звезд с учетом вышеуказанной особенности его характеристики обнаружения и того факта, что ориентация прибора относительно звездного неба заранее неизвестна.
Рассмотрение начнем в рамках пуассоновской модели звездного неба. Интенсивность
Pобн(m)
1
m1 m2 m Рис. 1. Типичный вид характеристики обнаружения звезд оптико-электронным прибором.

“Оптический журнал”, 78, 2, 2011

61

00 2 4 6 m
–2
–4
lg (m) Рис. 2. Зависимость среднего числа звезд в 1 кв. град небесной сферы от звездной величины.

наблюдаемого случайного поля в этой модели λ(m) представляет собой среднее число звезд в 1 кв. град небесной сферы, имеющих звездную величину от m до m + 1 (усреднение проводится по всей небесной сфере). На рис. 2 представлена зависимость λ(m), заимствованная из работы [1]. Обозначим телесный угол поля зрения прибора через Ω и выберем малый интервал звездных величин Δm таким, чтобы во всем интересующем нас диапазоне звезд вероятность Pпоп(m) попадания в поле зрения Ω звезды из диапазона m, m + Δm была мала, т. е.

Pпоп(m) = λ(m)ΔmΩ 0, точками mi, где i = 1, 2,…, k, на интервалы длиной Δm =

= mi + 1 – mi. Тогда полная вероятность выдачи прибором информации по крайней мере по

одной из звезд этого диапазона определяется

выражением

∏P1=1− k ⎡⎣1− P(m i )⎦⎤ =

i=1

∑=

1−

exp

⎡⎢⎣⎢

k i=1

ln(1−

P(mi

))⎥⎦⎤⎥.

(3)

Учитывая, что P(mi) 0.
Интегралу, стоящему в показателе экспо-
ненты выражения (4), можно дать наглядную
физическую интерпретацию. Если через Ω0 обозначить телесный угол, соответствующий
полной сфере, и положить Pобн(m) = 1 всюду в M, то величина

∫ ∫Ω0 Pîáí (m)λ(m)dm = Ω0 λ(m)dm = N
ÌÌ

будет полным количеством звезд в рабочем
диапазоне звездных величин прибора М. Если же Pобн(m) ≤ 1 в M, то величину Ω0∫Pобн(m)λ(m) dm можно назвать эффективным количеством звезд для прибора с характеристикой обнару-
жения Pобн(m) и обозначить Nэфф, где

∫Nýôô = Ω0 Pîáí (m)λ(m)dm.
Ì

(5)

Nэфф равно числу звезд, обнаруживаемых прибором, который имеет “прямоугольную” характеристику обнаружения: Pобн(m) = 1, m ≤ mэфф и Pобн(m) = 0, m > mэфф (здесь mэфф имеет такое значение, при котором число звезд с m ≤ mэфф равно Nэфф).
При “прямоугольной” характеристике об-
наружения количество рабочих звезд, которое
должно закладываться в приборный каталог, точно равно Nэфф. Если характеристика обнаружения отличается от “прямоугольной” (как
всегда бывает на практике), то для наиболее
полного использования обнаружительных воз-
можностей прибора число звезд в приборном
каталоге должно быть увеличено по сравнению с Nэфф так, чтобы все звезды , соответствующие зоне неуверенного обнаружения, вошли в этот
каталог.
С учетом сказанного, выражение для вероят-
ности выдачи информации прибором хотя бы
по одной звезде может быть записано в виде

62 “Оптический журнал”, 78, 2, 2011

P1

=

1



exp⎛⎜⎜⎜⎝−

Ω Ω0

Nýôô

⎟⎟⎞⎟⎠⎟.

(6)

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

ΩNýôô = Ω0ln(1− P1 )−1.

(7)

Аналогично можно получить выражение для вероятности выдачи информации по не менее, чем двум звездам

∏P2 =1− k ⎣⎡1− P(mi )⎦⎤ −

i=1

∑ ∏−

k l=1

P(ml

)i=k111−−

P(mi P(ml

)).

Выполняя преобразования, аналогичные пре-
дыдущему случаю, и учитывая, что 1 – P(mi) ≈ ≈ 1, можно получить следующее соотношение:

P2 =1−⎜⎝⎜⎛⎜1+

Ω Ω0

Nýôô

⎠⎞⎟⎟⎟⎟exp⎜⎜⎝⎛⎜−

Ω Ω0

Nýôô ⎟⎟⎞⎟⎠⎟.

(8)

Из выражений (6), (8) видно, что вероятности P1, P2 соответствуют пуассоновскому распределению с параметром Λ = (Ω/Ω0)Nэфф. Если же исходить из традиционного подхода к выбору
рабочих звезд, то в качестве параметра пуассоновского распределения следует брать Λ = = (Ω/Ω0)N1, где N1 – число звезд, соответствующее зоне уверенного обнаружения m ≤ m1, N1 < Nэфф. Таким образом при изложенном подходе за счет использования зоны неуверенно-
го обнаружения возможно увеличение числа рабочих звезд прибора с N1 до Nэфф с соответствующим увеличением вероятности выдачи
информации прибором. Для реализации тако-
го режима работы в приборный каталог в ка-
честве рабочих должны быть включены звезды с m ≤ m0, количество которых N0 > Nэфф > N1. При этом работа по каждой конкретной звезде
из зоны неуверенного обнаружения будет вы-
полняться не регулярно, но за счет большого
количества таких звезд вероятность выдачи
информации прибором повысится без измене-
ния его конструктивных параметров.
В общем случае для вероятности выдачи информации Pn по не менее, чем n звездам можно пользоваться соотношением

Pn =1− Ψ(n −1, Λ),

(9)

где Ψ(n – 1, Λ) – функция распределения вероятностей Пуассона с параметром Λ = (Ω/Ω0)Nэфф.

Функция распределения вероятностей Пуассона табулирована (см., например [3]).
Для реальной аппаратуры характеристика обнаружения Pобн(m) отличается от прямоугольной, причем форма ее в зоне неуверенного обнаружения трудно прогнозируема. Поэтому вычисление Nэфф по формуле (5) практически затруднительно, также как и определение соответствующей этому Nэфф предельной звездной величины mэфф. Поэтому в качестве Nэфф может быть принято то количество звезд, в котором для самой слабой звезды Pобн(mэфф) = 0,5. Пояснить это можно следующим рассуждением. При замене реальной характеристики обнаружения на “прямоугольную” с выбранным указанным способом mэфф отбрасываются звезды, обнаруживаемые с вероятностью менее 0,5. Звезды же, которые реально обнаруживаются с вероятностью более 0,5, но меньше 1, считаются обнаруживаемыми с вероятностью 1. Так как характер зависимости λ(m) экспоненциально возрастающий, то отбрасываемых звезд оказывается больше, чем звезд, у которых вероятность обнаружения завышается. Поэтому определенное таким образом Nэфф будет несколько меньше, чем в реальной аппаратуре. Другими словами, если за Nэфф принимать количество звезд, выбираемое из условия Pобн(mэфф) = 0,5, то фактическая вероятность выдачи информации прибором будет несколько больше, чем определяемая по формулам (6), (8), (9).
Изложенный качественный подход подтверждается результатами расчета вероятностей выдачи прибором информации, выполненного по формуле (6) для трех видов характеристик обнаружения, приведенных на рис. 3, и Ω = 320 кв. град, m1 = 3, m0 = 4, mэфф = 3,5 (в фотометрической системе прибора). Результаты расчета представлены в таблице.
Из данных расчета (второй столбец таблицы) видно, что, во-первых, для различных видов характеристики обнаружения с одним и тем же mэфф вероятности выдачи информации прибором отличаются мало, и, во-вторых, минимальное значение этой вероятности имеет место для прямоугольной характеристики обнаружения.
Для проверки степени приближения результатов, полученных при использовании пуассоновской модели звездного неба, к реальности было проведено математическое моделирование методом статистических испытаний по реальному расположению звезд на небесной сфере. Полученные при этом значения вероятно-

“Оптический журнал”, 78, 2, 2011

63

Pобн(m)
1

(а)

Pобн(m)
1

m1 (б)

m

Pобн(m)
1

m1 mэфф

m0 m

(в)

m1 mэфф

m0 m

Рис. 3. Три вида характеристик обнаружения

оптико-электронного прибора. а – прямо-

угольная, б – со спадом в – с линейным спадом.

вида

sin2

π 2

(m



m0

),

Вероятности выдачи информации по звезде

Вид Вероятности, Вероятности, характеристики вычисленные полученные при
обнаружения по формуле (6) моделировании

1

0,91

0,746

2

0,92

0,760

3

0,93

0,789

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

метра пуассоновского распределения. Но такая модификация лишила бы использование пуассоновской модели одного из достоинств – простоты. Поэтому необходимо иметь в виду, что простая пуассоновская модель дает оптимистические результаты и может использоваться только для ориентировочных оценочных расчетов. Для более точного определения характеристик приборов необходимо использовать статистическое моделирование с учетом реального расположения звезд на небе.
Из проведенного анализа можно сделать следующие выводы.
• Пуассоновская модель звездного неба является простым и наглядным средством для ориентировочной оценки характеристик оптикоэлектронных приборов при их работе по звездам (например, вероятностей выдачи информации прибором, формулы (6), (8), (9)).
• С помощью пуассоновской модели звездного неба продемонстрированы два интересных факта
– за счет использования зоны неуверенного обнаружения в характеристике обнаружения прибора возможно увеличение количества его рабочих звезд и вероятности выдачи им информации,
– для нахождения расчетных оценок вероятности выдачи информации прибором не требуется точного знания его характеристики обнаружения, достаточно знать величину звезды mэфф, для которой вероятность обнаружения равна 0,5.
• Пуассоновская модель звездного неба дает оптимистические оценки вероятностей выдачи информации оптико-электронным прибором, поэтому может использоваться только для ориентировочных расчетов. Для получения более точных количественных значений параметров приборов необходимо использовать статистическое моделирование с учетом реального распределения звезд на небесной сфере.
ЛИТЕРАТУРА
1. Астроследящие системы // Под ред. Чемоданова Б.К. М.: Машиностроение, 1980. 304 с.
2. Bezooijen R.W.H. Autonomous star referenced attitude determination // Guidance and control. 1989. AAS 89-003. P. 31–52.
3. Оуэн Д.Б. Сборник статистических таблиц. М.: ВЦ АН СССР, 1965. 586 с.

64 “Оптический журнал”, 78, 2, 2011