Сокращение порядка управляющего устройства при синтезе САУ полиномиальным методом Текст научной статьи по специальности «Электротехника, электронная техника, информационные технологии»

УДК 62-50-83

Сокращение порядка управляющего устройства при синтезе САУ полиномиальным методом1

В.В. Тютиков, Л.Г. Копылова, И.А. Тихомирова, Е.М. Шляцкая ФГБОУВО «Ивановский государственный энергетический университет имени В.И. Ленина»,

г. Иваново, Российская Федерация E-mail: tvv@ispu.ru

Авторское резюме

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

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

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

Ключевые слова: полиномиальный регулятор, синтез САУ, эквивалентные структурные преобразования, редукция управляющих устройств, электромеханическая система, передаточная функция регулятора, астатизм по управляющему воздействию, метод коэффициентов ошибок.

Reducing the dimension of control unit in the ACS synthesis by the polynomial method

V.V. Tyutikov, L.G. Kopylova, I.A. Tikhomirova, E.M. Shlyatskaya Ivanovo State Power Engineering University, Ivanovo, Russian Federation

E-mail: tvv@ispu.ru

Abstract

Background: Modern methods of regulator synthesis are often hard to realize technically as it is not easy to synthesize complex control devices. For example, input-output regulators are synthesized by the modal control method based on solving polynomial equations and employs pre-filters, the complexity of which is comparable with that of the regulators themselves. All this makes it quite urgent to develop an approach simplifying such control devices without affecting the qualitative parameters of the system as a whole.

Materials and methods: The regulators have been synthesized by the polynomial synthesis method ensuring the required quality of the processes in an automatic control system by measuring only the output point of the control object. The reduction of control device is based on the method of equivalent structural transformations, which guaranties the permanence of automatic control system transfer function. And the astatism of control influence is ensured by the error coefficient method.

Results: We propose a method of control device reduction based on implemention of the method of equivalent structural transformations, which essentially reduces the control device complexity. Unlike the known methods, this method does not require recalculating control device coefficients in order to preserve the performance indicators of the automatic control system.

Conclusions: The proposed method reduces the order of the control device transfer function twofold, which makes the device much easier to implement. It also reduces the execution time of the control program, which is important when developing modern automatic control systems composed of objects with a high processing speed and short execution

1 Работа выполнена в рамках проектной части государственного задания Минобрнауки.

time. The results of mathematical simulation have confirmed the identity of the performance indicators in the reduced system and the initial automatic control system in processing control and disturbance actions.

Key words: polynomial regulator, ACS synthesis, equivalent structural transformations, control device reduction, electromechanical system, regulator transfer function, astatism of control action, error index method.

DOI: 10.17588/2072-2672.2017.5.044-052

Введение. Синтез САУ методами полиномиального управления [1-8] предполагает получение передаточной функции регулятора на основе решения диофантова уравнения (или л-уравнения), неизвестными в котором являются полиномы. При этом степени полиномов передаточной функции регулятора (ПФР) задаются.

В задании степеней полиномов ПФР допустим некоторый произвол проектировщика. Например, для обеспечения произвольного задания некоторых коэффициентов ПФР необходимо соответствующим образом увеличить степень полинома числителя ПФР [8], получив в результате переопределенную совместную систему уравнений (количество переменных больше количества уравнений). Изменить фильтрующие свойства регулятора можно выбором степени полинома его знаменателя.

В зависимости от того, какие степени полиномов заданы, может быть получен ряд регуляторов:

- дифференцирующие (степень полинома числителя ПФР выше степени полинома знаменателя) и технически реализуемые (степень полинома числителя ПФР не выше степени полинома знаменателя);

- минимальные (степени полиномов ПФР минимально возможны для решения уравнения синтеза) и неминимальные (степени полиномов ПФР выше минимально необходимых для решения уравнения синтеза).

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

Использование неминимальных регуляторов дает проектировщику свободу в выборе необходимых ему свойств регулятора, например фильтрующих возможностей (при увеличении степени полинома знаменателя ПФР), или некоторых коэффициентов полиномов регулятора (при увеличении степени полинома числителя ПФР). В первом случае коэффициенты ПФР определяются однозначно (уравнение синтеза имеет единственное решение), во втором случае ряд коэффициентов при синтезе может быть задан произвольно.

