ISSN 0236-3933. Вестник МГТУ им. Н.Э. Баумана. Сер. «Приборостроение». 2012
196
2)
фильтрация в обратном времени — получение ПВ
1
,
N
k k k
p
Y λ α
и
1
;
N
k
k
p
Y α
3)
многоканальная интерполяция — получение ПВ
0
,
N
k
k
p
λ Y α
и
0
;
N
k
p
α Y
4)
итоговая оценка — получение ПВ
0
.
N
k
p
λ Y
Достоинством описанного выше подхода к синтезу предложенно-
го адаптивного алгоритма интерполяции на фиксированном интерва-
ле (применение в частности алгоритма двусторонней интерполяции)
является то, что для реализации текущей фильтрации в прямом вре-
мени и оценки ПВ
0
,
,
k
k
k
p
λ Y α
а также фильтрации в обратном
времени и оценки ПВ
1
,
N
k
k k
p
Y λ α
могут быть использованы из-
вестные и простые в вычислительном отношении алгоритмы квази-
оптимальной фильтрации на основе метода гауссовской аппроксима-
ции (в частности локальной гауссовской аппроксимации).
Предложенный алгоритм для задачи постобработки массива дис-
кретных выборок цифрового регистратора АПРС системы измерения
параметров траектории динамичных объектов с применением метода
широкополосной аналоговой ретрансляции может быть реализован
на базе алгоритмов комплексной обработки прямых и ретранслиро-
ванных сигналов ГНСС рассмотренных, например, в работах [6, 7,].
Результаты моделирования алгоритма.
Для проведения сравни-
тельного анализа предложенного адаптивного алгоритма интер-
поляции в среде MATLAB была разработана программно-матема-
тическая модель на базе метода статистических эквивалентов. При мо-
делировании использовалась траектория, характерная для реактивного
снаряда. Изменение во времени компонент вектора производной уско-
рения (рывок) и ускорения в геоцентрической связанной системе ко-
ординат для моделируемой траектории представлены на рис. 2.
Синтезированный алгоритм сравнивался с двумя известными ал-
горитмами:
1)
адаптивный алгоритм текущей фильтрации (одноэтапный ко-
герентный адаптивный алгоритм в форме многоканального измери-
теля [4, 8]);
2)
алгоритм интерполяции на фиксированном интервале
без
адаптации
(
одноэтапный когерентный алгоритм [4]).
Реализации во времени ошибок оценки координат и вектора ско-
рости объекта при наблюдении и обработке сигналов 8 НКА и ГФ = 2
для одноэтапного адаптивного алгоритма интерполяции приведены
на рис. 3. Для сравнения на рис. 4 показаны аналогичные ошибки для
одноэтапного алгоритма интерполяции
без адаптации
.