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

УДК 681.5.015

Марковские параметры и передаточные матрицы многомерных управляемых систем

А.Р. Гайдук1, И.А. Каляев1, С.Г. Капустян1, В.Н. Рябченко2 1ФГАОУВПО «Южный федеральный университет», г. Таганрог, Российская Федерация 2ОАО «Федеральная сетевая компания ЕЭС», г. Москва, Российская Федерация E-mail: [email protected], [email protected], [email protected]

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

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

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

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

Markov Parameters and Transfer Matrixes of Multidimensional Control Systems

A.R. Gaiduk1, I.A. Kalyaev1, S.G. Kapustyan1, V.N. Ryabchenko2 1 Southern Federal University, Taganrog, Russian Federation 2Public Corporation «Federal Grid Company of Unified Energy System», Moscow, Russian Federation E-mail: [email protected], [email protected]

Abstract

Background: The analytical correlations describing the connection of Markov parameters, in particular, of multidimensional systems with their order, transfer matrixes and measured input and output sizes are almost not considered in the known literature to the extent which is necessary for carrying out the identification. Thus, the urgent problem is to get analytical correlations, which demonstrate the relations between the quantative characteristics of the controllable system and measured values as well as to determine such quantative characteristics of indefinite system, which entirely have all its features and allow to synthesize the effective algorithms of controlling this system.

Materials and methods: The required mathematical correlations for continuous and discrete multidimensional systems are determined on the basis of decomposition in series of the attached matrix of a characteristic matrix and the expressions, determining values of output variables of discrete dynamic systems.

Results: The analytical correlations which connect Markov parameters of the full continuous or discrete multidimensional system with their order and with the factors of a transfer matrix are received. Connection between Markov parameters of discrete controlled systems with values of output variables is established at discrete control input. Conclusions: The received correlations can be applied to identify the controlled continuous or discrete dynamic systems. They allow to find both the system order and the transfer matrix coefficient by measurement results of scalar control input and reactions of the full multidimentional controlled discrete system.

Key words: controlled system, Markov parameters, system order, transfer matrix, identification, variable value.

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

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

тивные алгоритмы управления этой системой. К таким количественным характеристикам многомерных управляемых систем относятся, в частности, марковские параметры [2].

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

Ниже на основе формулы для матрицы, присоединенной к характеристической матрице [5, с. 88], устанавливается связь марковских параметров полных непрерывных и дискретных управляемых систем с их порядком и коэффициентами их передаточных матриц. Показана также связь марковских параметров полных дискретных управляемых систем с реакциями этих систем на дискретные управления.

Непрерывные управляемые системы. Рассмотрим линейную непрерывную динамическую управляемую систему с несколькими входами и несколькими выходами, уравнения которой в переменных состояния имеют следующий вид:

х = Ах + Ви, у = Сх + Ои, (1)

где х е Яп - вектор состояния размерности п; и е - вектор управлений; у е Я1 - вектор выходных переменных системы; А, В, С, О - числовые матрицы соответствующих размерностей.

Марковские параметры динамической

системы (1) - это величины, определяемые выражениями

ц0 = С„, цЦ; = С/Аи-1ВЦ, и = 1,2,3,..., (2)

где Сц - элементы матрицы О; С/ и В - строки и столбцы матриц С и В соответственно; / = 1,I;

Ц = Гд [2, 4].

Часто определяются матрицы марковских параметров Ми = [^ ] е Я1 хд по формулам

М0 = О, Ми = САи-1В, и = 1,2,3,.... (3)

Поставим задачу определения связи марковских параметров системы (1) с коэффициентами ее передаточной матрицы [6]:

(р) = С(рЕ - А)-1 В + Оо =

= А-1(р){Сас1](рЕ - А)В + Оо А(р)}, где А(р) = рЕ - А) или

А( р) = апр" +ап-1р"-1 +... + ар +а0, а п = 1. (5)

Применяя формулу для присоединенной матрицы общего вида [5, с. 88; 6, с. 231] к матрице (рЕ - А), получим следующее выражение:

(6)

(4)

а^(рЕ - А) = апЕр"-1 + (апА + ап-1Е)р"-2 + ... + (а пАп-1 + ап-1 Ап-2 + ... + а1Е). Подставляя соотношение (6) в равенство (4), получим

