Известия РАН. Теория и системы управления, 2020, № 2, стр. 17-26
Управление линейной mimo-системой по вектору измерения с использованием многоуровневой декомпозиции
Н. Е. Зубов a, *, Е. А. Микрин b, В. Н. Рябченко a
a МГТУ им. Н.Э. Баумана
Москва, Россия
b ПАО “РКК “Энергия”
Москва, Россия
* E-mail: nik.zubov@gmail.com
Поступила в редакцию 28.04.2019
После доработки 08.11.2019
Принята к публикации 25.11.2019
Аннотация
Разработан вариант метода синтеза управления линейной MIMO-системой по вектору выхода, обеспечивающий заданный спектр многомерной динамической системы в пространстве состояний. В основе лежит принцип дуальности задачи управления MIMO-системой по вектору состояния и задачи построения наблюдателя состояния. Алгоритм без каких-либо изменений применим как для непрерывного, так и для дискретного случаев описания математической модели MIMO-системы, не имеет ограничений по заданию элементов спектра замкнутой MIMO-системы, позволяет получать решения задач синтеза в аналитическом виде и осуществлять параметризацию множества эквивалентных законов управления с обратной связью. Приведен пример аналитического синтеза закона управления гипотетическим летательным аппаратом.
Введение. Задача управления спектром движения (полюсами, собственными значениями) линейной динамической MIMO-системой (multi input multi output) по вектору выхода системы относится в теории систем к разряду классических [1–11], однако до сих пор не имеет исчерпывающего решения и относится к задаче высокого уровня сложности [9]. В данной статье рассматривается метод синтеза закона управления с использованием принципа дуальности задачи управления по вектору состояния и задачи построения наблюдателя состояния по вектору измерения.
1. Декомпозиция динамической системы и подход к синтезу закона управления. Будем рассматривать MIMO-систему следующего вида:
где t – непрерывное $t \in {{\mathbb{R}}_{ + }} = [0, + \infty )$ или дискретное $t \in {{\mathbb{N}}_{0}} = \{ 0,1,2, \ldots \} $ время; $x \in {{\mathbb{R}}^{n}}$ – вектор состояния; $u \in {{\mathbb{R}}^{r}}$ – вектор входа; $y \in {{\mathbb{R}}^{m}}$ – вектор выхода; $\mathbb{R}$ – множество действительных чисел; $m + r > n$; σ – символ оператора дифференцирования $\sigma x(t)\mathop \doteq \limits^. \dot {x}(t)$ или оператора сдвига $\sigma x(t) \doteq x(t + 1)$.Предполагается, что спектр MIMO-системы (1.1) совпадает с множеством собственных значений матрицы A и равен
Рассматриваемый здесь метод предполагает выполнение алгоритма, содержащего два шага:
1) многоуровневая декомпозиция динамической системы с помощью техники ортогональных делителей нуля и псевдообратных матриц, т.е. декомпозиции матриц A, B и C (процесс “сверху–вниз”);
2) собственно синтез закона управления с использованием декомпозированных на соответствующих уровнях MIMO-систем, т.е. формировании (“сборки”) матрицы регулятора в законе управления u(t) = –Fy(t) (процесс “снизу–вверх”).
Указанное ранее применение принципа дуальности понимается в том смысле, что синтез закона управления по вектору выхода MIMO-системы (регулятора по выходу) равносилен синтезу наблюдателя Люэнбергера (наблюдателя состояния).
Введем в рассмотрение многоуровневую декомпозицию MIMO-системы (1.1) следующего вида:
нулевой уровень декомпозиции
первый уровень декомпозиции
(1.3)
${{A}_{1}} = {{(С_{0}^{{R \bot }})}^{ + }}{{A}_{0}}С_{0}^{{R \bot }},\quad {{B}_{1}} = {{(С_{0}^{{R \bot }})}^{ + }}{{A}_{0}}{{B}_{0}},\quad {{C}_{1}} = {{C}_{0}}{{A}_{0}}С_{0}^{{R \bot }},$k-й уровень декомпозиции ($1 < k < M$)
(1.4)
${{A}_{k}} = {{(С_{{k - 1}}^{{R \bot }})}^{ + }}{{A}_{{k - 1}}}С_{{k - 1}}^{{R \bot }},\quad {{B}_{k}} = {{(С_{{k - 1}}^{{R \bot }})}^{ + }}{{A}_{{k - 1}}}{{B}_{{k - 1}}},\quad {{C}_{k}} = {{C}_{{k - 1}}}{{A}_{{k - 1}}}С_{{k - 1}}^{{R \bot }},$M-й уровень декомпозиции
(1.5)
${{A}_{M}} = {{(С_{{M - 1}}^{{R \bot }})}^{ + }}{{A}_{{M - 1}}}С_{{M - 1}}^{{R \bot }},\quad {{B}_{M}} = {{(С_{{M - 1}}^{{R \bot }})}^{ + }}{{A}_{{M - 1}}}{{B}_{{M - 1}}},\quad {{C}_{M}} = {{C}_{{M - 1}}}{{A}_{{M - 1}}}С_{{M - 1}}^{{R \bot }}.$Здесь $M = {\text{ceil}}(n{\text{/}}m) - 1$, где ${\text{ceil}}\,( * )$ – операция округления числа “$ * $” в сторону большего значения, т.е. ceil(1.2) = 2, ceil(2.5) = 3 и т.д.; индексом “$R \bot $” обозначен правый делитель нуля [11–17], а индексом “$ + $” – псевдообратные матрицы Мура–Пенроуза [11–17].
Рассмотрим далее собственно синтез закона управления по вектору выхода на основе управления спектром декомпозированных MIMO-систем на соответствующих уровнях (в обратном порядке):
M-й уровень декомпозиции
k-й уровень декомпозиции ($k = \overline {0,М - 1} $)
(1.7)
${{F}_{k}} = B_{k}^{ + }(C_{k}^{ - }{{\Phi }_{k}} - {{A}_{k}}C_{k}^{ - }),\quad C_{k}^{ - } = C_{k}^{ + } - C_{k}^{{R \bot }}B_{k}^{{}}{{F}_{{k + 1}}}.$Следовательно, первая формула выражения (1.7) при k = 0 представляет собой матрицу коэффициентов обратной связи регулятора по выходу и собственно решение задачи синтеза управления по выходу.
2. Основной теоретический результат. Справедливо утверждение.
Теорема. Пусть следующие матрицы существуют и попарно полностью управляемые:
(2.1)
$G_{M}^{{}} = {{(B_{M}^{ \bot }C_{M}^{ + })}^{ + }}B_{M}^{ \bot }{{A}_{M}}C_{M}^{ + },\quad H_{M}^{{}} = {{(B_{M}^{ \bot }C_{M}^{ + })}^{{R \bot }}},$(2.2)
$G_{k}^{{}} = {{(B_{k}^{ \bot }C_{k}^{ - })}^{ + }}B_{k}^{ \bot }{{A}_{k}}C_{k}^{ - },\quad H_{k}^{{}} = {{(B_{k}^{ \bot }C_{k}^{ - })}^{{R \bot }}},$(2.3)
${{\Phi }_{i}} = {{G}_{i}} + {{H}_{i}}K_{i}^{{}} = {{(B_{i}^{ \bot }C_{i}^{ - })}^{ + }}B_{i}^{ \bot }{{A}_{i}}C_{i}^{ - } + {{(B_{i}^{ \bot }C_{i}^{ - })}^{{R \bot }}}K_{i}^{{}},$(2.4)
${\text{eig}}\left( {{{A}_{M}} + {{B}_{M}}{{F}_{M}}{{C}_{M}}} \right) = {\text{eig(}}{{\Phi }_{M}}{\text{),}}$(2.5)
${\text{eig}}\left( {{{A}_{k}} + {{B}_{k}}{{F}_{k}}{{C}_{k}}} \right) = \bigcup\limits_{i = k}^M {{\text{eig(}}{{\Phi }_{i}}{\text{)}}} ,\quad (k = \overline {0,М - 1} ),$(2.6)
${\text{eig}}\left( {{{A}_{0}} + {{B}_{0}}{{F}_{0}}{{C}_{0}}} \right) = {\text{eig}}\left( {A + BFC} \right) = \bigcup\limits_{i = 0}^M {{\text{eig(}}{{\Phi }_{i}}{\text{)}}} \doteq \Lambda .$Здесь индексом “$\, \bot $” обозначена ортогональная матрица левого делителя нуля [11], а при i = M имеет место равенство $C_{M}^{ - } = C_{M}^{ + }$.
Доказательство теоремы приведено в Приложении.
Как и в материалах, представленных в [12–17], в теореме при преобразованиях используются только ортогональные и псевдообратные матрицы, что позволяет, по крайней мере, не ухудшать обусловленность получаемых уравнений. Приведенный теоретический результат применим как для непрерывного, так и для дискретного случаев описания MIMO-системы (1.1), не имеет ограничений по заданию элементов спектра замкнутой MIMO-системы, позволяет получать решения задачи синтеза в аналитическом виде и осуществлять параметризацию множества эквивалентных законов управления с обратной связью. Это подтверждается, в том числе, математическим моделированием, результаты которого демонстрируют высокую точность управления спектром и практическое отсутствие ограничений на размерность MIMO-системы (1.1).
3. Аналитический синтез управления стабилизацией бокового движения маневренного летательного аппарата. Рассмотрим задачу управления движением маневренного летательного аппарата (ЛА) в горизонтальной плоскости (боковой канал движения и управления). В качестве органов управления будем использовать элероны, руль направления, переднее вертикальное оперение и отклоняемый вектор тяги двигателя. В этом случае линеаризованные уравнения бокового канала в отклонениях в соответствии с [18, 19] будут иметь вид
(3.1)
$\left[ {\begin{array}{*{20}{c}} {\dot {\beta }} \\ {{{{\dot {\omega }}}_{x}}} \\ {{{{\dot {\omega }}}_{y}}} \\ {\dot {\gamma }} \end{array}} \right] = \left[ {\begin{array}{*{20}{c}} {{{a}_{{11}}}}&{{{a}_{{12}}}}&{{{a}_{{13}}}}&{{{a}_{{14}}}} \\ {{{a}_{{21}}}}&{{{a}_{{22}}}}&{{{a}_{{23}}}}&0 \\ {{{a}_{{31}}}}&{{{a}_{{32}}}}&{{{a}_{{33}}}}&0 \\ 0&1&{{{a}_{{43}}}}&0 \end{array}} \right]\left[ {\begin{array}{*{20}{c}} \beta \\ {{{\omega }_{x}}} \\ {{{\omega }_{y}}} \\ \gamma \end{array}} \right] + \left[ {\begin{array}{*{20}{c}} {{{b}_{{11}}}}&0&{{{b}_{{13}}}}&0 \\ {{{b}_{{21}}}}&{{{b}_{{22}}}}&{{{b}_{{23}}}}&{{{b}_{{24}}}} \\ {{{b}_{{31}}}}&{{{b}_{{32}}}}&{{{b}_{{33}}}}&{{{b}_{{34}}}} \\ 0&0&0&0 \end{array}} \right]\left[ {\begin{array}{*{20}{c}} {{{\delta }_{э}}} \\ {{{\delta }_{{\text{н}}}}} \\ {{{\delta }_{{{\text{ОВТ}}}}}} \\ {{{\delta }_{{{\text{п}}{\text{.о}}}}}} \end{array}} \right]$.Здесь $\beta $ – угол скольжения; ${{\omega }_{x}}$ – угловая скорость крена; ${{\omega }_{y}}$ – угловая скорость рыскания; $\gamma $ – угол крена; ${{\alpha }_{0}}$ – угол атаки; ${{\vartheta }_{0}}$ – угол тангажа; ${{\delta }_{{\text{н}}}}$ – угол отклонения рулей направления; ${{\delta }_{э}}$ – угол отклонения элеронов; ${{\delta }_{{{\text{ОВТ}}}}}$ – угол отклонения вектора тяги; ${{\delta }_{{{\text{п}}{\text{.о}}}}}$ – угол отклонения переднего вертикального оперения; ${{a}_{{ij}}}$, ${{b}_{{ij}}}$ – коэффициенты линеаризации [19, 20].
Как видно, особенностью модели (3.1) является неполнота ранга матрицы управления, обусловленное “избыточностью” органов управления в боковой плоскости
В обозначениях (1.1) матрицы А и В с учетом (3.1) запишутся так:
(3.2)
$A = \left[ {\begin{array}{*{20}{c}} {{{a}_{{11}}}}&{{{a}_{{12}}}}&{{{a}_{{13}}}}&{{{a}_{{14}}}} \\ {{{a}_{{21}}}}&{{{a}_{{22}}}}&{{{a}_{{23}}}}&0 \\ {{{a}_{{31}}}}&{{{a}_{{32}}}}&{{{a}_{{33}}}}&0 \\ 0&1&{{{a}_{{43}}}}&0 \end{array}} \right],\quad B = \left[ {\begin{array}{*{20}{c}} {{{b}_{{11}}}}&0&{{{b}_{{13}}}}&0 \\ {{{b}_{{21}}}}&{{{b}_{{22}}}}&{{{b}_{{23}}}}&{{{b}_{{24}}}} \\ {{{b}_{{31}}}}&{{{b}_{{32}}}}&{{{b}_{{33}}}}&{{{b}_{{34}}}} \\ 0&0&0&0 \end{array}} \right].$Зададим матрицу выхода следующим образом:
что соответствует отсутствию информации об угловых скоростях движения ЛА.Пусть требуется найти в явном виде формулу закона управления (регулятора), обеспечивающего замкнутой системе “ЛА + система управления” следующий (вообще говоря, произвольный) спектр:
(3.4)
$\Lambda = \{ {{\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\frown}$}}{\lambda } }_{1}},{{\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\frown}$}}{\lambda } }_{2}},{{\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\frown}$}}{\lambda } }_{3}},{{\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\frown}$}}{\lambda } }_{4}}\} \subset \mathbb{C}.$Следует заметить, что аналитическое решение указанной задачи синтеза нельзя получить ранее опубликованными авторами методами в силу неполноты ранга матрицы $B$ (3.2).
Для системы (1.1) с матрицами (3.2), (3.3) рассмотрим многоуровневую декомпозицию, определенную в разд. 1, имеющую в данном случае два уровня декомпозиции: нулевой (1.2) и первый (1.3), для которых
(3.5)
$B_{0}^{ \bot } = \left[ {\begin{array}{*{20}{c}} 0&0&0&1 \end{array}} \right],\quad {{C}^{{R \bot }}} = \left[ {\begin{array}{*{20}{c}} 0&0 \\ 0&{ - 1} \\ 1&0 \\ 0&0 \end{array}} \right],\quad B_{0}^{ + } = \left[ {\begin{array}{*{20}{c}} {b_{{11}}^{ + }}&{b_{{12}}^{ + }}&{b_{{13}}^{ + }}&{b_{{14}}^{ + }} \\ {b_{{21}}^{ + }}&{b_{{22}}^{ + }}&{b_{{23}}^{ + }}&{b_{{24}}^{ + }} \\ {b_{{31}}^{ + }}&{b_{{32}}^{ + }}&{b_{{33}}^{ + }}&{b_{{34}}^{ + }} \\ {b_{{41}}^{ + }}&{b_{{42}}^{ + }}&{b_{{43}}^{ + }}&{b_{{44}}^{ + }} \end{array}} \right],\quad {{C}_{p}} = \left[ {\begin{array}{*{20}{c}} 1&0 \\ 0&0 \\ 0&0 \\ 0&1 \end{array}} \right],$Здесь элементы матрицы $B_{0}^{ + }$ из (3.5) в силу неполноты ее ранга вычислены в соответствии с [1] на основании следующей оригинальной формулы псевдообратной матрицы [21]:
Для проверки условий управляемости на первом уровне декомпозиции, определенных в теореме, вычислим следующие матрицы:
(3.6)
$H_{1}^{{}} = {{(B_{1}^{ \bot }C_{1}^{ + })}^{{R \bot }}} = I_{2}^{{}},\quad H_{1}^{ + } = I_{2}^{{}},$Вычислим для первого уровня декомпозиции ранг следующей блочной матрицы:
что совпадает с числом измеряемых компонент вектора состояния m = 2. Следовательно, первый уровень удовлетворяет критерию управляемости, определенному в теореме.Далее в соответствии с введенной формой регулятора (1.6), (1.7) зададим матрицу, собственные значения которой будут “приписаны” MIMO-системе на первом уровне декомпозиции. В силу очевидной обратимости матриц (3.6) для первого уровня декомпозиции мы вправе выбрать любую матрицу ${{{\text{Ф}}}_{1}}$, обладающую заданным спектром. Для простоты положим
На основании формулы (1.7) при k = 1 для первого уровня будем иметь
Для вычисления матрицы $C_{0}^{ - }$, необходимой для определения регулятора нулевого уровня, воспользуемся второй формулой в (1.7). В результате получим
В соответствии с теоремой проанализируем соблюдение условия управляемости MIMO-системы на нулевом уровне декомпозиции, для чего определим матрицы
Вычислим для нулевого уровня декомпозиции ранг следующей блочной матрицы:
что совпадает с числом измеряемых компонент вектора состояния, т.е. m = 2. Следовательно, нулевой уровень также удовлетворяет критерию управляемости, определенному в теореме.Теперь следует найти Φ0 для нулевого уровня декомпозиции. Для этого осуществим декомпозицию матриц $H_{0}^{{}}$, ${{G}_{0}}$ нулевого уровня на два подуровня и вычислим соответствующие матрицы. В результате получим
Назначим собственное значение ${{\left( {{{{\text{Ф}}}_{0}}} \right)}_{1}} = {{\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\frown}$}}{\lambda } }_{3}}$ и вычислим матрицу коэффициентов обратной связи первого подуровня нулевого уровня декомпозиции. Следовательно, имеем
Вычислим матрицу
Далее назначим собственное значение ${{\left( {{{{\text{Ф}}}_{0}}} \right)}_{0}} = {{\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\frown}$}}{\lambda } }_{4}}$ и найдем матрицу K0 по формуле
В результате с использованием выражения (2.3) получим искомую матрицу ${{{\text{Ф}}}_{0}}$ с собственными значениями ${{\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\frown}$}}{\lambda } }_{3}}$, ${{\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\frown}$}}{\lambda } }_{4}}$, а именно
Используя полученные выше выражения для $C_{0}^{ + }$, $C_{0}^{{R \bot }}$, ${{F}_{1}}$ вычисляем
и, применяя первую формулу из (1.7), в конечном итоге будем иметь следующую формулу регулятора в законе управления по вектору выхода:(3.7)
$F = {{F}_{0}} = B_{0}^{ + }(C_{0}^{ - }{{\Phi }_{0}} - {{A}_{0}}C_{0}^{ - }) = \left[ {\begin{array}{*{20}{c}} {{{F}_{{11}}}}&{{{F}_{{12}}}} \\ {{{F}_{{21}}}}&{{{F}_{{22}}}} \\ {{{F}_{{31}}}}&{{{F}_{{32}}}} \\ {{{F}_{{41}}}}&{{{F}_{{42}}}} \end{array}} \right],$Регулятор (3.7) представляет собой аналитическое решение рассматриваемой задачи синтеза и обеспечивает боковому управляемому движению гипотетического ЛА спектр (3.4).
Аналогичным образом можно получить решение в численном виде. Допустим, требуется обеспечить замкнутой системе “ЛА + система управления” следующий спектр:
когда матрицы A и B имеют вид(3.9)
$A = \left[ {\begin{array}{*{20}{c}} { - 0.1720}&{0.0631}&{0.9980}&{0.0510} \\ { - 26.0500}&{ - 2.7490}&{ - 0.5330}&0 \\ { - 4.3370}&{ - 0.0060}&{ - 0.3010}&0 \\ 0&1&{ - 0.0632{\text{ }}}&0 \end{array}} \right],$(3.10)
$B = \left[ {\begin{array}{*{20}{c}} { - 0.0340}&0&{ - 0.0100}&0 \\ { - 4.7570}&{ - 18.6640}&{ - 1.5000}&{ - 20.0000} \\ { - 3.0700}&{0.6660}&{ - 1.0000}&{ - 1.0000} \\ 0&0&0&0 \end{array}} \right].$Пусть начальное значение вектора состояния в системе единиц СИ равно:
(3.11)
$\left[ {\begin{array}{*{20}{c}} {{{\beta }_{0}}} \\ {{{\omega }_{{x0}}}} \\ {{{\omega }_{{y0}}}} \\ {{{\gamma }_{0}}} \end{array}} \right]{\text{ = }}\left[ {\begin{array}{*{20}{c}} { - 0.0938} \\ { - 0.0243} \\ { - 0.0220} \\ {0.0209} \end{array}} \right].$Числовое значение матрицы F в соответствии с (3.7) будет следующим:
(3.12)
$F = \left[ {\begin{array}{*{20}{c}} { - {\text{5}}{\text{.1600}}}&{ - {\text{16}}{\text{.8687}}} \\ { - {\text{7}}{\text{.6761}}}&{ - {\text{30}}{\text{.7344}}} \\ { - {\text{0}}{\text{.6561}}}&{ - {\text{1}}{\text{.2635}}} \\ {{\text{7}}{\text{.1490}}}&{{\text{32}}{\text{.8486}}} \end{array}} \right].$Графики переходных процессов и величин управляющих воздействий для заданного спектра (3.8), модели ЛА (3.9)–(3.11) и матрицы регулятора (3.12) приведены соответственно на рис. 1 и 2. При этом вычисление спектра матрицы A + BFС дает
что, согласно (3.8), и требовалось получить.
Заключение. Представлен алгоритм синтеза закона управления линейной MIMO-системой с обратной связью по вектору выхода, обеспечивающий заданный спектр многомерной динамической системы в пространстве состояний. Сформулирована и доказана соответствующая теорема.
В основе предлагаемого подхода лежит принцип дуальности задачи управления MIMO-системой по вектору состояния и задачи построения наблюдателя состояния. Алгоритм без каких-либо изменений применим как для непрерывного, так и для дискретного случаев описания математической модели MIMO-системы, не имеет никаких ограничений по заданию элементов спектра замкнутой MIMO-системы, позволяет получать аналитические решения задачи синтеза и осуществлять параметризацию множества эквивалентных законов управления с обратной связью.
Возможности метода продемонстрированы на решении задачи управления боковым движением моделью гипотетического ЛА. Особенностью данной модели является неполнота ранга матрицы управления, обусловленная “избыточностью” органов управления в боковой плоскости.
Список литературы
Bhattachrya S. Sparsity Based Ceedback Design: A new Paradigm in Opportunistic Sensing // Proc. American Control Conf. St. Louis, 2011. P. 3704–3709.
Blumthaler I., Oberst U. Design, Parameterization, and Pole Placement of Stabilizing Output Feedback Compensators Via Injective Cogenerator Quotient Signal Modules // Linear Algebra Appl. 2012. V. 436 (5-2). P. 963–1000.
Eremenko A., Gabrielov A. Pole Placement by Static Output Feedback for Generic Linear Systems // SIAM J. Contr. Opt. 2002. V. 41 (1). P. 303–312.
Franke M. Eigenvalue Assignment by Static Output Feedback – on a new Solvability Condition and the Computation of Low Gain Feedback Matrices // Intern. J. Contr. 2014. V. 87 (1). P. 64–75.
Kaiyang Yanga, Orsi R. Generalized Pole Placement Via Static Output Feedback: A Methodology Based on Projections // Automatica. 2006. V. 42. P. 2143–2150.
Peretz Y. A Randomized Approximation Algorithm for the Minimal-Norm Static-Output-Feedback Problem // Automatica. 2016. V. 63. P. 221–234.
Shimjith S.R., Tiwari A.P., Bandyopadhyay B. Modeling and Control of a Large Nuclear Reactor. Berlin, Heidelberg: Springer-Verlag, 2013.
Wang X.A. On Linear Solutions of the Output Feedback Pole Assignment Problem // IEEE Trans. Autom. Contr. 2013. V. 58(9). P. 2354–2359.
Fu M. Pole Placement Via Static Output Feedback is NP-hard // IEEE Trans. Autom. Contr. 2004. V. 49 (5). P. 855–857.
Van der Woude J.W. A Note on Pole Placement by Static Output Feedback for Single. Input Systems // Systems & Control Letters. 1988. V. 11. P. 285–287.
Зубов Н.Е., Микрин Е.А., Мисриханов М.Ш., Рябченко В.Н. Управление по выходу спектром дескрипторной динамической системы // ДАН. 2016. Т. 468. № 2. С. 134–136.
Зубов Н.Е., Лапин А.В., Микрин Е.А., Рябченко В.Н. Управление по выходу спектром линейной динамической системы на основе подхода Ван-дер-Воуда // ДАН. 2017. Т. 476. № 3. С. 260–263.
Зубов Н.Е., Микрин Е.А., Олейник А.С. и др. Оценка угловой скорости космического аппарата в режиме орбитальной стабилизации по результатам измерений датчика местной вертикали // Вестн. МГТУ. Приборостроение. 2014. № 5. С. 3–15.
Зубов Н.Е., Микрин Е.А., Мисриханов М.Ш., Рябченко В.Н. Стабилизация взаимосвязанных движений летательного аппарата в каналах тангаж-рысканье при отсутствии информации об угле скольжения // Изв. РАН. ТиСУ. 2015. № 1. С. 95–105.
Зубов Н.Е., Микрин Е.А., Мисриханов М.Ш., Рябченко В.Н. Управление по выходу продольным движением летательного аппарата // Изв. РАН. ТиСУ. 2015. № 5. С. 164–175.
Зубов Н.Е., Микрин Е.А., Рябченко В.Н. и др. Синтез законов управления боковым движением летательного аппарата при отсутствии информации об угле скольжения. Аналитическое решение // Изв. вузов. Авиационная техника. 2017. № 1. С. 14–20.
Зубов Н.Е., Микрин Е.А., Рябченко В.Н. Матричные методы в теории и практике систем автоматического управления летательными аппаратами. М.: Изд-во МГТУ им. Н.Э. Баумана, 2016. 666 с.
Рябченко В.Н. Вложение систем. Нерегулярные законы управления // АиТ. 2001. № 7. С. 198–210.
Боднер В.А. Системы управления летательными аппаратами. М.: Машиностроение, 1973. 504 с.
Буков В.Н. Адаптивные прогнозирующие системы управления полетом. М.: Наука, 1987. 230 с.
Зубов Н.Е., Рябченко В.Н. О вычислении псевдообратной матрицы. Общий случай // Вестн. МГТУ. Естественные науки. 2018. № 1 (76). С. 16–25.
Дополнительные материалы отсутствуют.
Инструменты
Известия РАН. Теория и системы управления