способ определения угловых координат гравитационно- волнового источника

Формула изобретения

Способ определения угловых координат гравитационно-волнового источника, основанный на измерении разности фаз сигнала от источника излучения, принятого в двух пространственно разнесенных точках, отличающийся тем, что при годовом орбитальном обращении Земли с гравитационно-волновым детектором последовательно в произвольно фиксированные моменты времени осуществляют N циклов когерентного накопления детектируемого гравитационно-волнового сигнала, а в каждом из N циклов с учетом движения Земли выделяют детектируемый гравитационно-волновой сигнал путем когерентного накопления его n периодов и запоминают фазу накопленного гравитационно-волнового сигнала в каждом N-м цикле, а между запомненными по выделенным в i- и j-x циклах фазами сигналов через интервал времени, равный времени орбитального прохождения Землей расстояния, соизмеримого с длиной волны детектируемого гравитационно-волнового сигнала и равного целому числу m периодов этого сигнала, измеряют разность фаз и по измеренной величине разности фаз определяют угловое положение источника гравитационного излучения относительно линии, соединяющей положения гравитационно-волнового детектора в пространстве при орбитальном движении Земли в моменты i- и j-x циклов когерентного накопления гравитационно-волнового сигнала, повторяют М раз измерение разности фаз между накопленными гравитационно-волновыми сигналами при последующих положениях гравитационно-волнового детектора при орбитальном движении Земли и осуществляют М раз определение углового положения гравитационно-волнового источника относительно той же перемещающейся в плоскости орбитального движения Земли линии и занимающей в этой плоскости М положений и по всем значениям углового положения гравитационно-волнового источника, определенных при годовом орбитальном перемещении гравитационно-волнового детектора, вычисляют угловые координаты источника детектируемого гравитационно-волнового сигнала относительно плоскости эклиптики и точки весеннего равноденствия.

Описание изобретения к патенту

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

Известен способ [1] отождествления источника по его гравитационному излучению - извлечение из принятого сигнала позиционной информации об источнике - основанный на схеме совпадений при функционировании сети гравитационных антенн. Способ предполагает использование большебазовых лазерно-интерферометрических и твердотельных антенн, предназначенных для детектирования импульсного гравитационного излучения от вспышечных источников. Это и является ограничением возможности способа.

Известен способ [2] определения угловых координат источника излучения, основанный на измерении разности фаз сигнала от источника излучения, принятого в двух пространственно разнесенных точках, который наиболее близок к заявляемому способу и может быть принят в качестве ПРОТОТИПА. Способ по прототипу предполагает измерение фазы принимаемого сигнала от источника излучения в двух разнесенных в пространстве (на Земле) точках и по разности между измеренными фазами и известному расстоянию (базы) между этими точками определение направления на источник излучения относительно линии, соединяющей эти точки, - фазоразностный метод. Точность измерения угловой координаты зависит от соотношения размера базы Б и длины волны принимаемого сигнала от источника, которое должно удовлетворять условию: Б . Этот способ может быть применен для определения угловых координат гравитационно-волновых источников. Однако в случае низкочастотных периодических гравитационно-волновых источников типа двойных релятивистских астрофизических объектов, у которых 10¹¹ м, упомянутое требование в земных условиях не выполнимо. И это является существенным недостатком прототипа.

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

Сущность изобретения заключается в том, что при определении угловых координат гравитационно-волновых источников по прототипу, основанный на измерении разности фаз сигнала от источника излучения, принятого в двух пространственно разнесенных точках, при годовом орбитальном обращении Земли с гравитационно-волновым детектором последовательно в произвольно фиксированные моменты времени осуществляют N циклов когерентного накопления детектируемого гравитационно-волнового сигнала, а в каждом из N циклов, с учетом движения Земли, выделяют детектируемый гравитационно-волновой сигнал путем когерентного накопления его n периодов и запоминают фазу накопленного гравитационно-волнового сигнала в каждом N-м цикле, а между запомненными по выделенным в i- и j-х циклах фазами сигналов через интервал времени, равный времени орбитального прохождения Землей расстояния, соизмеримого с длиной волны детектируемого гравитационно-волнового сигнала, и равного целому числу m периодов этого сигнала, измеряют разность фаз и по измеренной величине разности фаз определяют угловое положение источника гравитационного излучения относительно линии, соединяющей положения гравитационно-волнового детектора в пространстве при орбитальном движении Земли в моменты i- и j-х циклов когерентного накопления гравитационно-волнового сигнала, повторяют M раз измерение разности фаз между накопленными гравитационно-волновыми сигналами при последующих положениях гравитационно-волнового детектора при орбитальном движении Земли и осуществляют M раз определение углового положения гравитационно-волнового источника относительно той же, перемещающейся в плоскости орбитального движения Земли, линии и занимающей в этой плоскости M положений и по всем значениям углового положения гравитационно-волнового источника, определенных при годовом орбитальном перемещении гравитационно-волнового детектора, вычисляют угловые координаты источника детектируемого гравитационно-волнового сигнала относительно плоскости эклиптики и точки весеннего равноденствия.

