CC BY
8ВЕСТН. САМАР. ГОС. ТЕХН. УН-ТА. СЕР. ТЕХНИЧЕСКИЕ НАУКИ. 2021. Т. 29. № 1_
(М: 10.14498ЛесЬ.2021.1.5
УДК 517.977.5
СТРУКТУРНО-ПАРАМЕТРИЧЕСКИЙ СИНТЕЗ СИСТЕМЫ УПРАВЛЕНИЯ НЕНАБЛЮДАЕМЫМ ВЫХОДОМ ОБЪЕКТА С РАСПРЕДЕЛЕННЫМИ ПАРАМЕТРАМИ
А.В. Тычинин1, Ю.А. Тычинина2, Д.А. Рагазин 2
1 ЗАО «ТМ-Сервис»
Россия, 443029, г. Самара, 6-я просека, 157
2 Самарский государственный технический университет Россия, 443100, г. Самара, ул. Молодогвардейская, 244
Аннотация. Рассматривается синтез системы управления объектом с распределенными параметрами. Обозначены проблемы существующих методов описания таких объектов и разработки систем управления ими. Предложен подход для нивелирования недостатков существующих методов, связанных с допущениями о постоянстве физических характеристик объекта и с необходимостью управления ненаблюдаемой пространственной координатой объекта. В качестве примера рассматривается нагрев массивной заготовки в конвекционной печи. Управление по ненаблюдаемой координате объекта реализовано путем синтеза цепной структуры, реализующей передаточную функцию, близкую к обратной передаточной функции объекта. Слабая чувствительность к параметрической нестабильности обеспечивается решением обратной задачи динамики. В результате получена система, обладающая слабой чувствительностью к параметрическим возмущениям объекта и обеспечивающая требуемое поведение ненаблюдаемой координаты пространственной области распределенного объекта.
Ключевые слова: объект с распределенными параметрами, граничные условия 1-го рода, граничные условия 3-го рода, компенсация динамических свойств, обратная передаточная функция объекта, периодические структуры, устройство компенсации, параметрическая нестабильность, параметрические возмущения, обратная задача динамики, эталонная модель, ненаблюдаемая координата, конвекционный нагрев, температура в центре заготовки, изменение теплофизических характеристик объекта.
Введение
Современная наука имеет мощный математический аппарат для описания объектов с распределенными параметрами. Существуют различные подходы к задачам синтеза систем управления такими объектами.
Большинство подходов обладает общим недостатком - в них не учитывается параметрическая нестабильность объектов. Все модели объектов и алгоритмы управления рассчитываются с учетом допущений о постоянных или кусочно-постоянных физических параметрах объекта. Например, в задаче нагрева массивных заготовок объект рассматривается как распределенный и необходимо учитывать
Тычинин Александр Викторович (к.т.н.), начальник отдела разработки АСУ ТП. Тычинина Юлия Александровна (к.т.н.), доцент кафедры «Автоматика и управление в технических системах».
Рагазин Дмитрий Александрович (к.т.н.), доцент кафедры «Автоматика и управление в технических системах». 58
распределение температуры по его толщине. Подобные допущения о постоянстве теплофизических параметров вносят погрешность в математическое описание объекта управления.
При разработке систем управления таким объектом прибегают, как правило, к аналитическим методам, которые также базируются на допущении постоянства теплофизических параметров.
Еще одной проблемой при построении таких систем может быть необходимость управления ненаблюдаемой координатой объекта - например, температурой в середине объекта, измерить которую не представляется возможным. Приходится определять ее косвенно, прибегая к аналитическим методам.
Проблема инвариантности в системах управления широко рассматривается в современной науке, например в [1-4].
Применительно к объектам с сосредоточенными параметрами известна задача синтеза системы управления методами обратной динамики и задача компенсации динамических свойств объекта нахождением близкой к обратной передаточной функции с помощью цепных (периодических) структур.
Применительно к задачам управления метод обратной динамики подразумевает определение закона управления путем нахождения дифференциального уравнения, описывающего поведение управляемой величины по заданным требованиям к показателям качества, последующего вычисления по этому уравнению старшей производной и подстановки ее вместо старшей производной в уравнение объекта.
Такой метод был предложен Л.М. Бойчуком [5] для определения структуры при синтезе нелинейных систем управления. Широкому распространению метода способствовали работы П.Д. Крутько [6, 7].
Задача нахождения обратной передаточной функции объекта относится к классу «некорректных» и не имеет точного решения. Проблема состоит в необходимости компенсации динамических свойств объекта и, как следствие, в выполнении операций точного дифференцирования.
Идеальное дифференцирующее звено реализовать невозможно, поэтому, как правило, прибегают к различным приближенным способам реализации оператора дифференцирования. Большинство из них состоит в использовании метода малого параметра, где наличие малого коэффициента позволяет задать системе некоторый запас устойчивости, снизив при этом точность реализации.
Задачам нахождения обратных передаточных функций и решению некорректных задач посвящены работы Б.М. Менского, А.Н.Тихонова [8-12].
Для решения задачи компенсации динамических свойств объекта возможно использовать цепную схему, где последовательно соединены одинаковые звенья и каждое последующее звено все больше приближает структуру к обратной функции объекта [13-19].
Постановка задачи
Для синтеза системы управления процессом нагрева массивной металлической заготовки относительно ненаблюдаемой координаты (нагрев центральной точки) предлагается последовательно решать следующие задачи:
1. Разработка математической модели распределенного объекта с граничным управлением первого рода.
2. Разработка динамических структур, близких к обратной передаточной функции объекта.
3. Определение эталонной траектории изменения регулируемого параметра (температура в центре заготовки) и на основе обратной передаточной функции определение управляющего воздействия (температура на поверхности заготовки), обеспечивающего требуемый характер поведения ненаблюдаемой фазовой координаты.
4. Математическое описание объекта управления в граничных условиях третьего рода (конвективный нагрев в газовой печи).
5. Структурно-параметрический синтез системы граничного управления с обратной связью по температуре поверхности заготовки методом обратной динамики.
В результате получим систему управления нагревом ненаблюдаемой координаты объекта с обратной связью по температуре поверхности заготовки и температурой печи в качестве управляющего воздействия.
Как будет показано ниже, такая система будет обладать слабой чувствительностью к параметрическим возмущениям объекта управления.
Структура предлагаемого подхода представлена на рис. 1.
Синтезированная система управления Температура печи Объект управления с ГУ третьего рода темпера!ура поверхности ^
* Объект управления с ГУ первого рода Температура внутренней координаты
^ Ц
Рис. 1. Структура системы управления внутренней координатой объекта
В качестве объекта управления рассматривается массивная неограниченная пластина из алюминиевого сплава Д16, нагреваемая по толщине. В таблице приводятся принятые в расчетах характеристики нагреваемой заготовки.
Характеристики нагреваемых заготовок сплава Д16 (АА2024)
Параметр Значение
2Я, толщина заготовки, мм 250
X, коэффициент теплопроводности, Вт/м 0С 130
р, плотность, кг/м3 2800
с, теплоемкость, Дж/кглС 922
а, коэффициент температуропроводности, м2/с 4910-6
Объект управления описывается сначала в граничных условиях первого рода (управление температурой внутренней координаты по температуре поверхности). Затем осуществляется переход к управлению граничной температурой объекта с управляющим воздействием по температуре конвекционной печи. Для этого объект описывается в граничных условиях третьего рода.
Решение поставленной задачи
Решением первой краевой задачи (1) построена модель объекта с управлением по температуре поверхности пластины [20]:
dQ(x, t) = д d2Q(x,/), dt dt2 '
Q(x,0) = Qo(x), 0<x<R;
dQ(0, t) n
' ' = 0, t > 0;
dt ' ;
Q(R, t) = Qn(t), t > 0.
В (1) a=const - заданный постоянный коэффициент.
Краевая задача (1) может быть представлена в стандартной форме с однородными граничными и нулевыми начальными условиями:
dQ( x,t) d 2Q(x,t) . .. - -
' = a—' + a(x,t), t >0, 0<x<R dt dx2
Q( x,0) = 0, 0 < x < R (2)
6Q(0, t) = 6Q(R, t) = 0 t > 0 dt dt '
где a>(x, t) - стандартизирующая функция [20]:
t) = ab'(x — R) ■ Qn (t) + Q0 (x) ■ S(t). (3)
Здесь S(t) и S(x—R) - дельта-функции временного и пространственного аргументов соответственно, сосредоточенные в точках t=0 и x=R, и S'(x — R) - производная по х от S(x—R).
Переход к изображениям по Лапласу в уравнениях (2) приводит к следующему выражению для передаточной функции объекта управления (1) с граничным управляющим воздействием Qn (t) в форме ее разложения в бесконечный ряд по собственным функциям первой краевой задачи [20]:
4R2
Wx (x p) = Qm = % £(—1)" (2, — 1)CoI l^^^, (4)
Qn (p) R2 n=1 I 2R J 4R2 +1
an2(2n — 1)2
где ((х, р) и (2П (р) - изображения (х,?) и Qп (0 .
Модель объекта построена в соответствии с передаточной функцией (4) с учетом характеристик, приведенных в таблице. При построении модели (рис. 2) ограничились пятью слагаемыми бесконечного ряда.
Задачу компенсации динамических свойств объекта и получения приближенной обратной передаточной функции выполним с помощью цепной структуры [19,
22] (рис. 3), где последовательно соединены одинаковые звенья и каждое последующее звено все больше приближает структуру к обратной функции объекта [1419, 21].
Рис. 2. Модель объекта управления для нагрева центра пластины при граничных условиях первого рода
Рис. 3. Система компенсации динамических свойств объекта
На рис. 3 ^мод(р) — передаточная функция модели объекта управления,
а ^к (р) — передаточная функция корректирующего звена, описываемая уравнением
Щк (р) =
1-С + Щмод(р)
(5)
Коэффициент С варьируется в диапазоне [0... 1). При С=1 система оказывается физически нереализуемой. В данном подходе необходимо подобрать коэффициент С таким образом, чтобы система была устойчивой и реализуемой (например, в части ограничений на величину управляющего воздействия). Увеличением количества ячеек периодической структуры в устройстве компенсации добиваются приемлемой точности. При этом чем ближе значение коэффициента С к 1, тем выше скорость сходимости.
Несложно показать, что при бесконечном увеличении количества ячеек периодической структуры и при значениях С, принадлежащих указанному диапазону, передаточная функция структуры стремится к обратной передаточной функции модели объекта:
Щ( р)=I (р) Щ (р),
\( р)=I (р) • (1 - К (р) Щмод( р)),
и1( р) = А1( р) К (р),
Л2 (р) = А1 (р) • (1 - Жк (р) • Кмод( р)),
и2( р) = Л2( р) К (р),
лп (р) = Лп_1 (р) • (1 - К (р) • Кмод( р)),
ип (р) = Лп (р) Щ (р).
(6)
Из (6) получаем:
и (р) = Ёи (р),
1=0
и (р) = Щк (р) • I (р) • (1 - щ (р) • Щмод (р)У,
п
и(р) = Щк(р) • 1(р) • (1 + Х(1 -К(р) • Кмод(р))).
(7)
1=1
Из (5) следует, что
Пт/Кк (р) • Кмод(р)) = 1,
(8)
где С принадлежит диапазону [0.1).
Тогда в (7) при количестве ячеек имеем сумму бесконечной геометрической прогрессии, знаменатель которой по абсолютной величине меньше единицы:
\1 - Щк (р) Кмод(р)\ < 1, приС < 1.
(9)
Известно, что сумма бесконечно убывающей геометрической прогрессии равна первому члену этой прогрессии, деленному на разность между единицей и знаменателем этой прогрессии.
Таким образом, при бесконечном количестве ячеек структуры, представленной на рис. 3, она с абсолютной точностью реализует обратную передаточную функцию объекта управления при любом значении коэффициента С 6 [0,1):
и(р) = Жк (р)- /(р) 1 = /(р) (10)
(р) 'Кмод (р) Кмод (р)
На рис. 3 устройство компенсации включено последовательно с моделью объекта управления. Такое соединение при идеальной реализации компенсирующего устройства имеет единичную передаточную функцию и должно повторить входное воздействие / ( р).
В качестве объекта управления рассматривалась модель поведения регулируемой температуры в центре пластины:
ЖЭт(р) =-300-. (11)
5000(р2 + 330р + 1
Построен динамический компенсатор центра пластины при нагреве в граничных условиях 1 -го рода (рис. 4) и найдена необходимая траектория изменения температуры поверхности заготовки.
Значение коэффициента С=0.7 выбрано из соображений устойчивости системы.
Рис. 4. Динамический компенсатор температуры центра пластины в процессе нагрева при граничных условиях 1 -го рода
На рис. 5 приведена переходная характеристика (зад(0,0 эталонной модели полученной реакции ((0,1) на выходе модели объекта; управляющее воздействие (П (О - температура поверхности пластины, а также относительная погрешность, полученная при значении коэффициента С=0,7 и количестве ячеек 1 и 3.
Получена передаточная функция изменения граничной температуры нагреваемой заготовки с учетом заданной эталонной траектории температуры в центральной точке объекта (12). Она соответствует модели системы, представленной на рис. 4 с моделью распределенного объекта управления, представленного на рис. 2.
12р4 +21ръ +25р2 + 13р + 38
К эт.гр(р)
эшгр 50р6 + 200р5 + 313рл + 258ръ + 122р2 + 32.1р + 0.107' (12)
в
д
г
Рис. 5. Переходная характеристика эталонной модели (а), синтезированной системы (б), найденный закон изменения граничной температуры (в), погрешность отработки сигнала эталонной модели при 1 ячейке компенсационной структуры (г), при 3 ячейках компенсационной структуры (д)
д((х, г)=а д( X, г)
дг
дг2
0 < х < Я, г > 0;
((х,0) ((х), 0 < х < Я, г = 0; (0, г)
дг
= 0, г > 0;
ЛЩМ. + а.((Я,г) = а.(с(г), г>0. дг
Решением третьей краевой задачи (13) построена модель объекта [6] для температуры поверхности той же самой заготовки с управлением по температуре печи (рис. 6).
В (13) а=еот1 - заданный постоянный коэффициент.
Стандартная форма уравнений объекта (13) опять представляется в виде (2), где теперь стандартизирующая функция отличается от (3) и записывается следующим образом [20]:
а
<х, г) =—8(£- Я) - ((с (г)+((х) - Ъ(г). су
(14)
Передаточная функция объекта опять представляется в форме бесконечного ряда [20]:
Я2
а у
уЯ"кцп+5т(цп)Со8(цп)
2^„С08(Ц„ )
-008
V Я у
Я1
-. (15)
аЛп
р + 1
При построении модели, как и в прошлый раз, ограничились пятью слагаемыми. Также учтены характеристики, приведенные в таблице.
Рис. 6. Модель объекта управления для температуры поверхности пластины при граничных условиях третьего рода
2
а
2
Для синтеза системы управления, адаптивной к параметрическим возмущениям объекта, воспользуемся методами обратной динамики [5-7, 13].
Задача сводится к минимизации квадратичного отклонения значений производных для выходов эталонной модели и объекта управления соответственно.
Для этого предполагается минимизировать функционал (16) по градиентной схеме (17).
G(Qc(t))=\(y(t) - Q(xh t)f; (16)
dQM=AdGQcl, A=const (17)
dt dQc
Здесь y(t) - сигнал на выходе эталонной модели. За эталонную модель принимается передаточная функция (12).
Из (16) и (17) получаем:
Qc(t) = k(y(t)-Q{xl,t)), к = const. (18)
На основе уравнения (18) выполним структурный синтез системы управления. Для этого проинтегрируем почленно уравнение эталонной модели (12), выделив ее выход y(t) . Подставим его в (18) и выполним замену y(t) = Q(xi, t), обеспечив обратную связь по управляемой величине.
Указанные действия выполняются непосредственно при построении структуры системы управления (рис. 7).
На рис. 7: Q(0,p) - температура центра пластины; Q(R,p)- температура ее
поверхности; W^ (R, p) - передаточная функция модели объекта управления в условиях третьей краевой задачи; Qc (p) - граничное управляющее воздействие по температуре печи; WfY1(0, p) - передаточная функция модели объекта управления от температуры поверхности к температуре центра в условиях первой краевой задачи.
Следует обратить внимание, что коэффициенты синтезированной структуры соответствуют коэффициентам эталонной модели и не зависят от параметров самого объекта. Благодаря этому система получает свойство слабой чувствительности к параметрическим возмущениям.
Для проверки работоспособности системы в условиях параметрической нестабильности объект управления подменяли аналогичным объектом, но с другими параметрами (рис. 8а).
Требуемая переходная характеристика температуры поверхности заготовки и полученная при замене объекта моделью с другими теплофизическими параметрами представлены на рис. 86 и 8в.
Представленные переходные характеристики подтверждают слабую чувствительности синтезированной системы к параметрическим возмущениям объекта управления.
На полученных моделях проверялась способность системы отрабатывать другие виды возмущений. Для примера рассмотрено возмущение по входу объекта (помехи по температуре печи), которое может быть вызвано притоком воздуха разной температуры в печь извне.
12 21 13 38
-пя
1,107 32 123
мШ
Щ^я.р)-
QtR.Pi
6(0,0
а
Рис. 7. Структурная схема системы управления по температуре в центре заготовки (а), заданная переходная характеристика для температуры в центре пластины (б), полученная переходная характеристика (в), относительная погрешность синтезированной
системы (г)
На рис. 9 представлено случайное воздействие на температуру печи и относительное отклонение от эталонной траектории температуры ненаблюдаемой координаты распределенного объекта (нагреваемой заготовки).
б
Рис. 8. Параметрические возмущения объекта управления (а), заданная переходная характеристика для температуры на поверхности пластины (б), полученная переходная характеристика при параметрических возмущениях (в)
а
в
Рис. 9. Случайное воздействие на входе объекта управления (а), относительное отклонение от эталонной траектории (б)
Выводы
В работе выполнен синтез системы управления ненаблюдаемой координатой объекта с распределенными параметрами. Синтезированная система обладает свойством адаптивности к параметрическим возмущениям объекта управления. Рассматривался процесс нагрева массивной неограниченной пластины по толщине. В качестве материала для нагреваемой заготовки использован сплав Д16.
Для модели объекта управления в граничных условиях первого рода с помощью периодических структур получена близкая к обратной передаточная функция и найдено такое поведение граничной температуры, которое обеспечивает требуемый нагрев ненаблюдаемой центральной точки объекта.
Это поведение граничной температуры принято за эталонную модель, на основе которой решена обратная задача динамики для того же объекта, но в граничных условиях третьего рода (управление температурой поверхности с помощью температуры конвекционной печи).
В ходе исследования выявлены также основные ограничения предложенного подхода:
- порядок эталонной модели должен быть не ниже порядка объекта управления;
- инерционность эталонной модели не должна быть ниже инерционности объекта управления.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Нестеров В.Н. Моделирование инвариантных и квазиинвариантных систем на основе принципа многоканальности // Прикладная физика и математика. - 2020. - № 5. - С. 32-35.
2. Аполонский В.В., Копытова Л.Г., Тарарыкин С.В. Разработка и исследование селективно-инвариантных электромеханических систем с адаптацией регуляторов к изменениям уровня скорости // Известия Российской академии наук. Сер. Теория и системы управления. - 2020. - № 5.
- С. 28-43.
3. Хиен Ву Ань, Ягодкина Т.В. Синтез инвариантной системы адаптивного управления на базе следящей системы // Фундаментальные исследования. - М., 2016. - № 6-1. - С. 52-57.
4. Харьков В.П. Адаптивное управление динамической системой на основе обратных задач динамики // Труды международного симпозиума «Надежность и качество». - Пенза, 2016. - Т. 1. -С. 176-177.
5. Бойчук Л.М. Метод структурного синтеза нелинейных систем автоматического управления. -М.: Энергия, 1971.
6. КрутькоП.Д. Обратные задачи динамики управляемых систем. Линейные модели. - М.: Наука, 1987. - 304 с.
7. Крутько П.Д. Обратные задачи динамики в теории автоматического управления. - М.: Машиностроение, 2004. - 576 с.
8. Менский Б.М. Принцип инвариантности в автоматическом регулировании. - М.: Машиностроение, 1972. - 235 с.
9. Тихонов А.Н., Самарский А.А. Уравнения математической физики. - М.: Наука, 1966. - 724 с.
10. Тихонов А.Н., Арсенин В.Я. Методы решения некорректных задач. - М.: Наука, 1979.
11. ТихоновА.Н., Гончаровский А.В., Степанов В.В., Ягола А.Г. Численные методы решения некорректных задач. - М.: Наука, 1990.
12. Тихонов А.Н., КальнерВ.Д., Гласко В.Б. Математическое моделирование технологических процессов и метод обратных задач в машиностроении. - М.: Машиностроение, 1990.
13. КимД.П. Теория автоматического управления. - М.: ФИЗМАТЛИТ, 2004. - 464 с.
14. Тычинин А.В. Методы обратных задач динамики в задачах синтеза систем управления распределенным объектом // Вестник Самарского государственного технического университета. Сер. : Технические науки. - 2005. - № 39. - С. 21-25.
15. Тычинин А.В. Адаптивность в задачах синтеза систем управления распределенным объектом // Современные техника и технологии СТТ'2006: Труды XII Международной научно-практической конференции студентов, аспирантов и молодых ученых. - Томск, 2006. - С. 188-189.
16. Тычинин А.В. Синтез систем управления распределенным объектом управления решением обратной задачи динамики // Математическое моделирование и краевые задачи: Труды V Всероссийской научной конференции с международным участием. - Самара, 2008. - С. 128130.
17. Тычинин А.В., Тычинина Ю.А. Структурно-параметрический синтез системы управления объектом с распределенными параметрами // Вестник Самарского государственного технического университета. Сер.: Технические науки. - 2009. - № 1(23). - С. 74-79.
18. Тян В.К. Теория периодических структур в некорректных задачах синтеза инвариантных и автономных систем управления // Вестник Самарского государственного технического университета. Сер.: Технические науки. - 2006. - Вып. 41. - С. 47-54.
19. Тян В.К. Структурное представление обратного оператора в банаховом пространстве // Вестник Самарского государственного технического университета. Сер.: Технические науки. - 2007. -№ 1(14). - С. 197-199.
20. Рапопорт Э.Я. Структурное моделирование объектов и систем управления с распределенными параметрами: Учеб. пособие. - М.: Высш.шк., 2003. - 299 с.
21. Рапопорт Э.Я., Тян В.К. Достижение заданной инвариантности в стохастических системах комбинированного управления // Куйбышевский политехнический институт. Деп. В ВИНИТИ 20.06.89, № 4089-В89.
22. Рапопорт Э.Я. Анализ и синтез систем автоматического управления с распределенными параметрами: Учеб. пособие. - М.: Высшая школа, 2005. - 292 с.
Статья поступила в редакцию 10 декабря 2020 г.
STRUCTURAL AND PARAMETRIC SYNTHESIS OF THE CONTROL SYSTEM OF THE UNOBSERVED OUTPUT FOR THE OBJECT WITH DISTRIBUTED PARAMETERS
A.A. Tychinin1, Y.A. Tychinina2, D.A. Ragazin2
1 «TM-service»
157, 6th proseka, Samara, 443029, Russian Federation
2 Samara State Technical University
244, Molodogvardeyskaya st., Samara, 443100, Russian Federation
Abstract. Method of synthesis of a control system for the object with the distributed parametres, providing compensation of dynamic properties of the object and giving to it ofprop-erties of standard model of a set. The system possesses the property of an adaptability to parametrical interference of the object. Control along the unobservable coordinate of the object is implemented by synthesizing a chain structure that implements a transfer function close to the inverse transfer function of the object. Weak sensitivity to parametric instability is provided by solving the inverse problem of dynamics. As an example, the heating of a massive billet in a convection oven is considered.
Keywords: object with distributed parameters, boundary conditions 1 type, boundary conditions 3 type, dynamic compensation, object inverse transfer function, periodic structures, compensation device, parametric instability, parametric noise, inverse problem of dynamics, reference model, unobservable coordinate, convection heating, temperature at the center of the workpiece, change in the thermophysical characteristics of the object.
Aleksandr A. Tychinin (Ph.D. (Techn.)), Head of Department. Yulia A. Tychinina (Ph.D. (Techn.)), Associate Professor. Dmitry A. Ragazin (Ph.D. (Techn.)), Associate Professor.
REFERENCES
1. Nesterov V. N. Modeling of invariant and quasi-invariant systems based on the multichannel principle. Applied Physics and Mathematics. 2020. Vol. 5. Pp. 32-35. (InRussian).
2. Apolonsky V.V., KopylovaL.G., TtararykinS.V. Development and research of selective-invariant electromechanical systems with adaptation of regulators to change of speed level. Bulletin of the russian academy of sciences. Ser. control theory and systems. 2020. Vol. 5. Pp. 28-43. (InRussian).
3. Hien Vu Anh, Yagodkina T. V. Synthesis adaptive modal regulator for invariant system on the basic of tracking system. Fundamental research. 2016. Vol. 6-1. Pp. 52-57. (InRussian).
4. Kharkov V.P. Adaptive control of a dynamic system based on inverse problems of dynamics. Proceedings of the International Symposium "Reliability and Quality", Penza, 2016, Vol. 2. Pp. 176-177. (In Russian).
5. BoychukL.M. Method of structural synthesis of nonlinear automatic control systems. M.: Energiya, 1971.
6. KrutkoP.D. Inverse problems of the dynamics of controlled systems. Linear models. Moscow: Nauka, 1987. 304 p.
7. Krutko P.D. Inverse problems of dynamics in the theory of automatic control. Mashinostroenie, 2004. 576 p.
8. MenskyB.M. The principle of invariance in automatic control. M.: Mechanical engineering, 1972, 235 p.
9. TikhonovA.N., SamarskyA.A. Equations of mathematical physics. M.: Nauka, 1966. 724 p.
10. Tikhonov A.N., Arsenin V.Ya. Methods for Solving Ill-Posed Problems. M.: Nauka, 1979.
11. TikhonovA.N., GoncharovskyA.V., Stepanov V.V., Yagola A.G. Numerical Methods for Solving Ill-Posed Problems. M.: Nauka, 1990.
12. Tikhonov A.N., Kalner V.D., Glasko V.B. Mathematical modeling of technological processes and the method of inverse problems in mechanical engineering. M.: Mechanical Engineering, 1990.
13. Kim D.P. Automatic control theory. M.: FIZMATLIT, 2004. 464 p.
14. Tychinin A. V. Methods of inverse problems of dynamics in problems of synthesis of control systems for a distributed object. - steel. Vestnik Samgtu. Ser. Technical Sciences. 2005. Vol. 39. Pp. 21-25. (In Russian).
15. Tychinin A. V. Adaptability in problems of synthesis of control systems for a distributed object. Modern techniques and technologies CTT'2006 // proceedings of the XII International scientific-practical conference of students, graduate students and young scientists, Tomsk 2006. Pp. 188-189.
16. Tychinin A. V. Synthesis of control systems for a distributed control object by solving the inverse problem of dynamics. Proceedings of the fifth All-Russian scientific conference with international participation "Mathematical modeling and boundary value problems". Samara, 2008. Pp. 128-130.
17. Tychinin A. V., Tychinina Yu.A. Structural and Parametric Synthesis of a Distributed Parameter Object Control System. Vestnik Samgtu. Ser. Technical Sciences. 2009. Vol. 1(23). Pp. 74-79. (In Russian).
18. Tyan V.K. Theory of periodic structures in ill-posed problems of synthesis of invariant and autonomous control systems. Vestnik Samgtu. Ser. Technical Sciences. 2006. Vol. 41. Pp. 47-54. (In Russian).
19. Tyan V.K. Structural representation of the inverse operator in a Banach space. Vestnik Samgtu. Ser. Technical Sciences. Ser. Physics and mathematics. 2007. Vol. 1(14). Pp. 197-199.
20. RapoportE.Ya. Structural modeling of objects and control systems with distributed parameters: Textbook. allowance. M.: Vyssh. shk., 2003. 299 p.
21. Rapoport E.Ya., Tyan V.K. Achieving a given invariance in stochastic combined control systems. Kuibyshev Polytechnic Institute. Dep. AT VINITI 20.06.89, No. 4089-B89.
22. RapoportE. Ya. Analysis and synthesis of automatic control systems with distributed parameters: Textbook. allowance. M.: Vyssh. shk., 2005. 292 p.
