Синтез алгоритма оценки параметров многопозиционной системы наблюдения и исследование эффекта расходимости процесса фильтрации
- Авторы: Ерохин В.В.1, Лежанкин Б.В.1, Урбанский Д.Ю.2
-
Учреждения:
- Московский государственный технический университет гражданской авиации (Иркутский филиал)
- Восточно-Сибирский филиал ФГУП «ВНИИФТРИ»
- Выпуск: № 1 (2024)
- Страницы: 78–89
- Раздел: Системы наблюдения и организации воздушного движения
- URL: https://journals.rcsi.science/2312-1327/article/view/254498
- DOI: https://doi.org/10.51955/2312-1327_2024_1_78
- ID: 254498
Цитировать
Полный текст
Аннотация
Точность измерения местоположения воздушных судов (ВС) напрямую влияет на безопасность полётов и является одной из важнейших тактических характеристик. Внедрение новых перспективных средств наблюдения, таких как многопозиционные системы наблюдения (МПСН), может значительно повысить уровень безопасности полётов, а также улучшить эффективность использования воздушного пространства. В статье рассматривается задача улучшения качества функционирования МПСН и повышения точности оценки координат воздушных судов (ВС). Точность определения местоположения определяется погрешностью измерения времени прихода сигнала в условиях влияния шумов и помех. Случайные возмущения необходимо учитывать для обеспечения качественной работы МПСН, что достигается путем применения методов Калмановской теории фильтрации. Поэтому для решения задачи оценивания переменных состояния МПСН предлагается использовать фильтр Калмана (ФК). Эффективность применения фильтра Калмана зависит от адекватности математических моделей и реальных процессов. Неточности моделей, связанные с функционированием навигационных систем, приводят к расходимости ФК. В работе приведены результаты теоретических исследований и имитационного моделирования процессов функционирования МПСН на основе реализации алгоритма ФК.
Полный текст
Введение
Гражданская авиация является активно развивающейся отраслью промышленности и играет важную роль в обеспечении экономического развития Российской Федерации. Увеличение интенсивности воздушного движения послужило началом модернизации Единой системы организации воздушного движения (ЕС ОрВД). Для обеспечения эффективного функционирования системы УВД предлагается использовать многопозиционную систему наблюдения (МПСН).
Системный анализ принципов построения и функционирования МПСН показал значительные преимущества МПСН по сравнению с существующими классическими средствами наблюдения – возможность использования системы в труднодоступных районах для радиолокационных средств, а также простота обслуживания [Исследование…, 2023; Лежанкин и др., 2019].
Точность измерения местоположения воздушных судов (ВС) напрямую влияет на безопасность полётов и является одной из важнейших тактических характеристик. Внедрение новых перспективных средств наблюдения, таких как МПСН, может значительно повысить уровень эффективности использования воздушного пространства [Синтез…, 2022]. При этом важное значение имеет исследование точностных характеристик МПСН в зависимости от конфигурации параметров системы.
Точность определения местоположения ВС характеризуется качеством измерения времени прихода сигнала в условиях влияния шумов и помех [Арефьев и др., 2021; Применение…, 2021; Туринцев и др., 2022]. Случайные возмущения необходимо учитывать для обеспечения качественной работы МПСН, что достигается путем применения методов Калмановской теории фильтрации [Формирование…, 2008; Using…, 2022]. При синтезе алгоритмов оценки параметров МПСН, необходимо учитывать случайные возмущения для обеспечения качественной работы в условиях действующих шумов и помех, что достигается путем применения различных методов, в том числе, методов Калмановской теории фильтрации. Эффективность применения фильтра Калмана (ФК) существенно зависит от степени соответствия математических моделей реальным процессам. Неточности моделей, связанные с функционированием радионавигационных систем, приводят к появлению эффекта расходимости ФК.
На основе анализа вышеизложенного можно сделать вывод о том, что исследование точностных характеристик МПСН и повышение эффективности определения координат ВС в условиях воздействия случайных возмущений является актуальной научно-исследовательской задачей.
Цель работы – разработка алгоритма определения переменных вектора состояния многопозиционной системы наблюдения и исследование расходимости процессов фильтрации.
Постановка задачи
В основе принципа работы МПСН при определении положения приёмоответчика ВС лежат измерения на приемных пунктах значений псевдозадержки принятых сигналов, которые можно представить в виде:
,
где – принимаемый i-ым ПП сигнал с амплитудой Аi,; ω0i – несущая частота; φ0i – случайная начальная фаза; – псевдозадержка; – псевдодоплеровское смещение частоты; – координатная часть в виде преамбулы и информационное сообщение, содержащее данные о воздушном судне, параметрах его движения; – белый гауссовский шум (БГШ).
Измеренную псевдодальность представим в виде:
, (1)
где – расстояние от ВС до i-й приемной станции; – координаты ВС; – измеренное значение времени прихода сигнала от ВС до приемного пункта; с – скорость распространения радиоволн; – погрешность измерения дальности.
Все псевдодальномерные измерения можно объединить в вектор, состоящий из сигналов n видимых приемных станций:
, (2)
где – дискретные БГШ (ДБГШ) с нулевыми математическими ожиданиями и дисперсией ; – шаг дискретизации.
Синтез алгоритма
Применительно к задаче дискретной фильтрации, в предположении Марковского характера процесса , алгоритм вычисления значений апостериорной плотности вероятности (АПВ) имеет вид:
, (3)
, (4)
. (5)
Уравнения (3) – (5) позволяют рекуррентно вычислять значение АПВ на ν-м шаге по соответствующему значению той же плотности на предыдущем шаге. Начальные условия для такой рекуррентной процедуры описываются выражением (5).
При использовании метода максимума правдоподобия для обработки аддитивной смеси полезного сигнала и шума напряжение на выходе дискриминатора запишем в виде
,
где – принимаемые сигналы на интервале , вектор состояния – функция правдоподобия.
Применительно к рассматриваемой задаче АПВ представим в виде:
, (6)
На основании выражений (3) – (6) вычисляются оценки времени задержки и доплеровского сдвига частоты , соответствующие максимуму АПВ:
.
Модель динамики параметров МПСН представим в виде:
(7)
(8)
где T – время наблюдения; – ДБГШ с известной дисперсией ; – скорость изменения времени задержки сигнала; ν – временной индекс.
Уравнение наблюдения представим в виде:
(9)
где – ДБГШ с дисперсией .
На основе моделей (7) – (9) синтезирован алгоритм оценки параметров МПСН в скалярном виде:
,
,
где – коэффициенты усиления фильтра Калмана в выражениях для оценки соответствующих переменных; Rij – элементы главной диагонали матрицы ошибок R [Assad et al., 2019].
Матрица R вычисляется на основе выражений:
где – экстраполированная матрица.
При реализации перехода от скалярного вида к векторно-матричному получаем следующую форму записи уравнений (7) – (9):
Алгоритм оптимальной фильтрации вектора состояния представим в векторно-матричном виде:
, (10)
,
,
, (11)
Алгоритм позволяет решить задачу фильтрации параметров МПСН, состоящую в нахождении наилучшей оценки времени прихода сигнала по результатам псевдодальномерных измерений при наличии ошибок, вызванных действием случайных возмущений.
Результаты моделирования и исследования
В работе выполнено исследование характеристик алгоритма оптимальной фильтрации параметров МПСН на основе специализированного программного обеспечения. В качестве исходных данных выбраны значения параметров, которые представлены в [Исследование…, 2023; Лежанкин и др., 2019].
Для оценки качества измерений разности времён прихода переданного приемоответчиком объекта сигнала на рис. 1 представлены ошибки оценки и СКО , на рис. 2 приведены графики ошибки оценки и среднеквадратического отклонения . Представленные на рисунках 1, 2 результаты показывает достаточно высокую точность оценки параметров сигналов МПСН при совпадении фактических и расчетных параметров фильтра Калмана.
Рисунок 1 – Ошибка оценки времени задержки
Рисунок 2 – Ошибка оценки скорости изменения времени задержки
В то же время при работе фильтра Калмана может наблюдаться расходимость процессов фильтрации, как показано на графиках, приведенных на рисунках 3 и 4, при несоответствии между принятыми моделями и реальными процессами. Такое несоответствие, как показывают результаты исследований, представленные на рисунках 3 и 4, приводит к расходимости процессов фильтрации, т. е. к отличию между расчетным значением дисперсии ошибки оценивания и действительной ошибкой.
Расходимость возникает из-за того, что вероятностные характеристики шума наблюдений неизвестны и не учитываются в выражениях (4), (5) алгоритма фильтрации, при этом возникают расхождения между реальными процессами и моделями, принятыми в фильтре.
Расходимость процессов фильтрации приводит к увеличению дисперсии ошибки параметров МПСН, что демонстрируется на рисунках 5, 6.
Рисунок 5 – Дисперсия ошибки оценки времени задержки: 1 – при корректной работе; 2 – при расходимости.
Рисунок 6 – Дисперсия ошибки оценки скорости изменения времени задержки: 1 – при корректной работе; 2 – при расходимости.
Анализ представленных результатов показывает необходимость исследования эффекта расходимости процессов фильтрации и разработки новых подходов к оценке параметров МПСН для устранения данной проблемы, что и предопределяет актуальность выбранного направления исследований.
Анализ проблемы расходимости процессов фильтрации
При использовании алгоритмов Калмановского типа при обработке информации в МПСН возникают ситуации, когда действительные ошибки фильтрации превышают теоретически рассчитанные значения. Это явление называется расходимостью (неустойчивостью) фильтра Калмана. Причины расходимости заключаются в неточном задании моделей процессов сообщения и наблюдения, отсутствии достоверной информации о реальных физических процессах. Расходимость вызывают ошибки моделирования шумов в условиях статистической неопределенности их вероятностных характеристик [Болелов, 2021].
Основная причина расходимости процесса фильтрации – быстрое стремление к нулю коэффициента передачи (усиления) К фильтра Калмана, как показано на рисунках 7, 8. При этом процесс оценивания перестает быть зависимым от обновляющей информации, содержащейся в последовательности наблюдений, и рост ошибки не влияет на формирования оценки.
Рисунок 7 – Коэффициент усиления фильтра
Рисунок 8 – Коэффициент усиления фильтра
Следует отметить, что эффект расходимости – явление сложное и противоречивое, требующее проведения глубоких исследований применительно к решению конкретной задачи повышения эффективности функционирования МПСН.
В работе [Сейдж и др., 1976] рассматриваются методы устранения расходимости применительно к различным задачам. В основе рассмотренных методов лежит идея ограничения коэффициента усиления, чтобы избежать «нечувствительности» процесса фильтрации к вновь поступающей последовательности наблюдений.
Для уменьшения матричного коэффициента усиления предлагается использовать процедуру, состоящую в ограничении элементов матрицы снизу некоторой заранее выбранной величиной, либо можно добавить к элементам матрицы некоторую величину [Марковские модели…, 2019]. Например, в работе [Schmidt et al., 1968] для расчета матричного коэффициента усиления предлагается использовать следующее выражение:
.
Предложенный подход можно использовать для ограничения элементов матрицы дисперсий ошибок . В работе [Andrews, 1968] был развит другой подход, позволяющий контролировать нестабильность, который состоит в вычислении квадратного корня из матрицы . При рассмотрении данного способа вносится предположение, что шум модели равен нулю. Уравнение для матрицы ковариаций при записывается в виде:
,
где – экстраполированная матрица ковариаций.
Показано, что корень из матрицы равен
;
,
где используется нижняя треугольная матрица корней из дисперсий. Для случая непрерывного времени аналогичные выражения были предложены в [Andrews, 1968].
Перспективным направлением устранения эффекта расходимости является применение адаптивных алгоритмов фильтрации. Применение адаптивных алгоритмов представляется более эффективным ввиду того, что из-за расходимости величина ошибки становится относительно большой, что даёт возможность использовать большее количество информации для адаптации [Воскобойников, 2015]. Разработка и исследование алгоритмов адаптивной фильтрации параметров МПСН рассматривается в качестве приоритетного направления дальнейших исследований.
Заключение
Разработан алгоритм оценки переменных вектора состояния МПСН при использовании АЗН-В: времени, скорости и ускорения изменения времени задержки прихода сигнала на основе Калмановской теории фильтрации. Проведены исследования точностных характеристик МПСН наблюдения с применением специализированного программного обеспечения. Анализ представленных результатов показывает высокую точность оценки параметров при совпадении фактических и расчетных параметров фильтра Калмана. В то же время при несоответствии между принятыми моделями фильтра Калмана и реальными процессами на выходе алгоритма наблюдается расходимость процессов фильтрации, которая проявляется в отличии между расчётной и действительной ошибкой оценивания. Показано, что расходимость процессов фильтрации приводит к увеличению дисперсии ошибки оценки параметров МПСН. Обоснована необходимость разработки новых подходов к оценке переменных вектора состояния для устранения расходимости. Рассмотрены методы устранения неустойчивости функционирования разработанного алгоритма, при этом разработка и исследование алгоритмов адаптивной фильтрации параметров МПСН выбраны в качестве приоритетного направления дальнейших исследований.
Об авторах
Вячеслав Владимирович Ерохин
Московский государственный технический университет гражданской авиации (Иркутский филиал)
Автор, ответственный за переписку.
Email: ww_erohin@mail.ru
ORCID iD: 0000-0002-5549-3952
доктор технических наук
Россия, 664047, Иркутск, ул. Коммунаров, 3Борис Валентинович Лежанкин
Московский государственный технический университет гражданской авиации (Иркутский филиал)
Email: lezhbor@mail.ru
ORCID iD: 0000-0001-5504-0884
кандидат технических наук, доцент
Россия, 664047, Иркутск, ул. Коммунаров, 3Дмитрий Юрьевич Урбанский
Восточно-Сибирский филиал ФГУП «ВНИИФТРИ»
Email: urbdim.ru@gmail.com
ORCID iD: 0009-0000-2574-512X
Россия, 664056, Иркутск, Бородина, 57
Список литературы
- Арефьев Р. О. Применение помехоустойчивого кодирования при обработке сообщений локальной корректирующей станции / Р. О. Арефьев, С. В. Туринцев, М. С. Туринцева // Актуальные проблемы и перспективы развития гражданской авиации: сборник трудов X Международной научно-практической конференции, Иркутск, 14–15 октября 2021 года. Том 2. Иркутск: Иркутский филиал федерального государственного бюджетного образовательного учреждения высшего образования «Московский государственный технический университет гражданской авиации», 2021. С. 22-32. EDN YMDPPC.
- Болелов Э. А. Методы и алгоритмы комплексной обработки метеоинформации при метеорологическом обеспечении полетов воздушных судов гражданской авиации: специальность 05.22.14 Эксплуатация воздушного транспорта: диссертация на соискание ученой степени доктора технических наук / Болелов Эдуард Анатольевич, 2021. 421 с. EDN RAOAVI.
- Воскобойников Ю. Е. Критерий расходимости и алгоритм адаптации рекуррентного алгоритма оценивания вектора состояния / Ю. Е. Воскобойников // Научный вестник Новосибирского государственного технического университета. 2015. № 3(60). С. 7-22. doi: 10.17212/1814-1196-2015-3-7-22. EDN UNEYLH.
- Исследование точностных характеристик широкозонной многопозиционной системы наблюдения Иркутского регионального центра организации воздушного движения / Э. А. Болелов, Б. В. Лежанкин, М. А. Межетов, В. В. Ерохин // Вестник Санкт-Петербургского государственного университета гражданской авиации. 2023. № 3(40). С. 89-101. EDN LKMFMX.
- Лежанкин Б. В. Системный анализ задачи определения местоположения воздушного судна в многопозиционной системе наблюдения / Б. В. Лежанкин, В. В. Ерохин, В. С. Марюхненко // Информационные технологии и математическое моделирование в управлении сложными системами. 2019. № 1(2). С. 46-61. EDN ZDOOGT.
- Марковские модели технического состояния перспективных бортовых радиолокационных систем воздушных судов / А. С. Богачев, Э. А. Болелов, А. Т. Кудинов [и др.] // Научный вестник ГосНИИ ГА. 2019. № 26. С. 113-125. EDN YVVCPC.
- Применение технологии LoRa в беспилотных авиационных системах / М. А. Межетов, А. И. Тихова, У. С. Вахрушева, А. В. Федоров // Актуальные проблемы и перспективы развития гражданской авиации: сборник трудов X Международной научно-практической конференции, Иркутск, 14–15 октября 2021 года. Том 2. Иркутск: Иркутский филиал федерального государственного бюджетного образовательного учреждения высшего образования «Московский государственный технический университет гражданской авиации», 2021. С. 180-185. EDN UZUCLT.
- Сейдж Э. П. Теория оценивания и ее применение в связи и управлении / Э. П. Сейдж, Дж. Мелс [пер. с англ. под ред. Б.Р. Левина]. М.: Связь, 1976. 496 с.
- Синтез многопозиционных радиолокационных систем на базе сети специализированных излучателей / А. В. Журавлев, В. В. Кирюшкин, А. В. Коровин, Д. И. Савин // Успехи современной радиоэлектроники. 2022. Т. 76, № 4. С. 47-55. doi: 10.18127/j00338486-201807-21. EDN OQAAER.
- Туринцев С. В. Программная реализация алгоритма кодирования и декодирования местоположения ВС в дискретно-адресном режиме вторичной радиолокации / С. В. Туринцев, М. С. Туринцева // Актуальные проблемы и перспективы развития гражданской авиации : Сборник трудов XI Международной научно-практической конференции. посвященной празднованию 100-летия конструкторского бюро "Туполев", 55-летия Иркутского филиала МГТУ ГА, 75-летия Иркутского авиационного технического колледжа, Иркутск, 13–14 октября 2022 года. Том 2. Иркутск: Иркутский филиал федерального государственного бюджетного образовательного учреждения высшего образования «Московский государственный технический университет гражданской авиации», 2022. С. 115-121. EDN YGUUUY.
- Формирование радиолокационной карты подстилающей поверхности путем фильтрации случайных полей / О. Н. Скрыпник, Б. В. Лежанкин, Б. М. Миронов, Н. П. Малисов // Научный вестник Московского государственного технического университета гражданской авиации. 2008. № 133. С. 60-66. EDN KVVEYT.
- Andrews A. A square root formulation of the Kalman covariance equations // AIAA Journal. 1968. Vol. 6. pp. 1165-1166. doi: 10.2514/3.4696.
- Assad A. Novel Adaptive Fuzzy Extended Kalman Filter for Attitude Estimation in Gps-Denied Environment / A. Assad, W. Khalaf, I. Chouaib // Gyroscopy Navig. 2019. Vol. 10. pp. 131-146. doi: 10.1134/S2075108719030027.
- Schmidt S. F. (1968) Case Study of Kalman Filtering in the C-5 Aircraft Navigation System / S. F. Schmidt, J. P. Weinberg, J. S. Lukesh // Joint Automatic Control Conf. Univ. of Michigan. 1968. pp. 59-109.
- Using a MLAT Surveillance System to Locate Unmanned Aerial Vehicles Flying as a Swarm / E. A. Bolelov, B. V. Lezhankin, V. V. Erokhin, S. A. Zyabkin // 2022 XIX Technical Scientific Conference on Aviation Dedicated to the Memory of N.E. Zhukovsky (TSCZh). 2022. pp. 67-70. doi: 10.1109/TSCZh55469.2022.9802475.