(р) = А-1( р){С[апЕрп-1 + (апА + ап-1Е )рп-2 +... ... + (апАп-1 + а п-1Ап-2 +... + а1Е )]В + О0 А( р)} или, раскрывая скобки,

(р) = А-1(р){апСВрп-1 + (апСАВ + ап-1СВ)рп-2 +...

... + (апСАп-1В + ап-1САп-2В +... + а1СВ) + М0 А( р)}.

Таким образом, с учетом обозначений (3) передаточную матрицу (4) многомерной системы (1) можно представить следующим образом:

(р) = А~\ р){апМ1рп-1 + (апМ2 +ап-1М1)рп-2 +... ... + (апМп-2 + ап-1Мп-3 +... + а1М1) + М0 А( р)}.

(7)

Из (7) следует, что передаточная матрица Wyu(p) многомерной системы (1) полностью определяется п матрицами М° (3) марковских параметров ци, и = 0, п -1 (2) и п коэффициен-

(8)

тами а/, / = 0, п -1 характеристического полинома (5) этой системы, так как ап = 1.

С другой стороны, учитывая, что максимальная степень р в (5) равна п, передаточную матрицу Wyu(p) (4) можно записать следующим образом:

Wyu (р) = А-1( р ){Впрп + Вп-1рп-1 +...

... + В1р + В0 + В-1р-1 + В-2р-2 +...},

где Вр= [врр] - I х д -матрицы постоянных коэффициентов.

Сравнивая выражения (7) и (8), заключаем, что матрицы Вр определяются выражениями:

Вп =апМ0, Вп-1 =ап-1М0 +апМ1, Вп-2 = ап-2М0 + ап-1М1 + апМ2, ... или

Вп-р=]Гап-р+М, (9)

Ц=0

где р = 0,1,2,..., п -1, п, п +1,...

В частности, так как ап = 1, то, например, матрица

В0 = а0М0 + а1М1 +... + ап-1Мп-1 + Мп.

Подчеркнем, что в общем случае, в соответствии с выражениями (4) и (5), элементами рр матриц Вр при р = 0, п фактически являются

коэффициенты числителей скалярных передаточных функций Wy.u.(р), которые, в свою

очередь, являются элементами передаточной матрицы Wyu(p) (4). Эти функции в общем случае имеют следующий вид:

в0 + р1р+р2р2 +...+РП-у-1 +ряр" Ш (Р) = р" -Р" ^ и-1 Т . (10)

ао + а1Р + а2Р -

-ап-1р"-1 + рл

Подчеркнем, что если ^у.и.(р) ф 0, то по

крайней мере рр ф 0, ре[0, л]. Другими словами, в равенстве (8) при ре [0, л] по крайней мере одна из матриц Вр ф0, а при р = -1,-2,...

все матрицы Вр = 0 .

Из выражений (9) при ал = 1 вытекает (л + 1)-мерная система равенств:

0 0 0 • 0 м0 " а0 " Вл '

0 0 0 ■ м0 м1 а1 Вл-1

0 0 0 м1 м2 а2 = Вл-2

0 м0 м1 . . мл-2 мл-1 ал-1 В1

м0 м1 м2 . . мл-1 мл _ _ 1 _ . В0 _

(11)

а при любом 5 = 0, 1, ... бесконечная система равенств:

м5+1 м5+2 .. м5+л ' " м5+л+1"

м5+2 м5+3 .. м5+л+1 а0 м5+л+2

а1

м5+л м5+л+1 .. м5+2л-1 м5+2л1

м5+л+1 м5+л+2 .. м5+2л+1 _ал-1_ м5+л+1

(12)

При этом последнюю можно рассматривать как систему уравнений относительно коэффициентов а,-, / = 0, л -1.

Рассмотрим матрицу, составленную из первых л строк матрицы в левой части систе-

получим равенства:

мл =

5 +1

в+2

м м

в+2

в+3

м5

м

5+л+1

СА5В СА5+1В СА5+1В СА5+2В

СА5+л-1В СА5+лВ ••■ СА5+2л-1В

М5+л

5+л+1

м

5 + 2л

м

СА5+л-1В СА5+лВ

(13)

или

мл =

С СА

СА

л-1

А51 В АВ ... Ал-1В I.

