УСТРОЙСТВО ОПРЕДЕЛЕНИЯ ТЕХНИЧЕСКОГО СОСТОЯНИЯ БОРТОВЫХ СИСТЕМ ЛЕТАТЕЛЬНЫХ АППАРАТОВ ПО РЕЗУЛЬТАТАМ ОЦЕНКИ ПАРАМЕТРОВ Российский патент 2015 года по МПК G05B23/00 

Изобретение относится к области цифровой вычислительной техники и может быть использовано в автоматических и автоматизированных системах различного назначения для определения технического состояния по результатам идентификации параметров бортовых систем летательного аппарата.

Известно устройство оценки диагностических признаков в основе которого лежит метод Калмановской фильтрации [1]. Его эффективное функционирование возможно в случае полной априорной информации и законов распределения внешних воздействий, что является существенным недостатком.

Наиболее близким по технической сущности к заявленному изобретению является устройство идентификации параметров динамических систем на основе вариационных принципов [2]. Однако его использование не позволяет получить требуемую точность при определении технического состояния бортовых систем летательных аппаратов.

Цель изобретения - повысить эффективность определения технического состояния по результатам идентификации параметров бортовых систем летательного аппарата. Указанный технический результат достигается за счет устройства определения технического состояния типовых бортовых систем летательных аппаратов по результатам оценки параметров, которое содержит следующие блоки: блок хранения констант; первый, второй, третий, четвертый, пятый, шестой, седьмой, восьмой, девятый, десятый блоки формирования произведения; блок формирования суммы; блок транспонирования; первый, второй блоки формирования разности; первый, второй, третий, четвертый блоки формирования функций $$\frac{\partial H^T}{\partial \widehat{x}}$$ , $$f(\widehat{x}(\widehat{m}),\widehat{m},t)$$ , $$H(\widehat{x}(\widehat{m}),\widehat{m},t)$$ , $$\frac{\partial f}{\partial \widehat{m}}$$ соответственно; первый, второй, третий, четвертый блоки интегрирования; при этом на информационный вход блока хранения констант, который является входом устройства, поступает значение наблюдаемой величины; первый информационный выход блока хранения констант соединен с первым информационным входом блока формирования суммы; второй информационный выход блока хранения констант соединен с первым информационным входом первого блока формирования произведения и информационным входом блока транспонирования; третий информационный выход блока хранения констант соединен со вторыми информационными входами первого блока формирования произведения и первым информационным входом четвертого блока формирования произведения; четвертый информационный выход блока хранения констант соединен с первым информационным входом третьего блока формирования произведения и вторыми информационными входами шестого и девятого блоков формирования произведения; пятый информационный выход блока хранения констант соединен с первым информационным входом восьмого блока формирования произведения; шестой информационный выход блока хранения констант соединен с первым информационным входом второго блока формирования разности; седьмой информационный выход блока хранения констант соединен с первыми информационными входами первого, второго, третьего, четвертого блоков формирования функций $$\frac{\partial H^T}{\partial \widehat{x}}$$ , $$f(\widehat{x}(\widehat{m}),\widehat{m},t)$$ , $$H(\widehat{x}(\widehat{m}),\widehat{m},t)$$ , $$\frac{\partial f}{\partial \widehat{m}}$$ соответственно; восьмой информационный выход блока хранения констант соединен со вторыми информационными входами первого, второго, третьего, четвертого блоков формирования функций $$\frac{\partial H^T}{\partial \widehat{x}}$$ , $$f(\widehat{x}(\widehat{m}),\widehat{m},t)$$ , $$H(\widehat{x}(\widehat{m}),\widehat{m},t)$$ , $$\frac{\partial f}{\partial \widehat{m}}$$ соответственно; девятый информационный выход блока хранения констант соединен с третьими информационными входами первого, второго, третьего, четвертого блоков формирования функций $$\frac{\partial H^T}{\partial \widehat{x}}$$ , $$f(\widehat{x}(\widehat{m}),\widehat{m},t)$$ , $$H(\widehat{x}(\widehat{m}),\widehat{m},t)$$ , $$\frac{\partial f}{\partial \widehat{m}}$$ соответственно; информационный выход первого блока формирования произведения соединен с первым информационным входом второго блока формирования произведения; информационный выход блока транспонирования соединен со вторым информационным входом четвертого блока формирования произведения; информационный выход четвертого блока формирования произведения соединен со вторым информационным входом третьего блока формирования произведения и первым информационным входом пятого блока формирования произведения; информационный выход третьего блока формирования произведения соединен со вторым информационным входом второго блока формирования произведения; информационный выход второго блока формирования произведения соединен со вторым информационным входом блока формирования суммы, информационный выход которого соединен с информационным входом первого блока интегрирования; информационный выход первого блока формирования функции $$\frac{\partial H^T}{\partial \widehat{x}}$$ соединен со вторыми информационными входами пятого и восьмого блоков формирования произведения; информационный выход восьмого блока формирования произведения соединен с первым информационным входом девятого блока формирования произведения, информационный выход которого соединен с первым информационным входом десятого блока формирования произведения; информационный выход второго блока формирования функции $$f(\widehat{x}(\widehat{m}),\widehat{m},t)$$ соединен с первым информационным входом первого блока формирования разности; информационный выход третьего блока формирования функции $$H(\widehat{x}(\widehat{m}),\widehat{m},t)$$ соединен со вторым информационным входом второго блока формирования разности, информационный выход которого соединен со вторыми информационными входами седьмого и десятого блоков формирования произведения; информационный выход седьмого блока формирования произведения соединен с информационным входом второго блока интегрирования; информационный выход десятого блока формирования произведения соединен со вторым информационным входом первого блока формирования разности, информационный выход которого соединен с информационным входом третьего блока интегрирования; информационный выход четвертого блока формирования функции $$\frac{\partial f}{\partial \widehat{m}}$$ соединен с информационным входом четвертого блока интегрирования; информационный выход первого блока интегрирования соединен со вторым информационным входом первого блока формирования произведения и первым информационным входом четвертого блока формирования произведения; информационный выход второго блока интегрирования является выходом устройства, с него же снимается искомое значение неизвестного приведенного управляющего момента и соединен со вторыми информационными входами первого, второго, третьего, четвертого блоков формирования функций $$\frac{\partial H^T}{\partial \widehat{x}}$$ , $$f(\widehat{x}(\widehat{m}),\widehat{m},t)$$ , $$H(\widehat{x}(\widehat{m}),\widehat{m},t)$$ , $$\frac{\partial f}{\partial \widehat{m}}$$ соответственно; информационный выход третьего блока интегрирования является выходом устройства, с него же снимается значение вектора состояния и соединен с третьими информационными входами первого, второго, третьего, четвертого блоков формирования функций $$\frac{\partial H^T}{\partial \widehat{x}}$$ , $$f(\widehat{x}(\widehat{m}),\widehat{m},t)$$ , $$H(\widehat{x}(\widehat{m}),\widehat{m},t)$$ , $$\frac{\partial f}{\partial \widehat{m}}$$ соответственно; информационный выход четвертого блока интегрирования соединен с первым информационным входом первого блока формирования произведения и информационным входом блока транспонирования.

