СПОСОБ ИЗМЕРЕНИЯ ЧАСТОТЫ СИГНАЛОВ ПОСЫЛОК РАДИОБУЕВ В КОСМИЧЕСКОЙ СИСТЕМЕ ПОИСКА И СПАСАНИЯ Российский патент 2016 года по МПК G01R23/00 

Описание патента на изобретение RU2592050C1

Область техники

Изобретение относится к области радиотехники и предназначено для использования в создаваемом в настоящее время международной организацией КОСПАС-САРСАТ среднеорбитальном сегменте космической системы поиска и спасения терпящих бедствия судов, летательных аппаратов, отдельных людей или групп [1, 2, 3].

Уровень техники

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

Вместе с тем, этот сегмент имеет и существенный недостаток - большие ошибки определения координат терпящих бедствие при использовании существующих аварийных радиобуев, которыми в настоящее время оснащены более 500 тысяч судов, самолетов, вертолетов и проводится бурное оснащение различных туристических групп, альпинистов, геологов, лесников, охотников и других людей, находящихся и работающих в сложных и опасных условиях. Замена этих буев на новые требует огромных затрат и длительного времени. Большие ошибки определения координат при использовании существующих буев связаны с их узкой полосой частот излучаемых буями сигналов которая не должна быть больше 1,6 кГц, что не позволяет получить хорошую точность измерения псевдодальностей.

Кратковременность (всего 160 мс) излучения буем чистой несущей и существенно худшие геометрические факторы в среднеорбитальном сегменте по равнению с низкоорбитальным не позволяют получить достаточную точность измерения частоты сигнала буя для того, чтобы получить достаточную точность измерения псевдоскорости, при использовании доплеровского метода определения координат.

Раскрытие изобретения

Предлагаемое изобретение позволяет устранить этот недостаток за счет существенного (в 4,5 раза) уменьшения ошибок измерения частоты сигналов посылок радиобуев.

Это позволяет, как показывают расчеты и моделирование, получить требуемую точность определения координат. Ошибка определения координат радиобуя не превышает 5 км с вероятностью 95%, в любой момент времени в любой точке земного шара.

При этом большинстве случаев точность определения координат будет значительно (в 2-3 раза) лучше.

Существенное повышение точности измерения частоты сигналов радиобуев в предлагаемом способе получается за счет использования для измерения частоты всей длительности сигнала посылки радиобуя (440 мс) вместо использования только немодулированной части этой посылки (160 мс).

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

Технический результат достигается тем, способ измерения частоты сигналов посылок радиобуев в космической системе поиска и спасания, отличающийся тем, что измерения частоты сигналов посылок излучаемых радиобуями, ретранслированными спутниками-ретрансляторами и принимаемых наземной станцией, производят с использованием всей длительности сигналов посылок и всей энергии их, для чего на наземной станции усиливают принятые сигналы и переносят их спектр в область нулевых частот, производят их преобразование в цифровой формат, полученные комплексные цифровые отсчеты принятого сигнала Z k = Z ( t k ) = S ( t k ) + N ( t k ) = A s k × e i 2 π f изм . + N k , - частота сигнала S ( t k ) ,k=0,±1,±2,...,±N, t k = k Δ t - моменты отсчетов сигнала; Δ t = 10 мкс; A s k = а s ( t k ) e i φ ;

N k = X k + i Y k ; a s ( t k ) - амплитуда сигнала буя; X k , Y k - независимые случайные величины с σ X 2 = σ Y 2 = N 0 × Δ f ш ; Δ f ш = 80 кГц - шумовая полоса приемника, подвергают трехэтапной обработке, при этом результаты предыдущего этапа используют на последующем этапе, на первом этапе обработки посылок сигналов используют быстрое преобразование Фурье для обнаружения сигналов посылок радиобуев в шумах, при котором грубо определяют несущие частоты f изм * принятых сигналов, далее вычисляют взаимную корреляционную функцию принятых сигналов с эталонным сигналом синхронизирующего слова, с которого начинаются все посылки сигналов всех радиобуев, при превышении полученного значения взаимной корреляционной функции над порогом, определяемым шумами, судят о факте обнаружения посылки сигнала радиобуя, а момент времени в который эта взаимная корреляционная функция сигналов достигает максимума равен моменту начала 25 бита модулированной части посылки сигналов, используя значения несущей частоты и момента начала 25 бита каждой обнаруженной посылки сигнала радиобуя, на втором этапе обработки производят демодуляцию сигналов, посылок содержащих информацию, которая определяет тип радиобуя, страну принадлежности радиобуя, номер и информацию о причине включения, декодируют посылку сигнала, которая закодирована кодом Боуза-Чоудхури-Хоквингейма (БЧХ кодом), при обнаружении в сигнале посылки не более 2-х ошибок, эти ошибки исправляют, получая достоверную информацию I(n), на третьем этапе обработки, используя информацию выделенную из достоверной посылки сигнала и определенные на первом этапе параметры сигнала, сигнала и определенные на первом этапе параметры сигнала, формируют опорный комплексно-сопряженный сигнал S 0 * ( t k ) = e i [ 2 π f и з м * t k + φ * + 1,1 M ( t k ) ] принимаемому сигналу S ( t k ) с амплитудой, равной единице, и подставленными в него грубыми значениями несущей частоты f изм * , грубым значением фазы посылки φ и двоичной цифровой последовательностью сигнала посылки радиобуя I(n), где

