Способ защиты пользователя радионавигационного приемника по отношению к аберрантным измерениям псевдорасстояний

Изобретение относится к области радиотехники, а именно к радионавигации, и может быть использовано в спутниковой радионавигационной системе. Технический результат заключается в обеспечении защиты пользователя радионавигационного приемника от аберрантных измерений псевдорасстояний. Для этого погрешность измерения детектируют при помощи статистического метода оценки на основании вычисления остатков измерений, что позволяет, в частности, автономно от любого наземного сегмента (то есть с использованием функции RAIM) повысить эффективность имеющегося в наличии приемника (называемого «первичным») без функции контроля целостности, детектировать возможные погрешности, искажающие входные измерения вычисления положения, за счет использования робастного статистического алгоритма оценки, то есть алгоритма, не подверженного влиянию погрешностей измерений, и с применением динамического критерия, и вычислять робастную коррекцию для положения, выдаваемого первичным приемником, с исключением любой такой детектируемой погрешности. 7 з.п. ф-лы, 1 ил.

 

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

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

Предшествующий уровень техники

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

Чтобы ограничить эти погрешности и обеспечить защиту пользователей, необходимо располагать средствами, позволяющими идентифицировать такие погрешности с целью их устранения и вычислять предельную погрешность положения в зависимости от имеющихся в наличии измерений с учетом требований целостности и непрерывности, связанных с контекстом использования (посадка летательного аппарата и т.д.). Эти средства лежат в основе устройства, содержащего функцию RAIM (Receiver Autonomous Integrity Monitoring - мониторинг целостности автономного приемника). Современные приемники GNSS (глобальная система спутниковой навигации), используемые в гражданской авиации, невозможно использовать без функции RAIM.

Современное оборудование с функцией RAIM характеризуется двумя проблемами:

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

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

Краткое описание существа изобретения

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