Сущность изобретения поясняется фиг.1, где представлено устройство определения технического состояния типовых бортовых систем летательного аппарата, которое содержит:

1 - блок хранения констант;

2.1 - первый блок формирования произведения;

2.2 - второй блок формирования произведения;

2.3 - третий блок формирования произведения;

2.4 - четвертый блок формирования произведения;

2.5 - пятый блок формирования произведения;

2.6 - шестой блок формирования произведения;

2.7 - седьмой блок формирования произведения;

2.8 - восьмой блок формирования произведения;

2.9 - девятый блок формирования произведения;

2.10 - десятый блок формирования произведения;

3 - блок формирования суммы;

4 - блок транспонирования;

5.1 - первый блок формирования разности;

5.2 - второй блок формирования разности;

6.1 - первый блок формирования функции $$\frac{\partial H^T}{\partial \widehat{x}}$$ ;

6.2 - второй блок формирования функции $$f(\widehat{x}(\widehat{m}),\widehat{m},t)$$ ;

6.3 - третий блок формирования функции $$H(\widehat{x}(\widehat{m}),\widehat{m},t)$$ ;

6.4 - первый блок формирования функции $$\frac{\partial f}{\partial \widehat{m}}$$ ;

7.1 - первый блок интегрирования;

7.2 - второй блок интегрирования;

7.3 - третий блок интегрирования;

7.4 - четвертый блок интегрирования.

Устройство определения технического состояния типовых бортовых систем летательного аппарата работает следующим образом: на информационный вход блока 1 хранения констант поступает значение наблюдаемой величины; первый информационный выход блока 1 хранения констант соединен с первым информационным входом блока 3 формирования суммы; второй и третий информационные выходы блока 1 хранения констант соединены соответственно с первым и вторым информационными входами первого блока 2.1 формирования произведения; второй информационный выход блока 1 хранения констант соединен с информационным входом блока 4 транспонирования; третий информационный выход блока 1 хранения констант соединен с первым информационным входом четвертого блока 2.4 формирования произведения; четвертый информационный выход блока 1 хранения констант соединен с первым информационным входом третьего блока 2.3 формирования произведения; четвертый информационный выход блока 1 хранения констант соединен со вторыми информационными входами шестого блока 2.6 и девятого блока 2.9 формирования произведения; пятый информационный выход блока 1 хранения констант соединен с первым информационным входом восьмого блока 2.8 формирования произведения; шестой информационный выход блока 1 хранения констант соединен с первым информационным входом второго блока 5.2 формирования разности; седьмой информационный выход блока 1 хранения констант соединен с первыми информационными входами первого блока 6.1, второго блока 6.2, третьего блока 6.3 и четвертого блока 6.4 формирования функций $$\frac{\partial H^T}{\partial \widehat{x}}$$ , $$f(\widehat{x}(\widehat{m}),\widehat{m},t)$$ , $$H(\widehat{x}(\widehat{m}),\widehat{m},t)$$ , $$\frac{\partial f}{\partial \widehat{m}}$$ соответственно; восьмой информационный выход блока 1 хранения констант соединен со вторыми информационными входами первого блока 6.1, второго блока 6.2, третьего блока 6.3 и четвертого блока 6.4 формирования функций $$\frac{\partial H^T}{\partial \widehat{x}}$$ , $$f(\widehat{x}(\widehat{m}),\widehat{m},t)$$ , $$H(\widehat{x}(\widehat{m}),\widehat{m},t)$$ , $$\frac{\partial f}{\partial \widehat{m}}$$ соответственно; девятый информационный выход блока 1 хранения констант соединен с третьими информационными входами первого блока 6.1, второго блока 6.2, третьего блока 6.3 и четвертого блока 6.4 формирования функций $$\frac{\partial H^T}{\partial \widehat{x}}$$ , $$f(\widehat{x}(\widehat{m}),\widehat{m},t)$$ , $$H(\widehat{x}(\widehat{m}),\widehat{m},t)$$ , $$\frac{\partial f}{\partial \widehat{m}}$$ соответственно; информационный выход первого блока 2.1 формирования произведения соединен с первым информационным входом второго блока 2.2 формирования произведения; информационный выход четвертого блока 2.4 формирования произведения соединен со вторым информационным входом третьего блока 2.3 формирования произведения и первым информационным входом пятого блока 2.5 формирования произведения; информационный выход третьего блока 2.3 формирования произведения соединен со вторым информационным входом второго блока 2.2 формирования произведения; информационный выход второго блока 2.2 формирования произведения соединен со вторым информационным входом блока 3 формирования суммы; информационный выход первого блока 6.1 формирования функции $$\frac{\partial H^T}{\partial \widehat{x}}$$ соединен со вторыми информационными входами пятого блока 2.5 и восьмого блока 2.8 формирования произведения; информационный выход пятого блока 2.5 формирования произведения соединен с первым информационным входом шестого блока 2.6 формирования произведения; информационный выход шестого блока 2.6 формирования произведения соединен с первым информационным входом седьмого блока 2.7 формирования произведения; информационный выход восьмого блока 2.8 формирования произведения соединен с первым информационным входом девятого блока 2.9 формирования произведения; информационный выход девятого блока 2.9 формирования произведения соединен с первым информационным входом десятого блока 2.10 формирования произведения; информационный выход второго блока 6.2 формирования функции $$f(\widehat{x}(\widehat{m}),\widehat{m},t)$$ соединен с первым информационным входом первого блока 5.1 формирования разности; информационный выход третьего блока 6.3 формирования функции $$H(\widehat{x}(\widehat{m}),\widehat{m},t)$$ соединен со вторым информационным входом второго блока 5.2 формирования разности; информационный выход второго блока 5.2 формирования разности соединен со вторыми информационными входами седьмого блока 2.7 и десятого блока 2.10 формирования произведения; информационный выход десятого блока 2.10 формирования произведения соединен со вторым информационным входом первого блока 5.1 формирования разности; информационный выход блока 3 формирования суммы соединен с информационным входом первого блока 7.1 интегрирования; информационный выход седьмого блока 2.7 формирования произведения соединен с информационным входом второго блока 7.2 интегрирования; информационный выход первого блока 5.1 формирования разности соединен с информационным входом третьего блока 7.3 интегрирования; информационный выход четвертого блока 6.4 формирования функции $$\frac{\partial f}{\partial \widehat{m}}$$ соединен с информационным входом четвертого блока 7.4 интегрирования; информационный выход первого блока 7.1 интегрирования соединен со вторым информационным входом первого блока 2.1 формирования произведения и первым информационным входом четвертого блока 2.4 формирования произведения; информационный выход второго блока 7.2 интегрирования соединен со вторыми информационными входами первого блока 6.1, второго блока 6.2, третьего блока 6.3 и четвертого блока 6.4 формирования функций $$\frac{\partial H^T}{\partial \widehat{x}}$$ , $$f(\widehat{x}(\widehat{m}),\widehat{m},t)$$ , $$H(\widehat{x}(\widehat{m}),\widehat{m},t)$$ , $$\frac{\partial f}{\partial \widehat{m}}$$ соответственно; информационный выход третьего блока 7.3 интегрирования соединен с третьими информационными входами первого блока 6.1, второго блока 6.2, третьего блока 6.3 и четвертого блока 6.4 формирования функций $$\frac{\partial H^T}{\partial \widehat{x}}$$ , $$f(\widehat{x}(\widehat{m}),\widehat{m},t)$$ , $$H(\widehat{x}(\widehat{m}),\widehat{m},t)$$ , $$\frac{\partial f}{\partial \widehat{m}}$$ соответственно; информационный выход четвертого блока 7.4 интегрирования соединен с первым информационным входом первого блока 2.1 формирования произведения и информационным входом блока 4 транспонирования.

Заявленное изобретение направлено на повышение эффективности определения технического состояния бортовых систем летательных аппаратов, что также весьма важно на всех этапах создания, экспериментальной обработки и эксплуатации объектов ракетно-космической, авиационной, корабельной и других видов техники.

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

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

Модель системы ориентации и стабилизации представим в виде дифференциального уравнения

$$\{\begin{array}{l}\frac{d\varphi }{dt}=\dot{\varphi }=\omega ,\\ \frac{d\omega }{dt}=\dot{\omega }=Um,\left(1\right)\end{array}$$

ω(0)=30, $$\dot{\omega }\left(0\right)=1$$ ,

применив следующие обозначения x₀=ω, $$x_1={\dot{x}}_0=\varphi$$ , получим

$$\dot{x}=\left[\begin{array}{l}{\dot{x}}_0\\ {\dot{x}}_1\end{array}\right]=f(x,m,t)=\left[\begin{array}{l}{x}_1\\ Um\end{array}\right].\left(2\right)$$

где U - функция управления,

m=(M_упр+M_возм)/I - приведенный управляющий момент, [град/с²],

M_упр - управляющий момент,

M_возм - возмущающий момент,

I - момент инерции.

В идеальном случае, когда внешние возмущения отсутствуют, (M_возм=0), m=M_упр/I.

Приведенный управляющий момент позволяет осуществить обратную связь в системе ориентации и стабилизации, т.к. его значение зависит от отклонения текущего состояния системы g(t) от некоторого заданного состояния g(t)=(φ(t), ω(t)). Значение m выбирается на основании информации о текущем состоянии летательного аппарата. При этом источниками информации являются чувствительные элементы: датчики углов и датчики угловых скоростей с соответствующими уравнениями измерений

$$y=H\left(x,t\right)+n\left(t\right)=\left[\begin{array}{cc}{k}_{\omega }& 0\\ 0& k_{\varphi }\end{array}\right]\left[\begin{array}{l}{x}_0\\ x_1\end{array}\right]+\left|\begin{array}{l}\psi \\ \upsilon \end{array}\right|,\left(3\right)$$

где k_φ, k_ω - передаточные коэффициенты,

ψ, υ - соответствующие погрешности измерений.

Функция управления U имеет вид

$$U=\{\begin{array}{ccc}-1,& при& y\ge y_{ср},\\ 0,& при& \left|y\right|<y_{ср},\\ 1,& при& y\le -y_{ср},\end{array}\left(4\right)$$

где y_ср - векторный порог срабатывания управляющих элементов;

некоторое фиксированное значение наблюдаемого векторного параметра,

y - текущее значение наблюдаемого векторного параметра.

При k_φ, k_ω=1 обобщенная функция управления имеет следующий вид

$$U=\{\begin{array}{ccc}-1& при& \phi \ge {\phi }_{ср},\left|\omega \right|<{\omega }_{ср},\\ 0& при& \phi \ge {\phi }_{ср},\omega \le -{\omega }_{ср,}\\ 1& при& \left|\phi \right|<{\phi }_{ср},\omega \le -{\omega }_{ср},\\ 1& при& \phi \le -{\phi }_{ср},\left|\omega \right|<-{\omega }_{ср},\\ 1& при& \omega \le -{\omega }_{ср},\omega \le -{\omega }_{ср},\\ 0& при& \omega \le -{\omega }_{ср},\omega \ge {\omega }_{ср},\\ -1& при& \left|\phi \right|<-{\phi }_{ср},\omega \ge {\omega }_{ср},\\ -1& при& \phi \ge {\phi }_{ср},\omega \ge {\omega }_{ср},\\ 0& при& \left|\phi \right|<{\phi }_{ср},\left|\omega \right|<{\omega }_{ср}.\end{array}\left(5\right)$$

Модель неизвестного неслучайного параметра

$$\dot{m}=\eta ,m\left(0\right)=0.\left(6\right)$$

Критерий качества идентификации имеет вид

$$J_1=\frac{1}{2}{\displaystyle \underset{0}{\overset{T}{\int }}{[y-H(\widehat{x}(\widehat{m}),t)]}^T{N}^{-1}}[y-H(\widehat{x}(\widehat{m}),t)]dt\to \min ,\left(7\right)$$

Таким образом, модель системы ориентации и стабилизации летательного аппарата при вращении вокруг одной из осей может быть задана следующим образом: вращение летательного аппарата вокруг центра масс описывается уравнением (2) с функцией управления вида (5) и уравнениями измерений датчиков углов и датчиков угловой скорости вида (3).

Требуется провести идентификацию приведенного управляющего момента m=1.29, при условии что φ_ср=-0.5 [град], а ω_ср=2,9 [град/с²].

Алгоритм идентификации на основе использования вариационных принципов и итерационной регуляризации [3, 4] имеет вид

$$\begin{array}{l}\dot{\widehat{m}}=P{G}^T\frac{\partial H^T}{\partial \widehat{x}}{N}^{-1}(y-H(\widehat{x}(\widehat{m}),\widehat{m},t)),\widehat{m}(t_0)=0,\\ \widehat{P}={\alpha }^{-1}+P{G}^T{N}^{-1}GP,P(t_0)-0.1,\left(8\right)\\ \widehat{G}=\frac{\partial f}{\partial \widehat{m}},G(t_0)={\left[\begin{array}{cc}0& 0.1\end{array}\right]}^T,\\ \dot{\widehat{x}}=f(\widehat{x}(\widehat{m}),\widehat{m},t)-{\mu }^{-1}\frac{\partial H^T}{\partial \widehat{x}}{N}^{-1}(y-H(\widehat{x}(\widehat{m}),\widehat{m},t),\widehat{x}(t_0)={\left[\begin{array}{cc}0.5& 3\end{array}\right]}^T,\end{array}$$

где µ^-1=6.7, α^-1=1.

Результаты идентификации приведенного управляющего момента разработанным алгоритмом (8) и фильтром Калмана [5-7] представлены на фиг.2.

Относительная погрешность оценок составляет соответственно δОПМиИР=0,067%, δФК=0,603% в установившемся режиме. Анализ полученных данных позволяет сделать вывод, что фильтр Калмана уступает разработанному алгоритму по точности оцениваемого параметра.

Следовательно, в качестве алгоритма идентификации для определения технического состояния системы ориентации и стабилизации используем разработанный алгоритм (8).

Однако, помимо процедуры идентификации, процесс определения технического состояния предполагает проведение целого ряда этапов: таких как получение телеметрической информации с борта летательного аппарата, предварительную и первичную обработку телеизмерений, принятие решения о техническом состоянии. Поэтому рассмотрим пример решения задачи определения технического состояния бортовой системы летательного аппарата на примере типовой активной системы ориентации и стабилизации летательного аппарата [1, 8]. В соответствии со структурой процесса технического состояния примем, что априорно задано:

1. Множество классов технических состояний

Пусть на основании анализа функционирования системы ориентации и стабилизации летательного аппарата заданы следующие виды технического состояния {A_j}, j=0, 1, 2;

A₀ - система ориентации и стабилизации работоспособна;

A₁ - система ориентации и стабилизации неработоспособна вследствие изменения тяги двигателя ориентации и стабилизации;

A₂ - система ориентации и стабилизации неработоспособна вследствие выхода из строя датчика угловых скоростей.

2. Вектор диагностических признаков Q

На основании выбранной модели система ориентации и стабилизации определен вектор диагностических признаков, обеспечивающий наблюдаемость технических состояний Q=[m, k_ω]. При этом неизвестный параметр m рассчитывается в соответствии с алгоритмом (8), а передаточный коэффициент k_ω определяется по результатам расчетных и телеметрических измерений.

3. Модель системы диагностирования R

Модель системы диагностирования задается в виде матрицы соответствия R между техническими состояниями A_j и вектором диагностических признаков Q, т.е. отображения R:A→Q

Таблица 1 Зависимость R от A_j и Q. A Вид технического состояния Q A₀ A₁ A₂ m в норме не в норме в норме k_ω в норме в норме не в норме (k_ω=0)

4. Модель информационно-телеметрической системы

Модель информационно-телеметрической системы представляет собой оператор преобразования L:{g}→{λ}, где g - вектор состояния, λ - вектор телеметрируемых параметров. Для система ориентации и стабилизации оператор L по аналогии с уравнениями (2) можно представить системой уравнений:

$$\{\begin{array}{l}{\lambda }_{\varphi }\left(t\right)=h_{\varphi }\varphi \left(t\right)+\psi ,\\ {\lambda }_{\omega }\left(t\right)=h_{\omega }\omega \left(t\right)+\upsilon ,\end{array}\left(9\right)$$

или в векторно-матричной форме

$$\left|\begin{array}{l}{\lambda }_{\varphi }\left(t\right)\\ {\lambda }_{\omega }\left(t\right)\end{array}\right|=\left|\begin{array}{cc}{h}_{\varphi }& 0\\ 0& h_{\omega }\end{array}\right|\left|\begin{array}{l}\varphi \left(t\right)\\ \omega \left(t\right)\end{array}\right|+\left|\begin{array}{l}\psi \\ \upsilon \end{array}\right|,\left(10\right)$$

или в общем виде

$$\lambda \left(t\right)=Hg\left(t\right)+w,\left(11\right)$$

где $$H=\left|\begin{array}{cc}{h}_{\varphi }& 0\\ 0& h_{\omega }\end{array}\right|$$ - матрица наблюдений,

$$w={\left|\begin{array}{cc}\psi & \upsilon \end{array}\right|}^T$$ - вектор погрешностей измерений.

5. Модель системы предварительной и первичной обработки результатов телеизмерений

Модель системы предварительной и первичной обработки представляет собой оператор преобразования G:{λ}→{x^∗}, где x^∗ - оценка вектора состояния, реализованного набором алгоритмов предварительной и первичной обработки {G_r}, r=1…N. Для системы ориентации и стабилизации летательного аппарата

$$G:\left|\begin{array}{l}{\lambda }_{\varphi }\left(t\right)\\ {\lambda }_{\omega }\left(t\right)\end{array}\right|\to \left|\begin{array}{l}{\varphi }^{\ast }\left(t\right)\\ {\omega }^{\ast }\left(t\right)\end{array}\right|.\left(12\right)$$

6. Модель системы вторичной обработки результатов телеизмерений

Модель системы вторичной обработки результатов телеизмерений представляет собой оператор преобразования Id:{x^∗}→{Q^∗}, где Q^∗ - оценка вектора диагностических признаков, реализованного алгоритмом вторичной обработки (в условиях постановки данной задачи алгоритмом многопараметрической идентификации на основе объединенного принципа максимума и итерационной регуляризации): Id_p, p=1…M.

7. Показатель эффективности решения задачи идентификации

В качестве показателя эффективности решения задачи принята норма разности реализации вектора диагностических признаков Q^∗, оцененного в результате вторичной обработки, с априорно заданными значениями Q_j соответствующими различным видам технических состояний A_j:

$$d_j=\Vert Q^{\ast }-Q_j\Vert$$ .

8. Критерий эффективности решения задачи идентификации

В качестве критерия эффективности выберем требование обеспечения минимума выбранного показателя

J₁→min.

Таким образом, задача оценивания технического состояния системы ориентации и стабилизации формулируется следующим образом:

Дано:

1. Класс моделей системы в виде векторных дифференциальных уравнений (2).

2. Модель наблюдения (3).

3. Модель управления (5).

4. Модель динамики идентифицируемого параметра (6).

5. Критерий оценивания (7).

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

9. Практические рекомендации по принятию решения о техническом состоянии системы ориентации и стабилизации

Целью функционирования подсистемы принятия решений является обеспечение высокой достоверности принятия решения о техническом состоянии системы ориентации и стабилизации летательного аппарата. На практике принятие решения осуществляется в условиях, когда действуют факторы различной природы (случайные, неопределенные), когда возможны различные варианты, и требуется выбрать предпочтительный. Существует несколько вариантов решений: метрический, вероятностный и возможный. В условиях наглядной демонстрации преимуществ разработанного метода рассматривается только метрический вариант решения, заключающийся в последовательном уточнении оцениваемых параметров системы ориентации и стабилизации на каждой итерации до момента, когда точность и вероятность принятия решения о техническом состоянии системы ориентации и стабилизации будет однозначной.

Устройство работает следующим образом. В исходном состоянии на вход устройства подается y, а в блоке 1 хранения констант записываются начальные значения α^-1, G, P, N^-1, µ^-1, y, t, $$\widehat{m}$$ , $$\widehat{x}$$ . Значение α^-1 с первого информационного выхода блока 1 хранения констант поступает на первый информационный вход блока 3 формирования суммы. Значение G со второго информационного выхода блока 1 хранения констант поступает на первый информационный вход первого блока 2.1 формирования произведения и на информационный вход блока 4 транспонирования. Значение P с третьего информационного выхода блока 1 хранения констант поступает на второй информационный вход первого блока 2.1 формирования произведения и первый информационный вход четвертого блока 2.4 формирования произведения. Значение N^-1 с четвертого информационного выхода блока 1 хранения констант поступает на первый информационный вход третьего блока 2.3 формирования произведения и на вторые информационные входы шестого блока 2.6 и девятого блока 2.9 формирования произведения. Значение µ^-1 с пятого информационного выхода блока 1 хранения констант поступает на первый информационный вход восьмого блока 2.8 формирования произведения. Значение y с шестого информационного выхода блока 1 хранения констант поступает на первый информационный вход второго блока 5.2 формирования разности. Значение t с седьмого информационного выхода блока 1 хранения констант поступает на первые информационные входы первого блока 6.1, второго блока 6.2, третьего блока 6.3, четвертого блока 6.4 формирования функций $$\frac{\partial H^T}{\partial \widehat{x}}$$ , $$f(\widehat{x}(\widehat{m}),\widehat{m},t)$$ , $$H(\widehat{x}(\widehat{m}),\widehat{m},t)$$ , $$\frac{\partial f}{\partial \widehat{m}}$$ соответственно. Значение $$\widehat{m}$$ с восьмого информационного выхода блока 1 хранения констант поступает на вторые информационные входы первого блока 6.1, второго блока 6.2, третьего блока 6.3, четвертого блока 6.4 формирования функций $$\frac{\partial H^T}{\partial \widehat{x}}$$ , $$f(\widehat{x}(\widehat{m}),\widehat{m},t)$$ , $$H(\widehat{x}(\widehat{m}),\widehat{m},t)$$ , $$\frac{\partial f}{\partial \widehat{m}}$$ соответственно. Значение $$\widehat{x}$$ с девятого информационного выхода блока 1 хранения констант поступает на третьи информационные входы первого блока 6.1, второго блока 6.2, третьего блока 6.3, четвертого блока 6.4 формирования функций $$\frac{\partial H^T}{\partial \widehat{x}}$$ , $$f(\widehat{x}(\widehat{m}),\widehat{m},t)$$ , $$H(\widehat{x}(\widehat{m}),\widehat{m},t)$$ , $$\frac{\partial f}{\partial \widehat{m}}$$ соответственно. Значение GP с информационного выхода первого блока 2.1 формирования произведения поступает на первый информационный вход второго блока 2.2 формирования произведения. Значение G^T с информационного выхода блока 4 транспонирования поступает на второй информационный вход четвертого блока 2.4 формирования произведения, с информационного выхода которого значение PG^T поступает на второй информационный вход третьего блока 2.3 формирования произведения и первый информационный вход пятого блока 2.5 формирования произведения. Значение PG^TN^-1 с информационного выхода третьего блока 2.3 формирования произведения поступает на второй информационный вход второго блока 2.2 формирования произведения, с информационного выхода которого значение PG^TN^-1GP поступает на второй информационный вход блока 3 формирования суммы. Значение α^-1+PG^TN^-1GP с информационного выхода блока 3 формирования суммы поступает на информационный вход первого блока 7.1 интегрирования. Значение $$\frac{\partial H^T}{\partial \widehat{x}}$$ с информационного выхода первого блока 6.1 формирования функции $$\frac{\partial H^T}{\partial \widehat{x}}$$ поступает на вторые информационные входы пятого блока 2.5 и восьмого блока 2.8 формирования произведения. Значение $$P{G}^T\frac{\partial H^T}{\partial \widehat{x}}$$ с информационного выхода пятого блока 2.5 формирования произведения поступает на первый информационный вход шестого блока 2.6 формирования произведения, с информационного выхода которого значение $$P{G}^T\frac{\partial H^T}{\partial \widehat{x}}{N}^{-1}$$ поступает на первый информационный вход седьмого блока 2.7 формирования произведения. Значение $${\mu }^{-1}\frac{\partial H^T}{\partial \widehat{x}}$$ с информационного выхода восьмого блока 2.8 формирования произведения поступает на первый информационный вход девятого блока 2.9 формирования произведения, с информационного выхода которого значение $${\mu }^{-1}\frac{\partial H^T}{\partial \widehat{x}}{N}^{-1}$$ поступает на первый информационный вход десятого блока 2.10 формирования произведения. Значение $$f(\widehat{x}(\widehat{m}),\widehat{m},t)$$ с информационного выхода второго блока 6.2 формирования функции $$f(\widehat{x}(\widehat{m}),\widehat{m},t)$$ поступает на первый информационный вход первого блока 5.1 формирования разности. Значение $$H(\widehat{x}(\widehat{m}),\widehat{m},t)$$ с информационного выхода третьего блока 6.3 формирования функции $$H(\widehat{x}(\widehat{m}),\widehat{m},t)$$ поступает на второй информационный вход второго блока 5.2 формирования разности, с информационного выхода которого значение y - $$H(\widehat{x}(\widehat{m}),\widehat{m},t)$$ поступает на вторые информационные входы седьмого блока 2.7 и десятого блока 2.10 формирования произведения. Значение $$P{G}^T\frac{\partial H^T}{\partial \widehat{x}}{N}^{-1}(y-H(\widehat{x}(\widehat{m}),\widehat{m},t))$$ с информационного выхода седьмого блока 2.7 формирования произведения поступает на информационный вход второго блока 7.2 интегрирования. Значение $${\mu }^{-1}\frac{\partial H^T}{\partial \widehat{x}}{N}^{-1}(y-H(\widehat{x}(\widehat{m}),\widehat{m},t))$$ с информационного выхода десятого блока 2.10 формирования произведения поступает на второй информационный вход первого блока 5.1 формирования разности, с информационного выхода которого значение $$f(\widehat{x}(\widehat{m}),\widehat{m},t)-{\mu }^{-1}\frac{\partial H^T}{\partial \widehat{x}}{N}^{-1}(y-H(\widehat{x}(\widehat{m}),\widehat{m},t))$$ поступает на информационный вход третьего блока 7.3 интегрирования. Значение $$\frac{\partial f}{\partial \widehat{m}}$$ с информационного выхода четвертого блока 6.4 формирования функции $$\frac{\partial f}{\partial \widehat{m}}$$ поступает на информационный вход четвертого блока 7.4 интегрирования. Значение $$\widehat{P}$$ с информационного выхода первого блока 7.1 интегрирования подается на второй информационный вход первого блока 2.1 формирования произведения и первый информационный вход четвертого блока 2.4 формирования произведения. Искомое значение неизвестного приведенного управляющего момента $$\widehat{m}$$ с информационного выхода второго блока 7.2 интегрирования снимается и подается на вторые информационные входы первого блока 6.1, второго блока 6.2, третьего блока 6.3, четвертого блока 6.4 формирования функции $$\frac{\partial H^T}{\partial \widehat{x}}$$ , $$f(\widehat{x}(\widehat{m}),\widehat{m},t)$$ , $$H(\widehat{x}(\widehat{m}),\widehat{m},t)$$ , $$\frac{\partial f}{\partial \widehat{m}}$$ соответственно. Значение вектора состояния $$\widehat{x}$$ с информационного выхода третьего блока 7.3 интегрирования снимается и подается на третьи информационные входы первого блока 6.1, второго блока 6.2, третьего блока 6.3, четвертого блока 6.4 формирования функции $$\frac{\partial H^T}{\partial \widehat{x}}$$ , $$f(\widehat{x}(\widehat{m}),\widehat{m},t)$$ , $$H(\widehat{x}(\widehat{m}),\widehat{m},t)$$ , $$\frac{\partial f}{\partial \widehat{m}}$$ соответственно. Значение $$\widehat{G}$$ с информационного выхода четвертого блока 7.4 интегрирования подается на первый информационный вход первого блока 2.1 формирования произведения и информационный выход блока 4 транспонирования.

ИСТОЧНИКИ ИНФОРМАЦИИ

1. Мороз В.И., Ромашенко Е.К., Потюпкин А.Ю. Решение задач математической обработки результатов телеизмерений. МО РФ, М., 2001 г., 163 с.

2. Патент РФ рег. №2464615 от 20.10.2012 г. Устройство идентификации параметров динамических систем на основе вариационных принципов. Андрашитов Д.С., Костоглотов А.А., Кузнецов А.А., Лазаренко С.В., Сметанникова Н.А.

3. Андрашитов Д.С., Костоглотов А.А., Лазаренко С.В., Дерябкин И.В. Вариационный метод многопараметрической идентификации динамических систем на основе итерационной регуляризации // Успехи современной радиоэлектроники. №6, 2012, с.67-72.

4. Лурье А.И. Аналитическая механика. М.: Гос. изд. физ.-мат. лит., 1961, с.642-648.

5. Сейдж Э.П., Меле Д.Л. Идентификация систем управления. М.: Наука, 1974. с.204-212.

6. Справочник по теории автоматического управления / Под ред. А.А. Красовского. М.: Наука, 1987, с.347-360.

7. Тихонов А.Н., Арсенин В.Я. Методы решения некорректных задач. М.: Наука, 1986, с.15-18.

8. Фарина А., Студер Ф. Цифровая обработка радиолокационной информации. Сопровождение целей. - М.: Радио и связь, 1993, с.38-47.

Изобретение относится к области цифровой вычислительной техники и может быть использовано в автоматических и автоматизированных системах различного назначения для определения технического состояния по результатам идентификации параметров бортовых систем летательного аппарата. Техническим результатом является повышение эффективности определения технического состояния по результатам идентификации параметров бортовых систем летательного аппарата. Устройство содержит: блок хранения констант; первый, второй, третий, четвертый, пятый, шестой, седьмой, восьмой, девятый, десятый блоки формирования произведения; блок формирования суммы; блок транспонирования; первый, второй блоки формирования разности; первый, второй, третий, четвертый блоки формирования функций , , , соответственно; первый, второй, третий, четвертый блоки интегрирования. 2 ил., 1 табл.

Устройство определения технического состояния бортовых систем летательных аппаратов по результатам оценки параметров, содержащее блок хранения констант; первый, второй, третий, четвертый, пятый, шестой, седьмой, восьмой, девятый, десятый блоки формирования произведения; блок формирования суммы; блок транспонирования; первый, второй блоки формирования разности; первый, второй, третий, четвертый блоки формирования функций $$\frac{\partial H^T}{\partial \hat{x}}$$ , $$f(\hat{x}(\hat{m}),\hat{m},t)$$ , $$H(\hat{x}(\hat{m}),\hat{m},t)$$ , $$\frac{\partial f}{\partial \hat{m}}$$ соответственно; первый, второй, третий, четвертый блоки интегрирования; при этом информационный вход блока хранения констант является входом устройства, а информационный выход второго блока интегрирования является выходом устройства, отличающееся тем, что из устройства исключены второй блок хранения констант, третий блок формирования функции (1-γ_i-1), первый, второй, третий, четвертый и пятый блок формирования производной, одиннадцатый блок формирования произведения, третий блок формирования разности, второй и третий блок формирования суммы, второй блок транспонирования, пятый и шестой блоки интегрирования, в устройство включены первый блок формирования функции $$\frac{\partial H^T}{\partial \hat{x}}$$ , второй блок формирования функции $$f(\hat{x}(\hat{m}),\hat{m},t)$$ , третий блок формирования функции $$H(\hat{x}(\hat{m}),\hat{m},t)$$ , четвертый блок формирования функции $$\frac{\partial f}{\partial \hat{m}}$$ , при этом седьмой информационный выход блока хранения констант соединен с первыми информационными входами первого, второго, третьего, четвертого блоков формирования функций $$\frac{\partial H^T}{\partial \hat{x}}$$ , $$f(\hat{x}(\hat{m}),\hat{m},t)$$ , $$H(\hat{x}(\hat{m}),\hat{m},t)$$ , $$\frac{\partial f}{\partial \hat{m}}$$ соответственно, восьмой информационный выход блока хранения констант соединен со вторыми информационными входами первого, второго, третьего, четвертого блоков формирования функций $$\frac{\partial H^T}{\partial \hat{x}}$$ , $$f(\hat{x}(\hat{m}),\hat{m},t)$$ , $$H(\hat{x}(\hat{m}),\hat{m},t)$$ , $$\frac{\partial f}{\partial \hat{m}}$$ соответственно; девятый информационный выход блока хранения констант соединен с третьими информационными входами первого, второго, третьего, четвертого блоков формирования функций $$\frac{\partial H^T}{\partial \hat{x}}$$ , $$f(\hat{x}(\hat{m}),\hat{m},t)$$ , $$H(\hat{x}(\hat{m}),\hat{m},t)$$ , $$\frac{\partial f}{\partial \hat{m}}$$ соответственно, информационный выход первого блока формирования функции $$\frac{\partial H^T}{\partial \hat{x}}$$ соединен со вторыми информационными входами пятого и восьмого блоков формирования произведения, информационный выход второго блока формирования функции $$f(\hat{x}(\hat{m}),\hat{m},t)$$ соединен с первым информационным входом первого блока формирования разности, информационный выход третьего блока формирования функции $$H(\hat{x}(\hat{m}),\hat{m},t)$$ соединен со вторым информационным входом второго блока формирования разности, информационный выход четвертого блока формирования функции $$\frac{\partial f}{\partial \hat{m}}$$ соединен с информационным входом четвертого блока интегрирования.