В отличие от известного способа, где разность фаз измеряют между сигналом от источника, одновременно принятым в двух разнесенных точках (на Земле) с неперемещаемой базой в пространстве между ними, в заявляемом изобретении разность фаз измеряют между сигналом от источника, принятым последовательно во времени (при последовательном размещении в пространстве при орбитальном обращении Земли) гравитационно-волновым (ГВ) детектором. Кроме того, в заявляемом способе база, сформированная между одним и тем же ГВ-детектором, находящимся в разных точках пространства в разное время, перемещается в плоскости эклиптики относительно гравитационно-волнового источника (ГВИ). Последнее обстоятельство позволяет не только определить угловые координаты ГВИ (эклиптические широту и долготу) при максимально возможной базе на орбите (с учетом длины волны детектируемой ГВ), но и с максимальной возможностью уточнить угловые координаты при годовом обращении Земли.

При разработке заявляемого способа по определению угловых координат двойных релятивистских астрофизических объектов ГИ имелось в виду, что безразмерная амплитуда h ГВ-сигнала от этих источников меньше, как минимум на 2-3 порядка, по сравнению с амплитудой импульсного ГВ-сигнала от вспышечных источников. Отсюда вытекает необходимость их достаточно длительного накопления до требуемого отношения сигнал/шум. В [3] в качестве ГВ-детектора низкочастотных периодических ГВ-сигналов от двойных релятивистских астрофизических объектов с временем существования порядка 10¹³ с использованы компактные двухрезонаторные с геометрически и (или) пространственно неэквивалентными контурами с общей активной средой лазерные системы. Такие ГВ-детекторы представляют возможность в качестве метода оптимальной обработки периодических сигналов использовать метод совместного длительного межпериодного и внутрипериодного накопления с помощью гребенчатых фильтров. Этот метод обеспечивает также запоминание (фиксацию) начальной фазы накапливаемого периодического ГВ-сигнала на момент включения фильтров.

Схема, поясняющая процесс определения угловых координат, представлена на чертеже.

Определение угловых координат рассмотрим на примере низкочастотного периодического ГВ-источника PSR 15377+1155 (двойной релятивистский астрофизический объект [4] ), у которого период гравитационного излучения T_g = 2,892110³ с, а длина волны _g = 8,67810¹¹ м.

Принимаемый ГВ-сигнал запишем в виде [2]:

в точке 1: s₁(t) = hcos(_gt+_g),

в точке 2:

где _g(t) - случайная начальная фаза ГВ-сигнала с частотой _g и длиной волны _g;

h - безразмерная амплитуда ГВ-сигнала;

- направление на ГВИ, относительно линии, соединяющей точки 1 и 2;

Б - база (расстояние между положениями ГВ-детектора на орбите в точках 1 и 2);

точки 1 и 2 - положения (в пространстве) ГВ-детектора в моменты начала i- и j-го циклов накопления соответственно, в которых происходит измерение и запоминание фазы сигнала в произвольно фиксированные моменты времени. Для выделения (когерентного накопления) детектируемого периодического ГВ-сигнала (с учетом движения Земли) на фоне помех и запоминания его фазы используют межпериодный и внутрипериодный гребенчатые фильтры накопления [3].

Измеренная разность фаз _изм между запомненными фазами сигнала, измеренными ГВ-детектором в точках 1 и 2, будет равна [2]



откуда



Базу Б определим по формуле



где R_ОЗ = 1,49597876210¹¹ м - радиус орбиты Земли (ОЗ);

T_обр = 3,1556910⁷ с - время обращения (обр) Земли (1 год);

T_g, _g - период и длина волны детектируемого ГВ-сигнала, угловые координаты источника которого определяются;

m = 1, 2, 3 ... - количество периодов ГВ-сигнала на временной базе T_Б = mT_g.

В фиксированный момент времени t₁ (точка 1 на чертеже) в i-м цикле включают блок обработки сигналов ГВ-детектора [3] на время накопления T_н = nT_g = 6212,892710³ с 20,8 суток (n = 621 - количество периодов накапливаемого ГВ-сигнала в одном цикле) для получения на выходе ГВ-детектора отношения сигнал/шум q_вых.н = 10 (для достоверного выделения детектируемого ГВ-сигнала на фоне помех) и выключают в момент времени t₁+T_н. Метод обработки сигналов [3] с учетом априорной информации о движении Земли и угловых координат ГВИ в радио или оптическом диапазонах обеспечивает когерентное накопление ГВ-сигнала с запоминанием его фазы на входе ГВ-детектора на произвольно фиксированный момент времени t₁. При обращении Земли по орбите в произвольно фиксированные моменты времени осуществляют N циклов накопления (также с запоминанием начальной фазы ГВ-сигнала в каждом цикле) с временем накопления T_н = nT_g в каждом цикле, причем NnT_g = T_обр - t, где 0 t < T_g.