Наиболее часто в исследованиях используют минимальные технически реализуемые регуляторы. Однако и в этом случае основной трудностью технической реализации САУ является сложность управляющего устройства.

Модели и методы исследования.

Структура САУ с полиномиальным регулятором может быть различной (рис. 1, здесь Дф и Б(в), О(в) и Я(в) - полиномы знаменателя и числителя передаточных функций объекта управления и регулятора соответственно; уоСв), Д(5), и(в), у(в) - изображения по Лапласу сигналов задающего воздействия, ошибки управления, управления и выхода).

yo(s ) 1 У 0(s )

R(s)

A(s)

R(s) C(s )

u(s)

B(s ) A(s)

У (s)

а)

У o(s) 1 С^ ) B(s) A(s ) У (s )

R(s) C(s)

б)

Уo(s) A(s)

1 u(s ) B(s ) У (s )

C(s ) A(s )

R(s)

Рис. 1

в)

Для синтеза минимальных технически реализуемых регуляторов оперируют следующими степенями полиномов [7, 8]: degR(s) = degA(s) -1; degC(s) = degR(s); deg D(s) = deg A(s) + degC(s). (1)

Диофантово уравнение синтеза регулятора при этом имеет вид: A(s)C(s) + B(s)R(s) = D(s), (2)

где D(s) - характеристический полином САУ, определяющий качество переходных процессов и задаваемый проектировщиком.

Необходимо иметь в виду, что при использовании регулятора только в контуре главной обратной связи в передаточной функции САУ H(s) могут появляться «дополнительные» нули, определяемые корнями полинома R(s) (рис. 1 ,а)

R(s)B(s)

H(s ) =

A(s)C(s) + B(s)R(s)

или C(s) (рис. 2,б): H(s) = ■ C(S)ß(S)

A(s)C(s) + B(s)R(s)

Система, приведенная на рис. 1,в, лишена данного недостатка, но требует вычисления «чистых» производных, что на практике невыполнимо.

Влияние таких нулей может проявляться в появлении перерегулирования в переходных процессах. Для исключения их влияния [8] вводят префильтр вида HF(s) = 1/R(s) или HF(s) = 1/C(s) соответственно (рис. 1,а,б). Безусловно, это значительно усложняет техническую реализацию управляющего устройства, требуя выполнения большого количества дополнительных операций.

Часто одним из требований к системам управления является обеспечение астатизма по управляющему воздействию (следящие системы). В этом случае наличие префильтра не позволяет ограничиться введением соответствующего числа интегралов в контур управления, требуя использования префильт-ра с передаточными функциями:

HF s) = ™ или HF s) = ™.

