CC BY
45
УДК 621.396.96.001(07)
Н. П. Богомолов
ДЕЦЕНТРАЛИЗОВАННЫЕ АЛГОРИТМЫ ОБРАБОТКИ ИНФОРМАЦИИ В ДВУХКАНАЛЬНЫХ ИЗМЕРИТЕЛЬНЫХ СИСТЕМАХ
Рассматривается децентрализованная обработка информации в двухканальных измерительных системах при косвенном измерении для различных алгоритмов фильтрации оценки вектора состояния в измерительных пунктах и пункте обработки информации. Проводится сравнительный анализ результатов имитационного моделирования синтезированных алгоритмов.
Задача обеспечения высокой точности оценивания координат и параметров траектории движения объекта может быть решена за счет применения многоканальных измерительных систем с оптимальной централизованной обработкой [1]. Алгоритмы оптимальной многоканальной калмановской фильтрации предполагают объединение в пункте обработки информации (ПОИ) оценок векторов наблюдаемых параметров, рассчитанных в измерительных пунктах.
Централизованной обработке присущи следующие недостатки [2]:
- необходимость передачи с измерительных пунктов в пункт обработки информации всего объема получаемой информации;
- большой объем вычислительных операций в пункте обработки информации для определения координат и параметров траектории цели;
- необходимость предварительного приведения оценок измерительных пунктов к единой системе координат и др.
Для устранения этих недостатков представляется целесообразным предусмотреть фильтрацию оценок текущих траекторных измерений в измерительных пунктах. Применение фильтрации координатной информации в измерительных пунктах позволит в случае неисправности системы обработки в одном или нескольких из этих пунктов производить оценивание координат объекта с относительно высокими точностными характеристиками.
В работе производится сравнительный анализ эффективности функционирования двух систем децентрализованной обработки информации. В измерительных пунктах применяются алгоритмы фильтрации Калмана.
В первом алгоритме обработки результаты оценок векторов состояния, рассчитанных в измерительных пунктах, комплексируются в ПОИ (Ф-О-алгоритм). Во втором алгоритме в ПОИ осуществляется их фильтрация (ф-ф-алгоритм).
Ввиду того что оценки векторов состояния первого (а1) и второго пунктов (а2) независимы, то для синтеза алгоритма, основанного на алгоритме комплексирования соответствующих оценок а1 и а2 в ПОИ, используем методику, приведенную в [3].
Оценки а1 и а2 на выходе измерительных пунктов имеют нормальный закон распределения [3]
п 1
р(а) = (2л) 2 |С,.|?ехр[- (а - а)ТС,. (а - а¡)/2 , (1) где а - оцениваемый случайный вектор состояния; С1, С2 -симметричные матрицы точности соответственно для первого (г= 1) и второго (г = 2) измерительных пунктов (г= 1, 2).
Найдем выражение логарифма совместной плотности вероятности:
lnа2) = —(1/2)(а - (Xj)rCj(a - ctj) -- (1/2)а — а2)тС2(а — а2) + const. (2)
Преобразуем (2) и, изменяя постоянную, получим
ln p((x1, а2) = -(1/2)ат (C1 + С2)а + (1/2)ат (С1(х1 +
+ С2а2) + (1/2)(C1(i1 + С2а2)а + const. (3)
Введем матрицу
Ср = С1 + С2 (4)
и вектор-столбец
С рар = С1а1 + С2а2 • (5) Как и матрицы С1 = С1т, С2 = С2Т, матрица (4) обладает свойствами симметрии: Ср = Срт. Транспонируя (5), получаем
а тСр = а1 тС1 + а2ТС2. (6)
В силу (4), (5), выражение (3) после изменения постоянной принимает вид
ln р(а1; а2) = -(1/2)(а - ар)тСр (а - ар) + const • (7)
Определяя отсюда плотность вероятности р(а1;а2) и постоянную (7) по условию нормировки, приходим к стандартной форме записи многомерного нормального закона:
n 1
р(а1, а2) = (2п) 2 |Ср| 2ехр^(1/2)(а - ар )ТСр (а - ар ) • (8)
Формально введенные величины Ср и ар приобретают отчетливый смысл результирующей матрицы точности и результирующей оценки вектора состояния на выходе ПОИ. Последняя, в силу (5), определяется выражением
а р = Ср-1 (С1а1 + С2 а2) • (9)
Оценка (Xр соответствует максимуму послеопытной плотности вероятности р(а1;а2). Она оптимальна по критерию минимума как при квадратичной, так и при простой функции стоимости.
Вычисление результирующей матрицы точности сводится, согласно (4), к суммированию матриц точности С1 и С2. Результирующая корреляционная матрица ошибок Ср1 получается путем обращения результирующей матрицы точности Ср .
Результирующая оценка вектора состояния ар сводится к весовой сумме оценок а1 и (i2.
Произведем синтез Ф-Ф-алгоритма, который основывается на фильтрации в ПОИ оценок a1k и a2k, полученных в результате вторичной обработки в измерительных пунктах (k - номер такта измерения).
Результирующую оценку вектора состояния а рк, рассчитанную в ПОИ, получим на основе алгоритмов фильтрации Калмана [3]:
-1 + Кк [ к [к (•* к / к-1)],
с-1 рк = с-1к / к-1 - кнс-
(10)
кНкС -к/к-1, (11)
Кк = С ’к/к-іНкТ (НкС Ік/к-іИк + С V) 1, (12)
с/ к-1 = вк-1с-‘к-1ВТк-1 + дк-1, (із)
где а к / к-1 = Вк-1а к-1 - прогнозированная оценка вектора
состояния; Кк - матричный коэффициент усиления; Xк = (¿Т1*, аТ2к )Т - оценка вектора наблюдаемых параметров на входе ПОИ; Нк (а к / к-1) - прогнозированная оценка вектора наблюдаемых параметров; Нк, Вк, Qk - известные матрицы.
Покажем схему алгоритма децентрализованной обработки информации в двухканальной измерительной системе (рис. 1).
Н =
с-к=
Кк = с-1
(I, I)
I
І-Ч у
с э1 (I, I)+
с1 к 0
0
измерений определяется как
2
Ск1 = сэк -X Ки1сэк = I -X к
Результирующая оценка вектора состояния а к с учетом (16) имеет вид
а,, = а„
+ (К11 к ,К12 к )
а„
а.
'I'
I
V /
а„
+ (К11 к, К12 к )
(19)
Рис. 1. Структурная схема децентрализованной системы обработки координатной информации
Элементы вектора состояния аік представляют собой оценки координат объекта и оценки скорости их изменения в декартовой системе координат. Фильтрация оценок вектора состояния в пункте обработки информации может также производиться в декартовой системе координат. В данном случае блочная матрица статического пересчета Н имеет следующую структуру:
(іМ (Т тГ )Т
= (, Ч , (14)
уа=(Хк/к-1
Эа ®
где I - единичная матрица.
Корреляционная матрица ошибок текущего оценивания с-ь имеет блочную структуру и в силу независимости и одновременности измерений в измерительных пунктах определяется выражением
(15)
1 '
с-к о
.0 с-к,
Матричный коэффициент усиления Кк после подстановки в формулу (12) соотношений (14), (15) преобразуется к виду
(16)
После проведенных математических выражение для матричного коэффициента усиления Кк принимает вид
Кк = (КП к, К12 к , (17)
где К11 к, К12 к - матричные коэффициенты усиления соответственно для первого и второго измерительных пунктов.
Результирующая корреляционная матрица ошибок
(18)
В более компактной форме выражение (19) можно представить в виде
ак = аэк + Х К1 1 к (1к - ^эк) . (20)
1=1
В экстраполяторе результирующая оценка вектора состояния, сформированная в ПОИ, прогнозируется на следующий такт измерения. Невязка V 1к =а1к -азк, рассчитанная в каждом из измерительных пунктов приема, с установленным матричным весом добавляется к прогнозированной оценке, что и дает результирующую оценку вектора состояния [3].
Алгоритм упрощенного фильтра Калмана (УФК), в оценке вектора состояния которого присутствуют оценка координат и скорости их изменении, называют УФК-алгоритмом, а алгоритм расширенного фильтра Калма-на (РФК), в оценке вектора состояния которого имеется и оценка ускорения, - РФК-алгоритмом.
Рассмотрим следующие модификации фильтров сопровождения в декартовой системе координат:
- фильтрация в вынесенном приемном пункте (ВПП), УФК-фильтр + фильтрация в ПОИ, УФК-фильтр (УФК-УФК-алгоритм);
- фильтрация в ВПП, УФК-фильтр + фильтрация в ПОИ, РФК-фильтр (УФК-РФК-алгоритм);
- фильтрация в ВПП, РФК-фильтр + фильтрация в ПОИ, УФК-фильтр (РФК-УФК-алгоритм);
- фильтрация в ВПП, РФК-фильтр + фильтрация в ПОИ, РФК-фильтр (РФК-РФК-алгоритм).
Результаты имитационного моделирования приведены в виде графической зависимости среднеквадратической ошибки Ох и ошибок оценивания А в метрах от номера такта измерения к для одной декартовой координаты X. Для двух других координат (У и Z) результаты аналогичны.
Для полного анализа достаточно исследовать качество фильтрации оценок координат объекта при двух траекториях его движения:
- равномерном прямолинейном движении объекта в направлении на радиолокационную службу для выявления потенциальных возможностей;
- движением по окружности с заданной перегрузкой для выявления возможностей сопровождения маневрирующих объектов.
Рассмотрим качество фильтрации каждого алгоритма на этапе равномерного прямолинейного движения. Представив зависимости среднеквадратичной ошибки (СКО) и ошибки оценивания координаты X в метрах от такта измерения для всех предлагаемых к исследованию алгоритмов сопровождения (рис. 2). УФК-УФК-алгоритму соответствует кривая 1, РФК-РФК-алгоритму - кривая 2, УФК-РФК-алгоритму - кривая 3, РФК-УФК-алгоритму - кривая 4. Анализ приведенных графиков показывает, что УФК-УФК-алгоритм при слежении за неманеврирующим объектом имеет преимущество в точности сопровождения, СКО
оценивания в 1,5 раза меньше СКО алгоритмов, основанных на расширенных фильтрах Калмана. РФК-РФК-алго-ритм имеет максимальные среднеквадратические ошибки. Это объясняется тем, что при сопровождении объекта, движущегося равномерно и прямолинейно в декартовой системе координат, учет второй производной (ускорения) в векторе состояния приводит к увеличению ошибок оценивания. УФК-РФК-алгоритм имеет результаты фильтрации оценок координат, сравнимые с УФК-УФК-алгорит-мом. Это указывает на то, что применение в приемных пунктах фильтров с упрощенным вектором состояния позволяет получать высокие точности оценивания вектора состояния при сопровождении объектов с прямолинейной равномерной траекторией движения. РФК-УФК-ал-горитм занимает промежуточное положение среди анализируемых алгоритмов.
координаты х для прямолинейного движения для
модификаций УФК и РФК алгоритмов сопровождения
Таким образом, среднеквадратическая ошибка для УФК-УФК-алгоритма к десятому шагу фильтрации уменьшается в 1,2 раза, к двадцатому шагу - в два раза. Для модификации РФК-РФК результаты следующие: к десятому шагу фильтрации СКО уменьшается в 1,2 раза, а к двадцатому - в 1,4 раза. Иными словами, до десятого шага фильтрации результаты практически одинаковые, а потом фильтры сопровождения с упрощенным вектором состояния имеют значительный выигрыш в точности оценивания вектора состояния.
Работа рассматриваемых алгоритмов сопровождения имеет ряд особенностей. При сопровождении объектов используется несколько фильтров Калмана (см. рис. 1). Каждый из фильтров можно настраивать, учитывая особенности внешних и внутренних факторов, влияющих на точность оценивания. Основным элементом управления в фильтре Калмана является матрица дискретного маневра Qk, физический смысл которой заключается в регулировании полосы пропускания фильтра. Необходимо выбрать такой нижний порог полосы пропускания фильтра, который позволит с минимальными ошибками оценивать координаты объекта и при этом не будет наблюдаться процесс расходимости фильтра. Если полоса пропус-
кания фильтра уменьшается, то фильтр при поступлении обновляющей информации начинает слабо реагировать на нее и в этом случае основной вес в результирующей оценке имеет прогнозированное значение оцениваемого параметра. Чрезмерное уменьшение полосы пропускания приводит к тому, что фильтр практически полностью доверяет прогнозированному значению, не принимая во внимание обновляющую информацию. В этом случае может наблюдаться процесс расходимости фильтра сопровождения, заключающийся в увеличении ошибок оценивания. В случае если ошибки оценивания значительно превышают ошибки первичных измерений, то речь идет о срыве сопровождения объекта.
Проведем анализ влияния полосы пропускания на качество фильтрации результирующей оценки вектора состояния на выходе ПОИ.
Данный анализ целесообразно провести при движении объекта по дуге окружности. Если полоса пропускания не адекватна протекающим процессам, то должно наблюдаться увеличение среднеквадратических ошибок оценивания координат объекта. Ширина полосы пропускания регулировалась путем изменения величины элементов матрицы дискретного маневра Qk на два порядка. Необходимо отметить, что при моделировании изменялась полоса пропускания только того фильтра сопровождения, который используется в пункте обработки, в измерительных же пунктах ширина полосы пропускания фильтров оставалась неизменной.
Зависимости СКО и ошибок оценивания координаты X объекта при движении по окружности с перегрузкой, равной двум, приведены на рис. 3. Кривая 1 соответствует УФК-УФК-алгоритму с узкой полосой пропускания, кривая 2 - УФК-УФК-алгоритму с широкой полосой пропускания, кривая 3 - РФК-РФК-алгоритму с узкой полосой пропускания, кривая 4 - РФК-РФК-алгоритму с широкой полосой пропускания. УФК-УФК-алгоритм с узкой полосой пропускания имеет большие ошибки оценивания в момент маневра объекта, причем эти ошибки почти в два раза выше ошибок измерения первичных координат. При увеличении полосы пропускания ошибки оценивания уменьшаются и срыва сопровождения объекта не наблю-
су .
Рис. 3. Зависимость СКО <Jx и ошибок оценивания А координаты x при движении по окружности с перегрузкой, равной двум
дается. В случае использования РФК-РФК-алгоритма резкого увеличения ошибок оценивания не бывает в обоих случаях. Следует отметить, что данный алгоритм имеет высокие точностные характеристики при изменении полосы пропускания фильтра в широком диапазоне.
Кривые 3 и 4 (см. рис. 3) показывают, что ширина полосы пропускания фильтра в пункте обработки информации практически не влияет на точностные характеристики при сопровождении равномерно прямолинейно движущегося объекта, в силу этого обстоятельства примем полосу пропускания фильтра Калмана в пункте обработки широкой для всех предлагаемых к исследованию алгоритмов.
Анализ точностных характеристик алгоритмов сопровождения при движении объекта по окружности производился с перегрузкой, равной шести (рис. 4). На нем представлены зависимости СКО и ошибок оценивания от номера шага фильтрации к для УФК-УФК-алгоритма (кривая 1), УФК-РФК-алгоритма (кривая 2), РФК-УФК-алгоритма (кривая 3) и РФК-РФК-алгоритма (кривая 4). РФК-УФК и РФК-РФК-алгоритмы имеют лучшие точностные характеристики на участке маневра, однако при прямолинейном равномерном движении их СКО намного больше, СКО УФК-УФК и УФК-РФК-алгоритмов. Ошибка оценивания координаты х РФК-РФК-алгоритма в момент маневра составляет всего 250...300 м, в то время как остальные алгоритмы сопровождения имеют ошибки оценивания координаты объекта порядка 800...1 000 м. После окончания маневра через 7...10 тактов измерения среднеквадратические ошибки приближаются к СКО в установившемся режиме.
Рис. 4. Зависимость СКО <Jx и ошибок оценивания А координаты x при движении по окружности с перегрузкой, равной шести
Рассмотрим результаты моделирования для случая, когда в измерительных пунктах осуществляется фильтрация оценок вектора состояния объекта РФК-фильтрами, а в пункте обработки проводится комплексирование полученной от измерительных пунктов информации. Параметры фильтров в измерительных пунктах идентичны. Представим зависимость СКО и ошибки оценивания координаты X объекта при прямолинейном равномерном движении от номера такта измерения (рис. 5), причем кривая 1
соответствует комплексированию информации в ПОИ, а кривая 2 - обработке оценок вектора состояния объекта на основе алгоритмов фильтрации Калмана в ПОИ.
координаты х при прямолинейном равномерном движении для Ф-О и Ф-Ф-алгоритмов (20)
Результаты фильтрации оценок координаты X практически одинаковы для обоих алгоритмов обработки. Среднеквадратическая ошибка к десятому шагу уменьшается в 1,22 раза, а к двадцатому - в 1,37 раза. Ошибки оценивания при фильтрации оценок вектора состояния объекта в пункте обработки информации изменяются более плавно, чем в случае их комплексирования.
Результаты фильтрации оценок координаты X при движении объекта по дуге окружности с перегрузкой, равной шести, представлены на рис. 6. Кривая 1 соответствует комплексированию информации, полученной от измерительных пунктов, кривая 2 - фильтрации в пункте обработки информации.
Рис. 6. Зависимость СКО Ох и ошибок оценивания Аx координаты x при движении по окружности для Ф-О и Ф-Ф-алгоритмов (20)
В заключение можно сделать следующие выводы:
- анализ результатов моделирования функционирования УФК-УФК, УФК-РФК, РФК-УФК, РФК-РФК-алго-ритмов показывает, что при прямолинейном равномер-
ном движении рекомендуется применять УФК-УФК-ал-горитм, при маневре - РФК-РФК-алгоритм;
- среднеквадратичные ошибки при фильтрации в ПОИ оценок вектора состояния маневрирующего объекта соизмеримы с соответствующими СКО измерений. Кроме того, ошибки оценивания при маневрировании объекта в результате фильтрации оценок в ПОИ в 1,2 раза больше, чем в случае комплексирования оценок, а Ф-О-алгоритм требует значительно меньшей производительности вычислительных средств.
Таким образом при децентрализованной вторичной обработке информации в двухканальной измерительной системе целесообразно применять алгоритм комплекси-рования в ПОИ оценок векторов состояния измеритель-
ных пунктов, рассчитанных с применением алгоритма фильтрации Калмана.
Библиографический список
1. Гришин, Б. П. Динамические системы, устойчивые к отказам / Б. П. Гришин, Ю. М. Казаринов. М.: Радио и связь, 1985. 176 с.
2. Петров, А. В. Анализ и синтез радиотехнических комплексов / А. В. Петров, А. А. Яковлев; под ред. В. Е. Думвича. М.: Радио и связь, 1984. 248 с.
3. Ширман, Я. Д. Теория и техника обработки радиолокационной информации на фоне помех / Я. Д. Ширман, В. Н. Манжос. М.: Радио и связь, 1981. 416 с.
N. P. Bogomolov
DECENTRALIZED ALGORITHMS OF INFORMATION PROCESSING IN TWO-CHANNEL MEASURE SYSTEMS
It is covered a decentralized algorithms of information processing in two-channel measure systems in case of an indirect measuring for different filtration algorithms of a condition vector estimation at the reception measure station and the station of information processing. Comparative analysis is carried out with a help of imitation modeling of synthesized algorithms.
