Идентификация параметров исполнительных устройств регуляторов паровой энергетической турбины с использованием аппарата матричных операторов Текст научной статьи по специальности «Математика»

УДК 681.5.015.23

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

© Ю.П. Корнюшин, Н.Д. Егупов, П.Ю. Корнюшин Калужский филиал МГТУ им. Н.Э. Баумана, Калуга, 248000, Россия

Предложен алгоритм идентификации параметров — постоянных времени турбины с использованием градиентного метода с настраиваемой моделью. Настраиваемая математическая модель имеет такую же структуру, как и объект идентификации. Критерий идентификации формируется на основе функции потерь, которая представляет собой невязку между левой и правой частями уравнения, описывающего настраиваемую модель. Тем самым удается избежать необходимости нахождения в явном виде решения нелинейного уравнения для настраиваемой модели. Вместо выходного сигнала в модели используется сигнал, наблюдаемый на выходе идентифицируемого объекта. Поскольку математические модели являются нелинейными, для решения задачи применены линеаризация Ньютона -Канторовича и аппарат матричных операторов. Рассмотрены особенности вычисления вектора градиента, алгоритм идентификации и его организация. Приведены результаты идентификации двух постоянных времени для математической модели турбины ПТ-12/15-35/10М.

Ключевые слова: идентификация, турбина, градиентный оператор, критерий, параметры, постоянные времени.

Введение. Эффективность функционирования любой системы определяется тем, как работают ее управляющие устройства. Их проектирование предполагает наличие достаточно полной информации о математических моделях объектов управления. Без знания математической модели объекта нельзя синтезировать оптимальные для данной системы регуляторы. В связи с этим одной из определяющих задач при проектировании любой системы управления является идентификация ее математической модели. Наиболее часто это пара-мет-рическая идентификация. Результаты решения задачи идентификации служат основой для проектирования систем управления, не располагая которыми часто нельзя осуществить ни оптимизацию, ни синтез регуляторов, ни анализ систем управления. Для одних объектов управления идентификацию можно выполнять с помощью тестовых исследований с последующей обработкой выходных сигналов на заданные испытательные сигналы и построением математической модели. Для других объектов тестовые испытания недопустимы и идентификацию можно провести только в режиме их нормальной эксплуатации. К таким объектам управления относятся энергетические паровые турбины.

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

В работе предложен алгоритм, позволяющий проводить идентификацию в режиме нормальной эксплуатации объекта.

Отметим также, что рассматриваемый подход применим не только к конкретному объекту — турбине ПТ-12/15-35/10М, но и к другим турбинам и турбоагрегатам, например К-8-0,65ПА, К-800-130/3000 [1-3], поскольку они имеют близкие по структуре математические модели.

Постановка задачи. Полагаем, что математическая модель объекта идентификации — турбины — и ее системы регулирования по частоте вращения ротора описывается следующей системой дифференциальных уравнений [4]:

Здесь ф — относительное изменение частоты вращения ротора турбины; £ — относительное изменение расхода пара; ц — относительное отклонение поршня сервомотора; г — относительное отклонение поршня золотника; 5 — степень неравномерности системы регулирования; Тф, Т£, Тц, Тг — соответствующие постоянные времени; X — относительное отклонение электрической нагрузки генератора; Лг0 — относительная высота треугольного профиля кромок поршня отсечного золотника сервомотора.

Структурная схема системы регулирования турбины по угловой скорости изображена на рис. 1.

тф ж ф/ж = (1 -е)(£-х) / 5-0ф,

Т й й = ц-£,

Тц йц/й = <

2

г , ,

0,5-г при г < Лг0,

Лг0

