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

Михаил Георгиевич Кудинов

Сергей Борисович Силантьев

Андрей Васильевич Степовой

Сведения об авторах Военно-космическая академия им. А. Ф. Можайского, кафедра автономных систем управления, Санкт-Петербург; мл. науч. сотрудник; E-mail: [email protected]

канд. техн. наук, доцент; Военно-космическая академия им. А. Ф. Можайского, кафедра автономных систем управления, Санкт-Петербург; E-mail: [email protected]

канд. техн. наук; КБ точного машиностроения им. А. Э. Нудельмана, Москва; E-mail: [email protected]

Рекомендована кафедрой автономных систем управления ВКА

Поступила в редакцию 02.03.09 г.

УДК 62.50

Е. Д. Лихолетов, А. В. Ушаков, А. Ю. Цвентарный

АНАЛИЗ ПЕРЕКРЕСТНЫХ СВЯЗЕЙ В ДИНАМИЧЕСКИХ СИСТЕМАХ КЛАССА „ДВУМЕРНЫЙ ВХОД — ВЫХОД"

С ОДНОТИПНЫМИ КАНАЛАМИ

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

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

Постановка задачи. Рассматривается проблема построения автоматических систем, встраиваемых в измерительные мониторинговые комплексы, использующие принцип следящего преобразования. Такие системы, как правило [1—5], являются двухканальными с однотипным исполнением каналов и характеризуются наличием перекрестных межканальных связей с матрицей типа „матрица вращения" (МВ). Наличие перекрестных связей при определенных значениях аргумента д МВ становится причиной потери работоспособности двухка-нальных систем при устойчивых сепаратных каналах.

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

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

Рассмотрим динамическую систему (ДС) класса „двумерный вход — выход", структурная схема которой приведена на рис. 1. Здесь gt, s¿, yi, i = 1,2, — внешнее воздействие, ошибка слежения и выход i-го сепаратного канала соответственно; д — аргумент матрицы межканальных связей (матрицы вращения) T = col {[cos д sin д], [-sin д cos д]}; v1, v2 — переменные, образующие двумерный выход матрицы T; W(s) — передаточная функция прямой цепи i-го сепаратного канала спроектированной системы, имеющая представление

W ( s) = M ( s) N _1( s), (1)

где М(5), N(5) — полиномы с вещественными коэффициентами степенью т и п (т <п ) соответственно, образующие при д = 0 характеристический полином П( 5) каждого сепаратного канала в замкнутом виде в силу соотношения П( 5) = N (5) +М (5).

Рис. 1

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

Утверждение 1. Двухканальная система (см. рис. 1) с помощью процедуры комплекси-рования векторных переменных может быть приведена к скалярному представлению с передаточной функцией Жэкв (5) прямой цепи вида

^экв (5) = е-(5). (2)

Доказательство. Для рассматриваемой схемы справедливы следующие соотношения:

У1( s)

_ У2( s )_

W(s) 0 0 W(s)

cos | sin | - sin | cos |

S1( s) 82(s )

y (s) = W (s)(cos |i)S1 (s)+W (s )(sin |)8 2 (s) , У2 (s) = -W(s)(sin |)S1 (s)+W(s)(cos |)S2 (s) .

Введем в рассмотрение комплексированные переменные двухканальной системы:

* * *

g = g1 + jg2, 8 = 81 +jг2, У = У + jy2-Умножив соотношение (5) на J = , в результате получим

j>2 (s) = -W(s)(J sin |)б1 (s) + W(s)(cos |)(j82 (s)) .

Просуммируем выражения (4) и (7) и воспользуемся представлениями (6):