М ( t k ) - функция, определяемая информацией, заложенной в сигнал посылки радиобуя, М ( t k ) = { 0 п р и 22 000 k < 6000 B ( t k ) × [ 2 I ( n ) 1 ] п р и 6000 k 22 000 ,

где B ( t k ) - меандр с частотой, равной примерно 400 Гц и периодом, равным длительности одного бита посылки 2,5 мс = 250 Δ t ;

B ( t k ) = { 1 п р и t b n t k t b n + 125 Δ t 1 п р и t b n + 125 Δ t t k t b n + 250 Δ t ,

где

t b n - момент начала n-го бита сигнала посылки радиобуя;

t b 25 = t 0 - момент времени середины сигнала посылки радиобуя, принятый за нулевой;

t b n = t 0 6250 Δ t + 250 n Δ t ,

производят комплексное перемножение цифровых отсчетов Z k входного сигнала на значения S 0 * ( t k ) опорного сигнала, производят суммирование действительных значений k = N + N Re [ Z k × S 0 * ( t k ) ] и взвешенное суммирование мнимых значений полученных произведений с умножением на линейную функцию номера отсчета k k = N + N k Im [ Z k × S 0 * ( t k ) ] , где N равняется 22000 и вычисляют частоту сигнала посылки f изм = f изм + Δ F , где

Δ F = 6 2 π T N k = N + N k Im [ Z k × S 0 * ( t k ) ] k = N + N Re [ Z k × S 0 * ( t k ) ] .

Также в способе измерения частот сигналов используют суммирование сигналов нескольких посылок того же радиобуя, для чего комплексные цифровые отсчеты обнаруженных, но не достоверных посылок сигналов суммируют с комплексными цифровыми отсчетами других обнаруженных посылок сигналов, частоты и времени прихода 25 бита которых находятся в диапазонах возможных значений для этого радиобуя, после суммирования двух, трех или более посылок сигналов суммарный сигнал используют для получения достоверной информации I(n)посылки радиобуя для формирования опорных сигналов.

Основная идея способа состоит в следующем:

1. Восстановление чистой несущей на всей длине посылки буя, посредством фазовой модуляции принятого сигнала посылки достоверной информации, выделенной из нее, взятой с обратным знаком (ре-модуляции).

2. Измерении частоты этой ре-модулированной посылки оптимальным алгоритмом, дающим потенциальную точность измерения [4]:

σ f = 6 2 π × T × E N o (1)

где

Т - длительность сигнала;

E N 0 = H × T - отношение энергии сигнала к спектральной плотности шума;

H = P c N 0 - энергетический потенциал принятого сигнала;

P c - мощность сигнала.

Для лучшего понимания сути и полезности предлагаемого способа приведем краткое описание схемы проведения измерений частоты сигнала посылок радиобуев и использования этих измерений для определения координат радиобуев среднеорбитальным космическим сегментом поиска и спасания. На фигуре 1изображена функциональная схема сегмента, где обозначено:

1. НКА № 1…№ М - космические аппараты глобальных навигационных спутниковых систем - ГЛОНАСС, GPS, Gallileo;

2. РТР - ретрансляторы сигналов, излучаемых радиобуями;

3. НК - навигационные комплексы - аппаратура, служащая для формирования и передачи сигналов, излучаемых НКА;

4. А1- приемные антенны сигналов радиобуев;

5. А2 - передающие антенны ретранслированных сигналов радиобуев;

6. Анс - антенны навигационных сигналов;

7. СПОИ-СО - станция приема и обработки информации среднеорбитального сегмента;

8. ИИК-1…ИИК-М - информационно-измерительные комплексы;

9. КОВИ - комплекс обработки и выдачи информации;

10. КВЦ - координационно-вычислительный центр;

11. КВЦ других стран;

12. ПСС - поисково-спасательные службы;

13. fб - частота сигнала радиобуя;

14. f1, f2, ... fм - частоты сигналов радиобуя, приходящих на входы ретрансляторов НКА №1, НКА №2, НКА №М соответственно;

15. FРТР1, FРТР2, FРТРМ - частоты сигналов на выходах ретрансляторов, размещенных на НКА;

16. FНС1, FНС2, FНСМ - частоты навигационных сигналов, излучаемых НКА;

17. F1, F2, FМ - частоты сигналов ретрансляторов НКА, принимаемые СПОИ-СО;

18. FНС1, FНС2, FНСМ - частоты навигационных сигналов НКА, принимаемые СПОИ-СО.

На фигуре 2 представлены графики функций I(n), B(t), и M(t).

На фигуре 3 представлена Функциональная схема алгоритма обработки комплексных цифровых отсчетов с выхода ИИКа СПОИ-СО, где обозначено:

19. Блок обнаружения посылок;

20. Блок быстрого преобразования Фурье (БПФ);

21. Коррелятор;

22. Блок сравнения с порогом;

23, 24. Блоки разрешения на выдачу параметров посылки;

25. Блок демодуляции и выделения достоверной посылки;

26. Демодулятор;

27. Декодер Боуза-Чоудхури-Хоквингейма;

28. Блок уточнения частоты: f изм ;

29. Блок вычислений опорного сигнала;

30. Блок комплексного умножения;

31, 32. Блоки суммирования действительной и мнимой части комплексных произведений;

33. Блок вычисления уточненного значения частоты сигнала посылки;

34. Блок обмена информацией и управления вычислительным процессом.

На фигуре 4 представлены отсчеты амплитуд спектра, полученные БПФ в области максимума его огибающей.

Частота сигнала посылки радиобуя приходящая на m-й НКА, описывается выражением:

f m ( t ) = f б × ( t Δ t р а с п ) × ( 1 D ˙ m c 2 + V m 2 2 c 2 ψ ψ c c 2 ) , m=1, 2,…, M (2)

где

f б × ( t Δ t р а с п . ) = 406,05 М Г ц + Δ f б . o + δ f б × ( t Δ t р а с п . ) - частота, излучаемая буем;

Δ t расп . = D m c - время распространения сигнала буя до спутника-ретранслятора;

D m - дальность от буя до m-го НКА;

Δ f б . o - номинальное («литерное») значение частоты буя;

δ f б × ( t Δ t р а с п ) - мгновенное отклонение частоты буя от номинала из-за погрешности установки и всех видов нестабильностей;

D ˙ m = d D m ( t ) d t - радиальная скорость удаления спутника от буя;

V m 2 2 c 2 ψ ψ c c 2 - релятивистская поправка, первый её член за счет специальной, а второй за счет общей теории относительности;

Vm- модуль относительной скорости спутника относительно буя;

ψ , ψ c - гравитационные потенциалы Земли на её поверхности и на НКА соответственно;

с - скорость света.

В формуле (2) только одна величина фактически зависит от положения радиобуя - это D ˙ . Хотя формально релятивистская поправка V 2 2 c 2 тоже зависит от положения буя на поверхности Земли, абсолютная величина её влияния на f m не превышает 0,01 Гц и ею можно пренебречь.

D ˙ = d D ( t ) d t = ( V с V б ) × ( R с R б ) D , (3)

где

R с , R б - радиус-векторы спутника и буя;

V с , V б - векторы скорости спутника и буя;

D = | R с R б | - дальность от буя до спутника;

R с и V с - берутся на момент времени t, а R б и V б - на момент времени t t р а с п .

К сожалению, вычислению D ˙ m из формулы (2) мешает очень неточное (грубое) знание фактической частоты буя. Величина δ f б ( t ) может достигать сотен Гц.

Для устранения этой неприятности, как всегда делают в навигационных радиосистемах, используют разности измерений f m . :

Δ f m ,1 = f m f 1 , m = 2,... M (4)

В качестве первого (опорного) спутника можно выбирать любой.

Подставляя в формулы (4) выражения для соответствующих D ˙ m по формуле (3), получим систему уравнений с двумя неизвестными широтой(B) и долготой(L)

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

Для неподвижного буя достаточно иметь измерения f m относительно трех спутников, чтобы образовать две разности Δ f , для движущегося буя необходимо не менее шести спутников.

Используя измерения f m относительно трех спутников, образуем две разности f m .

Δ f 2,1 = 1 λ б × ( D ˙ 1 D ˙ 2 + V 1 2 2 c V 2 2 2 c + ψ с 2 ψ с 1 c )   Δ f 3,1 = 1 λ б × ( D ˙ 1 D ˙ 3 + V 1 2 2 c V 3 2 2 c + ψ с 3 ψ с 1 c ) } (5)

Ввиду малости релятивистских поправок и их близости для трех спутниковых систем (GPS, Глонасс, Galileo) их разностями в формулах (4 и 5) можно пренебречь.

Тогда получим:

D ˙ 1 D ˙ 2 = λ б × Δ f 2 ,1 D ˙ 1 D ˙ 3 = λ б × Δ f 3,1 } (6)