(14)

Таким образом, матрица Мл (13) при в ф 0

равна произведению трех матриц: матрицы наблюдаемости, матрицы Л5 и матрицы управляемости системы (1) [7]. Следовательно, матрица МЧ имеет обратную матрицу при б ф 0, если только Се1 А ф 0 , а система (1) является полной [7]. Так как Се1 А = а0, то, если Се1 М5ф0 = 0, а

система (1) полная, следовательно, а0 = 0. В этом случае справедлива следующая «укороченная» система:

5+3

5 + 4

м м

5 + 4

5+5

м

м5+л

5 +л+1

м5

м

5 + л+1

м

5 + 2л

1л-1.

5+л+1

5 +л+2

м

5 + 2л + 1

(15)

Матрица этой системы имеет обратную, если управляемая система (1) полная. Если же полнота системы (1) неизвестна, то для определения коэффициентов а,, / = 0, л -1, очевидно, целесообразнее использовать систему, соответствующую первым л строкам системы (12) при в = 0, т.е. с матрицей м^ (13).

Подчеркнем, что уравнения (12), (15) не включают ни матрицу м0, равную матрице О, ни параметры ц0 = d/j, которые описывают

прямые безынерционные каналы связи между входом и выходом системы (1).

Важно также, что в соотношения (11), (12) и (15) вместо матриц мм и Вр можно подставлять соответствующие марковские параметры и коэффициенты р?-, где и, р = 0, л , при какой-либо (любой, но при одной и той же) паре индексов /, j, причем 1 < / < I и 1 < j < д . Например, для некоторой системы при л = 3, 5 = 0, / = 2, а j = 3 можно записать следующие равенства:

В результате " 0 0 0 0 1 ^23 а0 >¡31

0 0 0 ^23 1 Н-23 а1 р23

0 ^23 1 ^23 а2 Р123

0 >23 1 ^23 2 ^23 3 ^23 _ _ 1 _ &

(16)

^23 ^23

2 3 ^23 ^23

3 4 ^23 ^23

^23

4 ^23

5 ^23

4 ^23

5 ^23

6 ^23

Другими словами, соотношения (11), (12), (15) справедливы как для всей многомерной системы, так и для каждого ее канала вход-выход. Соотношения указанного типа для одномерного случая приведены в [2, 6].

Утверждение 1. Если передаточная функция (р) (10) канала и! ^ у, полной

системы (1) порядка л не равна тождественно нулю, то среди л + 1 первых марковских пара-

метров , ие [0, п], этого канала по крайней

мере один не равен нулю.

Доказательство. Равенство (11) имеет место для каждого канала Uj ^ у, системы (1)

порядка п . При этом вектор правой части каждого такого равенства, соответствующего паре /, Ц, включает все коэффициенты рР-, р = 0, п ,

числителя не равной нулю передаточной функции Wy.u. (р) (10), а матрица в его левой

части - все марковские параметры .¡и, и = 0, п , данного канала u¡ ^ у, . Предположим, все марковские параметры ЦЦ, и = 0, п , некоторого канала u¡ ^ у, системы (1) равны нулю. В

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

р = 0, п , будут равны нулю. Но по условию функция Wy¡Uj (р) не равна нулю, т.е. по крайней мере один из коэффициентов рР- не равен

нулю. Полученное противоречие доказывает утверждение 1.

Соотношения (11), (12) имеют большое значение при решении задачи идентификации полных динамических систем. Практически, по значениям (2п + 2) -г- (2п + 6) марковских параметров можно оценить как порядок системы, так и коэффициенты ее уравнения вход-выход. Допустим, для некоторого канала u¡ ^ у,

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

параметров .у,-, и = 0,1, 2,. , и величины

Уц (з, п) =

з+1 +2

з . /

э+2 з+3

Ц

ив+п . 5+п+1

в+п

мв+п+1 .Ц

мв+2п-1 .Ц

(17)

при некоторых з = 0,1,. и п = 1.2..... Тогда

справедливо следующее утверждение.

Утверждение 2. Если при некоторых з = 0,1,... и всех значениях п = 1.2,..., п * -1 величина Уц (з, п*) имеет некоторые (в том числе, и нулевые) значения, при п = п * Уц(з,п*) * 0, а при всех п = п * +1, п * +2,... Уц(з,п) = 0, то порядок рассматриваемой полной системы равен п *, а числа а0, а1,..., ап*-1, определяемые одной из систем (12) при п = п*, являются коэффициентами ее характеристического полинома

(5). При этом числа р0, р1у,..., Рп , определяемые одной из систем (11), являются коэффициентами

числителя передаточной функции (10) ее канала

uj ^ у -/.

Покажем справедливость утверждения 2 на численном примере.

Пример. Система, уравнения которой имеют вид

" 0 1 0 " "0"

х = 0 0 1 х + 0 u1,

-1 0 -2 1 (18)

"3 0 0" "2"

у = х +

2 1 0_ _0_

согласно (2), имеет следующие марковские параметры. По каналу u1 ^ у1: = 2 , м11 = 0 ,

Ц31 = ^ М141 =-6, м51 = 12, м61 =-27,

.1 = 0 , м71 = 60,

каналу

м21 =0, м21 = 9,

м81 = -132, м91 = 291; по U1 ^ у2 : .21 = 0, м121 = 0 , .21 =1, .21 = 0 , .21 = -1, .21 = 2 , .21 = -4 , м21 = -20 .

Предположим, ни порядок, ни параметры системы (18) неизвестны, а значения указанных марковских параметров известны. В этой ситуации примем в выражении (17) / = 1, Ц = 1, 5 = 0, а п = 1, 2, ..., 5. Тогда имеем: Уц(0; 1) = .11 = 0, Уц(0; 2) = 0, Уц(0; 3) = -27, У11(0;4) = У11(0;5) = 0. Следовательно, согласно утверждению 2, порядок системы п = 3. При этом решение системы (15) с заменой М° на .и при п = 3 дает значения коэффициентов а0 = 1, а1 = 0, а2 = 2 .

Нетрудно убедиться, что аналогичные результаты получим, полагая = 1, Ц = 2, з = 2, а п = 1, 2, ..., 5 по-прежнему.

Решение системы (11) в форме матриц Вр при п = 3 имеет вид

В3 =

Подставляя полученные матрицы в формулу (8), будем иметь

1

"2" , В2 = "4" , В1 = "0" , В0 = "5"

_0_ _0_ _ 1_ _2_

Wyu (р) = — р -

или

Wy1 u1( р) = 2-

Wy2 u1(р) =

2р2 +1

3

2 3 4 2 0 5

р' + р + 1 р +

0 0 2

р3 + 2р р + 2

р3 + 2р2 +1

Наконец, пользуясь известными соотношениями перехода от передаточных функций к уравнениям в переменных состояния в канонической управляемой форме [6], нетрудно записать уравнения рассматриваемой управляемой системы (18) на основе полученных выражений для Wy р) и Wy2щ(р). Конечно, в общем слу-

/2 u1

чае, найденные рассмотренным способом урав-

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

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

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

Вик, ук = Охк + Ьик, (19)

к, ик, ук - векторы дискретных переменных состояния, управления и управляемых переменных.

Как видно, структура уравнений (19) полностью аналогична уравнениям (1). Поэтому для полных дискретных систем (19) справедливы все приведенные выше соотношения (2) -(17), а также утверждения 1 и 2 при соответствующих заменах матриц А, В, С, й, Ми на матрицы А, В, С, й, ММи и чисел ц}} на Д)у.

Например, марковские параметры канала и, ^ уI некоторой системы (19) опреде-

у10" 0 ц 11 Д 02 "и10"

У20 _ ~ 0 .Д 21 ~0 Д22 _ _и20 _

хк+1 = Ахк где хк, и

-^11 " " ~ 1 ц 11 ~1 " Д12 "и10" + " ~ 0 Д11 Д 02 ~и11~

у21. ~ 1 _Д 21 ~1 Д22 . _и20. ~ 0 .Д21 ~0 Д22 . и21.

"1

ляются выражениями

дЛ = Аи-1В1, и = 1,2,3,.... (20)

Для этой же системы третьего порядка

Д ° = ¿и

у к = ^ ММк_рир = ММ ки0 + ММк+... + ММ и_1 + М0ик.

р=0

(24)

Например, при I = 2, д = 2 и к = 0, 1, ... из этого выражения вытекают следующие равенства:

(25)

С точки зрения задачи идентификации управляемых систем необходимо, чтобы по результатам измерения входных воздействий и выходных переменных системы можно было бы найти значения параметров системы. В уравнениях (25) параметрами системы являются ее марковские параметры Д)}. Однако из выражений

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

В то же время уравнения (25) принимают совершенно другой вид, если на входы системы подавать скалярные воздействия сначала

вида и1 =[и1к 0]г, а затем иI =[0 и2к] . В

этом случае при тех же значениях I = 2, д = 2 и к = 0, 1, ... будем иметь:

0 0 0 0 0 0 Д 23 ~0 1 Д 23 ~1 Д23 х 0 хх 1 "р 23" р 23 -^10 у20 = ~ 0 Д11 ~ 0 Д21 и10 , у11 у21 = ~ 1 Д11 ~ 1 Д21 и10 + ~ 0 Д11 ~0 Д21

0 Д 23 ~ 0 ~ 1 23 I-123 1 Д 23 Д 23 ~2 Д23 Д23 . хХ 2 _ 1 . й3 р 23 (21) -^120 _ У20. = - 0 Д12 - 0 _Д 22 . и20 , У121 _У21_ = -1 Д12 -1 _Д 22 и20 + - 0 Д ^ - 0 _Д 22

'11 -

21 >

1 Д 23 2 Д 23 -3 1 Д 23 Х 0 Г~ 4 1 Д 43

2 Д 23 3 Д 23 ~4 Д 43 ХХ 1 = _ ~ 5 Д 23

3 Д 23 ~4 Д 23 ~5 Д 23 ХХ 2 ~ 6 Д 23

(26)

где

Ук

где ар, рр - коэффициенты передаточной функции

Р0 + Рк + Р2 +Р3^3

- векторы выходных переменных исследуемой системы при действии на ее входах скалярных управлений ик , } = 1,2. Из выражений (26) вытекают следующие рекуррентные соотношения:

И/У,} (z) =

Х г + а, г2 + г3

(22)

*0 + хх 1г + хх2г ' ^

Весьма существенным свойством дискретных систем является то, что их реакции на дискретные воздействия полностью выражаются через их марковские параметры. Действительно, как показано в работах [7, 8], реакция системы (19) на дискретное векторное воздействие ик определяется выражением

~0

Д у

0

Д 2 }

у10

у2 0

/ и

} 0'

-1 1 Д 1} у11 Г ~ 0 Д V

~1 Д2 }. К _ у21. ~ 0 Д 2 }

7у1

/и

} 0

к_1

(23)

ук = САкХ0 + Х ОАк_р-1ЕВир + йик.

р=0

При нулевых начальных условиях с учетом обозначений ММ0 = й, Ми= (СААи_1ЕЕ, аналогичных (3), из выражения (23) выводим

(27)

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

На основе выражений (24) и (27) можно сформулировать утверждение 3, которое имеет существенное значение для идентификации дискретных управляемых систем.

Утверждение 3. Для экспериментального определения на основе соотношения (24) всех марковских параметров Д- некоторого канала

иу ^ уI многомерной дискретной системы достаточно, чтобы хотя бы при одном значении к > 0 выполнялось условие и'к ф 0 , у е [1, д].

Для доказательства этого утверждения достаточно записать выражение (24) в раскрытой форме, начиная со значения к, при котором выполняется условие и]к ф 0 , у е [1, д].

Заключение

Полученные соотношения связывают марковские параметры полных, линейных, непрерывных и дискретных многомерных управляемых систем с порядком и коэффициентами передаточных функций отдельных каналов этих систем. Эти соотношения обратимы и позволяют по (2л + 1) марковским параметрам полной многомерной системы найти ее порядок, а также коэффициенты характеристического полинома. По значениям (л + 1) первых марковских параметров каждого канала этой системы можно найти коэффициенты числителя передаточной функции этого канала.

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

Исследование выполнено при поддержке РФФИ (грант № 12-08-90050-Бел_а).

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

1. Красовский А.А., Наумов А.И. Аналитическая теория самоорганизующихся систем управления с высоким уровнем интеллекта // Известия РАН. Теория и системы управления. - 2001. - № 1. - С. 69-75.