У* (s) = y (s) + jy2 (s) = W(s)(cos |)81 (s) + W(s)(sin |)82 (s)--W (s)(j sin |)81 (s) + W (s)(cos |)( j8 2 (s)) = = W (s )(cos |)81 (s )-W (s )(j sin |)(j 82 (s ))-W (s )(j sin |)81 (s) + W (s)(cos |)(j8 2 (s )) = = W(s) {cos j sin 81 (s) + W(s) {cos j sin (j82 (s)) =

= W (s )e- j| (81 (s)+j82 (s)) = e- j| W (s)8* (s) = W^ (s)8* (s).

(3)

(4)

(5)

(6) (7)

(8) ■

Доказанное утверждение делает справедливыми положения утверждения 2.

Утверждение 2. Двухканальная система (см. рис. 1) с симметричными каналами, описываемая передаточной функцией (2) относительно комплексированных переменных в*(^), _у*(£), оказывается на границе устойчивости, если аргумент д матрицы Т удовлетворяет условию

| = Аф, (9)

и становится неустойчивой при

' |>Аф. (10)

Доказательство. В силу критерия устойчивости Найквиста если выполняется условие

(9), то передаточная функция (2) характеризуется нулевым запасом устойчивости, при этом ее частотный годограф проходит через критическую точку (-1; 0), а при выполнении условия

(10) он охватывает эту точку. ■

Примечание. Если аргумент д оказывается интервальной величиной [д] = Д, Д согласно работам [9, 10], то в неравенстве (10) следует положить д = max { д , д }.

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

W(s) = M(S) =1--. (11)

W D(s) -M (s) s 1 + -1 v 1-. V 7

i=1

В замкнутом виде уравнение (11) формирует передаточную функцию „вход — выход":

0(s)=MS>=_W(£L=M(£)=_Vn__(12)

W gi (s) 1+W(s) D(s) " v '

s + / Уis i=1

Оценка запаса устойчивости сепаратного канала (12) осуществляется согласно следующей процедуре.

1. Задать ПДМ с характеристическим полиномом D(s), имеющим известные (Баттер-ворта, Ньютона) распределения мод или их модификации.

2. На основе выбранного в п. 1 представления ПДМ вычислить передаточную функцию разомкнутой системы вида (11).

3. Оценить запас устойчивости по фазе с помощью критерия устойчивости Найквиста, примененного к передаточной функции (11), с использованием процедуры „BODE" в программе MatLab Simulink.

При модификации распределения мод Баттерворта (МРМБ) изменяемым параметром является величина /2; 0] их локализации в левой полуплоскости на единичной окруж-

ности в секторе раскрывом 2у.

На рис. 2 приведен график зависимости Лф(у) для МРМБ с первого порядка по пятый (п = 1, 5). Для всех порядков ПДМ обнаруживается общая тенденция увеличения запаса устойчивости с уменьшением параметра у раскрыва сектора.

Лф,

90 80 70 60 50 40 30 20

10

X—к—X—X—х—

-0- — п=1

— п=2

— п=3

— п=4

— п=5

0

10

20

30

40

Рис. 2

50

60

70

У,

В табл. 1 приведены значения Лф для ПДМ с биномиальным распределением мод Ньютона (БРМН) при п = 1,5, здесь в аналитическом представлении полинома В(5, ©0) ©0 — характеристическая частота БРМН. Сравнение величин Лф, представленных на рис. 2 и в табл. 1, показывает, что запас устойчивости по фазе с БРМН превышает запасы устойчивости с МРМБ для равных порядков ПДМ при у^ 0, а при у = 0 эти величины совпадают. Таким образом, основное направление дальнейших модификаций распределений мод ПДМ будем связывать с возможностями модифицируемости БРМН.

Таблица 1

п Аналитическое представление полинома В(5, ©0) Аф,

1 5 + ©0 90

2 (+©0 )2 76,34

3 (+©0 )3 71,25

4 ( + ©0 )4 68,58

5 ( + ©0 )5 66,93

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

п-1

В(Х)=В(Х, ©0, V) = П(+©0 (1+/V)), (13)

1=0

где X — корень характеристического полинома.

Очевидно, что при v = 0 модифицированное биномиальное распределение принимает вид канонического биномиального.

При построении модифицированного БРМН в виде (13) учитывается доминирование одной моды над остальными, в соответствии с которым обеспечивается выбор параметра V. Следует ожидать, что при v>10 система с характеристическим полиномом вида (13) будет максимально приближена к апериодическому звену 1-го порядка.

Анализ перекрестных связей в динамических системах класса „ двумерный вход — выход " 39 На рис. 3 приведен график зависимости Аф(у) для ПДМ с первого порядка по пятый (п = 1, 5). Для всех порядков ПДМ обнаруживается общая тенденция увеличения запаса устойчивости с увеличением параметра V.

Аф, ...'

90

85

80

75 70 65

2

Рис. 3

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

Утверждение 3. Если ПДМ сепаратного канала системы имеет характеристический полином вида (13), то добротность = вг- уст (где в, уст — установившаяся ошибка в /-м сепаратном канале) канала по скорости определяется в соотношением

-1

а м-!+

1

1

1

(14)

ш0 (1+^ ш0 (l+2v)ш0 (1+(п-1^)со0

Доказательство строится на использовании аналитического представления добротности по скорости сепаратного канала на основе выражения (11) в виде = vn (сш )/^п-1 (ш0 ) . Если с помощью соотношения (13) сформировать вид коэффициентов vn-l(шо) и vn (Ш0), то их отношение приводит к уравнению (14). ■

Построение банка моделей ПДМ с динамическими показателями, параметризованными характеристической частотой Ш0 и аргументом д МВ. Полученные результаты по выбору полиномиальной модели с максимальным запасом устойчивости по фазе следует дополнить таблицей показателей {ву, у = 1, р} качества процессов в переходном и установившемся режимах сепаратного канала, вмонтированного в структуру двухканальной системы с перекрестными связями, характеризующимися аргументом д. Таким образом, таблица показателей будет содержать их значения в у (V, Ш0, д) .

В качестве примера в табл. 2 приведены значения показателей в у (V, Ш0, д) для ПДМ третьего порядка (п = 3) с модифицированным БРМН при v = 5, v = 10, v = 20 и де[0, 60°] Уш0; в таблице: а — перерегулирование; ^ — момент достижения переходной характеристикой ее максимального значения; 11 — время первого, после tа, пересечения переходной характеристикой границы допустимой ошибки А=5 %; t2 — время последнего пересечения

0

1

3

4

V

переходной характеристикой границы допустимой ошибки Д=5 %. Длительность переходного процесса ^ определяется выражением ^ = тах t2}.

_Таблица 2

v = 5 v = 10 v = 20

ст, % 'ст '1,c c ю0 ст, % 'ст '1,c c Ю0 ст, % 'ст '1,c c Ю0

0 0 — — — 0,769 0 — — — 0,87 0 — — — 0,930

6 0 — — — 0,774 0 — — — 0,874 0 — — — 0,935

12 0 — — — 0,786 0 — — — 0,889 0 — — — 0,951

18 0 — — — 0,809 0 — — — 0,914 0 — — — 0,978

24 1 3,54 Ю0 — — 0,842 0 — — — 0,952 0 — — — 1,018

30 1,5 3,07 Ю0 — — 0,888 1,06 3,03 ю0 — — 1,004 0,7 3,06 Ю0 — — 1,074

36 3,8 2,78 Ю0 — — 0,951 2,7 2,65 ю0 — — 1,075 2 2,64 Ю0 — — 1,15

42 7,7 2,58 Ю0 3,35 Ю0 — 1,035 5,6 2,4 ю0 3,17 ю0 5,28 Ю0 1,17 4,2 2,36 Ю0 — — 1,252

48 13,6 2,45 Ю0 3,24 Ю0 — 1,149 10,28 2,23 ю0 3,03 ю0 5,98 Ю0 1,3 8 2,15 Ю0 2,92 Ю0 5,52 Ю0 1,39

54 22,4 2,36 Ю0 3,16 Ю0 — 1,307 17,4 2,11 ю0 2,93 ю0 6,25 Ю0 1,479 13,8 2,01 Ю0 2,79 Ю0 5,86 Ю0 1,583

60 37,6 2,31 Ю0 3,11 Ю0 10,28 Ю0 1,539 27,8 2,03 ю0 2,85 ю0 9,1 Ю0 1,739 22,4 1,9 Ю0 2,7 Ю0 5,96 Ю0 1,861

Использование данных табл. 2 позволяет предложить следующую оформленную в виде алгоритма процедуру поканального синтеза систем класса „двумерный вход — выход", основанную на использовании возможностей модального управления.

Алгоритм.

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

2. Произвести оценку возможного диапазона вариаций параметра | — аргумента матрицы T — с последующим его представлением в виде интервального числа [|д] = | .

3. В силу доказанных утверждений принять | = max { |, |} .

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

5. Построить векторно-матричное ( A, B, C )-представление исходной динамической модели сепаратного канала, рассматриваемого в процедуре синтеза как объект управления.

6. Осуществить выбор типа ПДМ (параметризованной коэффициентом V, характеристической частотой Ш0 и аргументом д), доставляющей системе в условиях наличия перекрестных связей динамические показатели, сформированные в п. 1.

7. Построить модальную модель на основе векторно-матричного ( Г, Н )-представления ПДМ, выбранной при выполнении п. 6.

8. Осуществить канонический синтез [4—9, 10] сепаратных каналов двухканальной системы с использованием возможностей модельного управления.

9. Провести комплексное экспериментальное исследование в программе МайаЬ 81шиНпк динамических показателей спроектированной двухканальной ДС.

Результаты компьютерного эксперимента. Для иллюстрации полученных результатов проведено исследование двух версий двухканальной динамической системы (см. рис. 1), каналы которых построены с использованием ПДМ третьего порядка с модифицируемым биномиальным распределением для значений v = 1 и v = 10 при Ш0 =10 с-1 и д = 0, 30 и 75°.

Графики процессов в пространстве выходов двумерной системы при входном векторном скачкообразном единичном воздействии приведены на рис. 4, а—е.

а)

v=1, д=0

б)

у=10, д=0

у2 1,5 1

в)

