CC BY
127
Заметим, что при обоих значениях к1 в системе гарантируется устойчивость и требуемая точность функционирования при реакции на линейно-нарастающее входное воздействие g(t), рис. 4. Однако регулятор с коэффициентами, определенными из условия П1*=2,29 обеспечивает большее быстродействие, а также меньшее перерегулирование в системе по сравнению с регулятором, обеспечивающим 77 *,=0,325, что видно из рис. 5.
Выводы
1. Представлена методика параметрического синтеза ПИ-регулятора, обеспечивающего квазимакси-мальную степень устойчивости и заданную добротность системы автоматического управления.
2. Устойчивость системы автоматического управления достигается путем расположения ее полюсов левее максимальной оценки снизу степени устойчивости системы, что обеспечивает максимальное быстродействие и снижает колебательность.
3. Алгоритм синтеза использует коэффициентные оценки показателей качества системы, на основе которых формируются упрощенные алгебраические соотношения между оценкой снизу степени устойчивости системы и коэффициентами характеристического полинома.
4. Полученные результаты рассмотрены и апробированы на численном примере.
СПИСОК ЛИТЕРАТУРЫ
1. Петров Б.Н., Соколов Н.И., Липатов А.В. и др. Системы автоматического управления объектами с переменными параметрами: Инженерные методы анализа и синтеза. - М.: Машиностроение, 1986. - 256 с.: ил.
2. Волков А.Н., Загашвили Ю.В. Метод синтеза систем автоматического управления с максимальной степенью устойчивости при наличии ограничений // Известия РАН. Сер. Теория и системы управления. - 1997. - № 3. - С. 12-19.
3. Шубладзе А.М. Способы синтеза систем управления максимальной степени устойчивости // Автоматика и телемеханика. - 1980. - №1. - С. 28-37.
4. Шубладзе А.М. Методика расчета оптимальных по степени устойчивости ПИ-законов. I // Автоматика и телемеханика. -1987. - № 4. - С. 16-25.
5. Татаринов А.В., Цирлин А.М. Задачи математического программирования, содержащие комплексные переменные, и предельная степень устойчивости линейных динамических систем // Известия РАН. Сер. Теория и системы управления. -1995. - № 1. - С. 28-33.
6. Воронина Н.О., Татаринов А.В., Цирлин А.М. Предельная степень апериодической устойчивости и соответствующие ей настройки для типовых систем регулирования // Известия вузов. Сер. Приборостроение. - 1989. - № 3. - С. 2б-32.
Поступила 28.03.2012 г.
УДК 681.5.015
ИДЕНТИФИКАЦИЯ ОБЪЕКТОВ УПРАВЛЕНИЯ В ФОРМЕ ДИСКРЕТНЫХ ПЕРЕДАТОЧНЫХ ФУНКЦИЙ НА ОСНОВЕ ВЕЩЕСТВЕННОГО ИНТЕРПОЛЯЦИОННОГО МЕТОДА
В.А. Рудницкий, А.С. Алексеев, В.В. Курганкин
Томский политехнический университет E-mail: [email protected]
Для решения задачи параметрической идентификации линеаризуемых объектов использован вещественный интерполяционный метод. На основе аппарата чисел обусловленности рассмотрена возможность повышения точности определения структуры модели в виде дискретной передаточной функции.
Ключевые слова:
Идентификация, объект управления, числа обусловленности, вещественный интерполяционный метод.
Key words:
Identification, controlled object, condition numbers, real interpolation method.
Введение
В современной технике управления широко применяют различные цифровые устройства. Прогресс в развитии микроэлектроники привел к улучшению основных качественных и количественных показателей микроконтроллеров, что сделало возможным эффективно решать задачи адаптации,
оптимизации и диагностики систем управления. В то же время успешное их решение невозможно без получения и использования адекватного математического описания объектов и систем.
Методы идентификации цифровых объектов во многих случаях удобно классифицировать по области описания математической модели:
в области времени [1, 2] ив области изображений, в частности, в области ^-изображений [4] и частотной области [3]. Методы временной группы обладают сравнительно невысокой помехоустойчивостью и в общем случае требуют значительных вычислительных затрат. В случае представления модели системы с помощью частотных характеристик число операций оказывается сравнительно высоким из-за необходимости вычисления амплитудных и фазовых частотных характеристик объектов по довольно сложному алгоритму. Наиболее широкое распространение получило использование аппарата ^-преобразования для идентификации моделей цифровых объектов.
В связи с этим целью настоящего исследования является разработка эффективного алгоритма идентификации объектов управления, поведение которых можно достаточно адекватно описать с помощью линейных математических моделей, представленных в дискретном времени.
Задачу идентификации линеаризуемых объектов управления в форме дискретных передаточных функций (ПФ) можно решать двумя путями - на основе непрерывной и цифровой модели. В первом случае используется косвенный путь: по экспериментальным отсчетам входного X и выходного У сигналов находится ПФ W0(p), которая в последующем переводится в класс дискретных. Во втором случае определяется непосредственно дискретная ПФ W0(z) идентифицируемого объекта, структурная схема которого представлена на рис. 1.
Рис. 1. Модель объекта управления
При использовании любого из этих подходов задача идентификации объекта формулируется следующим образом. Известны входной Х(кТ0) и выходной У(кТ0) сигналы, где Т0 - период квантования; к - номер отсчета. Необходимо определить параметры ПФ вида
У(^) Ь zm + Ь ,хт-1 +... + Ь0 W0(z) = —— = Ьт--------т-1--;^, п > т. (1)
0( ) X(х) anzn + ап_^п-1 +... +1
Обычно структурные показатели т, п известны или заданы, что соответствует параметрической постановке задачи.
Общей проблемой различных алгоритмов идентификации остается довольно сложная процедура получения математических моделей исследуемых объектов, требующая для своей реализации существенных вычислительных затрат. По этой причине разработка более эффективных численных способов решения задачи идентификации дискретных объектов и систем остается актуальной. Определенные перспективы в реализации этих возможностей дает применение вещественного интерполяционного метода (ВИМ).
1. Вещественный интерполяционный метод
В основе использования ВИМ для решения задачи идентификации цифровых объектов лежит вещественное дискретное преобразование
да
^*(8) = £ /(кТ0)е-кТ°8, с < 8 < да, с > 0, (2)
к=0
в котором функции-оригиналу/(кТ0) соответствует изображение ^‘(8) [5]. В частном случае, когда рассматриваются устойчивые объекты и используются их импульсные переходные характеристики, параметр с может быть принят равным нулю [5, 6].
Для использования ВИМ в реальных задачах рекомендуется переходить от преобразования (2) к более простому в вычислительном отношении о-преобразованию. С этой целью достаточно выполнить подстановку о=ехр (Т08) в выражении (2), что приводит к соотношению
Яо)=£/(кТо)о-‘, о є [С, да], С > 1.
(3)
В формуле (3) переменная и не может принимать значения меньшие единицы, что соответствует области устойчивости на плоскости ^ Основные свойства дискретного и-преобразования во многом близки свойствам z-преобразования [5].
Понятие вещественного изображения распространяется на описание линейных импульсных систем и протекающих в них процессов. Если на некоторую линейную импульсную систему действует входной сигнал Х(кТ0) и известна реакция системы У(кТ0), то на основании (2) и (3) можно записать уравнение «вход-выход объекта» в вещественной области
У (и) = W (и) • X (и). (4)
Для выполнения численных действий над вещественными изображениями необходимо перейти от непрерывных функций Ди) к их дискретным аналогам в виде численных характеристик. Под ними, как и в случае дискретного вещественного преобразования, понимают совокупность значений
{Р(и)}п = {^(и1). Р^г),..., Яи,)}. П = т + п + 1, (5)
заданных на сетке 1<и1<и2<_<и1). Значения узлов интерполирования и, /еЩ, могут назначаться по известным рекомендациям, изложенным в [5, 6].
Рассмотрим приложение аппарата вещественных функций-изображений и численных характеристик к решению задачи идентификации. Положим в (4): Х(и) - единичный входной сигнал; У(и) - реакция объекта на это входное воздействие. По формуле (3) для принятого входного сигнала можно найти изображение Х(и)=и/(и-1), где ие[1,да). Выходной сигнал реакции объекта У(и) по определению будет представлять собой переходную функцию к(и). В этом случае при нулевых начальных условиях ур. (4) примет вид
к (о) = Ж (о) ———, о є [С, да].
о-1
(6)
Положим далее, что рассматривается класс объектов управления с самовыравниванием, что
соответствует абсолютному большинству практических задач. Тогда в соотношении (6) можно принять С=1 [7]. Приведенные сведения позволяют перейти непосредственно к решению задачи параметрической идентификации.
2. Алгоритм параметрической идентификации цифровых объектов
Задача заключается в определении коэффициентов вещественной ПФ вида
э =
.и т + ъ и т- +... + 1
пи п + апп" +... +1
(7)
нала
Чи) = Ёк(кТо) и-*> и е[1>да].
(8)
Параметр N в (8) определяется временем наблюдения переходного процесса. Далее на основании (6) определяется численной характеристикой искомой ПФ:
и- 1 —
Ш (и) = —-к(ц), г = 1,п, п = т + п +1.
и
Теперь имеются все необходимые данные для вычисления коэффициентов функции (7). Задача свелась к решению СЛАУ
Ш (и,) = т°і + "т -1^‘
ит-+...+-
’о иі -1
°пи" + ап-М” + ... + 1 Ц'
г = 1, п , П = т + п +1. (9)
Систему уравнений (9) при решении практических задач целесообразно представлять в матричной форме
Б ■ К = Ш, (10)
где
К = [Ьт , Ьт-1, ..., Ь0, ап , ап-1, ..., а1] ,
Ш = [Ш(01)01, Ш(о2)о2,..., Ш(о„)о„]г,
иГ‘ -и1т и1т -и1т-1 .. и,-1 -Ш(и1)ип+1 . -Ш(и1)и1
и"1 -и2т и2т -ит- .. и2-1 -Ш (и2)иП,+1 . -Ш(и2)и2
П1 и и," 1 -и"-1 .. и,_. -1 -Ш(и,_ 1)иП+11 . -Пи,-.и.
иГ-и,т и;-и;- .. и,_! -1 -Щи, )и"+1 . -Ш(и,)и,
по переходной характеристике, представленной своим вещественным изображением Н(и). Для решения задачи используются численная характеристика переходной функции [к(и)}п и входного сиг-
I . Размерность п характеристик (5)
и - Чп
принимается равной числу неизвестных коэффициентов ПФ (7). В этом случае равенство (6) можно развернуть в систему п линейных алгебраических уравнений (СЛАУ), которая имеет решение и оно единственно [5].
Техника формирования СЛАУ включает в себя получение численной характеристики единичного ступенчатого сигнала и переходной функции
{к(и)}„. Элементы X(и,) = , и,. > 1 находятся
и -1
расчетным путем, а получение значения {к(и)}п можно осуществить по формуле (3), зная последовательность к(кТ0), где к=0,1,2,...,Ж:
Поиск решения может осуществляться известными методами, т. к. СЛАУ не содержит каких-либо особенностей. Трудности могут возникнуть на этапе формирования СЛАУ, в частности, с получением первого уравнения - для и=и1. Рассмотрим особенности этого уравнения. Ранее принято ограничение и1>1, как условие существования и-преобразования и сходимости ряда (8), что однако не указывает на конкретное значение и1. Так, в случае неудачного выбора значения и1 матрица Б может оказаться плохо обусловленной. Для устранения этого недостатка существует подход, который снимает трудности выбора первого узла интерполирования [5].
С этой целью рассмотрим взаимную связь конечных значений оригиналов с их изображениями. Для них справедливо равенство
и — 1
Иш / (кТо) = Иш------/ (о),
к^да и^1 и
если эти пределы существуют. Узел интерполирования и1 и значение Дц) обеспечивают сохранение информации о поведении функции ЛкТ0) преимущественно на бесконечности. Эту информацию можно задать также и в области изображений по ПФ Ж(и) (7). Поэтому можно записать
и — 1
Иш к(кТ0) = Иш -— к(о) =
к^ю и^1 и
= Иш —1Ш (и) — = Ш(1).
и^1 и и-1
В большинстве случаев значение к(нТ0) в установившемся режиме известно из исходных данных, поэтому преобразуем предыдущее выражение к виду
Н(кГо)к
ап + ап-1 + ... + 1
(12)
обеспечивающему совпадение установившихся значений выходов идентифицируемого объекта и его модели. Включение последнего уравнения в состав СЛАУ (9) снимает проблему выбора первого узла интерполирования и1. Узел и2 не оказывает критического влияния на обусловленность матрицы СЛАУ даже при равномерном распределении узлов интерполирования [5].
Для моделей объектов с самовыравниванием в форме (7), полученных при начальных нулевых условиях, возможно осуществить переход к г-форме Щ>(г) (1) путем формальной замены о^г., что является важной практической чертой ВИМ и вещественного дискретного преобразования.
Последним этапом в решении задачи идентификации цифрового объекта или системы является
оценивание погрешности полученного решения и, если это необходимо, его коррекция, обычно сводящаяся к минимизации погрешности. В этой связи необходимо сформировать критерий близости в виде
Ah = max |Ah( kT0)| =
- max
k
100 • ho(kTo) - h„ (kTo)
ho(kTo)
(11)
IK
Оно может быть легко получено, если применить общую теорию обусловленности к частному случаю СЛАУ.
Число обусловленности Эйрда-Линча определено зависимостью
K2( D) =J
IId-1w||
W
где Нм(кТ0) - переходная функция идентифицированной модели, а также описать процедуру минимизации АН.
При решении задачи идентификации цифровых объектов на основе ВИМ может возникнуть ситуация, когда отсутствует достаточная априорная информация о порядке ПФ исследуемого объекта. С целью формализации выбора порядка модели возможно привлечение аппарата обусловленности матрицы Б. В виду важности и практической значимости проблемы поиска подходящей структуры модели рассмотрим этот подход подробнее.
3. Применение чисел обусловленности для повышения точности решения задачи идентификации
Изменение значений степеней полиномов числителя т и знаменателя п ПФ модели объекта означает изменение размерности соответствующих матриц, описывающих модели сигналов и их взаимосвязи. Высокая размерность матриц может привести к результатам, которые могут не соответствовать физическому содержанию задачи. Свойства идентифицированной модели могут весьма значительно отличаться от свойств самого объекта. Для минимизации негативного влияния этого обстоятельства и получения результатов, имеющих практическое значение, применяются различные методы регуляризации. Эти методы нуждаются в наличии определенного объема априорных данных о свойствах объекта, в целом ряде практических задач такие сведения отсутствуют [8]. В связи с этим возникает вопрос о необходимости введения в алгоритм решения задачи идентификации некоторого механизма, который позволил бы контролировать корректность задачи.
В практических задачах целесообразно использовать математический аппарат чисел обусловленности, применяемый для апостериорной оценки ошибок при решении СЛАУ (10) [9, 10]. Этот подход рассматривает числа обусловленности как количественную оценку меры обусловленности задачи, т. е. оценку того, насколько она чувствительна к ошибкам исходных данных.
Выделяют три типа чисел обусловленности [9]. Естественное число обусловленности К(П) определено выражением
||б -
КД Б) = -'
Важное свойство числа К2(Б) состоит в том, что позволяет получать верхнюю и нижнюю границы ошибки при обращении матрицы Б.
Числа обусловленности К(П), К2(Б) непосредственно зависят от нормы |Д, а само точное решение в задаче идентификации не известно. Поэтому перспективы применения чисел К(Б) и К2(Б) в рассматриваемой задаче весьма ограничены. С этой точки зрения более перспективно использование стандартного числа обусловленности А'3(Б)=||Б-1||^||Б||, которое может быть найдено на основе вычисления векторных норм для матриц Б-‘иБ [9]:
И =
max
D = max
ZV,
i=1
n ,
ZN
где V, V— - элементы матриц D и D-1 соответ-
4/’ '4/
ственно
Формулы для ||Б|| и ||Б-1|| свободны от отмеченных недостатков первых двух стандартных чисел обусловленности и потому использование числа К3(Б) более перспективно.
Однако общим недостатком каждой из вышеприведенных оценок остается необходимость нахождения нормы ||Б-1||, поскольку найти норму матрицы без вычисления самой обратной матрицы Б- удается чрезвычайно редко. Вычисление элементов обратной матрицы Б-1 сзаранее гарантированной достаточной точностью в некорректных задачах, к которым относится задача идентификации, практически невозможно. В этой связи большие преимущества для применения идентификации в рамках ВИМ может дать использование подхода, предложенного в [8].
Известно, что любая норма матрицы не меньше ее наибольшего собственного значения, т. е. ||Б||>шах|ЯБ|. Кроме того, собственные значения матриц Б и Б1 взаимно обратны и поэтому справедливо соотношение:
D-1 > max
1
1
N min\VD
Отсюда выражение для нахождения стандартного числа обусловленности K3(D) принимает вид:
max |ЯЛ I
K3(D) = —^.
min \Ad I
Отслеживая изменения величины стандартного числа обусловленности в каждом конкретном слу-
i=1
чае можно найти наилучшее соотношение параметров структуры модели. Значение К3(Б)=1 свидетельствует о наиболее благоприятном случае, когда оценка относительной погрешности решения задачи идентификации совпадает с оценкой относительной погрешности исходных данных. К сожалению такая ситуация редко возникает в практических задачах [9, 10]. Зачастую К3(Б)>>1 и имеетпо-рядок сотен, тысяч и более.
Числа обусловленности К3(Б) позволяют делать некоторые выводы относительно структуры модели объекта идентификации. При первоначальных малых значениях параметров модели т и п значение К3(Б) будет наиболее близким к наилучшему случаю К3(Б)=1, хотя величины относительной ошибки АН(кТ0) временных динамических характеристик объекта идентификации и его модели могут быть значительными. С целью уменьшения ошибки АН(кТ0) следует увеличивать т и п. При этом будет происходить рост величины К3(Б). На некотором этапе может сложиться ситуация, когда дальнейшее уменьшение относительной ошибки АН(кТ0) за счет усложнения структуры модели окажется не в состоянии компенсировать погрешность, вносимую погрешностями вычислений и некорректностью задачи. Это проявляется в виде резкого, иногда лавинообразного нарастания величины К3(Б), что показывает возрастание ошибки при решении СЛАУ. Полученная информация позволяет обоснованно выбрать значения т и п, при которых модель объекта или системы будет иметь наименьшую погрешность в соответствии с критерием (11). Для подтверждения сказанного рассмотрим численный пример.
4. Численный пример
Пусть известна переходная функция объекта управления, аналитическое выражение которой имеет вид
Н0(/)=1-2,3-0,368'+1,5-0,135'-0,2-(1,832-10-2)'.
Известно время квантования Т0=0,5 с и время наблюдения /„=12 с. Найдем модель объекта в форме дискретной ПФ
Таблица. Результаты идентификации объекта управления
О г + ь
№ (г) =—-----------т-1
а г" + <:
т-1 . . і
г +... + О
+... +1
п > т
т п Ш ПФ идентифицированных моделей АЛ, %
0 1 1,67 №(г) = -0’75 -1,75 г +1 21
1 1 44,41 „7, . 9,81-10-3г -0,69 № ( г) = — ; -1,69 г + 1 14
1 2 294,29 , -0,15г - 0,53 Ж (г) = - -1,39г2 -0,28г +1 7
2 2 8,32-103 „гг , -2,19-10-3г2 -0,16г-0,59 Ж (г) = 2 -1,61г2 - 0,14 г +1 5
2 3 6,78 -108 Ж -12,11г2 -0,88г + 2,41 1
-34,28г3 + 32,73г2 -10,04г +1
для различных значений т и п. Для каждой из полученных моделей найдем огибающие их переходных характеристик, оценим близость полученных решений к Н0(пТ0) по критерию (11). Значения числа К3(Б) для различных т и п используем для выбора модели, имеющей наибольшую точность по (11). Расчеты по приведенному ранее алгоритму идентификации дискретных объектов выполнялись в среде программного пакета MathCAD. Результаты расчетов сведены в таблицу.
На рис. 2 построены графики переходных функций объекта Н0(/) и модели Нм(кТ0), для которой получено наименьшее значение критерия (11).
Рис. 2. Графики переходных характеристик
Из таблицы видно, что при возрастании параметров т и п точность идентификации увеличивается. При этом число К3(В) имеет минимальное значение 1,67 при т=0, п=1. Последующее увеличение т и п приводит к снижению относительной ошибки АН и росту числа К3(В), например, при т=2, п=3 значение К(В) возрастает до 6,78-108, а относительная ошибка АН не превышает 1 %. Дальнейшее увеличение значений т и п модели даже на единицу приводит к тому, что получить модель объекта в классе устойчивых систем уже не удается. Наиболее вероятная причина - резкий рост вычислительной погрешности, обусловленной некорректностью задачи. На основании полученных результатов можно определить структуру модели с параметрами т=2 и п=3. Дальнейшее снижение погрешности идентификации возможно на основе применения методов регуляризации, но эти вопросы в данной статье не рассматриваются.
Выводы
Применение вещественного интерполяционного метода для идентификации дискретных объектов и систем управления позволяет получить эффективный в вычислительном плане алгоритм идентификации. Для получения моделей, в наибольшей степени удовлетворяющих условиям кон-
кретной задачи по точности, сложности математического описания и т. д., целесообразно использовать математический аппарат стандартных чисел обусловленности, который достаточно просто вписывается в структуру алгоритма идентификации. Эти обстоятельства позволяют рассматривать ве-
щественный интерполяционный метод как эффективный инструмент для решения практических задач адекватного математического описания цифровых объектов и систем.
Работа выполнена в рамках ФЦП «Научные и научно-педагогические кадры инновационной России» на 2009-2013 гг.
СПИСОК ЛИТЕРАТУРЫ
1. Балакирев В.С., Дудников Е.Г, Цирлин А.М. Экспериментальное определение динамических характеристик промышленных объектов управления. - М.: Энергия, 1967. - 232 с.
2. Ордынцев В.М. Математическое описание объектов автоматизации. - М.: Машиностроение, 1965. - 360 с.
3. Гроп Д. Методы идентификации систем. - М.: Мир, 1979. -302 с.
4. Цыпкин Я.З. Основы информационной теории идентификации. - М.: Наука, 1984. - 320 с.
5. Гончаров В.И. Синтез электромеханических исполнительных систем промышленных роботов. - Томск: Изд-во ТПУ, 2002. -100 с.
6. Алексеев А.С., Антропов А.А., Гончаров В.И., Замятин С.В., Рудницкий В.А. Вещественный интерполяционный метод
в задачах автоматического управления. - Томск: Изд-во ТПУ, 2008. - 217 с.
7. Дёч Г. Руководство к практическому применению преобразования Лапласа и /-преобразования. - М.: Наука, 1971. - 288 с.
8. Бахвалов Н.С., Жидков Н.П., Кобельков ГМ. Численные методы. - М.: Бином. Лаборатория Знаний, 2003. - 630 с.
9. Райс Дж. Матричные вычисления и математическое обеспечение. - М.: Мир, 1984. - 265 с.
10. Трофимов А.И., Егупов Н.Д., Дмитриев А.Н. Методы теории автоматического управления, ориентированные на применение ЭВМ. Линейные стационарные и нестационарные модели. - М.: Энергоатомиздат, 1997. - 651 с.
Поступила 06.04.2012 г.
УДК 62-533.65
СИНТЕЗ МАТЕМАТИЧЕСКОЙ МОДЕЛИ ПОДСИСТЕМЫ УПРАВЛЕНИЯ ПРОГРЕВОМ ВАКУУМНОЙ КАМЕРЫ КАЗАХСТАНСКОГО МАТЕРИАЛОВЕДЧЕСКОГО ТОКАМАКА КТМ
А.Г. Коровиков, В.М. Павлов*, Д.А. Ольховик
Институт атомной энергии НЯЦ Республики Казахстан, г. Курчатов *Томский политехнический университет E-mail: [email protected]
Исследован процесс прогрева вакуумной камеры токамака КТМ как объект управления, разработана математическая модель процесса как объекта управления. Сделан вывод о том, что для управления технологическим процессом прогрева необходимо идостаточно использовать систему управления, построенную на основе ПИ-регуляторов. Определены коэффициенты параметров настройки регулятора для каждой зоны прогрева.
Ключевые слова:
Система управления, математическая модель, вакуумная камера, нагрев.
Key words:
Control system, mathematical model, vacuum chamber, heat.
В настоящее время в г. Курчатов, Республика Казахстан, ведутся работы по реализации проекта строительства термоядерной установки Казахстанский материаловедческий токамак КТМ.
На этапе подготовки установки к эксперименту вакуумная камера токамака КТМ и патрубки должны прогреваться до температуры 200 °С с целью обезгаживания и уплотнения вакуумных конструктивов. Процесс должен проходить с обеспечением контроля равномерности прогрева. Техническая структура системы управления подробно рассмотрена в статье [1].
Синтез систем автоматического регулирования (САР) является главной задачей теории управле-
ния, в результате определяются состав, структура САР и параметры всех ее устройств из условия соответствия заданному комплексу технических требований [2, 3]. Основной задачей, которая должна быть решена на данном этапе, является получение модели объекта управления, его математического описания, которое представляет собой, систему обычных линейных/нелинейных дифференциальных или алгебраических уравнений.
При построении математической модели подсистемы омического и индукционного нагрева воспользуемся структурным методом. На первом этапе исходная принципиальная схема заменяется упрощенной функциональной схемой. Для