2 - 0,5Лг0в1§п г при |г| > Л(

0

(1)

Рис. 1. Структурная схема системы регулирования турбины

Требуется определить постоянные времени сервомотора Т^ и отсечного золотника Т2.

Решение задачи. В основу положена известная концепция идентификации с настраиваемой моделью [5]. Алгоритм идентификации можно пояснить схемой, представленной на рис. 2.

Рис. 2. Схема алгоритма идентификации

Полагаем, что настраиваемая модель имеет такую же структуру, как и объект идентификации. Идентификации подлежат только параметры объекта Р = [р1 р2 ], р1 = Т^, р2 = Тг. Соответственно настраиваемая модель описывается системой уравнений вида

Тф * =-е*,+Х2 -1 *, ф йг 5 2 5

Т

Р1

йх2

йг

йх3

йг

■ = - х2 + xз,

0,5-г при |г| < Аг0,

А^

2 - 0,5Аг081§и г при |г| > Аг0,

2

Р2

&х4 1 & ~ 5

Х1 — Х3 — Х4.

(2)

Под выходами объекта идентификации и настраиваемой модели будем понимать соответствующие их математическим моделям векторы фазовых координат: Хк (г) = [ф(г )£(г )д(г) 1(1 )]т — вектор фазовых координат объекта, Хм(г) = [х1(^)х2(г)х3(г)х4(г)]т — вектор фазовых координат модели. Поскольку вектор Хм (г) косвенно зависит от вектора настраиваемых параметров Р, можно записать: Хм (г, Р). Необходимо также отметить, что все элементы вектора Хк (г) — измеряемые величины.

В ненормированной нормальной форме Коши уравнения можно записать следующим образом:

где

Р • Хм (г, Р) = Р (X м (, Р)) + ВА(),

т 0 0 0

0 0 0

0 0 Р1 0

0 0 0 Р2

Р (Хм (г, Р)) =

—0х+ х2

о

—х2 + х3

2

2 | | 0, 5-2 при 2 < Д0

А

2 — 0,5 Д^ВР 2 при 2 > Д 1

--х1 — хз — х4

(3)

В

—1 о о о 5

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

Е(г, Р) = ХЛ (г) — Хм (г, Р).

Для оценки соответствия настраиваемой модели идентифицируемому объекту используем следующий критерий качества идентификации:

3(Р) = ^ (Е(г, Р)) . (4)

Здесь ¥ (Е(', Р)) — функция потерь:

¥ (Е(', Р)) = | (Xл (') - Xм (', Р))т (Xл (') - Хм (', Р))Л,

'г

А' = ('г, ) — интервал времени, на котором выполняется один шаг

итерационного процесса идентификации.

Таким образом, критерий качества идентификации принимает вид

3(Р) = ¥ (Хл ('), Хм (', Р)). (5)

Чем меньше значение функции потерь, тем выше качество идентификации. Однако данная функция потерь требует нахождения в явном виде вектора Xм (', Р) фазовых координат модели.

Более эффективной с точки зрения реализации градиентного метода идентификации является функция потерь, формируемая как невязка между левой и правой частями уравнения (3). В данном случае необходимо вместо выходного сигнала Xм (') модели использовать сигнал Xл ), наблюдаемый на выходе идентифицируемого объекта. Тогда функцию потерь можно представить в виде

¥ ( X л ('), X м (', Р ) )= |

Р • X л (') - Р (X л ('))- )

г _ .

Р • X л (') - Р (X л ('))- )

Ж.

(6)

Критерий качества идентификации по-прежнему имеет форму (5).

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

(Р) =

д3 (Р) д3 (Р)

дР1 дР2

= 0.

(7)

УУТ 3 (Р) Ч^)

> 0.

(8)

}г, ]=1

Полагаем также, что достаточное условие оптимальности (8) выполнено. Тогда можно построить следующий простейший градиентный алгоритм идентификации:

Рг+1 = Р-РгУЗ (Р), (9)

где рг — коэффициент, изменяющий скорость движения по градиенту.

Ясно, что задача сводится через искомые параметры настраиваемой модели к выражению для вектора градиента V/(Р) .

Вычисление вектора градиента V/(Р). Недостаток рассмотренной функции потерь (6) заключается в том, что в подынтегральное выражение входит производная Хя (г) измеряемого сигнала. Выходной сигнал объекта, как правило, измеряется с использованием цифровых датчиков, поэтому при вычислении его производных возможны большие погрешности, которые вносят ошибки в результаты идентификации. Эффект усиливается, поскольку вектор измеряемых фазовых координат содержит четыре элемента.

Избежать операции дифференцирования сигнала Хя (г) можно, если преобразовать функцию потерь (6). Для этого необходимо воспользоваться линеаризацией Ньютона - Канторовича и методом матричных операторов.

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