-0,5 0 0,5 1 1,5 у1 v=1, д=30°

1

0,5 0

-0,5

г)

д)

0 0,5 1 1,5 у1 у=1, д=75°

-0,5 0 0,5 1 1,5 у1 у=10, д=30°

1

0,5 0

-0,5

у2 2

1

0

е)

0 0,5 1 1,5 у1 v=10, д=75°

1

0,5 0

-0,5

0 0,5

1,5 у1

Рис. 4

ПДМ третьего порядка с модифицированным БРМН при v = 1 обладает запасом устойчивости Аф = 73,277°, а при v = 10 — Аф = 83,711°. Как и следовало ожидать, при д = 75° и v = 1 (см. рис. 4, д) система оказывается неработоспособной, в то время как при д = 75° и v = 10 (см. рис. 4, е) система остается работоспособной.

1

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

1. Акунов Т. А., Сударчиков С. А., Ушаков А. В. Синтез фотоэлектрической следящей системы на основе интервальных модельных представлений. Часть I. Построение интервальной модели компонентов системы // Изв. вузов. Приборостроение. 2004. Т. 47, № 1.

2. Акунов Т. А., Сударчиков С. А., Ушаков А. В. Синтез фотоэлектрической следящей системы на основе интервальных модельных представлений. Часть II. Синтез управления, обеспечивающего стабильные эллипсоидные показатели качества системы // Там же. 2004. Т. 47, № 2.