Разность фаз _изм между сигналом, принятым в точке 1 в момент времени t₁ в i-м цикле, и тем же сигналом, принятым через временную базу T_Б = mT_g в j-м цикле в момент времени t₂ (точка 2 на чертеже), зависит только от разности t₂-t₁ = mT_g. С учетом требования выполнения соизмеримости Б и _g получаем, что m = 1735. Тогда база Б между точками 1 и 2, вычисленная по вышеприведенной формуле, равна Б = 1,42410¹¹ м, при этом mT_g = 1,74210³ с 58,3 суток, а = arccos(0,97 _изм). Ошибка в измерении угла равна



где _изм - аппаратурная погрешность при измерении разности фаз. При /2 = 0,01 (что достижимо, поскольку инструментальная точность порядка 3% может быть реализована) получаем = 0,01_g/Б = 60 мрад. За время T_обр = 1 год может быть получено M = 17 раз независимых измерений разности фаз, что приведет к точности измерений в раз, т.е. _год = 60 мрад/ = 14,5 мрад.

Процедура вычисления координат из множества измерений состоит из двух этапов: определение предварительных значений координат и уточнение предварительных значений.

Определим предварительные значения координат. Пусть p, q и r - единичные векторы направлений, определенные M раз между попарными положениями ГВ-детектора в различные моменты времени t_i, t_j (см. чертеж). Если эти векторы некомпланарны, можно однозначно определить направление на источник x по измерениям углов _p, _q и _r, решив систему уравнений:

px = cos _p, qx = cos _q, rx = cos _r,

где _p, _q и _r - углы между вектором x и векторами p, q и r соответственно.

Поскольку детектор находится на Земле, то все его положения практически лежат в одной плоскости - плоскости эклиптики. В этом случае вместо единственного решения получается пара решений, симметрично расположенных относительно плоскости эклиптики, т.е. у них одинаковая эклиптическая долгота , а эклиптические широты различаются только знаком. Для нахождения этих решений достаточно иметь только два измерения - _p и _q.

Сориентируем прямоугольную систему координат так, чтобы первая координатная ось была направлена в точку весны, а третья - в северный полюс эклиптики. Для нахождения предварительных значений координат направления на источник принимают, что третьи компоненты векторов p и q равны нулю. Тогда система уравнений примет вид:

p₁x₁ + p₂x₂ = cos _p ;

q₁x₁ + q₂x₂ = cos _q.

Поскольку вектор x единичный, эту систему дополняют уравнением:



где знаки "+" и "-" относятся к решениям для северного и южного положений (относительно плоскости эклиптики) ГВ-источника соответственно.

Предварительные значения угловых эклиптических координат вычисляют по формулам



Для получения уточненных координат , оценим значения поправок , : = ₀+; = ₀+.

Связь между угловыми эклиптическими координатами направления на ГВ-источник , , координатами _i, _i, характеризующими направление i-го вектора, соединяющего два положения ГВ-детектора, и измеренным между этими направлениями углом _i, определяют из соотношения cos_i = sinsin_i+coscos_icos(-_i).

Линеаризация этого уравнения приводит к следующему условному уравнению: a_i+b_i+l_i = 0, где коэффициенты a_i, b_i и l_i имеют следующий вид:

a_i = -cos₀cos_isin(₀-_i),

b_i = cos₀sin_i-sin₀cos_icos(_o-_i),

l_i = sin_ii.

Невязка _i находится по схеме "измеренное минус вычисленное"



Решение системы условных уравнений по методу наименьших квадратов позволяет получить наиболее вероятные значения поправок , , оценки точности этих поправок и измерений .

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

Таким образом, заявляемый способ выгодно отличается от прототипа возможностью определения угловых координат низкочастотных (длинноволновых) периодических ГВ-источников путем формирования разности фаз между сигналами, принятыми от ГВ-источника последовательно во времени при последовательном размещении в пространстве ГВ-детектора при орбитальном перемещении Земли.

Источники информации

1. Милюков В.К., Руденко В.Н., // Итоги науки и техники ВИНИТИ АН СССР, серия Астрономия, 1991, т. 41, с. 147-193.

2. Пестряков В. Б. Фазовые радиотехнические системы. - М.: Сов. радио, 1968, 460 с. (прототип).

3. Balakin A.B., Murzakhanov Z.G., Skochilov A.F. Gravitation & Cosmology, Moscow, Vol. 5, N 3 (19), 1999

4. G.H. Taylor, R.N. Manchester, A.G. Lyne. //Astrophysical Journal Supplement, 1993, Vol. 88, p. 529.