где

λ б = c f б .

В этой формуле f б достаточно знать с погрешностью не большей, чем 1 к Г ц . При этом ошибки вычисления ( D ˙ 1 D ˙ 2 ) и ( D ˙ 1 D ˙ 3 ) не превысят 5 мм/с.

Из изложенного следует, что для определения координат радиобуя необходимо знать только разности f m ( Δ f m ,1 ), а не сами f m . По этой причине ошибки измерения
f m , вызванные нестабильностью частоты буя ( δ f б ( t Δ t р а с п ) ) полностью исключаются.

Осуществление изобретения

Суть данного способа измерения частоты сигналов посылок радиобуев в космической системе поиска и спасания КОСПАС-САРСАТ на станции приема и обработки информации со среднеорбитальных ИСЗ систем «Глонасс», GPS и Gallileo которые оснащены ретрансляторами сигналов радиобуев заключается в следующем. Для уменьшения ошибок измерения частот принимаемых сигналов, вызванных собственными шумами приемников спутников-ретрансляторов и наземных станций, эти измерения производятся с использованием всей длительности сигнала посылки радиобуя (440 мс), а не только по участку длительностью 160 мс - участку излучения чистой несущей частоты радиобуя, и соответственно всей энергии сигнала.

Для этого производится модуляция принятых наземной станцией (станцией приема и обработки информации со среднеорбитальных ИСЗ систем «Глонасс», GPS и Gallileo) сигналов аварийных радиобуев достоверной цифровой информацией, заложенной в сигналы передаваемые тем же самым аварийным радиобуем и выделенной из принятого сигнала в процессе его демодуляции и декодирования, взятой с обратным знаком (ремодуляция сигнала). Это преобразует весь принятый сигнал посылки этого радиобуя в немодулированную синусоиду, чем и обеспечивается получение минимально возможной ошибки измерения его частоты.

