Алгоритм оптимального оценивания ошибок динамических объектов контроля в дискретном времени




Конструкторский раздел

Структура системы контроля на основе оптимального оценивающего фильтра

При практической реализации диагностических моделей на основе наблюдателей параметров состояния необходимо учитывать следующие особенности реальных объектов контроля.

Нелинейная взаимосвязь входных и выходных сигналов объекта контроля;

Нелинейная взаимосвязь сигналов объекта контроля и наблюдений;

Влияние возмущений на входные и выходные сигналы объекта контроля, а также на наблюдателя.

Поэтому возникает необходимость построения стохастических диагностических моделей. Технология формирования таких моделей опирается на математическое описание функционирования эталонного и реального объектов контроля. Таким ОК ставятся в соответствие идеальный и реальный векторы параметров состояния, которые описываются следующими дифференциальными уравнениями:

для идеального ОК: ; (%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)

 

 

Свойства оптимального фильтра Калмана позволяют реализовать на его основе диагностические модели как сигнального, так и параметрического типа. Ядром таких моделей являются уравнения наблюдений и ошибок объекта контроля. При этом должны учитываться следующие особенности включения ОФК в структуру системы контроля:

- необходимо одновременно решать нелинейные динамические уравнения ОК и линеаризованные уравнения его ошибок;

- ОФК должен функционировать в условиях априорной неопределённости относительно структуры и параметров диагностических моделей ОК, а также статистических характеристик возмущений и шумов.

В настоящее время ОФК составляет ядро программно-математического обеспечения современных навигационных комплексов. Это создает необходимую основу для расширения функциональных возможностей ОФК на решение задач контроля и диагностики.



Поделиться:




Поиск по сайту

©2015-2024 poisk-ru.ru
Все права принадлежать их авторам. Данный сайт не претендует на авторства, а предоставляет бесплатное использование.
Дата создания страницы: 2019-07-23 Нарушение авторских прав и Нарушение персональных данных


Поиск по сайту: