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

А

нализ и синтез систем управления

УДК 681.5.01

0 СВОЙСТВАХ РЕГУЛЯРНЫХ ДИНАМИЧЕСКИХ СИСТЕМ1

Л.А. Мироновский, Т.Н. Соловьева

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

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

ВВЕДЕНИЕ

В теории управления известны различные типы линейных динамических систем, отличающиеся расположением нулей и полюсов, видом частотных характеристик и другими свойствами [1—5]. В настоящей статье вводится в рассмотрение еще один тип систем — так называемые регулярные системы. Они характеризуются специальным видом матриц сбалансированного представления и регулярной структурной реализацией. Коэффициенты передаточной функции (ПФ) регулярной системы полностью определяются значениями ее ганкелевых сингулярных чисел (ГСЧ).

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

1 Работа выполнена при финансовой поддержке РФФИ (грант 11-08-00240).

1. ОПРЕДЕЛЕНИЕ РЕГУЛЯРНЫХ СИСТЕМ И УСЛОВИЯ РЕГУЛЯРНОСТИ

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

Определение регулярных систем удобно дать, воспользовавшись понятием сбалансированного представления. Напомним, что реализация устойчивой системы, заданной описанием в пространстве состояний

X = АХ + Ьы, у = еХ, (1)

где X е Яп — вектор переменных состояния; ы, у е Я1 — входной и выходной сигналы; А, Ь и е — постоянные матрицы соответствующих размеров; называется сбалансированной, если ее грамианы

управляемости и наблюдаемости диагональны и равны между собой [1—3]: = Wo = Ша§(ар ап).

Такие реализации существуют для каждой устойчивой линейной системы. Числа аг- 1 0 называются ганкелевыми сингулярными числами системы (сингулярными числами ее ганкелева оператора). Для систем с одним входом и одним выходом они совпадают с собственными числами vi кросс-грамиана системы (т. е. собственными числами ее ганкелева оператора), взятыми

го

по абсолютной величине: = | е А%се А*Ж.

о

Сбалансированную реализацию можно задать набором ганкелевых сингулярных чисел и компонент Ь вектора Ь описания (1). Если все ГСЧ различны, то элементы матриц А и с определяются формулами:

а0

ь, ь,

s I Sj а I + аj•

с

(2)

где s/ = ±1; I = 1, п — знаки собственных чисел

кросс-грамиана (^ = s/а/).

В частности, для диагональных элементов матрицы А имеем

а,, = —

£ 2 а '

(3)

Определение. Система (1) с попарно различными ГСЧ называется регулярной, если все диагональные элементы матрицы А ее сбалансированного представления одинаковы. ♦

Такой специальный вид системной матрицы обусловливает простоту структурной реализации

Рис. 1. Пример реализации регулярной системы

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

А =

сбалансированного представления.

-1 2 -3 2 -1 1 3 1 -1

Характеристическое уравнение этой матрицы

3 2

имеет вид р + 3р + 7р + 5 = 0.

Заметим, что масштабированием времени диагональные элементы матрицы А всегда можно сделать единичными: аи = —1. В этом случае из равенства (3) следует, что элементы вектора Ь регулярной системы определяются формулами

Ь = 72а-, * = 1, п.

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

£ = s2, ..., sn].

Перейдем к рассмотрению регулярных систем третьего порядка. Будем считать, что ГСЧ системы упорядочены по убыванию: а1 > а2 > а3 > 0.

Передаточная функция регулярной системы

третьего порядка имеет вид [1]: Щ р) =

= 2к

р2 2(У 1 (у 1 - У2) , У2(У2 - У3) У3(У 1 - Уз)л|р , а р + т 1-+---)Р + а0

к1 V1 + V 2 V 2 + V 3 V1 + V 3 '

23

3

32 р + 3 р + а 1 р + а0

(4)

где а1 = 2 + 2 + 2,

1 1 + V 1 + V 3' 2 + V 3'

= Г - У2)(У1 - У3)(У2 - У 3 ^2 , к = V . + У2 + У3 .

0 ^ 1 + V 2 XV1 + V 3 XV 2 + V 3 У 1 2 3

В соответствии с формулой (4) одним из условий того, что произвольная ПФ третьего порядка

Щ р)

= в 2- + в1 р + в 0

3 2

р + а 2р + а 1р + а 0

(5)

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

А =

3 в2 в0 а2 а0

0 ^ 3а0в2 = а2в0.

(6)

Это равенство представляет собой первое необходимое условие регулярности системы третьего порядка.

Условие (6) учитывает только соотношение коэффициентов а0, а2, р0 и р2. При более подробном рассмотрении связей между коэффициентами формулы (4) получаем второе условие регулярности:

Б2 — 8 а 2 р1р2(а0 — а1а2) + + а2 р 1 (3ао — а^ + 2 а 2) + + 4 в2 (9 а0 - 9а

оа^2 + 2а0а2 + 2а 1 а2) — 0. (7)

Одновременное выполнение равенств (6) и (7) является необходимым и достаточным условием регулярности системы.

Для получения графического условия регулярности систем третьего порядка удобно использовать их представление на плоскости Вышнеград-ского.

2. ПРЕДСТАВЛЕНИЕ СИСТЕМ НА ПЛОСКОСТИ ВЫШНЕГРАДСКОГО

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

требованной. Наиболее характерным примером служит известный метод ^-разбиения [4].

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

Пусть система третьего порядка задана передаточной функцией (5).

Делая замену переменной р — ^00 д, получаем нормированное характеристическое уравнение

д3 + Лд2 + Бд + 1 — 0,

А — ^ (з/00) , Б — ^ З/00 , а 0 а 0

(8)

Рис. 2. Классическая диаграмма Вышнеградского

где коэффициенты Л и Б называются параметрами Вышнеградского.

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

На классической диаграмме Вышнеградского [4] выделены области с разными типами корней характеристического полинома (рис. 2). Границей области устойчивости на диаграмме служит гипербола Б — 1/А. Область, лежащая выше этой кривой, разбивается на три части I, II и III, соответствующие разным вариантам расположения корней характеристического уравнения. Например, часть II, имеющая вид остроконечного клюва, отвечает случаю вещественных корней.

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

3. ГРАНИЦА РЕГУЛЯРНЫХ СИСТЕМ НА ПЛОСКОСТИ ВЫШНЕГРАДСКОГО

Для определения области, соответствующей регулярным системам на плоскости Вышнеградского, была проведена серия компьютерных экспери-

Рис. 3. Расположение регулярных систем на плоскости Вышне-градского

ментов. Для каждого из 1000 случайных наборов ГСЧ, взятых из интервала от 1 до 10, и каждого из четырех вариантов сигнатурных векторов S, строилась матрица А сбалансированного представления регулярной системы. Элементы матрицы определялись по формуле (2). Ее характеристический полином приводился к форме Вышнеградского (8) и соответствующая точка наносилась на плоскость. Компьютерные эксперименты показали, что точки, соответствующие регулярным системам, лежат в пределах строго определенной области на плоскости Вышнеградского (рис. 3).

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

Преобразуем знаменатель ПФ (4) к форме Вышнеградского (8). Тогда получим следующие значения параметров А и Б:

А —

3

V 2Л 2/3 Г V1 - V 3Л 2/3 Г V 2 - V 3Л 2/3 '

V л ) ^ 1 + V ч) ^ л + V ч )

V1 + V 2

Б —

( V1 - V 2 ) 2 + (^Ч^^З^ + ( V1 + V 2 ) 2 ( V1 + V 3 ) 2 ( 2 + V 3 ) 2

П^Л 4/3 4/3 ( V 2 - V 3Л 4/3 '

V1 + V 2 ' ^ 1 + V 3' ^ 2 + V 3'

(9)

Для сокращения числа переменных перейдем от собственных чисел кросс-грамиана к их отно-

шениям д1 — v1/v2, д2 — v2/v3 и получим систему из двух уравнений вида:

Л — /1(д1, д2), Б — ^(др д2).

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

6 Л1 5 . 2 4 Л 3 3 к 2 2

д 1 - 6 -- д1 + 3 -2 д1 - 4 -3 д1 + 3 -2 д1 -

а2

а2

к1

- 6 -1 дх + 1 — 0,

где к1 — 2Л3 + 3ЛБ - 9, к2 — 20Л6 + 12Л4Б -

- 27Л2Б2 - 612Л3 + 162ЛБ - 234,

к3 — 20Л6 + 36Л4Б + 34Л2Б2 + 756Л3 - 486ЛБ + 36,

а — 2Л3 - 9ЛБ + 27.

Это уравнение является возвратным, поэтому заменой переменных х — д1 + 1/д1 его порядок понижается вдвое, в результате чего получаем кубическое уравнение:

- 61 х2 + 12 - х - 8 - — 0

к4

к5

(10)

где к4 — 4Л6 + 12Л4Б - 27Л2Б2 - 180Л3 + 162ЛБ - 234, к5 — 4Л6 + 36Л4Б + 81Л2Б2 + 324Л3 - 486ЛБ + 729.

Выполняя обратную замену переменных, получаем, что каждому решению уравнения (10) соответствует пара взаимно обратных корней уравнения д 1 - д1х + 1 — 0. Выбирая один из них

(больший по модулю, так как ГСЧ упорядочены по убыванию), получаем три решения для отношений собственных чисел кросс-грамиана д1, д2 и д3 — v3/v1 — 1/(д1д2).

Нас интересуют только вещественные д1 и д2,

а, следовательно, и вещественные корни уравнения (10). Очевидно, что критическим случаем их существования будет равенство нулю дискриминанта [6] этого уравнения:

А — 4Л6 - 12Л4Б - 2Л3Б3 + 20Л3 + 15Л2Б2 -

- 36ЛБ + 27 — 0.

(11)

Уравнение (11) описывает границу области существования регулярных систем на плоскости Вы-шнеградского (кривая 1 на рис. 3). Из него вытекает третье необходимое условие регулярности,

3

х

согласно которому дискриминант Б3 должен быть отрицательным.

Если при подстановке параметров Вышнеград-ского А и В для данной системы в формулу (11) получено положительное число Б3 > 0, система не является регулярной. Этим условием удобно пользоваться графически, как и показано далее в § 6, поэтому в дальнейшем будем называть его графическим условием регулярности.

Случаю Б3 = 0 соответствует наличие у уравнения (10) пары кратных корней, т. е. равенство отношений q1 = q2 = к. Это означает, что точкам кривой 1 (см. рис. 3) отвечают регулярные системы, у которых ганкелевы сингулярные числа расположены в геометрической прогрессии: а2 = hа1,

а3 = к2аг

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

Передаточная функция таких систем получается из формулы (4) при q1 = q2 = h и имеет вид:

W(p) = 2vx

( h2 + h + 1 ) (h2 + 1) p + (h - 1 )2

х ( к + 1 ) 2р2 + 2 ( /г2 + 1 ) ( к2 - 1 ) 2р + ( к + 1 ) 6 (к2 + 1)(к + 1 )2р2 + 2(к2 + к + 1)(к + 1 )2р + (к - 1 )4 '

Приводя знаменатель к форме Вышнеградско-го и исключая переменную к из выражений для его коэффициентов А и В, снова получаем уравнение (11).

Каждой точке в области регулярных систем на плоскости Вышнеградского соответствует в общем случае пара регулярных систем (с точностью до умножения их ПФ на число). Две указанные системы отличаются друг от друга перестановкой q1 и q2. На границе области q1 = q2, и системы становятся равными, т. е. каждой граничной точке отвечает только одна регулярная система, и ее ГСЧ образуют геометрическую прогрессию.

4. УЧЕТ ЗНАКОВ ГАНКЕЛЕВЫХ СОБСТВЕННЫХ ЗНАЧЕНИЙ

Системе третьего порядка с одним и тем же набором ГСЧ может соответствовать 2 = 8 различных сигнатурных векторов £ (сочетаний знаков). Легко показать, что при умножении ПФ на —1

собственные числа кросс-грамиана меняют свои знаки на противоположные. При этом характеристический полином системы не изменится. Это позволяет сократить число рассматриваемых сочетаний знаков до четырех, а именно (+ + +), (— + +), (+---+), (+ + —). Компьютерные эксперименты показали, что каждому из таких сочетаний будут отвечать свои области на плоскости Вышнеградского (см. рис. 3).

На рис. 3 выделены три внутренние области регулярных систем: (+ +) — системы с одинаковыми знаками V,.; (—) — системы, в которых v2 отличается по знаку от v1 и v3; (+—), (—+) — остальные системы. Эти обозначения соответствуют знакам

отношении

4l = v1/v2' 42 = v2/v3'

Найдем уравнение штриховой кривой, разделяющей три указанные области.

Очевидно, что пограничным случаем описанного расположения корней характеристического полинома кросс-грамиана является равенство одного из них нулю. При этом порядок ПФ понижается. Например, при v1 = 0 передаточная функция (4) принимает вид:

= 2(v2 + v3)

W(p) =

(v2 + v3)2P + (v2 - v3)2

(v2 + v3)2P2 + 2(v2 + v3)2P + (v2 - v3)2

При этом в дроби (4) числитель и знаменатель сокращаются на v2 ^ (р + 1), формулы (9) упрощаются и уравнение (11) принимает вид:

2A3 - 9 AB + 27 = О.

(12)

Ему отвечает кривая 2 на рис. 3. Она касается кривоИ 1 в точке с координатами (3; 3) и делит область существования регулярных систем на 3 части. В первоИ части (++), справа от кривоИ, лежат точки регулярных систем, у которых q1 и q2 положительны (т. е. все три v; имеют одинаковый знак). Во второй части (—), слева от кривоИ 2, находятся точки систем с отрицательными q1 и q2. Это системы, в которых v2 отличается по знаку от v1 и v3. Третья часть (+-), (-+) расположена над кривоИ, еИ соответствует точки всех остальных регулярных систем. Точка (3; 3) является особоИ, в неИ все три части соприкасаются.

Отметим, что кривая (12) уже имеется на клас-сическоИ плоскости Вышнеградского. Она составляет левую границу области III (см. рис. 2) и соответствует расположению корнеИ характеристического полинома в арифметическоИ прогрессии. Однако ранее смысл этоИ границы влево от точ-

ки (3; 3) оставался неясным, и обычно эта часть кривой просто не изображалась.

В случае кратных ганкелевых сингулярных чисел (у1 = у2) точки систем асимптотически приближаются к кривой 2, и система так же, как и в предыдущем случае, оказывается неминимальной, т. е. в числителе и знаменателе ПФ системы появляются одинаковые корни, и после их сокращения порядок системы понижается. Коэффициенты А и В для таких систем становятся бесконечно большими. Это означает, что ни одна из точек на кривой 2 не отвечает регулярной системе. В подтверждение этого факта заметим, что левая часть выражения (12) стоит в знаменателях членов уравнения (10). Следовательно, не существует решений qx и д2, для которых выполнялось бы условие (12).

5. ОБУСЛОВЛЕННОСТЬ РЕГУЛЯРНЫХ СИСТЕМ

Рис. 4. Кривые равной обусловленности

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

2 _ тах(X¡) тт(X'

ц

/ = 1, п ,

где Хг. — собственные числа матрицы АА .

Чем больше число ц, тем хуже обусловлена матрица, и тем больше уровень вычислительных погрешностей. Значение ц = 1 соответствует идеальной обусловленности, которая достигается в случае ортогональных матриц.

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

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

Точки, лежащие на одной кривой, отвечают регулярным системам с одинаковым числом обусловленности матрицы А. Чем ниже расположена кривая на плоскости, тем меньше число обусловленности. Докажем, что регулярных систем с ц = 1 не существует. Действительно, для выполнения условия ц = 1 требуется, чтобы матрица А была ортогональной. Решая матричное уравнение ААТ = Е,

где Е — единичная матрица, для матрицы А сбалансированного представления регулярных систем, получаем, что в системе должна быть пара кратных ГСЧ (а1 = а2), что противоречит определению регулярных систем.

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

Таким образом, положение точки с координатами (А; В) регулярной системы на плоскости Вы-шнеградского позволяет оценить обусловленность матрицы А сбалансированного представления и определить соотношения знаков собственных чисел кросс-грамиана.

6. ПРИМЕРЫ РЕГУЛЯРНЫХ СИСТЕМ

Пример 1. Пассивная электрическая схема. Рассмотрим электрическую ^С-цепь лестничного типа, схема которой приведена на рис. 5.

Входным сигналом цепи является ток /, выходным — напряжение и на конденсаторе Сг Известны значения пяти параметров цепи: С1 = 0,141 Ф, С2 = 0,634 Ф, С3 = 18,043 Ф, Ях = 7,723 Ом, Я2 = 4,699 Ом. Требуется

Рис. 5. Схема лестничной ^С-цепи

Рис. 6. Оценка параметров систем на плоскости Вышнеградского

определить, можно ли подбором значения сопротивления Я3 сделать эту систему регулярной.

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

W( p) =

C2 Cз R ! R 2 Язр2 + ( C2 R ! ( R 2 + R3 ) +

C1 C2 C3 R1R2 R3p +

+ C3 R 3 ( R1 + R2 ) ) p + R1 + R2 + R 3

+ £6р + р + 1

где = С1С2^1(^2 + Я3) + С3^3(С1(^1 + Я2) + С2Я2), = С1(^1 + Я2 + Я3) + С2(Я2 + Я3) + С3Я2. Подстановка заданных значений параметров дает:

Щ р) =

= 415,13 Я3р2 + (229 Я3 + 23) р + 12,42 + Я3 58,53 Я3р3 + (86 Я3 + 3,24) р2 + (18,82 Я3 + 4,73) р + 1'

Выписывая первое необходимое условие регулярности (6), получим уравнение для

86 Я3 - 173,29Я3 + 40,29 = 0.

Оно имеет два решения: = 1,732 и = 0,286. Для первого из них ПФ принимает вид:

= 7,078pz + 4,131 p + 0,1393

2

Щ р) = 3

р3 + 1,5 р" + 0,3673р + 0,0098

При подстановке полученных коэффициентов ПФ в условие (7) получаем, что второе необходимое условие регулярности = 0 выполняется. Следовательно, система является регулярной, так как параметры системы удовлетворяют условиям (6) и (7) одновременно.

Приводя характеристический полином к форме Вы-шнеградского (8), находим значения коэффициентов А

и В: А = 7, В = 8. Положение соответствующей точки на плоскости Вышнеградского показано на рис. 6. Из графика видно, что точка, соответствующая исследуемой электрической цепи, лежит в области регулярных систем. Это означает, что для системы выполняется графическое необходимое условие регулярности.

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

A =

-0,5 0,36 -0,24 0,36 -0,5 0,45 -0,24 0,45 -0,5_

Диагональные элементы матрицы A одинаковы, следовательно, система является регулярной. Проверка второго значения R3 дала отрицательный результат.

Оценим знаки собственных чисел кросс-грамиана системы. Точка с координатами (A; B) попала в область систем с одинаковыми знаками vr Действительно, прямые вычисления показывают, что значения v для данной цепи положительны: vx = 5,704, v2 = 1, v3 = 0,3737. Такое сочетание знаков означает, что ПФ системы можно представить в виде:

W( p) =

_ r 1

+

r2

+

r3

р - р 1 р - р 2 р - Р 3

где г1 — некоторые положительные числа. Следовательно, весовая функция будет иметь вид суммы трех затухающих экспонент с положительными коэффициентами перед ними, и переходные процессы в системе будут иметь монотонный характер. Этот вывод справедлив для любой регулярной пассивной ^С- или ^Х-цепи.

Число обусловленности ц матрицы А в соответствии с построенными ранее кривыми равной обусловленности (см. рис. 4) составляет приблизительно 39. Это означает, что матрица сбалансированного представления системы достаточно хорошо обусловлена. Точное значение ц = 39,363, рассчитанное в пакете МаИаЪ, отличается примерно на 1 %.

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

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

Рис. 7. Модель двухмассовой системы

может быть описана следующими дифференциальными уравнениями:

ш1 х х + — &2) + (ИI = И,

ш2 х2 — — &2) + = 0,

где ( — коэффициент демпфирования.

В качестве выходного сигнала выберем скорость у1 первого тела. В этом случае порядок системы понижается до трех, и ПФ имеет вид:

2 d Е p + —p + _Ъ.

W( p) =

m2

m2

p3 + d( m2 + m l )p2 + (d +mxE + km2) + 2+Е

.(13)

m1 m

l m2

m1m

lm2

m1 m

lm2

Выясним, при каких сочетаниях параметров ПФ (13) система будет регулярной. Для этого обратимся к первому необходимому условию регулярности. Выписывая условие (6) для ПФ (13), получаем:

3«0Р2 - «2Р(

= +Е

"■0^2 "-2^0 -2 (5m2 mi)

ml m2

Отсюда следует, что массы должны отличаться в пять раз: ш1 = 5ш2. Пусть ш1 = 0,03 кг, тогда ш2 = 0,006 кг. При анализе второго условия (7) было установлено, что система будет регулярной при значениях параметров = 0,1 Н/м, ( = 0,015.

Подставив эти численные значения в ПФ (13), получим:

W(p)

= +2 + 2,498p + 15,63 p3 + 3p2 + 20,01p + 15,63

Для проверки графического условия приведем характеристический полином к форме Вышнеградского (8) и вычислим значения коэффициентов А и В: А = 1,2, В = 3,2. Нанесем полученную точку на плоскость Вышнеградского (см. рис. 6). По расположению точки можно сказать, что для системы выполняется графическое необходимое условие регулярности. Построим матрицу А сбалансированного представления объекта:

A =

-1 0,56 0,51 0,56 -1 -4,14 -0,51 4,14 -1

Диагональные элементы матрицы А одинаковы, следовательно, данная система является регулярной.

Оценим знаки собственных чисел кросс-грамиана системы. Точка с координатами (А; В) на рис. 6 попала в область систем, у которых большее или меньшее собственное число кросс-грамиана отличается по знаку от двух остальных. Действительно, находя собственные значения кросс-грамиана, получаем: = 0,483, у2 = 0,045, у3 = -0,028.

Число обусловленности матрицы А в соответствии с построенными кривыми равной обусловленности (см. рис. 4) составило примерно пять. Точное значение числа обусловленности, рассчитанное в МаИаЪ, несколько выше: р = 5,315. Это означает, что матрица сбалансированного представления хорошо обусловлена.

ЗАКЛЮЧЕНИЕ

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

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

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

Исследована обусловленность системной матрицы сбалансированного представления регулярных объектов третьего порядка и построены линии равной обусловленности на плоскости Вышнеградского.

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

ЛИТЕРАТУРА

1. Мироновский Л.А., Курмаев И.Р. Синтез трисингулярных динамических систем // Информационно-управляющие системы. - 2010. - № 6. - С. 77-85.

2. Ober R.J. Balanced parameterization of classes of linear systems // SIAM J. Control and Optimization. - 1991. - Vol. 29, N 6. -P. 1251-1287.

3. Мироновский Л.А. Линейные системы с кратными сингулярными числами // Автоматика и телемеханика. - 2009. -№ 1. - С. 51-73.

4. Бесекерский В.А, Попов Е.П. Теория систем автоматического управления. - СПб.: Профессия, 2003. - 752 с.

5. Щепетов А.Г. Об оптимальных формах переходного процесса и амплитудно-частотной характеристики линейной динамической системы // Проблемы управления. - 2008. -№ 3. - С. 30-36.

6. Фаддеев Д.К. Лекции по алгебре. - М.: Физматлит, 1984. -416 с.

Статья представлена к публикации членом редколлегии

В.Ю. Рутковским.

Мироновский Леонид Алексеевич - д-р техн. наук, профессор,

И [email protected],

Соловьева Татьяна Николаевна — аспирант,

И [email protected],

Санкт-Петербургский государственный университет

аэрокосмического приборостроения, S (812) 494-70-44.