3. КвакернаакХ., Сиван Р. Линейные оптимальные системы управления: Пер с англ. М: Мир, 1977.

4. Синтез дискретных регуляторов при помощи ЭВМ / В. В. Григорьев, В. Н. Дроздов, В. В. Лаврентьев, А. В. Ушаков. Л.: Машиностроение, 1983.

5. Николаев П. В., Сабинин Ю. А. Фотоэлектрические следящие системы. Л.: Энергия, 1969.

6. МирошникИ. В. Теория автоматического управления: Линейные системы. СПб.: Питер, 2005.

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

8. Ушаков А. В. Обобщенное модальное управление // Изв. вузов. Приборостроение. 2000. Т. 43, № 3.

9. Никифоров В. О., Ушаков А. В. Управление в условиях неопределенности: чувствительность, адаптация, робастность. СПб.: СПбГИТМО (ТУ), 2002.

10. Дударенко Н. А., Слита О. В., Ушаков А. В. Математические основы современной теории управления: аппарат метода пространства состояний: Учеб. пособие / Под ред. А. В. Ушакова. СПб.: СПбГУ ИТМО, 2008.

Сведения об авторах

Евгений Дмитриевич Лихолетов — студент; Санкт-Петербургский государственный университет информационных технологий, механики и оптики, кафедра систем управления и информатики; E-mail: [email protected] Анатолий Владимирович Ушаков — д-р. техн. наук, профессор; Санкт-Петербургский государственный

университет информационных технологий, механики и оптики, кафедра систем управления и информатики; E-mail: [email protected] Артем Юрьевич Цвентарный — аспирант; Санкт-Петербургский государственный университет ин-

формационных технологий, механики и оптики, кафедра систем управления и информатики; E-mail: [email protected]

Рекомендована кафедрой Поступила в редакцию

систем управления и информатики 25.12.08 г.