Конструкторский раздел
Структура системы контроля на основе оптимального оценивающего фильтра
При практической реализации диагностических моделей на основе наблюдателей параметров состояния необходимо учитывать следующие особенности реальных объектов контроля.
Нелинейная взаимосвязь входных и выходных сигналов объекта контроля;
Нелинейная взаимосвязь сигналов объекта контроля и наблюдений;
Влияние возмущений на входные и выходные сигналы объекта контроля, а также на наблюдателя.
Поэтому возникает необходимость построения стохастических диагностических моделей. Технология формирования таких моделей опирается на математическое описание функционирования эталонного и реального объектов контроля. Таким ОК ставятся в соответствие идеальный
и реальный
векторы параметров состояния, которые описываются следующими дифференциальными уравнениями:
для идеального ОК:
; (%1)
для идеального ОК:
(%2)
где
– вектор погрешностей действующий на ОК, характеризующийся ковариационной матрицей;
Параметры реального и идеального ОК связаны через уравнение ошибок. Для получения уравнения ошибок вычтем из соотношения (%2) выражение (%1) и разложим функцию
в ряд Тейлора относительно параметров идеального ОК, тогда
(%3)
Обозначим
– вектор ошибок ОК;
– матрица коэффициентов, характеризующих динамику изменения ошибок ОК (фундаментальная матрица динамической системы).
Тогда из соотношения (%3) можно получить следующее уравнение ошибок:
(%4)
Аналогичным образом можно определить взаимосвязь наблюдаемых параметров с ошибками ОК. Наблюдаемые параметры реального объекта контроля (ОКР) описываются уравнением
. (%5)
Модель внешнего по отношению к ОКР наблюдающего устройства имеет вид
, (%6)
где
– вектор возмущений в канале наблюдения.
– вектор наблюдений, формируемых эталонной по отношению к ОКР системой.
Вычитая из соотношения (%5) выражение (%6) и разложив функцию
в ряд Тейлора относительно параметров идеального ОК, получим
(%8)
Обозначим
– вектор наблюдений;
– матрица связи наблюдаемых параметров с вектором ошибок ОКР.
Тогда из уравнения (%8) можно получить модель наблюдения ошибок ОКР (модель сигналов)
(%10)
С учетом уравнений (%4) и (%10) сформируем алгоритм оптимального оценивания ошибок ОКР.
Алгоритм оптимального оценивания ошибок динамических объектов контроля в дискретном времени
На практике наблюдения
, отражающие состояние объекта контроля, формируются в дискретные моменты времени
и обозначается
. Для получения оптимальных оценок вектора ошибок ОКР по наблюдениям
применяют дискретный оптимальный фильтр Калмана. Для реализации такого фильтра необходимо определить динамику изменения оценок
и ковариационных матриц
за время
между моментами наблюдений
и
. Переходные матрицы
и
для прогнозирования указанных параметров находятся из решения дифференциального уравнения (%4). Для однородной части этого уравнения при неизменной матрице
на интервале прогноза решение определяется следующими отношениями:
(%32)
(%33)
, (%34)
где 
– переходная матрица для вектора ошибок, удовлетворяющая следующему дифференциальному уравнению.
(%35)
При
. (%36)
Для определения переходной матрицы
для вектора возмущений
покажем сначала, что решение неоднородного уравнения (%4) имеет вид
. (%37)
Дифференцируя соотношение (%37) по аргументу
и применяя правило Лейбница к интегралу свертки, а именно:
(%38)
Получим
. (%39)
С учетом равенства (%35) и условия (%36) соотношение (%39) можно привести к следующему виду:
(%40)
Можно видеть, что подстановка решения (%37) в соотношение (%40) приводит к уравнению, эквивалентному (4). В выражении (%37) можно ввести следующее обозначение:
(%42)
Тогда, дифференцируя интеграл свертки (%42) по правилу (%38) и учитывая соотношения (%35), (%36), получим дифференциальное уравнение для переходной матрицы вектора возмущений ОКР.
при
. (%43)
Таким образом, уравнение ошибок ОК (4) будет иметь вид
, (%44)
С учетом решения (%37) значения оценок
вектора ошибок ОК
и его ковариационной матрицы
на интервале прогноза
будут определяться следующими соотношениями:
(%45)

. (%46)
Соотношение (%46) сформировано с учетом некоррелированности начальных ошибок оценивания и шумов ОК, т.е.
для
.
Для того, чтобы гарантировать попадание точки, в которой наблюдается всплеск дельта-функции, в область интегрирования, необходимо вначале интегрировать по переменной, имеющей больший диапазон. В частности, при
необходимо выражение (%46) сначала проинтегрировать по
. В этом случае соотношение (%46) с учетом свойств дельта-функции (%27) можно представить в следующем виде:
(%47)
Свойства оптимального фильтра Калмана позволяют реализовать на его основе диагностические модели как сигнального, так и параметрического типа. Ядром таких моделей являются уравнения наблюдений и ошибок объекта контроля. При этом должны учитываться следующие особенности включения ОФК в структуру системы контроля:
- необходимо одновременно решать нелинейные динамические уравнения ОК и линеаризованные уравнения его ошибок;
- ОФК должен функционировать в условиях априорной неопределённости относительно структуры и параметров диагностических моделей ОК, а также статистических характеристик возмущений и шумов.
В настоящее время ОФК составляет ядро программно-математического обеспечения современных навигационных комплексов. Это создает необходимую основу для расширения функциональных возможностей ОФК на решение задач контроля и диагностики.