CC BY
13
УДК 519.246; 524
КВАЗИОПТИМАЛЬНЫЕ АЛГОРИТМЫ СТАТИСТИЧЕСКОГО АНАЛИЗА СОВПАДЕНИЙ В ГРАВИТАЦИОННО-ВОЛНОВОМ
ЭКСПЕРИМЕНТЕ
А. В. Гусев, В. Н. Руденко
(.ГАИШ)
E-mail: avg@sai.msu.ru, rvn@sai.msu.ru
В работе рассматриваются квазиоптимальные байесовские алгоритмы статистического анализа информации, полученной при эксплуатации пространственно разнесенных резонансных гравитационных антенн по схеме совпадений. Первичная обработка (накопление) гравитационных данных осуществляется по звездному времени.
1. Обработка информации, полученной с помощью двух пространственно разнесенных криогенных резонансных гравитационных антенн (РГА) с перекрывающимися диаграммами направленности, традиционно проводится по схеме совпадений [1]. Принцип действия схемы совпадений при обнаружении сильных и редких гравитационных импульсов основан на селекции (отборе) совпадающих во времени локальных максимумов выходных сигналов отдельных гравитационных антенн (предполагается, что амплитуды таких максимумов превышают высокий пороговый уровень). Фоновые совпадения рассматриваются как стационарный поток событий. В условиях априорной неопределенности средняя частота фоновых совпадений определяется по обучающей выборке.
Обработка информации по совпадениям [1] в гравитационно-волновом эксперименте носит чисто эвристический характер; статистика экстремальных значений пуассоновских потоков событий не учитывается. В предлагаемой работе статистический анализ совпадений в гравитационно-волновом эксперименте рассматривается с позиции байесовской задачи обнаружения и оценивания параметров ква-зидетерминированных сигналов [2]. Такой подход позволяет предложить специалистам, работающим в области гравитационно-волнового эксперимента, новые и эффективные алгоритмы статистической обработки информации. На этапе первичной обработки осуществляется накопление совпадений по звездному времени. Длительность интервала наблюдения — одни звездные сутки. Звездные сутки равны времени одного оборота Земли относительно точки весеннего равноденствия, момент верхней кульминации (прохождения через небесный меридиан при видимом суточном вращении небесной сферы) которой считается началом звездных суток (звездные сутки подразделяются на звездные часы, минуты и секунды). Применение звездного времени обеспечивает принципиальную возможность накопления гравитационных импульсов [1].
2. В режиме «быстрой фильтрации» («Fast filtering» [3]) выходной сигнал антенной решетки,
образованной двумя пространственно разнесенными РГА, можно представить в виде векторного случайного процесса у(t) = [yi(t)y2(t)]T, где ух(í) и |/2(í) — низкочастотные случайные процессы на выходах отдельных каналов. Пусть тр — время разрешения схемы совпадений. Тогда случайное событие
MC)-Í2¿(C)|<TP, ¿,¿ = 1,2,...,
где tn(C') и %(С) — моменты локальных максимумов, высоты которых превышают пороговый уровень С, можно рассматривать как совпадение случайных процессов yi(t) и При высоком
пороге С случайное число n(t) таких совпадений на интервале (0,t) представляет пуассоновский поток событий.
Как и в работе [1], будем предполагать, что интервал наблюдения (О, Т), длительность которого выбирается равной одним звездным суткам, разбивается на М = 24 отдельных подынтервалов. Длительность отдельного подынтервала равна одному звездному часу. Реалистическая модель такого потока n(t) определяется следующим выражением:
{пой при к = О, (М- 1), кф$,
(1)
Ans + при к =
где rife — случайное число совпадений на отдельном подынтервале (ífe, ífc+i), к = О, (М — 1), — фоновые совпадения, ns — сигнальные совпадения, А = (0,1) — параметр обнаружения, — неизвестный информативный параметр, подлежащий оценке.
Фоновые и сигнальные совпадения итальянской группой рассматриваются как статистически независимые случайные величины, распределенные по закону Пуассона:
Р{пок < m} = Fi(m,¿¿o), P{ns < m} = Fi(m,/xs),
(2)
где цо = (пofc) и fj,s = (ns) — среднее число фоновых и сигнальных совпадений на отдельном подынтерва-
ле ¿Ь ^+1) >
= У]р(г,ц) = -^-Г(т + 1, ц), *—4 т\
м
г=О
а'
р{т,ц) = —-ехр {-/л}, т\
(3)
Г(т + 1,/х) = J гт ехр
ч
— неполная гамма-функция [4]. В условиях априорной неопределенности для оценивания неизвестного параметра ¿¿о используется обучающая выборка.
Из (1), (2) и (3) находим условное отношение правдоподобия Л[п|/х8,1?] при обнаружении сигнальных совпадений пя по выборке п = (по, • • •, «м) [5]:
дм,,., *]==ехр {-».у ■
/V
(4)
3. При небайесовском алгоритме обнаружения неизвестные номер $ подынтервала, содержащего полезный сигнал, $ = О, М и среднее число сигнальных совпадений рассматриваются как неслучайные параметры. Оптимальный по обобщенному критерию максимального правдоподобия [2, 5] приемник (обнаружитель) формирует статистику
Л*[п] = тахЛ[п|^8,0] = -ехр{-^},
(5)
где в* и ц* — максимально-правдоподобные оценки неизвестных параметров 1? и . Решающее правило при небайесовском подходе можно представить в виде [2, 5]
1 при Л*[п]>са (1пЛ*[п] > 1пса), А=<{ (6)
О при Л*[п]^са (1пЛ*[п] < 1пса),
где са — пороговый уровень, зависящий от вероятности ложной тревоги а. С учетом (5) правило обнаружения (6) можно преобразовать к эквивалентной форме:
(по, • • •, пм) ^ Щ* = шах п п0
(7)
где па = [1пса + /х*]/[1п(/хо + /4) —^Мо] — пороговое число совпадений, зависящее от вероятности ложной тревоги а.
Пусть тр — квантиль одномерного распределения Пуассона ^(т,/хо): -Рх(тр,¿¿о) = Р- Тогда пороговое число совпадений па при небайесовском алгоритме обнаружения определяется следующей формулой:
па = т7, 7= (1 -а)1/м и 1 - -^ск.
Вероятность пропуска сигнальных совпадений при обнаружении по критерию максимального правдоподобия равна
/3 = Р{п^па\Х=1} =
При дополнительном условии А = 1 («сигнал есть») выражение (7) определяет оптимальный по критерию максимального правдоподобия алгоритм оценивания неизвестного информативного параметра Вероятность (Зу ошибочного решения при таком подходе равна:
/% = Р{Г = Л = 1} =
= Р {щ ^ щ = тах п|А = 1} = -Рх(%, \*>о +
4. При байесовском алгоритме обнаружения неизвестные параметры и 1? в формуле (4) рассматриваются как случайные величины с известными априорными распределениями [2, 5]. Решающее правило при байесовском обнаружении определяется выражением (6) при замене Л*[п] на безусловное отношение правдоподобия
м
Л[п] = <[п|^, 0])^ = (8)
к=О
Здесь
Ф(щ) = J к[п\ц,к]Ш9Г{ця)йц = о
1ш < М + Г и (9)
рк = Р{# = к}, к = 0,М,
где — априорная плотность вероятности
неизвестного неинформативного параметра . При частично байесовском подходе в качестве неизвестных распределений случайных параметров гравитационных импульсов выбираются равномерные распределения (постулат Байеса [2]):
1
ИргЫ = д->
А (10)
ц,€(0,А), р0 = ...=рм = (1/М).
Из (8), (9) и (10) в условиях параметрической априорной неопределенности находим:
ФЫ =
1
Ад
Пк
{//о} [Г(п*+1, цо+А) - Т(пк+1, щ>)] ■
5. При получившей на практике широкое распространение квазиоптимальной байесовской схеме раздельного обнаружения и оценивания параметров квазидетерминированного сигнала измерение информационного параметра 0 е (0, М) производится
после того, как на первом этапе принято решение А = 1 о наличии полезного сигнала. Оптимальный алгоритм оценивания зависит от выбранной функции потерь (платежа) П($, ё), где ё — оптимальная оценка неизвестного параметра рассматриваемого как дискретная случайная величина. При простой функции потерь вида П($, ё) = с\ — — ■&), с\ > О оптимальный по критерию максимума апостериорной вероятности алгоритм оценивания определяется следующим выражением:
ё = 2 при р,Ф(п,) = тах {рйФ(п^)}
(11)
р,1пФ(п,) = тах {рй1пФ(п^)}.
Так как функция Ф(п^) является монотонно возрастающей (Ф(п^1) < Ф(щ2) при Пк1 < щ2), то при равномерном априорном распределении р0 = ... = рм оптимальная оценка (11) оказывается оценкой максимального правдоподобия: ■&* = ■&. Оптимальный алгоритм обнаружения при равномерном априорном распределении определяется схемой
^ м
(п0,..., пм) ^ с>
к=О
где с — пороговый уровень, зависящий от выбранного критерия качества, существенно отличается от небайесовского алгоритма обнаружения по обобщенному критерию максимального правдоподобия (7). При квадратической функции потерь вида П($, ё) = ($ — ё)2 оптимальная по критерию минимума среднеквадратической ошибки оценка равна:
м м
Е кркф(пк) Е №(пк)
а _ к = 0_ _ к= 0_
" м м
Е РкЪ(пк) Е
к= 0 к= О
при р1 = ...=рм = (1/М).
6. Основные результаты и выводы. I. Ква-
зиоптимальная обработка выходных сигналов РГА
по схеме совпадений в соответствии с байесовской
задачей обнаружения и оценивания параметров ква-
зидетерминированных сигналов [2] осуществляется в два этапа. На первом этапе производится оптимальное обнаружение гравитационных импульсов. Априорные распределения неизвестных параметров ця (среднее число сигнальных совпадений на подынтервале и $ е (О, М) выбираются равномерными (постулат Байеса).
II. Оценивание информативного параметра § производится на следующем этапе после того, как будет принята гипотеза А = 1 «сигнал есть». Оптимальный (байесовский) алгоритм оценивания определяется выбранной функцией потерь П($, ё).
III. Реализация байесовских схем оптимального обнаружения и оценивания (при А = 1) представляет собой достаточно сложную в техническом отноше-
нии задачу. Наиболее простая схема обработки выборки п = (по,..., пм) соответствует небайесовскому подходу, основанному на применении обобщенного критерия максимального правдоподобия. При равномерном априорном распределении информативного параметра ê максимально-правдоподобная оценка ê* совпадает с оптимальной по критерию максимума апостериорной вероятности байесовской оценкой ё, ê* = ё.
IV. Схема обработки совпадений выходных сигналов РГА «Explorer» и «Nautilus», которая была использована итальянской группой [1], совпадает с оптимальным по критерию максимального правдоподобия алгоритмом совместного обнаружения гравитационных импульсов и оценивания (измерения) информативного параметра ê. Основной недостаток [1] связан с неправильной оценкой порогового числа совпадений па:
Fi(na, ß) = 1 — а (итальянская группа),
-Р\М(«а! ß) = l — a (критерий максимального *
правдоподобия),^
где ßa — полученная по обучающей выборке оценка среднего числа фоновых совпадений на отдельном подынтервале (tk,tk+1) (подобное замечание содержится в работе [6]).
V. Рассмотренные выше схемы обработки выходных сигналов РГА по совпадениям можно рассматривать как раздельное обнаружение сигнальных совпадений и оценивание информативного параметра ê (при А = 1) по оптимальным байесовским алгоритмам. Общая схема совместного (одновременного) байесовского обнаружения и оценивания параметров рассматривается в монографии [2]. При этом выигрыш оптимальной «в узком смысле» системы достигается за счет улучшения качества оценивания информативного параметра ê при недостоверном присутствии полезного сигнала (характеристики обнаружения сохраняются).
7. В заключение остановимся на анализе относительной эффективности байесовского обнаружителя по отношению к обнаружителю максимального правдоподобия. С целью упрощения при таком анализе ограничимся гауссовым приближением (асимптотика Муавра-Лапласа [4]). В гауссовом приближении характеристики обнаружения (вероятности ложной тревоги а и пропуска сигнала ß однозначно определяются отношением сигнал-шум
р и nsln\j2 (ns и fig » 1). Тогда, предполагая, что 1) ра = ... =рм = 1 /М и 2) W-pT(fis) =S(fj,g—a), при оптимальном (байесовском) алгоритме обнаружения имеем [7]:
р2« ^ylni - 1.4 + фп^ - 1.4 j +0.5 In M, (12) где а < 0.1 и /3 ^ 0.1.
При небайесовском подходе, основанном на применении обобщенного критерия максимального правдоподобия, характеристики обнаружения а и (3 определяются следующим выражением [7]:
—\ 2
In М + In -
1.4 •
а
In
1 1
1.4
(13)
Из (12) и (13) при дополнительном условии In ^ » 1.4, In ^ + In М » In jj получим
а
М
а
где а' и а" — вероятности ложной тревоги при байесовском и небайесовском подходах соответственно. Следовательно, при оптимальном (байесовском) алгоритме обнаружения вероятность ложной тревоги по отношению к алгоритму максимального правдоподобия заметно падает.
Литература
1. Astone Р., Babusci D., Bassati М. et al. // Class. Quantum Grav. 2002. 19. P. 5449.
2. Сосулин Ю.Г. Теоретические основы радиолокации и радионавигации. М., 1992.
3. Astone Р., Buttiglione S., Frasca S. et al. // Nuovo Cimento. 20C, N 1. P. 9.
4. Левин Б.Р. Статистическая радиотехника. М., 1994.
5. Тихонов В.И. Оптимальный прием сигналов. М., 1983.
6. Finn L.S. Ц gr-gc/0301092 v. 1. 23 Jan. 2003.
7. Гуткин Л.С. Оптимальные методы радиоприема на флюкту-ационных помехах. М., 1972.
Поступила в редакцию 17.11.03
