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

Тип работы:
Реферат
Предмет:
Физико-математические науки


Узнать стоимость

Детальная информация о работе

Выдержка из работы

Математическое моделирование. Оптимальное управление Вестник Нижегородского университета им. Н. И. Лобачевского, 2014, № 1 (1), с. 257−261
УДК 519. 711. 2
О ВЛИЯНИИ ВЫБОРА ПАРАМЕТРОВ КВАЗИИНВАРИАНТНОГО УПРАВЛЕНИЯ НА ТОЧНОСТЬ УСТАНОВИВШЕГОСЯ РЕШЕНИЯ И ВРЕМЯ ПЕРЕХОДНОГО ПРОЦЕССА
© 2014 г.
И.С. Гельфер
НИИ прикладной математики и кибернетики Нижегородского госуниверситета им. Н.И. Лобачевского
neymark@pmk. unn. ru
Пнступила в редакцию 06. 11. 2013
Работа посвящена дальнейшему исследованию свойств квазиинвариантного управления, а именно, даны оценки для точности установившегося решения и времени переходного процесса. Полученные аналитические результаты сопровождаются численными экспериментами.
Ключевок слнва: динамическая система, устойчивость, квазиинвариантное управление, численное моделирование.
С помощью синтеза квазиинвариантного управления объектом динамической системы, осуществлённого Ю. И. Неймарком, была разрешена проблема идеального регулятора Г. В. Щи-панова. При этом наметился ряд вопросов, связанных с возможностью и улучшением качества управления объектом, требующих дальнейшего изучения.
Система квазиинвариантного управления объектом описывается [1] уравнениями вида
Ап (р)х = -Вт (Р)(и + ^(Х)Х
ци = Dq (р)(х где Ап (р), Вт (р), Dq (р) — действительные полиномы от р = d|dt степеней п, т и q, ц -малый параметр, х — управляемая действительная переменная, и — управление, ?,(Х) — внешнее ограниченное воздействие. В [1] показано, что при любом Ап (р), гурвицевом Вт (р) и надлежащим образом выбираемых ц и Dq (р), q = п — т -1, возможно обеспечить устойчивость системы управления и спустя некоторое время Т & gt- 0 сколь угодную малость |х (Х)|. Но в
переходном процессе и может достигать очень больших значений, а время установления — быть длительным, в то время как |и| должен быть
ограниченным, а сам переходный процесс достаточно быстрым. Требование ограниченности |и| можно удовлетворить путём изменения вида управления, варианты которого, предложенные Ю. И. Неймарком, рассмотрены в [2]. Способ добиться малости времени установления Т, также предложенный Ю. И. Неймарком, кратко
описан в [3]. В данной работе подробно показывается, как выбор корней полинома Dq (р) и значений параметра ц влияет на время переходного процесса Т, а также даётся оценка точности установившегося решения |х (Х)|.
О быстроте переходного процесса
Общее решение уравнения (цАп (р) +
Вт (р^ (р))х = ЦВт (р^, которое получается из системы (1), если выразить и через х, есть сумма общего решения соответствующего однородного уравнения, х^Х), и частного решения неоднородного с нулевыми начальными условиями, х2(Х):
х (Х) =? с{* +ц? К | е^ (Х-^(тМт, (2)
i=1 i= 1 о
х1(х)=X/
где
се
определяет переходный
процесс, а х2(Х) = цХК1е*4(Х х) Е,(т^т — устано-
i=1 о
вившийся режим. Коэффициенты с. определяются из заданных начальных условий х (0) = х0, х (0) = х1, …, хп-1(0) = хп-1, а
К =
Вт (*,)
(Ц • Ап (р) + Вт (р^с (р)) р=
где Х1 — корни характеристического полинома цАп (р) + Вт (p)Dq (р), имеющие отрицательные действительные части в силу требуемой его устойчивости. Имеет место
Утверждение. Для любых сколь угодно малых е & gt- 0 и Т & gt- 0 при ограниченности возмущения Е,(Х) и начальных условий при т = 0 возможен такой выбор ц и корней полинома Dq (р), с = п -1, для которых при всех Х & gt- Т |х (Х)| & lt-е.
Доказательство. По условию Ы & lt- N & lt-да,
с, =
К… х,-і х х,+1… К
Л п-1 Л п-1 Л п-1 Л п-1
Кі … К- і X і X •. 1… К
1 ,-1 п-1, + 1 п • 1
і = 1, п.
1 1 1
К2 Кп
х1-- 1л п-1 X 2 … л п-1 Кп
с =
п-1 п-,
і + Хп-3 і Кі
і'- ,і *& lt-'-, і
V Л
-… + (-1)п-1 ХоПК- П (Кі-К,)
(3)
і = 0, п -1, |Е,(?) & lt-М & lt- да. Выберем корни многочлена Dn-1(р) действительными, отрицательными и различными. При ц — 0 п -1 корней полинома цАп (р) + Dn-1 (р) стремятся к корням Dn-1 (р), а ещё один, Кп, должен стремиться к
, т.к. Д,-1(Кп) * 0, а Мп (Кп) + Dn-1(Кn) = 0. Для устойчивости цАп (р) + Dn-1(р) требуется, чтобы Кп — -да. Поделив равенство цАп (Xп) +
+А,-1(Хп) = 0 на, ПолУчаем МОДп = -(^п-1 + + цап-1) -… — (й0 + ма0) Кпп& quot-1, откуда следует, что
при ц — 0 Кп «- -^п-1, где ап, йп1 — коэффи-апЦ
циенты при старших степенях полиномов Ап (р) и Dn-1 (р), и, чтобы Кп — -да, надо выбрать
ПЫ = ^К-^ ап).
Учитывая, что х2(?) берётся с нулевыми начальными условиями, имеем:
К1с1 + К 2С2 + … + К пСп = Х1,
Л п-1, Л п-1 «. Л п-1 «
К1 с1 + К2 С2 + … + Кп сп = Хп-1
откуда определяются
1…1×01… 1
Если 1 = п, то Xn входит в знаменатель выражения (3) в степени п -1, а в числителе X п отсутствует, и, т.к. при ц -- 0 Хп -- -да, а X у,
] = 1, п -1, конечны, то |сп | - 0.
Пусть 1 ф п. Хп входит в числитель и знаменатель выражения (3) в первой степени. Разделив числитель и знаменатель на Хп, получим
с =
2 х і+… +(-1)п х0 Пх і +
я (Х--1-
П (К і -К і) —
___ ч -1
Ї (х--1- ОЛ
Х0 ={Х0,…, Хп-1}, К1 — {К1,… Кп-1} ,
где Я = Хп-1 — Хп-2 2 Кі + … + (- 1) п 1 Х0 П Кі и
і *і, п jФi, n
/ = К, х X (Кі - К,) не зависят от Кп, а зна-
і *і, п Я (
чит, -, — стремятся к нулю при ц -- 0 и
х» к
~ Хп-3 2Кі + … + (-1)п Х0 ПК
п п
ц-0
і *і'-, п
0± 1'- V
і *і, п
П (Кі-Кі)
і'-*і, п
откуда следует ограниченность с., 1 = 1, п -1.
Теперь заменим X. на, где & quot-л — некоторое положительное число. Тогда
ц-0
(
л
^/лп-2 — Хп-3 і I
V
+… + (-1)пХ0ПК- лп-2П (К- - Кі)
і *і, п у
Л-1
л
V і *і'-, п
Знаменатель с1 — определитель Вандермон-да, равный П (Х. — X j), не равный нулю в си-
1& lt- ^& lt-. <-п
лу того, что X. Ф Х^. при 1 Ф у. После соответствующих преобразований, а именно, приведя числитель у с. к треугольному виду и сократив общие множители, получаем
При неограниченном возрастании ^ в числителе все члены, кроме последнего, стремятся к
х0 Пу
у Ф1, п
нулю и с,
ц-0,л-да
П (Кі -к)
і ^і, п
значит, с.
ограничен1 при сколь угодно больших, откуда следует, что |с.| & lt- С & lt-да для всех 1 = 1, п. Далее, задавая достаточно малую величину
+
К
К
п
п
с1 + с2 + … + сп = Х0
с
& gt-
+
с
е& gt- 0, можно при ц — 0 сделать |х1(Х)| & lt-
& lt- СX е1'-'-1 & lt- (п — 1) Се Хш"Т & lt--е, выбрав (за счёт
выбора п) максимальное из Xi, 1 = 1, п -1,
(

X шах & lt- -1П
Т у 2(п — 1) С Х Рассмотрим второй член в (2): К
г (Х)| & lt- цМX
X,.
К
где К1 =
-1
цап П ^ -X у)
уФ1
, 1 = 1, п. Если 1 = п.
то
К
ца
Л пп П (1 -X у/ X п) п-1
у Фп
и при ц — 0 стремится к величине
ц п-Чи-1
dn
п-1
К
1
и при ц — 0 стремится к
цап
ца
п Xtdn-1П (X1 -X у)
у Ф1, п
& lt- К & lt- да.
чаем |х2 (Х)| & lt- цМ
М = шах|(^(Х))|. Член XKiXi1еXit неограни-
1=1
ченно убывает с ростом Х в силу выбора Xi,
значит, |х2 (Х)| & lt- цМ
X КД-'
. Малость |х2(Х)|
при заданном М определяется величиной пара-
п
метра ц и значением выражения X KiX-1, где
1=1
к. = -1 цап П (X ,¦ - X у), 1 = 1, п, как было по/ уф1
казано выше, ограничены.
Следует заметить, что оценка для установившегося решения х2(Х) справедлива и при т Ф 0. В этом случае при 1 = п
Ът X пт П|1 —
X
Л
ЦanXпп П| 1
у=1
X.
Последнее выражение в силу того, что ц — 0 и Xn — -да, стремится к
Ъ Ъ цп-т-1аи-т-1
т т^ п
цОЛ п
d:
п-т-1
стремящейся к нулю вместе с ц. Если же 1 Ф п, то
Отсюда следует, что |К,. /X.| при ц — 0 стремится к нулю при т & lt- п -1 и ограничено при т = п -1 (X* - корни полинома Вт (р)). При
1 Ф п оценка для |К. /X.| получается аналогично случаю с т = 0 и отличается только конечным
множителем Ът П^ * ^).
Значит, х2 (Х) & lt- цМ (п — 1) К + о (ц), и при достаточно малом ц может быть меньше заданной малой величины е/2, а |х (Х)| & lt-е.
Оценка установившегося решения
Можно найти более точную оценку для решения в установившемся режиме х2(Х). После соответствующих преобразований полу-
где
Численные эксперименты
Для полинома А4(р) = р4 — р3 + 2р2 + 3р -2, постоянно действующего ограниченного случайного возмущения Е,(Х) = 25 sin (3t), при начальных условиях х (0) = 0. 5, х (0) = 0. 2, х (0) = 0.2 и
х (0) = 0. 1, е = 103 при надлежащем выборе полинома D3(р) и параметра ц определяется время Т, начиная с которого |х (Х)| & lt- е. Оказалось, что
3
при ц = 0. 001 для полинома D3(p) = П (р -X.)
1=1
при X1 = -1, X2 = -2, X3 = -3 время установления Т = 7. 44, для X1 = -10, X2 = -20, X3 = -30 Т = 0. 73 и для X1 = -100, X2 = -200, X3 = -300
Т = 0. 07 (рис. 1, соответственно а), б), в)).
В табл. 1 приведены значения времени переходного процесса Т и точности установившегося решения при различных X и ц.
Как видно из таблицы, время переходного процесса Т убывает с убыванием X1, X2, X3, а параметр ц существенного влияния на Т не оказывает, в то время как точность установив-
х
X
X
п
п
п
1
X
п
п-т
X
Таблица 1
X ц T X уст.
3 — & lt-N — - 0. 001 7. 440 5. 174e-04
3 — & lt-N — - 0. 0005 7. 505 2. 586e-04
3 — & lt-N — - 0. 0001 7. 605 5. 169e-05
-10, -20, -30 0. 001 0. 730 3. 927e-06
-10, -20, -30 0. 0005 0. 730 1. 964e-06
-10, -20, -30 0. 0001 0. 730 3. 927e-07
-100, -200, -300 0. 001 0. 070 4. 164e-09
-100, -200, -300 0. 0005 0. 070 2. 082e-09
-100, -200, -300 0. 0001 0. 070 4. 164e-10
-2, -3
-10, -15
-20, -30
-100, -150
T
7. 505
6. 460
6. 335
6. 235
Таблица 2
Х уст. 1
1. 64e-04
1. 56e-05
4. 10e-06
1. 66e-07
|x2(t)| не превосходит величины цМ
шегося решения хуСт = max|x2(t)| улучшается с убыванием X. и ц. В установившемся режиме
которая достигается при ^) = М. При Х1 = -1, Х2 = -2, X3 =-3, ц = 0. 001 и Щ) = 25 |хуст.| = =0. 0042, в то время как при ^) = 25sin (3t) хуст |, как видно из табл. 1, не превосходит 5. 174е-04.
Рассмотрим случай т = 1, Вт (р) = р +1, Ц),) = 25 эт (3^), начальные условия те же самые. В табл. 2 приведены значения Т и |хуст | при ц = 0. 0001 и корнях полинома В2(р) =
= (р -Х1)(р -Х2).
При убывании X. | хуст | становится сколь
угодно малым, в то время как Т стремится к конечной величине. Последнее происходит оттого, что максимальный корень полинома
4 6
Рис. 1
Вт (р)(р) нельзя сделать меньше максимального из заданных корней Вт (р).
Таким образом, точность установившегося решения и время переходного процесса целиком и полностью зависят от выбора параметров управления X и ц.
Работа выполнена при финансовой поддержке гранта РФФИ, проект № 11−01−379.
Список литературы
1. Неймарк Ю. И. Синтез и функциональные возможности квазиинвариантного управления // АиТ. 2008. № 10. С. 48−56.
2. Гельфер И. С., Неймарк Ю. И. Нелинейность как расширение возможностей квазиинвариантного управления // Нелинейные колебания механических систем. Труды VIII Всероссийской научной конференции. Т. 1. Н. Новгород, 2008. С. 132−136.
3. Неймарк Ю. И., Гельфер И. С. Квазиинвариантное управление с быстрым установлением // Модели, методы и программные средства: Тез. докл. итоговой научной конференции учебно-научного инновационного комплекса, Н. Новгород, 27−30 ноября 2007 г.
Н. Новгород, 2007. С. 27−30.
THE INFLUENCE OF THE CHOICE OF QUASI-INVARIANT CONTROL PARAMETERS ON THE ACCURACY OF A STEADY-STATE SOLUTION AND THE TRANSIENT PROCESS TIME
I.S. Gelfer
A further study is presented of quasi-invariant control properties- in particular, the estimates are given for the steady-state solution accuracy and the time of a transition process. The analytical results obtained are supported by numerical experiments.
Keywords: dynamic system, stability, quasi-invariant control, numerical simulation.
References
1. Nejmark Ju.I. Sintez i funkcional'-nye vozmozhnosti kvaziinvariantnogo upravlenija // AiT. 2008. № 10. S. 4856.
2. Gel'-fer I.S., Nejmark Ju.I. Nelinejnost'- kak rasshirenie vozmozhnostej kvaziinvariantnogo upravlenija // Neli-nejnye kolebanija mehanicheskih sistem. Trudy VIII Vse-
rossijskoj nauchnoj konferencii. T. 1. N. Novgorod, 2008. S. 132−136.
3. Nejmark Ju.I., Gel'-fer I.S. Kvaziinvariantnoe upravlenie s bystrym ustanovleniem // Modeli, metody i programmnye sredstva: Tez. dokl. itogovoj nauchnoj kon-ferencii uchebno-nauchnogo innovacionnogo kompleksa, N. Novgorod, 27−30 nojabrja 2007 g. N. Novgorod, 2007. S. 27−30.

ПоказатьСвернуть
Заполнить форму текущей работой