Для рассматриваемой нелинейности с параметром Дг0 = 0,3 при I = 5 коэффициенты ^ имеют следующие значения: с0 = 0,0; с1 = 0,475; с2 = 0,0; с3 = 1,295; с4 = 0,0; с5 =—0,832.

Полагаем, что в настраиваемой модели нелинейная характеристика является аналитической (10). Тогда настраиваемая модель описывается системой уравнений (3) с правой частью:

(10)

i=0

Р (Хм (г, Р))

—х2 + х3 5

i=0

1

— х1 — х3 — х4

о

Линеаризация Ньютона - Канторовича достаточно подробно рассмотрена в работе [6]. С использованием данной процедуры можно перейти от уравнения (3) к последовательности уравнений:

Р • Хм к+1 (г) = Ак (г )Хмк+1 (г) + БХ(г) + Zk (г),

Хмк+1(гд = к = 0, 1 2, ...,

(11)

где Ак (г) = РХ( Хмк(г));

Zк (г) = — Ак (г )Х мк(г) + Р (Хмк(г)), РХ (Хмк(г)) = ■

I 1i ,1=1

Уравнения (11) являются линейными на каждом шаге итерационного процесса по отношению к искомому решению Хмк(г) .

Для настраиваемой модели матрица Ак (г) и вектор Zk (г) имеют

дхк

вид:

Ак (г) =

—е 0

0

(1—) 0 5

—1 1

0 0

0

0 2 ^х4к (г) i=1

— 0

—1

—1

Zk (г) =

0 0 2 (1—0^4* (г) 0

i=0

С учетом уравнений (11) критерий качества идентификации можно записать следующим образом:

3 (Р) = |

Р • Хя (г) — Ак (г )Х я (г) — БХ (г) — Zk (г)

Р • Хя (г)— Ак (г)ХК (г)(г) — БХ(г) — Zk (г)

ёг.

(12)

Уравнения (11), а также критерий (12) можно параметризовать и тем самым избежать процедуры дифференцирования сигнала, наблюдаемого на выходе идентифицируемого объекта. Для этого воспользуемся аппаратом матричных операторов [7].

Сигналы Xк (1), Х(1), функцию Zk (1) и начальные условия Хк (^) разложим по заданному ортонормированному базису

Ф = {фу (0: j = t , ti+1 ]} :

хй ()=Бг4(/)Сх, х(о=БжОСх, zk (о = БТ4(о^ , х&)=бт4('Сх°.

Здесь ) = Жав (¥Т\г) ¥Т2(/)¥Т3(/) ¥Т\г)); БТ\(/) = ¥Т\г),

¥Т1 (I) =... = ¥Т4 (1) = [ф1 (1) ф2 (1) ... фN (1)].

х Л z Х

Спектральные характеристики С , С", Ск представим в следующем виде:

(СХ =

С С ^2 С ? С Х4

где Сх' =

С? С?'... С

С2' =

Сz = с С 22 С С2

С2' С2' С2' С1 С2 .

, С Х =

СХГХ ГХ 1 С2 .

*'+1 _ _ *'+1 _

С? = | X* (Офу (№, у = 1, Ы, к = 1,4; СХ = {Х(0фу (¿)Л, у = 1, N ,

'' ''

_ _

Ск = | 2* (1) ф у (1) Л, у = 1, N, к = 1,4.

В этом случае уравнение (11) и критерий (12) запишем как

Р • (СХ - СХ° ) = кАСХ+кХСХ +к1 Сz, (13)

3(Р) = (Р • (СХ - СХ°) - кАСХ - кХСХ - к7Сz)т X

X (Р • (СХ - СХ°) - кАСХ - кХССХ - к7 Сz), (14)

где Р =

Т Аф 0 0 0

0 0 0

0 0 р1 0

0 0 0 Р 2

0 = (°)4у=1; тф = (тФ тФ тФ тФ),

Т= = ё1а§ (Т Т Т Т4), Р1 = (р р р р),

Р 2 = (Р2 Р2 Р2 Р2 ) .

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

/ (Р) = ( Р • (Х —К А( Х — К} (} — К7 ( Z

X (Р • (Х - К А(Х — -К } — -К 7( Z)

Г X

(15)

А } 7

В выражения (13) и (15) входят матрицы К , К , К , представляющие собой спектральные характеристики ядер

К х (г, т) = {кХ (г, т) = 1(г, т)4 (х)}4==1, К х (г, т) = {к} (г, т) = 1(г—т)Ъ} , К 7 (г, т) = {к7 (г, т) = 1(г, т)}4=1

интегрального уравнения Фредгольма второго рода

Р •( Хмк+1(г) — Хмк+1(0) ) =

т т т

= | К X (г, т) Хмк+1 (т)ёт +1К} (г, т)}(т)ёт +1К 7 (г, т) Zk (т)ёт,

0 0 0

следующего из дифференциального уравнения (11) при его интегри-

А } 7

ровании. Соответственно матрицы К, К , К имеют вид: К А

{Ах =■ К ЧБ'1,, К7 = £i }м. А х={ар }1 „=,• "М* }1 „=,• i С „=,-

Элементы матриц Ах, В} и ZI определяются следующими фор-

мулами:

ТТ тт

Пу = Г Г Vх (1 ^^а™ (т\йтж ^i = Г

\\кХ (г, т)фя (г)ф р (т) ётёг, = Цк} (г, т)ф5 (г)Ф р (т) ётёг,

тт

2*р =ЯК7 (г, т)фя (г)Фр (т) ётёг .

г0 г0

Найденные спектральные характеристики однозначно определяют критерий качества идентификации. Выполняя операцию вычисления вектора V/ (Р), имеем

V/ (Р) = ЦР (Р • (СХ — к АсХ — К } — К7 (Z).

Таким образом, все элементы, необходимые для выполнения идентификации, определены.

Особенности алгоритма идентификации и его организация.

Назначение системы регулирования — стабилизация частоты вращения ротора турбины. Если электрическая нагрузка генератора не будет изменяться, система регулирования через определенное время выполнит свою задачу, т. е. частота вращения ротора будет стабилизирована. Поэтому даже при «ненастроенных» параметрах модели может возникнуть такая ситуация, когда выходные сигналы объекта и модели будут не только постоянными, но и равными. В этом случае функция потерь принимает нулевое значение и возникает эффект достижения поставленной цели, хотя задача идентификации не решена. Следовательно, обязательное условие правильно решенной задачи состоит в том, что электрическая нагрузка генератора в процессе идентификации должна изменяться. Реализацию алгоритма осуществляют следующим образом. Выбирают интервал времени, на котором выполняется один шаг итерационного процесса идентификации Лг = (гi, г1+1). Величину Дг можно выбирать экспериментально либо задавать соизмеримой или немного большей, чем максимальная идентифицируемая постоянная времени. Это связано с тем, что в заданный промежуток времени не должны произойти переходные процессы в случае установления постоянного значения электрической нагрузки. В течение заданного интервала времени Лг вычисляются элементы градиента по формуле (15) и выполняется один шаг итерационного градиентного алгоритма идентификации, определяемого формулой (9). Далее процесс повторяется. Условием окончания итерационного процесса является достижение заданной точности. Кроме того, организовывать итерационный процесс, предписанный линеаризацией Ньютона - Канторовича, нет необходимости, поскольку в выражение для критерия идентификации входит выходной сигнал объекта, а он всегда представляет собой точное решение системы уравнений (2).

Пример. Приведем результаты идентификации указанных выше двух постоянных времени: Т^ и Т2. В качестве базиса используются функции Уолша ^ = 32).

Начальные значения постоянных времени: Тц = 0,04 с и Т = 0,01 с. Известные параметры объекта — турбины ПТ-12/15-35/10М: Т^= 0,12 с, Тф = 6 с, 5 = 0,04, Лг0 = 0,3, коэффициент самовыравнивания е = 0,05.

Интервал времени, на котором проводится идентификация, Лг = 0,1с. Параметры итерационного процесса: р1 = 0,0085, р2 = 0,005.

На рис. 3, 4 приведены зависимости идентифицируемых параметров — постоянных времени — от числа ^итер итераций.

Рис. 3. Изменение постоянной времени Г^ в процессе итераций (1), точное значение (2)

Рис. 4. Изменение постоянной времени Т2 в процессе итераций (1), точное значение (2)

Выводы. Из рис. 3, 4 следует, что предложенный алгоритм обеспечивает необходимое качество идентификации. Потребовалось лишь Ыитер = 2° итераций для получения заданной точности. Традиционно используемые в системах регулирования данного класса турбин постоянные времени имеют значения, близкие к следующим: Т^= °,1 с, Т2 = °, °3 с. Ясно, что идентификация выполнена вполне адекватно текущей ситуации работы турбины.

Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований и Правительства Калужской области (гранты №№ 14-41-03071, 14-48-03013).

ЛИТЕРАТУРА

[1] Мельников Д.В., Фишер М.Р. Математическая модель контура регулирования частоты вращения ротора паровой турбины К-8°°-13°/3°°°. Вестник МГТУ им. Н.Э. Баумана. Сер. Машиностроение, 2°11, спец. вып. «Энергетическое и транспортное машиностроение», с. 197-215.

[2] Корнюшин Ю.П., Мельников Д.В., Егупов Н.Д., Корнюшин П.Ю. Исследование и расчет параметров элементов системы регулирования частоты вращения ротора турбины с учетом параметрической неопределенности математической модели. Вестник МГТУ им. Н.Э. Баумана. Сер. Естественные науки, 2014, № 1, с. 78-93.

[3] Мельников Д.В., Егупов Н.Д. Синтез систем регулирования энергетических турбин в условиях параметрической неопределенности. Известия Тульского государственного университета. Технические науки, 2°11, вып. 5, ч. 1, с. 108-113.

[4] Калашников А.А. Динамика регулирования турбин. Москва, Энергоатом-издат, 1999, 328 с.

[5] Цыпкин Я.З. Основы информационной теории идентификации. Москва, Наука, 1984, 320 с.

[6] Гайский В.А., Егупов Н.Д., Корнюшин Ю.П. Применение функций Уолша в системах автоматизации научных исследований. Киев, Наукова думка, 1993, 212 с.

[7] Пупков К.А., Егупов Н.Д., ред. Методы инженерного синтеза сложных систем управления: аналитический аппарат, алгоритмы приложения в технике. Часть II. Вычислительно-аналитический эксперимент: аппарат матричных операторов и вычислительные технологии. Москва, Изд-во МГТУ им. Н.Э. Баумана, 2012, 416 с.

Статья поступила в редакцию 04.12.2014

Ссылку на эту статью просим оформлять следующим образом: Корнюшин Ю.П., Егупов Н.Д., Корнюшин П.Ю. Идентификация параметров исполнительных устройств регуляторов паровой энергетической турбины с использованием аппарата матричных операторов. Математическое моделирование и численные методы, 2015, № 2(6), с. 73-86.

Корнюшин Юрий Петрович — д-р техн. наук, профессор, заведующий кафедрой «Системы автоматического управления» Калужского филиала МГТУ им. Н.Э. Баумана. Автор более 100 печатных работ в области энергетики и электротехники, оптимальных систем, моделирования технических систем. Область научных интересов: нелинейные, оптимальные, робастные системы управления. e-mail: theroland@yandex.ru.

Егупов Николай Дмитриевич — д-р техн. наук, профессор кафедры «Системы автоматического управления» Калужского филиала МГТУ им. Н.Э. Баумана. Автор более 200 печатных работ в области теории матричных операторов, энергетики, оптимальных систем. Область научных интересов: теория матричных операторов, нелинейные, оптимальные, робастные системы управления. e-mail: theroland@yandex.ru.

Корнюшин Петр Юрьевич — аспирант кафедры «Системы автоматического управления» Калужского филиала МГТУ им. Н.Э. Баумана. Автор 23 печатных работ в области энергетики и электротехники, моделирования технических систем, управления техническими системами. Сфера научных интересов: нелинейные, оптимальные, робастные системы управления. e-mail: kornyushin.petr@gmail.com.

Identification of parameters of regulator actuators for steam power turbines using matrix operator apparatus

© Yu.P. Kornyushin, N.D. Egupov, P.Yu. Kornyushin

Bauman Moscow State Technical University, Kaluga Branch, Kaluga, 248000, Russia

In the paper we propose an algorithm ofparameters (time constants of the turbine) identification using the gradient method with an adaptive model. The adaptive mathematical model has the same structure as the identification object. The identification criterion is based on the loss function, which is the misalignment between the left and right sides of the equation, which describes the adaptive model. Thus it is avoided the need of finding the solution of a nonlinear equation for the adaptive model in an explicit form. In the model the signal observed at the output of the identified object is used instead of the output signal. Since mathematical models are nonlinear, the Newton - Kantorovich linearization and the matrix operator apparatus are applied to solve the problem. The features of gradient vector computation and features of the identification algorithm and its organization are considered. The results of the two time constants identification for the mathematical model of the turbine PT-12/15-35/10M are presented.

Keywords: identification, turbine, gradient operator, criterion, parameters, time constants.

REFERENCES

[1] Melnikov D.V., Fisher M.R. Vestnik MGTU im. N.E. Baumana. Seria Mashinostroenie - Herald of the Bauman Moscow State Technical University. Series: Mechanical Engineering, 2011, spetsialnyi vypusk "Energeticheskoe i transportnoe machinostroenie" [special issue "Power and transport mechanical engineering"], pp. 197-215.

[2] Kornyushin Yu.P., Melnikov D.V., Egupov N.D., Kornyushin P.Yu. Vestnic MGTU im. N.E. Baumana. Seria Estestvennye nauki - Herald of the Bauman Moscow State Technical University. Series: Natural Sciences, 2014, no. 1, pp. 78-93.

[3] Melnikov D.V., Egupov N.D. Izvestiya Tulskogo gosudarstvennogo universiteta. Tekhnicheskie nauki - Proceedings of the Tula State University. Engineering sciences, 2011, issue. 5, part 1, pp. 108-113.

[4] Kalashnikov A. A. Dinamika regulirovaniya turbin [Turbine control dynamics]. Moscow, Energoatomizdat Publ., 1999, 328 p.

[5] Tsypkin Ya.Z. Osnovy informatsionnoy teorii identificatsii [Basics of information theory of identification]. Moscow, Nauka Publ., 1984, 320 p.

[6] Gayskiy V.A., Egupov N.D., Kornyushin Yu.P. Primenenie funktsiy Walsha v sistemakh avtomatizatsii nauchnykh issledovaniy [Application of the Walsh functions in research automation]. Kiev, Naukova Dumka Publ., 1993, 212 p.

[7] Pupkov K.A., Egupov N.D., eds. Metody inzhenernogo sinteza slozhnykh system upravleniya: analiticheskiy apparat, algoritmy prilozheniya v tekhnike. Chast II. Vychislitelno-analiticheskiy eksperiment: apparat matrichnykh operatorov i vychislitelnye tekhnologii [Methods for the engineering synthesis of complex control systems: an analytical apparatus, algorithms of applications in engineering. Part II. Computational and analytical experiment: The matrix operators apparatus and the computational technology]. Moscow, BMSTU Publ., 2012, 416 p.

M.n. KopnwrnuH, H.ff. Ezynoe, n.M. KopnmwuH

Kornyushin Yu.P., Dr. Sci. (Eng.), professor, head of the Automatic Control Systems Department at Bauman Moscow State Technical University, Kaluga Branch. The author of about 100 publications in the field of power engineering and electrical engineering, optimal systems, simulation of engineering systems. Research interests: nonlinear, optimal, robust control system. e-mail: theroland@yandex.ru

Egupov N.D., Dr. Sci. (Eng.), professor of the Automatic Control Systems Department at Bauman Moscow State Technical University, Kaluga Branch. The author of about 200 publications in the field in the theory of matrix operators, power engineering, optimal systems. Research interests: theory of matrix operators, nonlinear, optimal, robust control system. e-mail: theroland@yandex.ru

Kornyushin P.Yu., post-graduate of the Automatic Control Systems Department at Bauman Moscow State Technical University, Kaluga Branch. The author of 23 publications in the field of power engineering and electrical engineering, simulation of engineering systems, control of engineering systems. Research interests: nonlinear, optimal, robust control system. e-mail: kornyushin.petr@gmail.com