Способ защиты в соответствии с настоящим изобретением в самом общем случае отличается тем, что взвешивают N остатков от оценки вектора состояния во время оценки навигационного решения, при этом взвешивание осуществляют по оцениваемой статистике погрешностей этих измерений, используя робастные статистические оценки (такие как Least Trimmed Squares Estimator (оценка по методу наименьших квадратов), или Minimum Covariance Determinant Estimator (по минимальному ковариантному детерминанту), или M-оценка, или A-, D-, GM-, L-, MM-P, R-, S- или W-оценка, или MSTD, причем это взвешивание может быть полным или частичным. Если оно является полным, учет остатков является полным с весовым коэффициентом, равным 1 или равным 0 для исключения, и если оно является частичным, весовой коэффициент находится между 0 и 1. Случай полного взвешивания охватывает, кроме всего прочего, робастные методы “RAIM-MSTD” и “RAIM-LTS”, а случай частичного взвешивания охватывает робастные методы типа “RAIM-M-оценок”. Робастные статистические оценки сами по себе хорошо известны и описаны, например, в статье «Википедия», которую можно найти по следующему адресу: , а также в ссылках, указанных в конце этой статьи.

Согласно другому отличительному признаку изобретения, выбирают подгруппу из h остатков оценки вектора состояния с наименьшей дисперсией, применяют эту дисперсию, умноженную на коэффициент адаптации, а также среднее значение этой же подгруппы, чтобы усреднить и затем нормализовать N входов, которые после возведения в квадрат сравнивают со статистическим порогом и отбрасывают, если они превышают этот порог, при h=N-1 и/или N-2, при этом N является числом остатков. Сравнение со статистическим порогом предназначено для того, чтобы систематически отделять h остатков от оценки вектора состояния. Это и есть метод “RAIM-MSTD”.

Согласно еще одному отличительному признаку изобретения, выбирают подгруппу из N-h остатков оценки вектора состояния, при которой сумма этих измерений в квадрате является минимальной, оценивают среднее значение и дисперсию этой подгруппы и используют это среднее значение и эту дисперсию (предпочтительно последнюю умножают на коэффициент адаптации, превышающий 1), чтобы усреднить и затем нормализовать N входов, которые после возведения в квадрат сравнивают со статистическим порогом и отбрасывают, если они превышают этот порог. Это является методом “RAIM-LTS”.

Согласно еще одному отличительному признаку изобретения, N остатков оценки вектора состояния взвешивают при помощи итеративного метода, который при каждой итерации состоит в оценке среднего значения и дисперсии остатков оценки вектора состояния, в вычислении весового коэффициента при помощи функции взвешивания, входными данными которой являются усредненные и нормализованные по дисперсии остатки, и остатки умножают на эти весовые коэффициенты, при этом итеративный процесс прекращается, если сумма квадратов разностей между двумя последовательными остатками меньше порогового значения, при этом среднее значение и дисперсию (последнюю умножают на коэффициент адаптации), полученные в результате этого процесса, используют, чтобы усреднить и затем нормализовать N входов, которые после возведения в квадрат сравнивают со статистическим порогом и отбрасывают, если они превышают этот порог. Это является методом “RAIM-М-оценок”.

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

- формируют все подгруппы остатков с кардиналом h=N-1 и/или N-2, где N является числом остатков,

- вычисляют типовое отклонение для каждой подгруппы и определяют при наименьшем типовом отклонении σ_min, используемом для взвешивания остатков, h_min значение кардинала, при котором получают σ_min, Y_min вектор соответствующего кардинала h_min, и m_min среднее значение Y_min,

- вычисляют первый вектор, определяемый как: $$r_{vec\mathrm{,1}}=\frac{\sqrt{{(Y-m_{\min })}^2}}{{\sigma }_{\min }}$$

- векторы r_vec,1 сортируют в порядке возрастания, чтобы получить: $$r_{vec\mathrm{,1}}^{сорт}$$

- вычисляют коэффициент f, определяемый как: $$f={({\chi }^2)}^{-1}\left(\frac{h_{\min }}{N}\mathrm{,1}\right)$$ , где (χ²)(.,1) является обратной величиной распределения χ² при одной степени свободы,

- после этого осуществляют повторное взвешивание σ_min следующим образом:

σ_min,2= $$\sqrt{\frac{r_{vec\mathrm{,1}}^{сорт}(h_{\min })}{f}}*{\sigma }_{\min }$$

- вычисляют новый вектор остатков: $$r_{vec\mathrm{,2}}=\frac{\sqrt{{(Y-m_{\min })}^2}}{{\sigma }_{\min \mathrm{,2}}}$$

- определяют порог Т: $$T={({\chi }^2)}^{-1}(P\mathrm{,1})$$ , где Р является вероятностью, оптимальное значение которой необходимо определить вероятностным моделированием,

- каждый элемент вектора r_vec,2 сравнивают с Т, и если $$r_{vec\mathrm{,2}}(i)>T$$ , детектируют неисправность, и спутник i исключают.

Краткое описание чертежей

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

Подробное описание предпочтительных вариантов воплощения настоящего изобретения

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

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

Вычисление информации целостности далее содержит следующие этапы, начиная с вычисления остатков (5) измерений псевдорасстояний (обозначаемых Yi, где i является индексом спутника, а N является числом этих остатков), на которых:

1. Формируют все подгруппы остатков с кардиналом h=N-1 и/или N-2,

2. Вычисляют типовое отклонение каждой подгруппы (6). Наименьшее типовое отклонение обозначают σ_min. Его используют в качестве опорной величины (7). После этого определяют:

а. h_min значение кардинала, при котором получают σ_min,

b. Y_min вектор соответствующего кардинала h_min,

с. m_min среднее значение Y_min.

3. σ_min используют для взвешивания остатков. Вычисляют первый вектор: $$r_{vec\mathrm{,1}}=\frac{\sqrt{{(Y-m_{\min })}^2}}{{\sigma }_{\min }}$$

4. Для разных спутников первые векторы r_vec,1 сортируют в порядке возрастания, чтобы получить: $$r_{vec\mathrm{,1}}^{сорт}$$ . Кроме того, вычисляют коэффициент f, определяемый как: $$f={({\chi }^2)}^{-1}\left(\frac{h_{\min }}{N}\mathrm{,1}\right)$$ , где (χ²)(.,1) является обратной величиной распределения χ² при одной степени свободы.

5. После этого осуществляют повторное взвешивание σ_min следующим образом: σ_min,2= $$\sqrt{\frac{r_{vec\mathrm{,1}}^{сорт}(h_{\min })}{f}}*{\sigma }_{\min }$$

6. Вычисляют новый вектор остатков: $$r_{vec\mathrm{,2}}=\frac{\sqrt{{(Y-m_{\min })}^2}}{{\sigma }_{\min \mathrm{,2}}}$$

7. Определяют порог Т: $$T={({\chi }^2)}^{-1}(P\mathrm{,1})$$ , где Р является вероятностью, оптимальное значение которой необходимо определить методом вероятностного моделирования типа «Монте-Карло».

8. Каждый элемент вектора r_vec,2 сравнивают с Т. Если $$r_{vec\mathrm{,2}}(i)>T$$ , отмечают детектирование неисправности, и спутник i исключают (8).

9. После осуществления детектирования возможных ошибочных измерений, можно произвести вычисление районов защиты (10) обычным путем с учетом числа произведенных измерений.

Следует отметить, что Р, представляющее собой вероятность исключения неошибочного спутника, вместе с тем нельзя напрямую соотнести с вероятностью исключения, определенной нормами ИКАО. С другой стороны, действие по пункту 8 осуществляют отдельно для каждого спутника: поэтому нет детектирования без исключения.

Таким образом, выходы этого устройства обеспечивают:

- решение положения, в случае необходимости, скорректированное (11), если детектирована погрешность во входных данных измерений,

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

Способ в соответствии с настоящим изобретением позволяет, автономно от любого наземного сегмента (то есть с использованием функции RAIM):

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

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

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

- вычислять предельные погрешности положения в зависимости от имеющихся в наличии измерений и с учетом требований целостности и непрерывности, соответствующих контексту использования (например, в фазе посадки летательного аппарата). Этими предельными погрешностями являются пороговые значения, которые не должны быть превышены более одного раза для N измерений расстояний, например, при N=от 10⁵ до 10⁷,

- эффективность превышает эффективность стандартного алгоритма RAIM (использующего метод наименьших квадратов).

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

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

- Способ RAIM, который благодаря изобретению становится робастным, позволяет повысить надежность оценки точки местонахождения (географического положения) за счет повышения характеристик контроля целостности (детектирование и исключение) по сравнению со стандартным способом RAIM.

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

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

1. Способ защиты пользователя радионавигационного приемника по отношению к аберрантным измерениям псевдорасстояний, отличающийся тем, что он содержит следующие этапы вычисления информации целостности путем вычисления остатков (5) измерений псевдорасстояний, полученных на основании географического положения и временного сдвига (2), выдаваемых первичным радионавигационным приемником, и на основании измерений (3), используемых первичным радионавигационным приемником для получений этих значений, причем упомянутые остатки обозначены Y_i, где i является индексом спутника, на которых:
a) формируют подгруппы остатков с кардиналом h=N-1 и/или N-2, где N является числом остатков,
b) вычисляют стандартное отклонение для каждой подгруппы (6) и определяют при наименьшем стандартном отклонении σ_min, используемом в качестве опорного значения (7) для взвешивания остатков, h_min - значение кардинала, при котором было получено σ_min, Y_min - вектор соответствующего кардинала h_min и m_min - среднее значение Y_min,
c) вычисляют первый вектор остатков, определяемый как
$$r_{vec\mathrm{,1}}=\frac{\sqrt{{(Y-m_{\min })}^2}}{{\sigma }_{\min }};$$
d) векторы r_vec,1 сортируют в порядке возрастания, чтобы получить
$$r_{vec\mathrm{,1}}^{copm};$$
е) вычисляют коэффициент f, определяемый как
$$f={({\chi }^2)}^{-1}\left(\frac{h_{\min }}{N}\mathrm{,1}\right),$$
где (χ²)(.,1) является обратной величиной хи-квадратичного статистического распределения χ² при одной степени свободы;
f) после этого осуществляют повторное взвешивание σ_min следующим образом:
$${\sigma }_{\min \mathrm{,2}}=\sqrt{\frac{r_{vec\mathrm{,1}}^{сорт}(h_{\min })}{f}}\cdot {\sigma }_{\min };$$
g) вычисляют новый вектор остатков $$r_{vec\mathrm{,1}}=\frac{\sqrt{{(Y-m_{\min })}^2}}{{\sigma }_{\min }};$$
h) определяют порог Т: Т=(χ²)^-1(P,1), где Р является вероятностью, оптимальное значение которой определяют вероятностным моделированием;
i) каждый элемент вектора r_vec,2 сравнивают с Т и, если r_vec,2(i)>Т, отмечают детектирование неисправности и спутник i исключают.

2. Способ по п.1, отличающийся тем, что метод вероятностного моделирования, применяемый для определения порога Т, является методом типа «Монте-Карло».

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

4. Способ по п.1, отличающийся тем, что при детектировании погрешности во входных данных измерений вычисляют скорректированное положение (9).

5. Способ по п.1, отличающийся тем, что при детектировании погрешности во входных данных измерений вычисляют скорректированное значение времени (11).

6. Способ по п.1, отличающийся тем, что районы защиты (10), связанные с измерениями, вычисляют известным способом.

7. Способ по п.1, отличающийся тем, что при детектировании погрешности во входных данных измерений передают тревожный сигнал (12).

8. Способ по п.1, отличающийся тем, что способ применяют для радионавигационного приемника, изначально не имеющего источника робастной защиты от погрешностей измерения.