Однако здесь возникает ряд проблем вычислительного порядка, связанных с увеличением размерности задачи оценивания, что ведет к значительному увеличению объема вычислений, и как следствие - к неустойчивости машинных алгоритмов.
В таких условиях чрезвычайно важное, если не решающее, значение приобретает исследование возможностей сведения исходной задачи совместной идентификации, имеющей большую размерность, к задачам меньшей размерности, т.е. возможностей декомпозиции задачи оценивания расширенного вектора состояния (сведение ее к задачам меньшей размерности).
Поэтому в данной работе в качестве одного из путей построения алгоритма оценивания в СМСН ДПЛА предлагается сначала провести идентификацию параметров
, , и оценить параметры , , , a затем, уже с учетом этой информации, осуществить идентификацию параметров , , . В свете этих предложений идентификацию параметров целесообразно осуществлять в два этапа. Таким образом, идея этого метода интересна, так как позволяет существенно упростить алгоритмы обработки информации СМСН ДПЛА. Следуя намеченной схеме, решение задачи нелинейной фильтрации при поиске оценок параметров относительного движения ЛА ищется с помощью модифицированного нелинейного фильтра.При решении практических задач идентификации возникает проблема обеспечения устойчивости процессов обработки информации. В некоторых случаях реальные ошибки оценивания могут во много раз превосходить расчетные значения, характеризуемые матрицей ковариации ошибок оценивания. Наиболее часто встречающимся источником расходимости являются неточности моделей состояния и наблюдения, неточность задания характеристик шумов и т.д. Так как синтез алгоритмов обычно проводится в условиях априорной неопределенности, то для ее преодоления используют различные адаптивные алгоритмы.
В данной работе с целью устранения в определенных пределах незнания о характеристиках шумов измерений используется идея J-адаптивного фильтра. Такое название он получил по первой букве фамилии автора “J” предложивший алгоритм J-адаптивного фильтра для дискретных систем в работе [7].
В работе [7] J-адаптивный фильтр использован для оценивания параметров движения баллистического объекта: положения и скорости, а также для оценивания вектора аэродинамических параметров (коэффициентов лобового сопротивления и подъемной силы) или составляющих вектора немоделируемых ускорений по осям декартовой СК. Как указано в этой работе, J-адаптивный фильтр идентифицирует 99% немоделируемых ускорений.
Пусть динамика системы описывается следующим нелинейным уравнением:
(18) |
В дискретной форме это уравнение (с точностью до членов второго порядка) будет иметь вид [7]:
, | (19) |
где
вектор состояния размерности , непрерывная векторная функция, которая по крайней мере дважды дифференцируема по своим аргументам,, | (20) |
матрица Якоби вектор-функции
, a отражает погрешность дискретизации и неопределенности моделирования в уравнении динамики. Предполагается, что является белой гауссовской последовательностью с нулевым средним и матрицей ковариации :, | (21) |
где
дельта-функция Кронекера.Начальное значение вектора состояния
принимается гауссовской случайной переменной с известными средним значением и матрицей ковариации вида:(22) |
Процесс измерения описывается уравнением:
, | (23) |
, |
где
вектор погрешностей измерений типа белой последовательности с нулевым средним и матрицей ковариации , вектор размерности , непрерывная векторная функция, дважды дифференцируемая по своим аргументам.Тогда обобщенный калмановский фильтр описывается следующими уравнениями:
, | (24) |
, | (25) |
где
весовая матрица (матрица коэффициентов усиления фильтра).Уравнения экстраполяции матрицы ковариации:
Уравнения коррекции матрицы ковариации:
, | (27) |
. |
Далее на основе [8] рассмотрим метод получения J-адаптивного нелинейного фильтра.
Предполагается, что в уравнения состояния и измерений входят одни и те же определенные параметры, величины которых известны неточно. Такие параметры можно представить как случайные переменные с заранее известной статистикой
, При этом предполагается, что , , , некоррелированы, и матрица ковариации остается неизменной на всем временном интервале процесса оценивания.При синтезе фильтра в вектор состояния
включается вектор . Вводя для вектора параметров уравнения вида, | (28) |
можно получить расширенную систему уравнений с вектором состояния:
, | (29) |
и уравнениями измерений: