Вестник Военного инновационного технополиса «ЭРА», 2022, T. 3, № 4, стр. 387-391
ОБ ОПТИМИЗАЦИИ РЕГУЛЯТОРОВ СУДОВОГО ВЫЧИСЛИТЕЛЯ, ОБЕСПЕЧИВАЮЩИХ АВТОНОМНОСТЬ УПРАВЛЯЮЩЕЙ СИСТЕМЕ
Г. Е. Панамарев 1, *, В. В. Попов 2, А. П. Троеклазов 2
1 Военный инновационный технополис “ЭРА”
Анапа, Россия
2 Государственный морской университет им. адмирала Ф.Ф. Ушакова”
Новороссийск, Россия
* E-mail: era_lab5@mil.ru
Поступила в редакцию 15.03.2022
После доработки 20.03.2022
Принята к публикации 20.03.2022
Аннотация
Показано построение обратной связи для обеспечения системе свойства автономности при работе судового вычислителя в режиме функционирования ввода управляющих воздействий при сложной обстановке движения судна на курсе.
ВВЕДЕНИЕ
В современном регламенте управления судном на курсе широко применяются оптимизированные системы обработки информационных кортежей для помощи в принятии решений судоводителем и непосредственной работы с системными пакетами программного продукта с элементами искусственного интеллекта [1].
Постановка задачи. Пусть поведение объекта регулирования как сегмента управления судового вычислителя описывается системой дифференциальных уравнений вида
где x – n-мерный вектор фазовых координат системы, u – r-мерный вектор управления, y – r‑мерный вектор постановки управляющей задачи, φ – n-мерная вектор-функция, C – матрица $r \times n$.ПРИКЛАДНЫЕ ЗАДАЧИ УПРАВЛЯЮЩИХ ВОЗДЕЙСТВИЙ
В ряде прикладных задач имеется возможность развязывать управляющие воздействия $\left( {{{u}_{1}},{{u}_{1}},...,{{u}_{r}}} \right)$ = u. Относительно координат вектора $y = \left( {{{y}_{1}},{{y}_{1}},...,{{y}_{r}}} \right)$ построение системы предполагается так, чтобы управление ${{u}_{i}}$ действовало лишь на координату ${{y}_{i}}$, не возмущая координат ${{y}_{j}}\left( {j \ne i} \right)$. В исследовании систему (1), обладающую такими свойствами, примем как автономно управляемую, где значение координаты ${{y}_{i}}$ в любой момент времени $t$ будет функционалом системы вида
(2)
${{y}_{i}}\left( t \right) = {{y}_{i}}\left( {x\left( 0 \right),u\left( \tau \right),{{{\bar {u}}}_{i}}\left( \tau \right)} \right);\quad \left( {0 \leqslant \tau \leqslant t} \right),$Итак, система (1) автономно управляема по $i$‑й координате при функционале (2), не зависящем от ${{\bar {u}}_{i}}$. Тогда система полностью автономна при управляемости по всем координатам. Примем, что объект регулирования описывается системой уравнений
где A, B, C – постоянные матрицы размера $n \times n$, $n \times r$, $r \times n$. Система (3) свойствами автономности не обладает [2].Проведем исследования вариации применимости обработки данных в различных конфигурациях постановки задач законов автономной системы и асимптотического регулирования. Проанализируем множество законов обратной связи, обеспечивающих системе заданные свойства. Примем входное воздействие в виде двух слагаемых
где $f\left( t \right)$ – внешнее управляющее воздействие, $\vartheta \left( x \right)$ – определяющий закон обратной связи.Условия автономного регулирования. Исследуем скалярные произведения вида
где ${{d}^{i}}$ – $i$-вектор столбец, ${{d}_{j}}$ – $j$-я вектор-строка произвольной матрицы D.Пусть ${{n}_{{ij}}}$ – наименьшее значение si, при котором скалярное произведение (5) отличается от нуля, тогда
Чтобы система обладала обратной связью $\vartheta \left( x \right)$, обеспечивающей системе (3), (4) автономность по i-й координате, необходимо и достаточно выполнение условий
Пусть существует $\vartheta \left( x \right)$, при котором система (3) и (4) полностью автономна и условия (6) не выполнимы, тогда хотя бы для одного ${{n}_{{ij}}}$:
(7)
$\begin{gathered} ({{c}_{i}},A_{{ij}}^{n}{{b}^{j}}) = e_{i}^{j} \ne 0;\quad ({{c}_{i}},{{A}^{s}}{{b}^{j}}) = 0e_{i}^{j}; \\ \left( {s = 0,1,...,n - 1;{{n}_{{ij}}} - 1} \right). \\ \end{gathered} $Проведем в уравнениях (3) модификацию относительно замены переменных:
В исследовании принято, что каждая из вектор-строк ${{l}_{1}},...,{{l}_{p}}$ матрицы D линейно независима с остальными вектор-строками. Тогда векторы ${{c}_{i}},A{\text{*}}{{c}_{i}}...A{{{\text{*}}}^{{nij}}}{{c}_{i}}$ линейно независимы, следовательно, матрица D, не особая.
Пусть существуют числа ${{p}_{i}}$, такие что
где ${{p}_{i}} = 0$ на все величины. Тогда исследуем равенство, определяющее регулирование уровня воздействия видаТогда в новых переменных модифицированное уравнение принимает вид
(8)
$\begin{gathered} {{{\dot {z}}}_{{\left( {{{n}_{{ij}}} + 1} \right)}}} = \sum\limits_{p = 1}^n {{{h}_{p}}{{z}_{p}}} + e_{i}^{j}{{u}_{j}} + k{{u}_{i}}, \\ {{{\dot {z}}}_{m}} = \sum\limits_{p = 1}^n {{{g}_{{mp}}}{{z}_{p}} + \sum\limits_{p = 1}^r {{{h}_{{mp}}}{{u}_{p}}} } , \\ \end{gathered} $Пусть в (8) выбран произвольный непрерывно дифференцируемый закон обратной связи $\vartheta \left( z \right)$ при условии $\left( {u = \vartheta \left( z \right) + f\left( t \right)} \right)$, тогда вариации фазовых координат при варьировании управления описываются системой уравнений
(9)
$\begin{gathered} (\dot {\delta }{{{\dot {z}}}_{{({{r}_{{ij}}} + 1)}}}) = \sum\limits_{p = 1}^n {{{h}_{p}}\delta {{z}_{p}}} + e_{i}^{j}\sum\limits_{p = 1}^n {\frac{{\partial {{\vartheta }_{i}}}}{{\partial {{z}_{p}}}}} \delta {{z}_{p}} + \\ \, + k\sum\limits_{p = 1}^n {\frac{{\partial {{\vartheta }_{i}}}}{{\partial {{z}_{p}}}}} \delta {{z}_{p}} + e_{i}^{j}\delta {{f}_{j}} + k\delta {{f}_{j}}, \\ \left( {\dot {\delta }{{z}_{m}}} \right) = \sum\limits_{p = 1}^n {{{g}_{{mp}}}\delta {{z}_{p}}} + \sum\limits_{p = 1}^r {{{h}_{{mp}}}} \left( {\sum\limits_{p = 1}^n {\frac{{\partial {{\vartheta }_{p}}}}{{\partial {{z}_{s}}}}} + \delta {{f}_{s}}} \right), \\ \end{gathered} $РЕГУЛЯТОРЫ УПРАВЛЕНИЯ
Если система (3), (4) автономно управляема по i-й координате и равенство при произвольных $\delta {{f}_{i}}\left( {j \ne i} \right)$:
(10)
$\delta {{f}_{i}}\left( t \right) \equiv 0,\quad {\text{то}}\quad \delta {{y}_{i}}\left( t \right) \equiv 0.$Тогда из (9), (10) следует, что δz1(t) = $ = \delta {{z}_{2}}(t) = ... = \dot {\delta }{{z}_{{({{n}_{{ij}}} + 1)}}} \equiv 0$, и уравнение системы (9) дает необходимое условие автономности:
(11)
$\begin{gathered} \sum\limits_{p = \left( {{{n}_{{ij}}} + 2} \right)}^n {{{h}_{p}}\delta {{z}_{p}}} + \sum\limits_{p = \left( {{{n}_{{ij}}} + 2} \right)}^n {\left( {e_{i}^{j}\frac{{\partial {{\vartheta }_{j}}}}{{\partial {{z}_{p}}}} + k\frac{{\partial {{\vartheta }_{j}}}}{{\partial {{z}_{p}}}}} \right)\delta {{z}_{p}} + } \\ \, + e_{i}^{j}\delta {{f}_{i}} = 0. \\ \end{gathered} $Поскольку $\delta {{z}_{i}}\left( t \right)$ непрерывная функция, а $\delta {{f}_{i}}\left( t \right)$ независимы и произвольны, то условие (11) может быть выполнено лишь при $e_{i}^{j} = 0$, что противоречит (7) и удовлетворяет условию оптимального управления [4].
Пусть условия (6) выполнены с оптимальным положением, тогда
(12)
$({{c}_{i}}{{A}^{s}}{{b}^{j}}) = 0\quad {\text{при}}\quad s \leqslant {{n}_{i}}\quad {\text{и}}\quad i \ne j.$Проведем в уравнениях (3) модификацию путем замены переменных, тогда
В выражении каждая из вектор-строк ${{l}_{1}}...,{{l}_{p}}$ матрицы D линейно независима с остальными вектор-строками. Исследуем и покажем, что векторы
(13)
${{c}_{1}}...A{{{\text{*}}}^{{{{n}_{{11}}}}}}{{c}_{1}},{{c}_{2}},...,A{{{\text{*}}}^{{{{n}_{{22}}}}}}{{c}_{2}},...,{{c}_{r}},...A{{{\text{*}}}^{{{{n}_{{rr}}}}}}{{c}_{r}}$Пусть существуют числа ${{p}_{{is}}}$, такие что
Тогда имеем
Отсюда при j = 1 с учетом (12) находим ${{p}_{{1n11}}}e_{1}^{1}0$ и ${{p}_{{1n11}}} = 0$. Если при ${{n}_{{11}}} = {{n}_{{22}}}$ j = 2, получаем ${{p}_{{2n22}}}e_{2}^{2}$ = 0 и ${{p}_{{2n22}}} = 0$.
Исследуем выражения
В новых условиях переменные величины уравнения объекта примут вид
Выражения представим как совокупность связанных подсистем уравнений
(15)
$\begin{gathered} {{{\dot {z}}}^{i}} = {{{\bar {A}}}_{i}}z + {{{\bar {B}}}^{i}}{{u}_{i}};\quad {{{\dot {z}}}^{0}} = {{{\bar {A}}}_{0}}z + {{{\bar {B}}}_{0}}{{u}_{i}}; \\ {{y}^{i}} = z_{i}^{i},\quad {{z}^{i}} = (z_{1}^{i},z_{2}^{i},...,z_{{\left( {{{n}_{{ii}}} + 1} \right)}}^{i});\quad \left( {i = 1,2,...,r} \right). \\ \end{gathered} $Здесь первая подсистема состоит из первых ${{n}_{{ii}}} + 1$ уравнений (14), а вторая из следующих ${{n}_{{22}}} + 1$ и так далее, последняя – из оставшихся $s = n - \sum\limits_{i = 1}^r {{{n}_{{ii}}}} - r$ уравнений (14).
Матрицы ${{\bar {A}}_{i}}$ и векторы ${{\bar {b}}_{i}}$ имеют вид
Последняя строка матрицы $\bar {A}$ состоит из коэффициентов разложения вектора $A{{{\text{*}}}^{{\left( {{{n}_{{ii}}} + 1} \right)}}}$ по вектор-строкам матрицы D [5]. Выберем для объекта исследования (15) законы обратной связи вида
(16)
$\begin{gathered} {{\vartheta }_{i}} = - \frac{1}{{e_{i}^{i}}}\sum\limits_{j = 1}^n {{{{\bar {a}}}_{{{{k}_{i}}j}}}} {{z}_{j}} + {{{\bar {\vartheta }}}_{i}}({{z}^{i}}) = \\ \, = - \frac{1}{{e_{i}^{i}}}({{c}_{i}},A{\kern 1pt} {{*}^{{\left( {{{n}_{{ii}}} + 1} \right)}}}x) + {{{\bar {\vartheta }}}_{i}}\left( {{{y}_{i}},...,\frac{{{{d}^{{{{n}_{{ii}}}}}}{{y}_{i}}}}{{d{{t}^{{{{n}_{{ii}}}}}}}}} \right), \\ \end{gathered} $(17)
${{\dot {z}}^{0}} = {{\bar {A}}^{0}}z + {{\bar {B}}^{0}}\bar {\vartheta } + {{\bar {B}}_{0}}f;\quad \left( {i = 1,2,...,r} \right),$Функционалы (2) для системы (17) имеют вид
(18)
$\begin{gathered} {{y}_{i}}\left( t \right) = z_{1}^{i}\left( t \right) = z_{i}^{1}({{z}^{i}}(0),{{f}_{i}}(\tau )) = {{y}_{i}}\left( {x(0),{{f}_{i}}(\tau )} \right); \\ \left( {0 \leqslant \tau \leqslant t;i = 1,2,...,r} \right), \\ \end{gathered} $АВТОНОМНЫЕ РЕГУЛЯТОРЫ И ОПТИМАЛЬНОСТЬ
Пусть выполнены условия (6), тогда класс регуляторов будет вида
(19)
${{\vartheta }_{i}} = \frac{1}{{({{c}_{i}},{{A}^{{{{n}_{{ii}}}}}}{{b}^{i}})}}({{c}_{i}},{{A}^{{\left( {{{n}_{{ii}}} + 1} \right)}}}x) + {{\bar {\vartheta }}_{i}};\quad \left( {i = 1,2,...,r} \right),$Асимптотическая устойчивость и автономность. В прикладных задачах требуется асимптотическая устойчивость системы (3)–(4) при $f\left( t \right) \equiv 0$. Автономные регуляторы (19) не обладают данными свойствами, хотя устойчивы по координатам ${{y}_{1}},{{y}_{2}},...,{{y}_{r}}$ при выполнении условий (6). Исследуем при отсутствии внешнего воздействия $f\left( t \right) \equiv 0$, тогда в системе (17) получим
(20)
$\begin{gathered} {{{\dot {z}}}^{0}} = {{A}_{0}}{{z}^{0}} + {{A}_{{01}}}\bar {z} + {{{\bar {B}}}_{0}}\bar {\vartheta }\left( {\bar {z}} \right), \\ {\text{при}}\quad \bar {z} = ({{z}^{1}},{{z}^{2}},...,{{z}^{r}}). \\ \end{gathered} $Выбор линейных регуляторов ${{\bar {\vartheta }}_{i}}({{z}^{i}})$. Первые r блоков системы (17) можно принять асимптотически устойчивыми при любом значении $\bar {z}\left( 0 \right)$, модифицировав
Тогда система асимптотически устойчива при устойчивости уравнения
Следовательно, для существования регуляторов, обеспечивающих системе (14) и (16) и одновременно системе (3) и (4) полную автономность и асимптотическую устойчивость при $f\left( t \right) \equiv 0$, достаточным условием можно считать модерацию (21) [7].
Находим уравнения в переменных значениях уравнения управления:
Модифицируем и получаем систему уравнений
(22)
$\begin{gathered} \dot {x} = \left( {A + BK} \right)x + B\bar {\vartheta } \\ K = \left| {\begin{array}{*{20}{c}} {A{\kern 1pt} {{*}^{{\left( {{{n}_{{ii}}} + 1} \right)}}}{{c}_{1}}{\text{/}}e_{1}^{1}} \\ { \cdots \cdots \cdots \cdots \cdots } \\ { \cdots \cdots \cdots \cdots \cdots } \\ {A{\kern 1pt} {{*}^{{\left( {{{n}_{{rr}}} + 1} \right)}}}{{c}_{r}}{\text{/}}e_{r}^{r}} \end{array}} \right|. \\ \end{gathered} $Выбором ${{\bar {\vartheta }}_{i}}\left( {{{y}_{i}},...,\frac{{{{d}^{{{{n}_{{ii}}}}}}{{y}_{i}}}}{{d{{t}^{{{{n}_{{ii}}}}}}}}} \right)$ в виде линейной функции можно обеспечить устойчивость первых r блоков в (17) в равновесном состоянии:
(23)
$\frac{{d{{y}_{i}}}}{{d{{t}^{{\left( {{{n}_{{ii}}} + 1} \right)}}}}} = (A{\kern 1pt} {\text{*}}{{c}_{i}},x) = 0,$Исследуем систему алгебраических уравнений (23), поскольку система $rank\_P = s = \sum\limits_{i = 1}^r {{{n}_{{ii}}}} + r$, тогда система может быть разрешена относительно s координат вектора $x = \left( {{{x}_{1}},{{x}_{2}},...,{{x}_{n}}} \right)$ с помощью обозначения совокупности координат вектором $\bar {x}$, остальные ${{x}_{0}}$ получим:
В исследовании примем $\bar {\vartheta } = 0$, уберем строки, соответствующие $\bar {x}$, и перейдем к системе дифференциальных уравнений порядка $\left( {n - s} \right)$:
Регуляторы класса (16) обладают свойством устойчивости при ${{\bar {\vartheta }}_{i}}$, как линейный регулятор, обеспечивающий асимптотическую устойчивость дифференциальному уравнению
Групповая автономность. Пусть поведение объекта описывается системой дифференциальных уравнений
Примем матрицы $A,B$ постоянными размером $n \times n$ и $n \times r$, наблюдаемые координаты разобьем на r групп:
где ${{C}_{i}}$– постоянные матрицы размером ${{m}_{i}} \times n$ и ${{y}^{i}} = \left( {y_{1}^{i},y_{2}^{i},...,y_{{{{m}_{i}}}}^{i}} \right)$, тогда функционалы (2) для системы (26)–(27) имеют вид(28)
$\begin{gathered} y_{k}^{i} = y_{k}^{i}\left( {x\left( 0 \right),{{u}_{i}}\left( \tau \right),{{{\bar {u}}}_{i}}\left( \tau \right)} \right); \\ \left( {0 \leqslant \tau \leqslant t;\;k = 1,2,...,{{m}_{i}};\;i = 1,2,...,r} \right). \\ \end{gathered} $Для существования регуляторов, обеспечивающих системе полную групповую автономность, достаточно выполнить условия:
Тогда законы управления примут вид
ЗАКЛЮЧЕНИЕ
Таким образом, с помощью достаточных условий оптимального управленческого воздействия на исполнительные системы выведены дифференциальные уравнения для нахождения синтезирующего функционала и оптимального управления с обратной связью. Оптимальное управление динамической частью системы реализуется линейным регулятором по Вознесенскому, а оптимальное управление логической частью определяется рекуррентным уравнением. Исследована применимость оптимального процесса управления судовым вычислителем при управлении судном на курсе, в том числе со счетным множеством переключений логической части системы, которые происходят в фиксированный момент времени.
Список литературы
Лицкевич А.П., Старжинская Н.В., Попов В.В. Математические методы в электродинамике. Новороссийск: ГМУ им. адм. Ф.Ф. Ушакова, 2019. 212 с.
Долматов Б.М., Попов В.В. Информатика. Новороссийск: ГМУ им. адм. Ф.Ф. Ушакова, 2017. 60 с.
Гантмахер Ф.Р. Теория матриц. М.: Наука, 2019. 576 с.
Колмогоров А.Н., Фомин С.В. Элементы теории функций и функционального анализа. 3-е изд. М.: Наука, 2018. 496 с.
Стренг Г. Линейная алгебра и ее применение. М.: Мир, 2020. 456 с.
Райфа Х. Анализ решения. М.: Наука, 2017. 405 с.
Аблязов К.А., Катрюк И.С., Попов В.В. Основы теории надежности и диагностики. Новороссийск: ГМУ им. адм. Ф.Ф. Ушакова, 2018. 212 с.
Малкин И.Г. Теория устойчивости движения. М.: Наука, 2019. 345 с.
Дополнительные материалы отсутствуют.
Инструменты
Вестник Военного инновационного технополиса «ЭРА»