F ( ; R(s) F ( ' C(s)

для систем на рис. 1,а,б соответственно.

В этом случае требуемый порядок аста-тизма по управлению обеспечивается совокупным использованием интеграторов в контуре управления и соответствующих коэффициентов полинома Q(s) [8, 9].

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

Результаты исследования. Примем за основу структуру на рис. 1,а. Для структур на рис. 1,б,в результат будет идентичным.

Очевидно,что

C(s)u(s) = R(s)A(s), u(s) =

R(s ) C(s ) 1

A(s),

A(s) = y0 (s) - y (s), y0 (s) =—— y0 (s).

R(s y

Тогда можно записать:

u(s) = Rs ( yS (s) - y (s )) 1 yo(s) - Rêy(s),

R(s) f 1

C(s) { R(s )

-y0(s) - у (s) =

C(s)

C(s)

C(s)u(s) = R(s) (y;(s) - y(s)):

= R(s) Ус (s) - y (s ) j = y0 (s) - R(s)y (s).

Если учесть (1) и принять, что deg^(s) = n, а также то, что полиномы ^(s) и D(s) при реше-

нии уравнения (2) нормированы (коэффициенты при высших степенях равны единице), получим

2° 2 + ■■■ + + Сд I

(s" 1 + cn_2sn 2 +... + qs + c0^)u(s) = = Уо (s) - ( rn-1s"-1 + Гп-2S"- 2 + ... + r1s + r0 ) y (s)

или

s"-1 u(s) = Уо (s ) - cn-2sn u(s ) -... - csu{s) - C;u{s) -

- Гп-1°п У + гп-2°п У(5) + ■■■ + гл°У + Го у(5). Выразив и(в), будем иметь

11 11

1 1 1

-Сои(в ) - Гп _1У (5 )--Гп-2У -...--— Г1У (5 )--— ГоУ (5),

° °п 2 °п I

или, в более удобном виде,

1 1

и° = Уо(°)-Сои(в)-ГоУ+ х 1

х(-С^и(5)- ^(5)) + ■■■ + -(-Сп-2и(°) - Гп-2У- (3)

-Гп- М5У

Соответствующая структурная схема приведена на рис. 2,а.

Таким образом, для реализации регулятора, не вносящего в передаточную функцию САУ дополнительных нулей, вместо 2п - 2 интеграторов потребуется п - 1. Более того, метод не предполагает выполнения дополнительных расчетов, поскольку оперирует только коэффициентами с, и г, полученными в результате синтеза регулятора.

Приняв выходы интеграторов управляющего устройства (рис. 2,а) за переменные состояния х,; уо(£) и у(£) - за входы, а и(в) - за его выход, на основе (3) можно записать для него следующие уравнения в форме Коши: 5Х ф = у о (5) - с0и(в) - г0 У(в),

вХ2 (5) = х1 (5) - си(в) - ГУ (в), 5Х3 (в) = Х2 (в) - С2и(в) - Г2У(в),

вХп-1 (в) = Хп-2 (в) - Сп-2и(в) - Гп-2У(в) ■

Если, согласно рис. 2,а, записать и(в) = хп ^(в)- гп _1 У(в) и переписать уравнения

5X1(5) = Уо (в) - Со хп_1(в) + СоГп_1У(в)-Го У (в), 8Х2(в) = Х1(в)-С1Хп_1(в) + С1Гп_1У(в)-Г1У(в), вХз (в) = Х2 (в) - С2Хп_1 (в) + С2Гп_1У(5) - ^(в),

°Хп-1 = Хп-2 (5°) - Сп-2Хп-1(5) + Сп-2Гп- 1У(5) - Гп-2У(5) то векторно-матричная форма регулятора примет вид канонической формы наблюдаемости [9]:

5хя (5) = Ая хя (5) + Вяия (5),

и(5) = хк (5) + йяия (5),

где

XR = [Xl(S) X2(S) ■■■ Xn-2(S) Xn-l(S)] ,

Bc

У o(s ) У (s )

c1rn—1 r

n-2 ( 0 0

-co —С

—С,

n—3

n—2

В этих условиях передаточные функции регулятора и префильтра для структуры САУ, изображенной на рис. 1,а, примут следующий вид:

Hr (S ) =

HF (s ) =

R(s) _ 59,45s2 + 441,4s + 997 C(s) " s2 — 12,45s + 295,7

1 1

CR =[0 0 ■■■ 0 1],

о« =[0 -гл_1].

Аналогичный результат получен в [7]. Пример. Синтезируем регулятор для объекта, описываемого передаточной функцией

4(s ) =

B(s) A(s )

s2 + 5s +100 s3 + 3s2 + 3s +1

R(s) 59,45s2 + 441,4s + 997 ' Согласно (3), уравнение управляющего устройства будет следующим: 1

u(s ) = 4 ( у0 (s) - 259,7u(s) - 997у (s)) + s2

1

+- (12,45u(s) - 441,4y(s)) - 59,45y(s)

или в векторно-матричной форме (4) с учетом (5):

Уo(s)" =г0 -295,7" у (s ) ' R 1 12.45

(5)

С учетом (1) degR(s) = 2, degC(s) = degR(s) = 2, degD(s) = 5.

Примем полином D(s) вида D(s) = (s + 5)5.

xr (s ) = Bc

X1(s ) X2(s ) 1 —997 0 —441,4

; Ur (s) =

; Cr =[0 1]; Dr =[0 —59,45].

a)

59,45

y,(s)

—

T

59,45

1 1

s s

441,4

y0(s) A(s)

997 *■

441,4

©<-—12,45

1

295,7 •*■

997 |-KB' u(s )

s2 + 5s +100

s3 + 3s2 + 3s+1

y(s)

б)

—295,7

У0(s )

12,45

1 kí

s

u(s )

s2 + 5s +100

s3 + 3s2 + 3s +1

У (s )

—997 —441,4 —59,45

Рис. 2

в)

На рис. 2,б приведена структура, аналогичная изображенной на рис. 1, а, реализованная на интеграторах (префильтр и регулятор представлены в управляемой (нормальной) канонической форме), а на рис. 2,в - структура, реализованная по предложенному методу согласно (5). Переходные процессы при единичном ступенчатом управляющем воздействии в обеих системах идентичны (рис. 3).

У(«)

Рис. 3

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

При использовании предложенного метода надо иметь в виду, что степень полинома 0(5) должна быть больше или равна степени полинома Я(в).

Астатические системы. В ряде случаев может быть необходимо обеспечение астатиз-ма по управляющему воздействию и возмущениям. Обычно это обеспечивают введением соответствующего количества интеграторов в контур управления [10]. Однако при таком подходе наличие префильтра (рис. 1) с передаточной функцией 1/Я(в) или 1/0(5), равно как и его кажущееся отсутствие (рис. 2, а) не позволяют получить астатизм по задающему воздействию.

В этом случае необходимо использование структур, приведенных на рис. 4. Наличие полинома 0(э) [11, 12] позволяет выбором соответствующих значений коэффициентов я, обеспечить порядок астатизма у< ЬедЯ(5)-1 по задающему воздействию. Степень полинома 0(5) при этом Ьед0(5) < ЬедЯ(5)-1.

Использование предложенного подхода упрощения структуры управляющего устройства возможно и в этом случае:

C(s)u(s) = R(s)A(s), u(s) =

R(s ) C(s )

Л(5),

A(s) = уС (s) - y (s), y0 (s) = Q(sl Ус (s),

R(s )

C(s)u(s) = R(s) f Q(s\ y0(s) - y(s) 1 = Q(s)y0(s) - R(s)y(s).

Я(5

yo(s ) Q(s) y0(s )

R(s )

R(s ) C(s )

u(s )

B(s ) A(s )

y(s )

а)

y0{s) -► Q(s ) y0(s ) —*(

N(s )

т.

А(5 )

У (s )

Я(5)

С(5)

Рис. 4

б)

Если принять degC(s) = degR(s) = n -1, degQ(s) = degR(s) -1 = n- 2, то

(sn-1 + cn_2sn-2 +... + C^s + c0 ) u(s) = = (Qn-2sn-2 +... + qs + Яс )Ус (s) --( rnS-1 + Гп-2sn-2 + ... + ris + Гс ) У (s), sn- 1u(s) = Яп- 2sn- 2 Ус (s ) + ... + Я^Ус (s) +

+Яс yc(s)- rn-isn-1 y (s )- rn-2sn-2 y (s ) -... - rsy (s) -

-Гс y (s) - Cn-2sn 2u(s ) -... - Cisu(s ) - Ccu(s),

1 1

u(s ) = Яп-2- yo(s ) +... + y0(s)+ s sn 2

1 1 1

+Я0—У0(s ) - гп-1У s - rn-2- y(s ) "... ~r1— y(s ) -

sn 1 s sn 2

1

1

■Го— У (s) -C—-2- u(s)-...- C u(s) -Co— u(s )

Или в более удобной форме: u(s ) = Гп (s ) + 1

+ -(Яп-2Уо (s) - Гп-2У(s) - Cn-2u(s)) + ... +

2 (ЯУ o(s )- Г У (s )- Cu(s)) +

(6)

^^пт (Яо У о (5) - Го У (5) - Сои(5)).

При этом V коэффициентов д,, рассчитанных из условия равенства нулю V коэффициентов ошибки (см., например, [11]), будут равны коэффициентам г, для любой из структур рис. 4:

д = г , I = о, V -1,

поэтому никаких дополнительных вычислений при синтезе регуляторов, обеспечивающих астатизм у-го порядка к задающему воздействию, вообще говоря, также не требуется.

Для переменных состояния х, входных сигналов у0(5) и у(в), выходного сигнала и(в) регулятора на основе (6) можно записать следующие уравнения: 5X1(5) = %Уо (5) - сои(8) - гоУ (5) -5Х2 (5) = 91Уо (5) + XI (5) - си(5) - /1^(5) -5Хз (5) = Ц2Уо (5) + Х2 (5) - С2и(5) - /^(5) -

5Х„-1 (5) = qn-2Уo (5) + Х„_2 (5) - Сп^и(5) - Гп^У(5) -или с учетом и(5) = ^^5) - Гп_У(5) (см. схему на рис. 2, а):

5X1 (5) = qoУo (5) - Сохп_1 (5) + СоГп_ 1У(5) - ^(5) -5X2 (5) = ^Уо (5) + X! (5) - С^- (5) + ^-^(5) - Г,У(5) -5Хз (в) = q2Уо(s) + X2(s) - 02X^^(3) + С2Гп-1У(в) - Г2У(в)-

5Xn_1 (5) = qn_2Уо (5) + Xn_2 (5) + Сn_2fn_1У(5) - Г^У^) и для системы (4):

X« =[ Х1(5) х2(5) ... 2(5)

" о о

xn- 1(s)]

Ur (s) =

Уо (s ) У (s)

1 о

0 0

0 1

-оп

-о.

-On

-о,

n-3

n-2

Bc

D,

Q0 Q1

Qn-3 Qn-2

-O0rn-1 - r0

Cr =[0 0 ... 0 1],

= [0 -rn-1].

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

В первом случае процедура синтеза будет полностью аналогична описанной выше. Второй рассмотрим подробнее.

Неминимальные астатические регуляторы. Использование неминимальных регуляторов [8] позволяет получить возможность назначения некоторых коэффициентов регулятора. Необходимым условием для этого является наличие полинома C(s) и выполнение условия degR(s) > degA(s) -1. Количество k назначаемых коэффициентов будет определяться выражением:

k = degR(s) - (deg A(s) -1). (7)

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

Пример. Пусть необходимо обеспечить астатизм по задающему воздействию и по возмущению объекта, структура которого приведена на рис. 5. Задающее воздействие является линейно нарастающим сигналом у0(t) = 1-1, возмущение - ступенчатым сигналом f(t) = 1(t). При приложении ступенчатого задающего воздействия необходим апериодический переходный процесс длительностью не выше 0,2 с.

u(s )

220

a 1

0,1s +1

f (s)

0,1s

У (s ) i—►

Рис. 5

Передаточная функция объекта имеет вид 2200

H (s ) = 2200-■

s2 + 10s +100

Используя возможности неминимальных регуляторов и учитывая высшую из требуемых степеней астатизма v = 2 (для обеспечения нулевой ошибки при линейно нарастающем задающем воздействии), примем

degR(s) = degA(s) — 1 + v = 3 , тогда, согласно (7), можно обеспечить произвольное задание двух коэффициентов регулятора: k = degR(s) — (deg A(s) — 1) = 3 — 1 = 2.

В нашем случае для обеспечения требуемых астатических свойств зададим c0 = c1 = 0.

Примем для технически реализуемого регулятора degC(s) = degR(s) = 3 , тогда, согласно (1), degD(s) = 5. Задавшись для обеспечения необходимого быстродействия значением среднегеометрического корня Q0 = 50

полинома Ньютона D(s) = (s + Q0)5 и решив уравнение синтеза, получим передаточную функцию регулятора:

R(s) 10,23s3 + 557,3s2 +14204,5s +142045

Hr (s ) =

C(s )

Матрицы векторно-матричного описания управляющего устройства будут иметь следующий вид:

1

xR (s) =

Be

^i(s)" x2(s) *a(s)_ 142045 14204,5 0

Ur (s) =

У o(s) У (s)

0 0 0 "

; Ar = 1 0 0

0 1 -240

-142045 -14204,5 -240 • (-10,23) - 557,3

Ск =[о о 1]; Ок =[о -Ю-23].

Структурная схема САУ приведена на рис. 6,а.

Введение ошибки регулирования Д(5) = У0 (5) - У(5) позволяет придать САУ окончательный вид (рис. 6, б). Переходные

процессы по ошибке регулирования в САУ при приложении задающего воздействия вида y0(t) = 1-1 в нулевой момент времени и возмущающего воздействия в виде ступенчатого сигнала f(t) = 1(t) в момент времени 1 с приведены на рис. 7.

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

У 0 (s)

Ф

-142045

-КБ

-14204

г 1

142045 14204 -240

1

— -н

s

-557,3

u(s)

220

-10,23

k 1

0,1s +1

f (s)

1

0,1s

У (s)

>—►

а)

14204

У 0(s)

a(s)

142045

1

— —и

s

-240

Xu(s)

-557,3

220

-10,23

k 1 kf

* 0,1s +1

f (s)

1

0,1s

У (s)

б)

Рис. 6

Д. рад

/

li t С

0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6

Рис. 7

Выводы

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

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

Использование предложенного способа позволяет проводить редукцию регулятора и в случае синтеза следящих систем, в которых требуемая степень астатизма по управлению обеспечивается совокупным использованием интеграторов в контуре управления и соответствующих коэффициентов полинома Q(s) пре-фильтра.

Список литературы

1. Волгин Л.Н. Элементы теории управляющих машин (метод полиномиальных уравнений в задачах синтеза систем автоматического управления с цифровыми вычислительными машинами). - М.: Советское радио, 1962. - 164 с.

2. Гайдук А.Р. Синтез систем управления по передаточным функциям // Автоматика и телемеханика. - 1980. -№ 1. - С. 11-16.

3. Kucera V. Discrete Linear Control. The Polynomial Equation Approach. - Prague: Academia, 1979.

4. Волгин Л.Н. Оптимальное дискретное управление динамическими системами. - М.: Наука. Гл. ред. физ.-мат. лит., 1986. - 240 с.

5. Скворцов Л.М. Синтез линейных систем методом полиномиальных уравнений // Автоматика и телемеханика. - 1991. - № 6. - С. 54-59.

6. Grimble M.J., Kucera V. Polynomial Methods for Control Systems Design. Springer-Verlag, 1996.

7. Гайдук А.Р. Теория и методы аналитического синтеза систем автоматического управления (полиномиальный подход) - М.: ФИЗМАТЛИТ, 2012. - 360 с.

8. Тютиков В.В., Тарарыкин С.В. Робастное модальное управление технологическими объектами. - Иваново, 2006. - 256 с.

9. Phillips C., Harbor R. Feedback Control systems. Prentice hall, Inc. 2000.

10. Мееров М.В. Синтез структур систем автоматического регулирования высокой точности. - М.: Наука, 1967. - 424 с.

11. Гайдук А.Р. Основы теории систем автоматического управления. - М.: УмиИЦ «Учебная литература», 2005. - 405 с.

12. Бесекерский В.А., Попов Е.П. Теория систем автоматического регулирования. - М.: Наука. Гл. ред. физ.-мат. лит., 1966. - 992 с.

References

1.Volgin, L.N. Elementy teorii upravlyayushchikh mashin (metod polinomial'nykh uravneniy v zadachakh sinteza sistem avtomaticheskogo upravleniya s tsifrovymi vychisli-tel'nymi mashinami) [Elements of machine control theory (method of polynomial equations in problems of automatic control system synthesis with digital computing machines]. Moscow, Sovetskoe radio, 1962. 164 p.

2. Gajduk, A.R. Sintez sistem upravleniya po pereda-tochnym funktsiyam [Transfer function synthesis of control machines]. Avtomatika i telemekhanika, 1980, no. 1, pp. 11-16.

3. Kucera, V. Discrete Linear Control. The Polynomial Equation Approach. Prague, Academia, 1979.

4.Volgin, L.N. Optimal'noe diskretnoe upravlenie di-namicheskimi sistemami [Optimal discrete control of dynamic machines]. Moscow, Nauka, Glavnaya redaktsiya fiziko-matematicheskoy literatury, 1986. 240 p.

5. Skvortsov, L.M. Sintez lineynykh sistem metodom polinomial'nykh uravneniy [Synthesis of linear systems by the polynomial equations method]. Avtomatika i telemekhanika, 1991, no. 6, pp. 54-59.

6. Grimble, M.J., Kucera, V. Poly-nomial Methods for Control Systems Design. Springer-Verlag, 1996.

7. Gayduk, A.R. Teoriya i metody analiticheskogo sin-teza sistem avtomaticheskogo upravleniya (polinomial'nyy podkhod) [Theory and methods of analytical systhesis of automatic control (polynomial approach)]. Moscow, FIZMATLIT, 2012. 360 p.

8. Tyutikov, V.V., Tararykin, S.V. Robastnoe modal'noe upravlenie tekhnologicheskimi ob"ektami [Robust modal control of industrial facilities]. Ivanovo, 2006. 256 p.

9. Phillips, C., Harbor, R. Feedback Control systems. Prentice hall, Inc. 2000.

10. Meerov, M.V. Sintez struktur sistem avtomati-cheskogo regulirovaniya vysokoy tochnosti [Synthesis of high precision automatic control system structures]. Moscow, Nau-ka, 1967. 424 p.

11. Gayduk, A.R. Osnovy teorii sistem avtomati-cheskogo upravleniya [Theoretical fundamentals of automatic control systems]. Moscow, UmilTs «Uchebnaya literature», 2005. 405 p.

12. Besekerskiy, V.A., Popov, E.P. Teoriya sistem av-tomaticheskogo regulirovaniya [Theory of automatic control systems]. Moscow, Nauka, Glavnaya redaktsiya fiziko-matematicheskoy literatury, 1966. 992 p.

Тютиков Владимир Валентинович,

ФГБОУВО «Ивановский государственный энергетический университет имени В.И. Ленина»,

доктор технических наук, профессор, проректор по научной работе,

телефон (4932) 41-50-24,

e-mail: tvv@ispu.ru

Tyutikov Vladimir Valentinovich,

Ivanovo State Power Engineering University,

Doctor of Engineering Sciences (Post-doctoral degree), Professor, Vice-Rector for Research, tel. (4932) 41-50-24, e-mail: tvv@ispu.ru

Копылова Лариса Гзннадьевна,

ФГБОВО «Ивановский государственный энергетический университет имени В.И. Ленина», кандидат технических наук, доцент кафедры электроники и микропроцессорных систем, телефон (4932)26-97-55, e-mail: klg@eims.ispu.ru

Kopylova Larisa Gennadievna,

Ivanovo State Power Engineering University,

Candidate of Engineering Sciences (PhD), Associate Professor of the Electronics and Microprocessor Systems Department, tel. (4932) 26-97-55, e-mail: klg@eims.ispu.ru

Тихомирова Ирина Александровна,

ФГБОУВО «Ивановский государственный энергетический университет имени В.И. Ленина»,

старший преподаватель кафедры электроники и микропроцессорных систем,

телефон (4932)26-97-55,

e-mail: tia@eims.ispu.ru

Tikhomirova Irina Aleksandrovna,

Ivanovo State Power Engineering University,

Senior Lecturer of the Electronics and Microprocessor Systems Department, tel. (4932) 26-97-55, e-mail: tia@eims.ispu.ru

Шляцкая Елена Михайловна,

ФГБОУВО «Ивановский государственный энергетический университет имени В.И. Ленина», магистрант,

телефон (4932)26-97-55,

e-mail: lenashlyatskaya@mail.ru

Shlyatskaya Elena Mikhailovna,

Ivanovo State Power Engineering University,

Master Course Student,

tel. (4932) 26-97-55,

e-mail: lenashlyatskaya@mail