CC BY
67
ВЫЧИСЛИТЕЛЬНАЯ ТЕХНИКА
УДК 62-50
DOI: 10.17586/0021-3454-2017-60-4-311-317
СВОЙСТВА ФИЛЬТРА КАЛМАНА В ПЕРЕХОДНОМ РЕЖИМЕ
А. Ю. Каплин1, М. Г. Степанов1, А. Г. Ярмолич12
1ОАО „Радиоавионика", 190005, Санкт-Петербург, Россия E-mail: [email protected] 2Балтийский государственный технический университет „ВОЕНМЕХ" им. Д. Ф. Устинова,
198005, Санкт-Петербург, Россия
Рассматриваются условия достижения высокой точности фильтра Калмана в переходном режиме. Определяются соотношения между соответствующими значениями параметров уравнений состояния и измерения. Сформулировано несложное правило предварительного определения характера переходного процесса. На примере практически важной задачи калмановской оценки экспоненциально-коррелированного случайного процесса расчетным путем иллюстрируется применение предложенного правила при различных исходных данных.
Ключевые слова: переходный процесс, фильтр Калмана, установившийся режим фильтрации, точность оценивания
Фильтр Калмана традиционно используется в радиолокационных и навигационных комплексах, системах управления и связи, центрах обработки траекторных измерений [1—6] и др. Не ослабевает интерес к калмановской фильтрации и в настоящее время [7—10].
При использовании калмановских алгоритмов в конкретных комплексах и системах разработчики справедливо считают, что в рамках принятых моделей (уравнений) состояния и измерения фильтр Калмана обеспечивает минимальную (в среднеквадратическом смысле) ошибку оценивания по сравнению с любой другой линейной системой фильтрации вне зависимости от соотношения значений параметров этих моделей. Однако при этом разработчики зачастую пренебрегают характером переходного процесса или предполагают, что при оценке стационарных случайных процессов, описываемых уравнениями состояния и измерения с постоянными параметрами, зависимость ошибки оценивания каждого параметра от объема измерительной выборки представляет собой монотонно убывающую функцию, асимптотически приближающуюся к установившемуся значению. В результате не учитывается соотношение исходного уровня точности априорных данных (начальных условий) об оцениваемом процессе с уровнем точности его апостериорной оценки в установившемся режиме фильтрации. Примером такого подхода являются работы [7, 9].
В то же время существует ряд практических задач, где переходный режим фильтра Калмана не менее важен, чем установившийся. Продолжительность и точность переходного процесса, установившийся уровень точности фильтрации являются важными факторами в таких случаях, как наведение боевых элементов кассетных боеприпасов с активным ударным ядром (время стабилизации и прицеливания, как правило, меньше времени переходного режима), автоматическое сопровождение беспилотных летательных аппаратов квадрокоптер-ного типа или комплексирование систем спутниковой и автономной навигации в информационной системе бойца (спонтанное резкое изменение параметров уравнения состояния
наблюдаемого объекта) и т.п. Предельным случаем является оценка нестационарных процессов, задаваемых уравнениями состояния и измерения с переменными параметрами. В этом случае фильтр Калмана также является линейной динамической системой с переменными параметрами, перманентно находящейся в режиме переходного процесса, характер которого определяет весь процесс фильтрации, и не учитывать который уже нельзя. На практике такая задача может иметь место при автоматическом сопровождении статически неустойчивых аэродинамических летательных аппаратов.
В разное время предпринимались отдельные попытки скорректировать переходный процесс фильтра Калмана, предлагались меры по его ускорению и повышению точности (см., например, [1, 11, 12]). Однако полученные результаты носили частный характер. Проанализируем возможные особенности фильтра Калмана и определим условия обеспечения высокой точности его переходного режима.
Цели работы — анализ условий, благоприятствующих высокой точности фильтра Калмана в переходном режиме; определение соотношений между соответствующими значениями параметров уравнений состояния и измерения; формулировка несложного правила предварительного определения характера переходного процесса.
Для анализа выберем базовый вариант фильтра [13, 14], являющийся основой большинства известных теоретических обобщений и практических модификаций процедуры калма-новской фильтрации. Подлежащий оценке стационарный векторный случайный процесс хд формируется в соответствии с уравнением состояния
хд+1 = Фхд + ™ д , д=1, 2 (1)
где Ф — постоянная переходная матрица; w д — векторный дискретный нормально распределенный (гауссов) белый шум возмущений с нулевым математическим ожиданием и ковариационной матрицей 0, w д -N(0,0). Математическое ожидание и ковариационная матрица заданного начального условия х0 соответственно равны х и Р0 .
Уравнение измерения имеет вид
Уд = Нхд + Уд , (2)
где Н — матрица измерения; Уд ~М(0,К) — векторный шум измерений.
В этих условиях оптимальной по критерию минимума среднеквадратической ошибки несмещенного оценивания процесса хд по измерениям Уд (минимума ковариационной матрицы ошибок Рд = М (хд -хд)(хд -хд)Т , где хд — оценка вектора хд, М[*] — математическое ожидание) является линейная динамическая система, задаваемая соотношениями при начальных условиях х и Ро :
хд = Фхд-1 + РдНТЯ-1 (уд - НФхд-1), (3)
Рд / д -1 = ФРд-1ФТ + 0,
1 НРд/д -1. (4)
Рд/1д-1 + нТ Я 1н
-1 _ Т = Рд/д-1 - Рд/д -1Н
НРд/д-1нТ + Я
Анализ показывает, что при определенных соотношениях между параметрами 0, Я и Ро в уравнениях (1), (2) переходный процесс и установившийся уровень точности фильтра
Калмана могут различаться.
Обратимся к следующему практически важному, наглядному и простому случаю. Оценке подлежит скалярный случайный процесс, описываемый уравнением состояния [3]
хд +1 = Фхд + wk , (5)
где Ф = ехр(-а ХТ) — постоянный коэффициент; — дискретный белый гауссов шум с ну-
2 2
левым математическим ожиданием и дисперсией Q = оХ(1 -Ф ), и О,0); начальное условие х0~Л"(Х,Ро); Т = гк +1 - гк — интервал дискретизации (шаг) по времени.
Рекуррентное соотношение (5) является для любых Т строгим дискретным аналогом стохастического дифференциального уравнения
дх (г) / дг = -аХх (г ) + ^/ 2ах а2Хп (г) (6)
(п (г) — стандартный белый гауссов шум с нулевым математическим ожиданием и единичной интенсивностью), описывающего непрерывный экспоненциально-коррелированный скалярный гауссов процесс х (г) с шириной спектра ах (обратно пропорциональна постоянной
времени корреляции Тх ) и стационарным значением дисперсии а1Х [3].
Модель состояния (6) широко используется для описания статистической динамики случайной составляющей параметров движения маневрирующих объектов [2], при калманов-ской оценке корректируемых погрешностей датчиков инерциальных навигационных систем [7] и в других задачах.
Измерения проводятся в соответствии с уравнением
Ук = хк + ук,
где Ук — скалярный шум измерений, Ук -N(0, оУ = Я ).
В этом случае уравнения фильтра Калмана (3), (4) принимают следующий вид при наЛ г>
чальных условиях Хо и Ро :
хк = ФХк-1 + Я (Ук- ФХк-1), (7)
Рк / к-1 = ф 2 Рк-1 + 0,
РыЬ-лЯ
( Ф2 Рк-1 + 0 ) Я
Р = -к / к= -(8)
Рк / к-1 + Я Ф2 Рк-1 + 0 + Я
Здесь Рк =ок — текущее значение дисперсии ошибки оценки Хк, Рк/к-1 — экстраполиро-
2
ванное на шаг значение этой дисперсии, Ро = Оо — начальное значение дисперсии Рк, для Хо = Х Ро =о2и =0. Возможны другие варианты задания начального условия Хо и соответствующего значения дисперсии Ро.
Очевидно, что если с увеличением номера шага фильтрации к дисперсия Рк становится меньше начального значения Ро и равномерно убывает (нисходящий переходный процесс), стремясь к установившемуся значению Рю (Рю < Ро), фильтр Калмана обеспечивает повышение точности оценивания в обоих — переходном и установившемся — режимах. При таком характере переходного процесса применение фильтра целесообразно. Если переходный процесс сопровождается ростом Рк и приводит к Рю > Ро (возрастающий характер переходного
процесса), вопрос о целесообразности калмановской фильтрации должен решаться в каждом конкретном случае исходя из целей и условий решаемой задачи.
Учитывая монотонный характер изменения Рк с ростом к, для определения характера переходного процесса калмановского фильтра достаточно ограничиться сравнением величин Р1 и Ро на первом шаге и проверить выполнение неравенства Р1 < Ро, которому в соответствии с соотношением (8) при к=1 эквивалентно условие
0(Я - Р0) < Ф2Р2 + роЯ(1" Ф2). (9)
Выполнение неравенства (9) свидетельствует о нисходящем характере переходного процесса фильтра, невыполнение — о возрастающем.
Поскольку в (9) при всех возможных значениях параметров ах и Т (см. уравнение (5))
0<Ф<1, правая часть в неравенстве (9) всегда положительна. Возможны два варианта выполнения (невыполнения) условия (9).
Вариант 1: Я< Ро. Неравенство (9) выполняется при любых значениях дисперсии шума состояния Q и любом способе задания начального условия Ро (начального условия хо в уравнении фильтра (7)), в том числе при Ро = Я. Последнее означает, что в качестве начального условия Хо выбирается первое измерение >1 с отсчетом последующих, начиная с >2 .
Вариант 2: Я > Ро. Из (9) видно, что при различных значениях Я, Q, Ро и Ф возможен как нисходящий (неравенство (9) выполняется), так и возрастающий характер переходного процесса (неравенство (9) нарушается). В последнем случае (9) приобретает вид
Q(Я - Ро) > Ф Ро + РоЯ(1 - Ф ), откуда легко получить пороговое значение Я = Япор :
Я = Ро(ф2 Ро + Q) (1о)
Япор=фо^рОТО-РО ■ (1о)
являющееся границей между областью повышения и областью понижения точности фильтрации в переходном режиме. Если при заданных условиями задачи значениях 0, Ро и Ф (заданном Япор ) измерения ук проводятся с ошибкой Ук , имеющей дисперсию Я < Япор, переходный процесс имеет нисходящий характер; при Я > Япор — восходящий. Заметим, что в отмеченном ранее практически важном случае Ро = 0 пороговое значение Япор = 0(Ф2 +1)/ Ф2.
Возможен иной способ предварительного определения характера переходного процесса фильтра, заключающийся в сравнении начального значения дисперсии Ро с установившимся значением дисперсии ошибки оценки р при к ^го . При Ро>— нисходящий характер переходного процесса, при Ро< — возрастающий.
Значение можно получить, подставив в уравнение (8) единственное неизвестное
2 2
Рк-1 = Рк = Р. В результате приходим к квадратному уравнению Ф + -0Я = 0, положительный корень которого есть искомое значение :
0+(1-Ф2)
Я
-0 - (1 - ф2 ) Я + ¡0+(1 - ф2 ) я ]2 + 4 Ф20Я
Р«=--^Г2-"-. (11)
2Ф 2
На практике предпочтителен более простой способ проверки выполнения неравенства (9) путем сравнения дисперсии ошибки измерений Я с Япор (10). В то же время более сложная процедура расчета дисперсии (11) позволяет установить конечную достижимую точность фильтрации.
Приведенные соотношения были использованы при расчете дисперсии ошибки оценивания Рк как функции номера шага фильтрации к для двух примеров. Первый (рис. 1) иллюстрирует переходный процесс калмановской оценки случайной составляющей некоторого
обобщенного параметра движения маневрирующего объекта. Типичными для подобных задач являются значения ширины спектра случайной составляющей а х=1—3 Гц и шага измерений 7=0,1—1 с. В примере были выбраны а х=1 Гц, 7=0,1 с, 0=0,362. Расчеты проводились
при фиксированной дисперсии случайной составляющей ох = 2 для двух значений дисперсии ошибки измерений Я=4 (а) и 6 (б). Использовались три варианта задания начального значения дисперсии Ро : как одного из двух значений дисперсии Я (1), равного дисперсии шума
состояния 0 (напомним, что 0 = ъ\ [1 _ехР (~2ахТ)]) (2) и равного 0,1 (3). Последний
вариант соответствует случаю высокоточных априорных данных о начальных условиях фильтрации.
а)
Рк
3 2
б)
( • \ * \ V * 1 Рк 5 4
1 3
■ 2
2 - * ■ ШЬ 1
3
О
Г \ . » : \ * \ * 1
2 ... ... . -------- ГА&КЛ
>«* V» -.--- ------—1
к
Рис. 1
Второй пример (рис. 2) соответствует задаче калмановской оценки погрешностей датчиков инерциальных навигационных систем. Типичными здесь являются значения а х =5—10 Гц,
7=0,01—0,05 с. В расчетах использовались величины а х =5 Гц, 7=0,05 с, 0=0,786. Остальные
данные аналогичны первому примеру.
а)
Рк \
3 2
\ » \ г \ » V 1 1
51
2 "г- ■ ■ ■ ■'VI Им
^ 3
* ... .... .....
б)
* * \ » : \ * \ 1
* N.
2 —
для
Рис. 2
Кривые на рис. 1, 2 подтверждают сделанный ранее вывод о зависимости характера переходного процесса фильтра Калмана (нисходящий или возрастающий) от сочетания исходных данных (параметров уравнений состояния и измерения), задаваемых условиями решаемой задачи. Продолжительность переходного процесса при этом может быть значительной.
В заключение важно отметить, что в настоящей статье речь идет о потенциальных точностных свойствах фильтра Калмана. Вопросы точности и сходимости (расходимости) фильтра, связанные с несоответствием реальной динамики случайных процессов принятым моделям (1), (2), погрешностями линеаризации при решении задач нелинейной фильтрации, вычислительными аспектами и другими факторами, здесь не рассматриваются.
1
0
1
2
3
4
5
к
0
1
2
3
4
5
1
к
0
1
2
3
4
5
к
0
1
2
3
4
5
список литературы
1. Кузьмин С. З. Основы теории цифровой обработки радиолокационной информации. М.: Сов. радио, 1974. 432 с.
2. Шебшаевич В. С., Дмитриев П. П., Иванцевич Н. В. и др. Сетевые спутниковые радионавигационные системы. М.: Радио и связь, 1993. 408 с.
3. ЯрлыковМ. С. Статистическая теория радионавигации. М.: Радио и связь, 1985. 344 с.
4. Кузовков Н. Т., Салычев О. С. Инерциальная навигация и оптимальная фильтрация. М.: Машиностроение, 1982. 216 с.
5. Сейдж Э., Мелс Дж. Теория оценивания и ее применение в связи и управлении. М.: Связь, 1976. 496 с.
6. ЖданюкБ. Ф. Основы статистической обработки траекторных измерений. М.: Сов. радио, 1978. 384 с.
7. Иванов А. В. Навигационные системы подвижных наземных объектов. Алгоритмы обработки информации в угломерном канале // Радиотехника. 2013. № 4. С. 23—29.
8. Каплин А. Ю., Степанов М. Г. Использование автономной навигационной системы высокоточного позиционирования пешехода на местности // Информационно-управляющие системы. 2015. № 6. С. 86—92.
9. Каплин А. Ю., Степанов М. Г. Модель и алгоритм комплексной обработки информации азимутального канала пешеходной навигационной системы // Изв. вузов. Приборостроение. 2016. Т. 59, № 3. С. 181—188.
10. Советов В. М. Дискретный фильтр Калмана в задачах оптимального приема сигналов // Радиотехника. 2016. № 3. С. 43—49.
11. Степанов М. Г. Использование редуцированных оценок для повышения устойчивости фильтра Калмана в переходном режиме // Изв. вузов МВ и ССО СССР. Приборостроение. 1984. Т. 27, № 1.
12. Степанов М. Г. Редуцированное оценивание вектора состояния линейной динамической системы // Изв. АН СССР. Техническая кибернетика. 1984. № 4.
13. Ли Р. Оптимальные оценки, определение характеристик и управление. М.: Наука, 1966. 176 с.
14. Браммер К., Зиффлинг Г. Фильтр Калмана—Бьюси. М.: Наука. Гл. ред. физ.-мат. лит., 1982. 200 с.
Сведения об авторах
Александр Юрьевич Каплин — канд. техн. наук; ОАО „Радиоавионика";
E-mail: [email protected]
Михаил Георгиевич Степанов — д-р техн. наук, профессор; ОАО „Радиоавионика";
E-mail: [email protected]
Алексей Григорьевич Ярмолич — ОАО „Радиоавионика"; ведущий инженер; БГТУ „ВОЕНМЕХ" им. Д. Ф. Устинова; аспирант; E-mail: [email protected]
Рекомендована Поступила в редакцию
БГТУ „ВОЕНМЕХ" им. Д. Ф. Устинова 01.09.16 г.
Ссылка для цитирования: Каплин А. Ю., Степанов М. Г., Ярмолич А. Г. Свойства фильтра Калмана в переходном режиме // Изв. вузов. Приборостроение. 2017. Т. 60, № 4. С. 311—317.
THE KALMAN FILTER PROPERTIES IN TRANSITION REGIME
A. Yu. Kaplin1, M. G. Stepanov1, A. G. Yarmolich1,2
1 Radioavionica JSC, 190005, St. Petersburg, Russia E-mail: [email protected] 2D. F. Ustinov Baltic State Technical University (Voenmech), 198005, St. Petersburg, Russia
The conditions providing high accuracy of the Kalman filter in transition regime are considered. Relations between the respective values of parameters of state equations and measurement are determined. A simple rule for preliminary evaluation of the transition process nature is formulated. By the example of practically important problem of Kalman estimation of exponentially-correlated random process, application of the proposed rules is illustrated by presented calculations with different initial data.
Keywords: transition regime, Kalman filter, steady-state filter mode, estimation accuracy
Data on authors
Alexander Yu. Kaplin — PhD; Radioavionica JSC; E-mail: [email protected] Mikhail G. Stepanov — Pr. Sci., Professor; Radioavionica JSC; E-mail: [email protected] Alexey G. Yarmolich — Radioavionica JSC; Leading Engineer; Post-Graduate Student; D. F. Ustinov Baltic State Technical University (Voenmech), E-mail: [email protected]
For citation: Kaplin A. Yu., Stepanov M. G., Yarmolich A. G. The Kalman filter properties in transition regime // Journal of Instrument Engineering. 2017. Vol. 60, N 4. P. 311—317 (in Russian).
DOI: 10.17586/0021-3454-2017-60-4-311-317