Хотя, как было показано выше для определения координат радиобуев необходимы только разности f m ( Δ f m ,1 ), технически удобнее в приемниках каждого ИИКа СПОИ-СО измерять f m , а разности их получать при обработке этих измерений в алгоритме определения координат. По этой причине мы опишем способ измерения f= f m ,однако, при оценке ошибок измерения будем учитывать только ошибки, возникающие из-за собственных шумов ретрансляторов и, в меньшей степени, из-за шумов СПОИ-СО.

Здесь и далее при описании способа измерения частоты и алгоритма обработки сигнала его реализующего, с целью упрощения формул, мы будем опускать индекс «m» у сигналов и их параметров (f, φ и т.п.) обозначающий номер КА-ретранслятора и соответствующего ему ИИК'а так, как алгоритм обработки сигналов принимаемого любым ИИК'ом совершенно одинаков.

Как указывалось выше, основная идея способа измерения FOA состоит в следующем:

- восстановлении чистой несущей на всей длине посылки буя, за счет фазовой модуляции принятого сигнала посылки достоверной информацией выделенной из посылки, с обратным знаком (ремодуляции) и измерении частоты этой ре-модулированной посылки с использованием алгоритма, дающего потенциальную точность.

Приведем подробное описание этого способа.

Входной информацией для алгоритма измерения частоты f являются комплексные цифровые отсчеты сформированные аналого-цифровым приемником ИИКа СПОИ-СО из принятого им сигнала, передаваемого спутником-ретранслятором, гетеродинированного в область нулевых частот:

Z k = Z ( t k ) = S ( t k ) + N ( t k ) = A s k × e i 2 π f изм . + N k , k=0±1,±2,...,±N (7)

где

t k = k × Δ t - моменты отсчетов сигнала;

Δ t = 10мкс;

A s k = а s ( t k ) × e i φ ;

N k = X k + i Y k ;

a s ( t k ) - амплитуда сигнала буя;