Гайдук Анатолий Романович, ФГАОУВПО «Южный федеральный университет», профессор кафедры систем автоматического управления, телефон (8634) 37-16-89, e-mail: [email protected]

Каляев Игорь Анатольевич, Научно-исследовательский институт многопроцессорных вычислительных систем имени академика А.В. Каляева ФГАОУВПО «Южного федерального университета», член корреспондент рАн, директор телефон (8634) 36-07-57, e-mail: [email protected]

Капустян Сергей Георгиевич,

Научно-исследовательский институт многопроцессорных вычислительных систем имени академика А.В. Каляева

ФГАОУВПО «Южного федерального университета»,

доктор технических наук, заведующий отделом,

телефон (8634) 31-54-94,

e-mail: [email protected]

Рябченко Владимир Николаевич,

ОАО «Федеральная Сетевая Компания ЕЭС»,

доктор технических наук, референт генерального директора,

телефон (495) 962-81-82,

е-mail: [email protected]

2. Гайдук А.Р. Алгоритмическое обеспечение самоорганизующихся регуляторов с экстраполяцией // Известия РАН. Теория и системы управления. - 2002. - № 3. - С. 56-63.

3. Гайдук А.Р., Капустян С.Г. Самоорганизующийся алгоритм действий интеллектуальных автономных роботов // Наука и образование на рубеже тысячелетий: сб. НИР. - М.: Учлитвуз, 2011. - Вып. 1. - C. 5-15.

4. Марковские параметры многомерных динамических систем управления / И.А. Каляев, А.Р. Гайдук, С.Г. Капустян, В.Н. Рябченко // Вестник ИГЭУ. - 2013. - № 1. -С. 45-50.

5. Гантмахер Ф.Р. Теория матриц. - М.: Наука, 1988.

6. Гайдук А.Р. Непрерывные и дискретные динамические системы. - М.: Учлитвуз, 2004.

7. Гайдук А.Р. Теория автоматического управления. -М.: Высш. шк., 2010.

8. Куо Б. Теория и проектирование цифровых систем управления: пер. с англ. - М.: Машиностроение, 1986.

References

1. Krasovskiy, A.A., Naumov, A.I. Analiticheskaya te-oriya samoorganizuyushchikhsya sistem upravleniya s vy-sokim urovnem intellekta [The Analytical Theory of Self-organizing Control Systems with High Intellect Level]. Izvestiya RAN. Teoriya i sistemy upravleniya, 2000, no. 6, pp. 37-41.

2. Gayduk, A.R. Algoritmicheskoe obespechenie samoorganizuyushchikhsya regulyatorov s ekstrapolyatsiey [Algorithmic Maintenance of Self-organizing Controllers with Extrapolation]. Izvestiya RAN. Teoriya i sistemy upravleniya, 2002, no. 3, pp. 56-63.

3. Gayduk, A.R., Kapustyan, S.G. Samoorganizu-yushchiysya algoritm deystviy intellektual'nykh avtonomnykh robotov [Self-organizing Actions Algorithm of Intellectual Independent Robots]. Nauka i obrazovanie na rubezhe tysyach-eletiy: sbornik NlR [Science and Education on Boundary of Millennia: Science and Research Collected Works]. Мoscow, Uchlitvuz, 2011, issue 1, pp. 5-15.

4. Kalyaev, I.A., Gayduk, A.R., Kapustyan, S.G., Ryab-chenko, V.N. Markovskie parametry mnogomernykh di-namicheskikh sistem upravleniya [Markov Parameters of Multidimensional Dynamic Control Systems]. Vestnik IGEU, 2013, issue 1, pp. 45-50.

5. Gantmakher, F.R. Teoriya matrits [The Theory of Matrixes]. Moscow, Nauka, 1988.

6. Gayduk, A.R. Nepreryvnye i diskretnye dinamiches-kie sistemy [Continuous and Discrete Dynamic Systems]. Мoscow, Uchlitvuz, 2004.

7. Gayduk, A.R. Teoriya avtomaticheskogo upravleniya [Automatic Control Theory]. Мoscow, Vysshaya shkola, 2010.

8. Kuo, B. Teoriya i proektirovanie tsifrovykh sistem upravleniya [Theory and Designing Digital Control Systems]. Мoscow, Mashinostroenie, 1986.