a s ( t k ) = { a п р и 2,2 10 5 < t k 2,2 10 5 м к с 0 п р и о с т а л ь н ы х з н а ч е н и я х t k ;

X k , Y k - независимые случайные величины с σ X 2 = σ Y 2 = N 0 × Δ f ш ;

Δ f ш = 80 кГц - шумовая полоса приемника.

В этой формуле присутствует сигнал S ( t k ) только одного буя. На самом деле одновременно в Z k могут входить сигналы нескольких буев, а также нескольких помех. В дальнейшем описании алгоритма обработки отсчетов Z k мы для упрощения изложения также будем говорить о сигнале только одного буя. Необходимые указания об изменениях или дополнениях к алгоритму, возникающие из-за наличия сигналов нескольких буев и помех, не меняющие сути описываемых алгоритмов, будут даваться в тексте, по мере необходимости.

Частота сигнала S ( t k ) в комплексных отсчетах на выходе аналого-цифрового приемника ИИКа СПОИ-СО в формуле (7) - f изм  равна:

f изм = F ртр ( 1 D . c + Δ r e l . ) f гет . спои (8)

где

F ртр - частота, излучаемая спутником ретранслятором (НКА);

f гет . спои - частота гетеродина приемника СПОИ-СО;

Δ r e l = V 2 2 c 2 Ψ c Ψ c 2 - релятивистская поправка;

F ртр = f + f гет . ртр ;

f гет . ртр - частота гетеродина ретранслятора НКА.

Частоты гетеродинов ретрансляторов, установленных на НКА, формируется из частот высокостабильных, высокоточных генераторов навигационных сигналов НКА, а частота гетеродина приемника СПОИ-СО формируется из принятой на СПОИ-СО частоты навигационного сигнала, излучаемого этим спутником, и равна:

f гет . спои = ( 406,05 10 6 Г ц + f г е т . р т р ) × ( 1 D ˙ c + Δ r e l . ) . (9)

Подставляя это значение в (8), получим:

f изм = f 406,05 10 6 Г ц . (10)

Подставляя f= f m из формулы (2), получим:

f изм  = Δ f б . o 5 10 4 Г ц + δ f б 1 λ б × ( D ˙ c Δ r e l . ) (11)Возможные значения f изм лежат в области ± 40 к Г ц .

Фаза сигнала буя φ k равна:

φ k = φ ( t k ) = φ o + 1,1 р а д × М ( t k ) , (12)

где

φ o - начальная фаза сигнала, то есть фаза сигнала в момент времени начала 25 бита цифровой информации сообщения посылки буя ( t 0 ),

М ( t k ) - функция, определяемая информацией, заложенной в посылку буя.

М ( t k ) = { 0 п р и 22 000 k < 6000 B ( t k ) × [ 2 I ( n ) 1 ] п р и 6000 k 22 000 (13)

где

I ( n ) = { 1 0 - двоичная цифровая последовательность посылки буя,

n - номер бита этой последовательности,

B ( t k ) - меандр с частотой F t 400 Г ц и периодом, равным длительности одного бита 2,5 мс = 250 Δ t ;

B ( t k ) = { 1 п р и t b n t k t b n + 125 Δ t 1 п р и t b n + 125 Δ t t k t b n + 250 Δ t (14)

где

t b n - момент начала n-го бита;

t b 25 = t 0 - момент времени середины посылки, принятый за нулевой;

t b n = t 0 6250 Δ t + 250 n Δ t (15)

На фигуре 2 представлены графики функций I(n), B(t) и M(t). На этих графиках моменты времени t k не показаны.

Алгоритм обработки комплексных цифровых отсчетов Z k с выхода приемника СПОИ-СО можно разделить на четыре взаимодействующих блока, каждый из которых получает входные данные - поток Z k и (кроме первого блока) данные, полученные другими блоками в результате их работы. Входной поток Z k поступает на все блоки пакетами. Начало каждого пакета привязано к бортовому времени, для чего используется время, полученное от приемника навигационного сигнала, принимаемого с того же спутника, что и сигналы ретранслированных им радиобуев, тактовая частота опроса (100 кГц) также формируется из принятого навигационного сигнала. Таким образом, любое значение Z k является точно привязанным к бортовому времени спутника-ретранслятора.

На фигуре 3 представлена функциональная схема алгоритма обработки отсчетов Z k , показывающая их взаимодействие.

Блок обнаружения посылок (блок 19) состоит из двух алгоритмов, решающих задачи: обнаружения и предварительного (приближенного) измерения их параметров:

- алгоритма быстрого преобразования Фурье (блок 20, БПФ), который обнаруживает несущую частоту и фазу посылки и определяет их значения
( f изм , φ ). Здесь и далее верхним индексом * мы будем обозначать приближенные значения параметров посылки;

- алгоритма корреляционного анализа наличия синхрослова (первые 24 бита цифровой части посылки, одинаковые для всех посылок, блок 21) и определения
F Τ * и t 0 * .

Для получения более точных значений f изм и t 0 * используется квадратичная аппроксимация дискретных значений основного пика результата БПФ и корреляционной функции (см. фигуру 4).

Используются три значения амплитуды спектральных отчетов БПФ:

- наибольшее - a 1 с частотой f 1 ;

- отстоящее от него на - Δ по частоте - a 2 ;

- отстоящее от a 1 на + Δ по частоте - a 3 .

Значение частоты максимум a в спектре вычисляется по формуле

f и з м * = f 1 + δ = f 1 a 2 a 3 a 1 a 2 + a 3 2 ; (16)

Аналогично вычисляется и t 0 * .

Кроме этих функций блок обнаружения посылки определяет значение тактовой частоты посылки ( F Τ * ), измеряет отношение ( P c / N 0 ) * - сигнал/шум посылки и проверяет посылку на достоверность её длительности и положение временной метки t 0 * в середине посылки.

Полученное в результате работы алгоритмов обнаружения максимальное значение корреляционной функции сравнивается с заданным порогом (блок 22) и в случае превышения его, посылка считается обнаруженной, а ее параметры ( f , φ , F Τ * , t 0 * , ( P c / N 0 ) * ) поступают в блок обмена информацией и управления вычислительным процессом (блок 34). На фигуре 3 разрешение на выдачу параметров посылки условно обозначается ключами (блок 23 и блок 24).

Полученные в результате работы блока обнаружения приближенные значения параметров посылки f изм , φ 0 , t 0 * , F Τ * и ( P c / N 0 ) * используются последующими блоками в качестве необходимых исходных данных.

Блок обнаружения должен обнаруживать все посылки от всех буев, имеющиеся в потоке входных данных Z k .

Блок демодуляции и выделения достоверной посылки(блок 25), используя параметры f изм , φ 0 , t 0 * , F Τ * и ( P c / N 0 ) * производит демодуляцию посылки (блок 26), выделяя из неё информацию I * ( n ) и затем, с помощью декодера БЧХ (блок 27), исправляет в I * ( n ) все неверные биты. Если это исправление произошло правильно, то есть исправлено не более 2-х ошибок и, кроме того, в синхрослове все 9 бит, начиная с 16 и по 24, верные - посылка считается достоверной, и её можно обозначать I ( n ) (без звездочки). Если же в результате процедур декодирования БЧХ и проверки синхрослова не удалось получить достоверную посылку, то включается процесс накопления сигнала посылки путём суммирования соответствующих значений Z k , полученных от двух или более соседних по времени посылок.

Распознавание посылок от одного и того же буя производится по параметрам f изм и t 0 * . Демодуляция посылки, полученной в результате суммирования дает более достоверную информацию ( I * ( n ) ), которая после декодирования БЧХ и верификации синхрослова дает правильное I ( n ) .

Блок уточнения f изм (блок 28), уточняет f изм , получая в результате более точное значение f изм .

Значение f изм вычисляется для каждой обнаруженной посылки, в том числе и для тех посылок, по которым не получено правильной информации ( I * ( n ) I ( n ) ), а достоверная информация I ( n ) получена после суммирования сигналов нескольких посылок.

Полученная информационная последовательность I ( n ) может быть использована для вычисления f изм любой посылки этого буя, причем, в том числе,
в каналах приема сигналов этого буя через другие спутники-ретрансляторы. Такое использование I ( n ) существенно увеличивает число посылок, параметры f изм и t o которых могут быть измерены и, следовательно, повышает надежность и точность определения координат этого буя.

Уточненное значение f изм вычисляется по формулам:

f изм = f изм + Δ F , (17)

где

Δ F = 3 × k = N + N Im [ Z k S 0 * ( t k ) ] × k 2 π × a * × N × ( N + 1 ) × ( 2 N + 1 ) × Δ t 6 2 π × T × N × k = N + N k Im [ Z k × S 0 * ( t k ) ] k = N + N Re [ Z k × S 0 * ( t k ) ] , (18)

В этих формулах

Δ F = f и з м f изм * - уточнение предварительной оценки частоты f и з м *
(блок 33);

a * - оценка амплитуды сигнала буя:

a * = 1 2 N × k = N + N Re [ Z k × S 0 * ( t k ) ] , (блок 31) (19)

Re [ ... ] , Im [ ... ] - действительная и мнимая части комплексных произведений Z k × S 0 * ( t k ) ,

Z k ̶ комплексные цифровые отсчеты сигнала с выхода приемника (см. формулы (7), (11) и (12)),

k-номер отсчета Z k .

За нулевой номер (k = 0) принят момент времени t o ,

Δ t = 10 м к с - интервал квантования,

N = T 2 Δ t = 22000,

T = 2 N Δ t - длительность посылки буя (440 мс);

S 0 * ( t k ) - комплексно-сопряженное значение сигнала (с амплитудой, равной 1, блок 29) с подставленными в него значениями частоты f изм * , фазы φ и I ( n ) , то есть:

S 0 * ( t k ) = e i [ 2 π × f и з м * × t k + φ * + 1,1 × M ( t k ) ] , (20)

Блок 32 производит взвешенное (умноженное на k) суммирование мнимых частей комплексных произведений.

Формулы (18) и (19) дают достаточно высокую точность определения f изм и a * , если ошибка предварительной оценки f и з м не превышает 0,5 Г ц , что обеспечивается алгоритмом обнаружения, описанным выше, при P c N 0 30 д Б Г ц . При больших ошибках предварительного определения f и з м можно сделать второе приближение для вычисления Δ F 2 , использовав в формулах (17) - (20) вместо
f и з м - f и з м первого приближения.

Среднеквадратическая ошибка измерения f и з м за счет шумов ретранслятора и СПОИ-СО равна:

σ f = 6 2 π T ( P c / N 0 ) × T , (21)

а относительная С.К.О. амплитуды посылки равна:

σ a a = 1 2 ( P c / N 0 ) × T . (22)

При P c / N 0 = 30 д Б Г ц , σ f = 0,042 Г ц , а σ a a = 3,37 %

При этих оценках не учтены ухудшения, вызванные не идеальностью формы моделирующего сигнала, а также возникающие из-за ошибок t 0 * и F Τ * . Поэтому реальные значения σ F O A и σ a a в ~ 1,5 раза хуже, то есть 0,063 Гц, а ошибка измерения разности f и з м ( Δ f m ,1 ) будет σ Δ f и з м = 2 σ f = 0,09 Г ц

При больших ( P c / N 0 ) ошибки σ F O A и σ a будут, естественно уменьшаться.

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

Проведенные оценки этих ошибок показали, что в периоды большой солнечной активности С.К.О. ионосферных ошибок Δ f m ,1 могут достигать величины:

σ f и о н . = 0,05 Г ц .

Тогда суммарная величина вычисления разности f m ( Δ f m ,1 ) при P c / N 0 = 30 д Б Г ц может достигать:

σ Δ f с у м = 0,1 Г ц .

Список использованной литературы

1. А.И. Балашов, Ю.Г. Зурабов, Л.С. Пчеляков, В.И. Рогальский, В.С. Шебшаевич «Международная космическая радиотехническая система обнаружения терпящих бедствие» под ред. В.С. Шебшаевича, Москва, Радио и связь, 1987,

2. Ю.М. Урличич, Ю.Ф. Макаров, А.С. Селиванов, И.В. Никушкин, В.И. Рогальский, Ю.Г. Зурабов «История создания и перспективы развития международной космической системы поиска и определения местоположения терпящих бедствие судов и самолетов Коспас-Сарсат», Телекоммуникации и транспорт, 4-2012.

3. «MEOSARImplementationPlan (MIP) C/SR.012 (План внедрения системы СССПС (MEOSAR) 406 МГц КОСПАС-САРСАТ, ISSUE 1 - ред.10, октябрь 2014 г.

4. В.А. Котельников «Теория потенциальной помехоустойчивости», Москва, Госэнергоиздат М.-Л. 1956.

Похожие патенты RU2592050C1

название год авторы номер документа
СИСТЕМА ОПРЕДЕЛЕНИЯ МЕСТОПОЛОЖЕНИЯ ПОДВИЖНОГО ОБЪЕКТА ПО СИГНАЛАМ ГЛОБАЛЬНЫХ НАВИГАЦИОННЫХ СПУТНИКОВЫХ СИСТЕМ 2012
  • Баринов Сергей Прокопьевич
  • Непогодин Дмитрий Сергеевич
RU2490665C1
СПОСОБ СИНХРОНИЗАЦИИ ЧАСОВ И УСТРОЙСТВО ДЛЯ ЕГО РЕАЛИЗАЦИИ 2013
  • Ипатов Александр Васильевич
  • Дикарев Виктор Иванович
  • Койнаш Борис Васильевич
  • Смоленцев Сергей Георгиевич
RU2539914C1
СПОСОБ СИНХРОНИЗАЦИИ ЧАСОВ И УСТРОЙСТВО ДЛЯ ЕГО РЕАЛИЗАЦИИ 2013
  • Варганов Михаил Евгеньевич
  • Дикарев Виктор Иванович
  • Койнаш Борис Васильевич
  • Смоленцев Сергей Георгиевич
  • Жуков Евгений Тимофеевич
RU2535653C1
СПОСОБ СИНХРОНИЗАЦИИ ЧАСОВ И УСТРОЙСТВО ДЛЯ ЕГО РЕАЛИЗАЦИИ 2013
  • Ипатов Александр Васильевич
  • Дикарев Виктор Иванович
  • Койнаш Борис Васильевич
  • Смоленцев Сергей Георгиевич
RU2528405C1
СПОСОБ ОПРЕДЕЛЕНИЯ ПОЛОЖЕНИЯ ПОТРЕБИТЕЛЕЙ НАВИГАЦИОННОЙ ИНФОРМАЦИИ СПУТНИКОВЫХ НАВИГАЦИОННЫХ СИСТЕМ 2013
  • Себряков Герман Георгиевич
  • Щербаков Владимир Иванович
RU2517176C1
УСТРОЙСТВО ПЕЛЕНГАЦИИ ИССКУСТВЕННЫХ ИОНОСФЕРНЫХ ОБРАЗОВАНИЙ 2013
  • Цимбал Владимир Анатольевич
  • Пашинцев Владимир Петрович
  • Коваль Станислав Андреевич
  • Стрекозов Владимир Иванович
  • Сивоплясов Дмитрий Владимирович
  • Зинкин Александр Александрович
  • Горжий Василий Александрович
  • Коваль Станислав Андреевич
  • Форсов Георгий Львович
  • Шиманов Сергей Николаевич
RU2523912C1
БУЙ ДЛЯ ОПРЕДЕЛЕНИЯ ХАРАКТЕРИСТИК МОРСКИХ ВЕТРОВЫХ ВОЛН 2011
  • Жильцов Николай Николаевич
  • Грязин Дмитрий Геннадьевич
  • Зеньков Андрей Федорович
  • Аносов Виктор Сергеевич
  • Чернявец Антон Владимирович
  • Чернявец Владимир Васильевич
RU2490679C1
СПОСОБ ИНЕРЦИАЛЬНО-СПУТНИКОВОЙ НАВИГАЦИИ ЛЕТАТЕЛЬНЫХ АППАРАТОВ 2013
  • Бабуров Владимир Иванович
  • Гальперин Теодор Борисович
  • Иванцевич Наталья Вячеславовна
  • Саута Олег Иванович
  • Соколов Алексей Иванович
  • Юрченко Юрий Семенович
RU2536768C1
СПОСОБ СИНХРОНИЗАЦИИ ЧАСОВ И УСТРОЙСТВО ДЛЯ ЕГО РЕАЛИЗАЦИИ 2013
  • Варганов Михаил Евгеньевич
  • Смоленцев Сергей Георгиевич
  • Дикарев Виктор Иванович
  • Койнаш Борис Васильевич
  • Жуков Евгений Тимофеевич
RU2583894C2
Станция (система) приёма и обработки информации от среднеорбитального сегмента космической системы поиска и спасания и способ управления наведением антенн этой станции 2015
  • Федосеев Андрей Викторович
  • Селезнев Владимир Васильевич
  • Антонов Дмитрий Валентинович
  • Белоглазова Надежда Юрьевна
  • Дедов Николай Вадимович
  • Суринов Анатолий Серафимович
  • Семин Виктор Иванович
  • Архангельский Вячеслав Андреевич
  • Литвин Анатолий Иванович
RU2622390C2

Иллюстрации к изобретению RU 2 592 050 C1

Реферат патента 2016 года СПОСОБ ИЗМЕРЕНИЯ ЧАСТОТЫ СИГНАЛОВ ПОСЫЛОК РАДИОБУЕВ В КОСМИЧЕСКОЙ СИСТЕМЕ ПОИСКА И СПАСАНИЯ

Изобретение относится к области радиотехники и предназначено для использования в среднеорбитальном сегменте космической системы поиска и спасения терпящих бедствия судов, летательных аппаратов, отдельных людей или групп. Согласно способу измерения производятся с использованием всей длительности сигнала посылки радиобуя (440 мс), а не только по участку длительностью 160 мс - участку излучения чистой несущей частоты радиобуя, и соответственно всей энергии сигнала. Для этого производится модуляция принятых наземной станцией (станцией приема и обработки информации со среднеорбитальных ИСЗ систем «Глонасс», GPS и Gallileo) сигналов аварийных радиобуев достоверной цифровой информацией, заложенной в сигналы, передаваемые тем же самым аварийным радиобуем и выделенной из принятого сигнала в процессе его демодуляции и декодирования, взятой с обратным знаком (ремодуляция сигнала). Это преобразует весь принятый сигнал посылки этого радиобуя в немодулированную синусоиду, чем и обеспечивается получение минимально возможной ошибки измерения его частоты. Технический результат заявленного изобретения заключается в повышении точности измерений частоты сигналов радиобуев. 1 з.п. ф-лы, 4 ил.

Формула изобретения RU 2 592 050 C1

1. Способ измерения частоты сигналов посылок радиобуев в космической системе поиска и спасания, отличающийся тем, что измерения частоты сигналов посылок излучаемых радиобуями, ретранслированными спутниками-ретрансляторами и принимаемыми наземной станцией, производят с использованием всей длительности сигналов посылок и всей энергии их, для чего на наземной станции усиливают принятые сигналы и переносят их спектр в область нулевых частот, производят их преобразование в цифровой формат, полученные комплексные цифровые отсчеты принятого сигнала Z k = Z ( t k ) = S ( t k ) + N ( t k ) = A s k × e i 2 π f и з м + N k , f и з м - частота сигнала S ( t k ) , k = 0, ± 1, ± 2,..., ± N , t k = k × Δ t - моменты отсчетов сигнала; Δ t = 10 мкс; A s k = a s ( t k ) × e i ϕ ; N k = X k + i Y k ; a s ( t k ) - амплитуда сигнала буя; X k , Y k - независимые случайные величины с σ X 2 = σ Y 2 = N 0 × Δ f ш ; Δ f ш = 80 кГц - шумовая полоса приемника, подвергают трехэтапной обработке, при этом результаты предыдущего этапа используют на последующем этапе, на первом этапе обработки посылок сигналов используют быстрое преобразование Фурье для обнаружения сигналов посылок радиобуев в шумах, при котором грубо определяют несущие частоты f и з м * принятых сигналов, далее вычисляют взаимную корреляционную функцию принятых сигналов с эталонным сигналом синхронизирующего слова, с которого начинаются все посылки сигналов всех радиобуев, при превышении полученного значения взаимной корреляционной функции над порогом, определяемым шумами, судят о факте обнаружения посылки сигнала радиобуя, а момент времени, в который эта взаимная корреляционная функция сигналов достигает максимума, равен моменту начала 25 бита модулированной части посылки сигналов, используя значения несущей частоты и момента начала 25 бита каждой обнаруженной посылки сигнала радиобуя, на втором этапе обработки производят демодуляцию сигналов, посылок, содержащих информацию, которая определяет тип радиобуя, страну принадлежности радиобуя, номер и информацию о причине включения, декодируют посылку сигнала, которая закодирована кодом Боуза-Чоудхури-Хоквингема (БЧХ кодом), при обнаружении в сигнале посылки не более 2-х ошибок, эти ошибки исправляют, получая достоверную информацию I(n), на третьем этапе обработки, используя информацию, выделенную из достоверной посылки сигнала и определенные на первом этапе параметры сигнала, формируют опорный комплексно-сопряженный сигнал S 0 * ( t k ) = e i [ 2 π f и з м * × t k + ϕ * + 1,1 M ( t k ) ] принимаемому сигналу S 0 ( t k ) с амплитудой, равной единице, и подставленными в него грубыми значениями несущей частоты f и з м * , грубым значением фазы посылки ϕ * и двоичной цифровой последовательностью сигнала посылки радиобуя I(n), где
M ( t k ) - функция, определяемая информацией, заложенной в сигнал посылки радиобуя,
M ( t k ) = { 0 п р и 22000 k < 6000 B ( t k ) [ 2 I ( n ) 1 ] п р и 6000 k 22000 ,
где B ( t k ) - меандр с частотой, равной примерно 400 Гц, и периодом, равным длительности одного бита посылки 2,5 мс = 250 Δ t ;
B ( t k ) = { 1 п р и t b n t k t b n + 125 Δ t 1 п р и t b n + 125 Δ t t k t b n + 250 Δ t , где
t b n - момент начала n-го бита сигнала посылки радиобуя;
t b 25 = t 0 - момент времени середины сигнала посылки радиобуя, принятый за нулевой;
t b n = t 0 6250 Δ t + 250 n Δ t ,
производят комплексное перемножение цифровых отсчетов Z k входного сигнала на значения S 0 * ( t k ) опорного сигнала, производят суммирование действительных значений k = N + N Re [ Z k × S 0 * ( t k ) ] и взвешенное суммирование мнимых значений полученных произведений с умножением на линейную функцию номера отсчета k k = N + N k × Im [ Z k × S 0 * ( t k ) ] , где N равняется 22000 и вычисляют частоту сигнала посылки f и з м = f и з м * + Δ F , где
Δ F = 6 2 π × T × N k = N + N k × Im [ Z k × S 0 * ( t k ) ] k = N + N Re [ Z k × S 0 * ( t k ) ] .

2. Способ измерения частот сигналов по п.1, в котором используют суммирование сигналов нескольких посылок того же радиобуя, для чего комплексные цифровые отсчеты обнаруженных, но не достоверных посылок сигналов суммируют с комплексными цифровыми отсчетами других обнаруженных посылок сигналов, частоты и времени прихода 25 бита которых находятся в диапазонах возможных значений для этого радиобуя, после суммирования двух, трех или более посылок сигналов суммарный сигнал используют для получения достоверной информации I(n) посылки радиобуя для формирования опорных сигналов.

Документы, цитированные в отчете о поиске Патент 2016 года RU2592050C1

RU 94039162 A1, 10.09.1996
СПОСОБ ОПРЕДЕЛЕНИЯ КООРДИНАТ РАДИОПЕРЕДАТЧИКОВ С ИСПОЛЬЗОВАНИЕМ ГЛОБАЛЬНЫХ НАВИГАЦИОННЫХ СИСТЕМ И УСТРОЙСТВО ДЛЯ ЕГО РЕАЛИЗАЦИИ (ВАРИАНТЫ) 2005
  • Емельянов Геннадий Саулович
  • Колобанов Павел Михайлович
  • Кригер Виталий Борисович
  • Скрипкин Алексей Анатольевич
  • Терехов Юрий Павлович
  • Щербачев Владимир Александрович
RU2321017C2
US 7881663 B2, 01.02.2011
US 2014266871 A1, 18.09.2014.

RU 2 592 050 C1

Авторы

Петушков Александр Михайлович

Романов Евгений Олегович

Селезнев Владимир Васильевич

Архангельский Вячеслав Андреевич

Дедов Николай Вадимович

Даты

2016-07-20Публикация

2015-05-26Подача