
книги из ГПНТБ / Зубов В.И. Математические методы исследования систем автоматического регулирования
.pdfчего можно всегда добиться при помощи замены искомых функций
Hi = 4 + ai |
= ! . • ■ ■ . «)• |
Предположим, что переходные процессы рассматриваемой си стемы можно описать при помощи системы дифференциальных уравнений
dys |
|
p ^ i + |
- j |
p :{ |
nn) |
Umi |
' ' ' |
y mn |
|
|
dt |
i=l |
|
•'l |
n |
t |
|
||||
|
mi+m2+ '''+mn>1 |
|
|
|
|
|
|
(5.5) |
||
|
|
|
|
|
|
|
|
|
|
|
правые части |
которых — голоморфные функции |
от |
у х, |
. |
. ., |
уп. |
||||
Предположим далее, |
что коэффициенты psi и |
P s(m‘..... тп) |
ве |
щественны и вещественные части корней характеристического уравнения | Р — КЕ J отрицательны. Построим для этой системы уравнений функции, о которых идет речь в теореме 13, для того чтобы с их помощью построить область асимптотической устойчи вости нулевого решения (5.5) или же аппроксимировать эту об ласть изнутри некоторыми вложенными областями. Для этого рас смотрим уравнение в частных производных 1
Е |
г- -1 ^ р (=У и - - .<, У п ) |
L |
1 + |
J |
Е |
* 1 |
|
/=i |
°У1 |
|
i=i |
|
|
||
Здесь функция ц>(уи |
уп) — положительно-определенная |
||||||
квадратичная форма у ъ . . ., |
уп. |
В качестве функции ф (уъ . . . |
. . ., уп) можно также брать любую голоморфную положительную функцию, наименьшие члены в разложении которой по целым поло
жительным степеням величины у ъ . |
. ., |
уп образуют положитель |
|||||||
но-определенную форму степени 2т\ |
m |
1. |
|
|
|
||||
|
|
П |
|
|
|
|
|
|
|
Множитель 1 + ^ Y}, стоящий в правой части уравнения (5.6), |
|||||||||
может |
быть |
i=i |
любой |
другой |
голоморфной |
функцией |
|||
заменен |
|||||||||
G (уъ |
. . ., уп), такой, что все интегральные кривые системы урав |
||||||||
нений |
|
|
|
|
|
|
|
|
|
|
|
|
dys |
Ys |
|
|
|
|
|
|
|
|
dx |
~ G |
|
|
|
|
|
продолжимы |
на все значения т; — оо <С т << + |
оо. |
|
|
|||||
Если таким свойством обладают интегральные кривые системы |
|||||||||
(5.5), |
то в качестве функции G (уъ |
. . ., |
уп) можно выбрать еди |
||||||
ницу. |
|
|
|
|
|
|
|
|
|
Как следует из вспомогательной теоремы Ляпунова |
[26], |
при |
|||||||
таком |
выборе функции |
ф (уг, |
. . ., |
уп) |
функция |
V {ух, . . ., |
уп) |
может быть получена единственным образом в виде степенного
ряда, сходящегося при достаточно |
малых |«/s| и обращающегося |
|||
в нуль |
при у г = |
у 2 = ■■• = |
уп = |
0. |
1(5.6) |
получается |
из условия |
4 теоремы 14 при замене V на In (1 + V). |
60
Будем искать решение V (уу, |
...,«/„) |
в виде ряда |
|
||||||
^ {Уъ • ■ •> Уп) |
— |
^ 2 |
(i/x> • • |
■>Уп) Ь П |
(У и ■ • • , Уп) ~\г ■■ ■ |
|
|||
|
|
• |
• |
• “Б Ут(*/]'ь • • |
•. Уп) + |
• • •» |
(5.7) |
||
где Vm (уI, |
. . ., |
уп) — однородная |
форма т-й степени относи |
||||||
тельно величин |
у х, . |
. ., |
уп. |
|
|
|
|
||
Подставляя |
(5.7) |
|
в |
(5.6), |
получим |
для определения |
форм |
||
Ут{Уъ ■• |
Уп) |
систему |
уравнений: |
|
|
2 |
*12 |
(5.8) |
|
|
|
«=1 |
w - *k=\ pikUk = Rm (Уъ • • • ’ |
= з - 4- • • |
где Rm является известной формой m-й степени, если найдены уже
Из системы (5.8) последовательно найдем формы V2, V3, . . .
Отыскание этих форм осуществляется следующим образом. Берем произвольную квадратичную форму V2 с коэффициентами, под лежащими определению. Подставляем ее в систему (5.8) и прирав ниваем слева и справа коэффициенты при одинаковых произведе ниях yyyj (i,j = 1, . . . , « ) . После этого получаем линейную си стему алгебраических уравнений, определитель которой отличен от нуля, поэтому коэффициенты формы V2 определяются однознач но через коэффициенты формы ср и параметры системы. После этого находим функцию R a, а затем функцию V3 в виде формы третьей степени с коэффициентами, подлежащими определению.
Подставляя V3в (5.8) при т = 3 и приравнивая справа и слева члены при одинаковых степенях у ъ . . ., уп, получаем систему линейных уравнений с отличным от нуля определителем. Продол жая этот процесс дальше, можно найти все члены ряда, определяю щего функцию V.
Известно, что этот ряд будет сходящимся во всяком случае в достаточно малой окрестности начала координат и квадратичная форма V2 будет определенно-отрицательной. Отметим, что выбор функции ф (уу, . . ., уп) в уравнении (5.6) влияет на область сходи-
п
мости ряда (5.7); например, при ф = 2 2j У1 уравнение
7 = ф(^ ’ •••’ ^ )(1+1/)
61
имеет решение
п
V (Ух, • • ■, Уп) = е
П
( = 1
а при <р (ух, • • •, Уп)
/ [
V (Ух, ■■■Уп) = е |
— 1. |
Если искать в этом примере V (уг, ■■•, Уп) в виДе ряда, то областью сходимости его в первом случае будет все пространство, тогда как во втором —■лишь его ограниченная часть.
Рассмотрим уравнение первого порядка
гДе f (У) — целая функция, не имеющая вещественных корней, отличных от у = 0. Уравнение для функции будет иметь вид
/( г /) |: = ф Ы (1 + / 2)(1 + ^ ) .
Нетрудно показать, что при ф — f2 V также получается целой функцией; больше того, если f (у) не имеет линейного члена в своем разложении по целым степеням у, выбором функции ф можно до биться того, чтобы и в этом случае V получилось в виде ряда.
То же самое следует выяснить для случая системы уравнений. Дадим теперь метод построения области, целиком погруженный
в область А. Рассмотрим |
семейство поверхностей V2 (Уъ ■■■ |
■■-,Ьп) = —Р 1см. (5.7)], где |
(0, +оо), р = const. Если поверх |
ность s — граница области А, то найдется такое значение р = ц, |
|
что поверхность V2 (уi, ■■■, |
Уп) — —Р будет касаться в некото |
рой точке поверхности s. Действительно, так как семейство поверх ностей V2 = —р заполняет все пространство, то найдется такое значение р, при котором У2 = —Р будет пересекать поверхность s.
Обозначим через —р наибольшее значение V2 (у2, ■■■, уп) на куске поверхности s, заключенном внутри поверхности V2 = —р.
Тогда V2 (у1, ■■•, уп) = —р будет касаться s.
62
Как было |
показано |
ранее, |
поверхность |
s имеет уравнение |
||
1 + V (X) = |
0. Тогда в точке касания X получаем |
|||||
|
|
дУ |
Y, = |
0; |
|
|
|
|
t=1 д х . |
|
|
|
|
|
дУ |
дУ |
|
|
дУ |
|
|
д х х |
д х 2 |
|
|
д х п |
|
|
д У 2 |
дУ 2 ~ ' |
~ |
д У 2 |
|
|
|
дх. |
д х а |
|
|
дхп |
|
откуда имеем |
|
|
|
|
|
|
|
|
у |
^ 2 |
|
|
|
|
|
1 |
dxt |
|
|
|
|
|
; = i |
|
|
|
|
(случай, когда s лежит целиком в области А, |
исключен). |
|||||
Заметим, что V2 (X) является функцией Ляпунова для системы |
||||||
(5.5). Найдем |
в силу системы (5.5) |
|
|
|||
|
|
|
дУ 2 |
|
ф (X) |
|
|
|
!=1 |
dyi Yi ~ |
|||
|
|
|
|
|
|
|
|
|
р {" |
пп) |
|
УпП |
|
|
|
■+тп > \ |
|
|
|
|
|
|
|
|
|
|
|
Функция W положительно-определенна. Найдем множество |
||||||
точек X, удовлетворяющих уравнению |
W {X) = 0, и обозначим |
|||||
его через W 0, считая, что точка X = |
0 не принадлежит множеству |
|||||
UV Найдем наибольшее значение функции |
V2 (X) на множестве |
|||||
W 0и обозначим его через —р 0. При этом р 0 > |
0, так как W (X) — |
положительно-определенная функция.
Имеет место следующее предложение: поверхность V2 (X) = = —Ро целиком содержится в области А.
Предположим противное, т. е. что эта поверхность частично содержится в области А. Это семейство поверхностей представляет
собой семейство |
эллипсоидов, заполняющих |
все пространство |
|
при 0 < р < + |
оо. Тогда на куске поверхности s существует |
||
точка, в которой функция |
V2 (X) достигает наибольшего значения |
||
—р 0; как было показано, |
в этой точке W = |
0, что невозможно, |
ибо —р 0 есть наибольшее значение V2 (X) на множестве W„. Сле довательно, значение —р 0 не может быть наибольшим для функции V2 на множестве W 0, откуда и следует наше утверждение.
Отметим, что в наших рассуждениях функция <р является произ вольной положительно-определенной квадратичной формой: таким образом,
К (X) - - Р о |
(5.9) |
63
есть семейство кривых, зависящих от параметров
ащ = —Ф"xtxh.
Пользуясь этим, для установления характера области устойчи вости можно сформулировать ряд признаков.
Для того чтобы область устойчивости была ограничена, необ ходимо, чтобы семейство поверхностей (5.9) было ограниченным, а для неограниченности области устойчивости достаточно, чтобы
это семейство было неограниченным. |
можно |
||
Для практического |
применения указанного метода |
||
с успехом пользоваться семейством поверхностей |
|
||
где |
sn = |
—Р, |
(5.Ю) |
|
|
|
|
sn ~ |
Уг W |
+ ' • ' + Vn {X), |
|
считая при этом, что sn = —р, ограничивает область sn > —р, содержащую начало координат X = 0.
Аналогичным образом определим для этого семейства функцию Wn (А) = ds и множество ее нулей W йп, причем X = 0 не при
надлежит Won. Как и раньше, можно показать, что поверхность s„ = —р 0„ целиком содержится в области устойчивости и сфор мулированные выше выводы можно перенести на семейство (5.10). Этим самым предполагается метод построения семейства областей, целиком принадлежащих области А.
Идея построения некоторой области, целиком принадлежа щей А, основанная на использовании функции Ляпунова, была применена самим Ляпуновым, а затем, в различных формах, дру гими авторами. В настоящем параграфе эта идея используется в несколько иной форме.
Вообще говоря, при увеличении я семейство sn = —р 0„ все более точно представляет границу области А.
Следует заметить, что особую роль играют те системы автома тического регулирования, у которых имеется единственный уста новившийся режим, асимптотически устойчивый в целом, т. е. такой установившийся режим, область притяжения которого совпадает со всем фазовым пространством.
Само собой разумеется, что общие теоремы об устойчивости в целом указывают лишь направление для исследования таких систем. Поэтому весьма важно иметь для исследования конкрет ных классов систем автоматического регулирования развитый математический аппарат, учитывающий особенности этих систем.
Способ построения специальных функций Ляпунова, обнару живающих устойчивость в целом, впервые был указан в работе [25] и более подробно рассмотрен в работах Е. А. Барбашина, Н. Н. Красовского, Н. П. Еругина, В. А. Плиса и ряда других авторов.
Г л а в а И
ПОСТРОЕНИЕ И СТАБИЛИЗАЦИЯ ПРОГРАММНЫХ ДВИЖЕНИЙ
6. Построение программных движений в управляемых системах
Рассмотрим сначала построение программных движений в ли нейных управляемых системах в простейшем случае. Пусть дана система п дифференциальных уравнений в векторной форме
|
|
■£ = />(*)x + Q(/)u + FW. |
(6.1) |
||
Будем |
считать, |
что элементы |
матриц |
Р (t) = \п X п), |
|
Q (t) = {п X г) |
и компоненты вектора |
F (t) заданы при t 5» О, |
|||
вещественны и |
непрерывны. |
|
— начальное |
||
Пусть |
заданы два |
постоянных вектора х 0 и |
и конечное положение системы (6.1). Требуется найти такую век
тор-функцию |
u = u (t), |
чтобы |
решение системы |
(6.1), |
начинаю |
||
щееся при t = 0 в точке х = |
х 0, попадало при |
t = Т |
в точку |
||||
х = x lt |
при |
этом |
т |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| и* (т) и (т) da |
|
|
|
|
|
|
|
о |
|
|
|
|
должен |
быть |
ограничен. |
|
|
называть |
программными. |
|
Такие управления и = и (t) будем |
|||||||
Рассмотрим вопросы |
существования |
программных |
управле |
ний и их построения. Обозначим через Y (t) фундаментальную си стему решений уравнений
~ = P(t)x- |
Y = P(t)Y |
(6.2) |
с начальными условиями Y (0) = |
Е, где Е — единичная матрица. |
Сделаем в системе (6.1) замену искомой функции по формуле
х = Y (t) z, |
(6.3) |
5 В. И. Зубов |
65 |
где вектор г представляет собой новую искомую вектор-функцию. Эта функция удовлетворяет уравнению
§ |
= |
в (0 11 + 0 (0 , |
(6.4) |
где |
Б (0 = У ~ г ( t ) Q ( t ) |
|
|
|
|
||
G { t ) = Y ~ ' { t ) F ( t ) . |
|
||
Произведем в системе (6.4) |
замену |
|
|
г |
— I |
-f Jt G (т) dr. |
(6.5) |
|
|
о |
|
Тогда векторная функция \(t) будет удовлетворять |
уравне |
||
нию |
|
|
|
|
§ |
= Б (0 и. |
(6.6) |
В результате этих преобразований над исходной системой краевые (начальное х0 и конечное хх) условия переходят в
I (0) = 10 = х0;
т
I (Т) = t1 = Y ' 1(Т) х, - j Y ' 1(т) F (т) (dx). |
(67) |
|
|
о |
|
Для простоты ПОЛОЖИМ |
Г — 1. |
|
Интегрирование от 0 до Т уравнения (6.6) дает |
|
|
|
т |
|
I i - S o = J b ( t ) u ( t ) d t . |
( 6 . 8 ) |
|
|
о |
|
Таким образом, для |
нахождения программного |
управления |
и (t) мы получили систему (6.8) линейных интегральных урав нений.
Решение уравнения (6.8) |
будем искать в виде |
|
и (0 = |
В*с + v (0, |
(6.9) |
где с — постоянный вектор, подлежащий определению, a v (t) — некоторая функция, такая, что
т |
|
jfcs (О V (о dt = 0 (s - 1, . . ., п). |
(6.10) |
о |
|
Здесь bs (t) (s = 1, . . ., п) означают компоненты вектора В, и равенство (6.10) выражает собой условие ортогональности функ ции v (t) ко всем компонентам вектора В.
66
|
Подставляя (6.9) |
в (6.8), |
находим |
|
|
|
||
|
А (Т) == Jг В (/) |
!i |
1>о — ^4 (Т) с, |
|
( 6 . 11) |
|||
где |
В* (t) dt. |
|
|
|
|
|
||
|
|
о |
|
|
|
если det А (Т) |
Ф 0 |
либо |
|
Уравнение (6.11) имеет решение, |
|||||||
ранг матрицы А (Т) совпадает с |
рангом расширенной |
матрицы |
||||||
А = {А (Т); Ь - Ы - |
|
|
|
|
|
|||
|
|
|
Теорема |
15 |
|
|
|
|
= |
Для |
того, чтобы существовало программное управление и = |
||||||
и (t), |
переводящее систему (6.1) из любого положения х 0 |
в лю |
||||||
бое другое положение х г за |
время Т, |
необходимо и достаточно, |
||||||
чтобы матрица |
г |
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
А - j ВВ* dt (В = |
F-!Q) |
|
|
|||
|
|
|
о |
|
|
|
|
|
была неособой. При этом все множество программных управлений дается формулой
и — В*с + v,
где v (/) — суммируемая с квадратом на [О, Т ] функция и
|
|
j В vdt — 0; |
|
||
|
|
о |
|
т |
|
|
|
|
|
|
|
с - А - 1 (Т) |
У - Д ^ х , — { У - 1 (т) F (т) dx — х 0 ]. |
||||
|
|
Определение |
14 |
||
Функции Ьг (t), |
. . ., bn (t) |
будем |
называть линейно незави |
||
симыми при t ^ |
0, если существует такое число Т > 0, что из ра |
||||
венства |
|
|
|
|
|
|
|
t c |
t b{ ( 0 - 0 |
(6.12) |
|
|
|
г=1 |
|
|
|
при t 6 [0, Т J |
следует с(- == |
0 |
(г == 1, |
. . ., п). |
Теорема 16
Для того чтобы функции Ьг (t), . . ., bn (t) были линейно-неза висимыми при t 0, необходимо и достаточно, чтобы существо вало такое число Т >> 0, чтобы
т
А (Т) = | ВВ* dt
о
была положительно-определенной матрицей,
5* |
67 |
Здесь В (t) — вектор с компонентами |
|
||||
|
|
bi (t), • • |
bn (t). |
|
|
Необходимость. |
Пусть функции b1 (t), . . ., |
bn (t) линейно-не |
|||
зависимы при t ^ |
0. |
Покажем, что матрица А (Т) положительно |
|||
определенна. |
из линейной независимости указанных функ |
||||
Действительно, |
|||||
ций следует, что существует такое Г |
> 0, что при любом выборе |
||||
постоянного вектора |
с, компоненты |
которого |
одновременно не |
равны нулю, будет существовать по крайней мере одна точка в про
межутке [0, Т], |
в которой В* с ф 0. |
Отсюда вытекает, что |
[В* с ]2 > 0 при |
£ £ [0, Т ] и вследствие |
непрерывности компо |
нент вектора В (£) будет |
|
|
|
т |
|
|
J [В*с]2 dt > 0 . |
(6.13) |
|
о |
|
Неравенство (6.13) можно записать в виде
г
с* j ВВ* dt-с = с*А (Т) с > 0
о
при любом с; это неравенство и означает положительную опреде ленность матрицы А (Т).
Заметим, что матрица А (Т) будет также положительно-опреде ленной при всех t ^ Т.
Достаточность. Пусть матрица А (Т) положительно-опреде ленная; тогда для любого вектора с = (сх, . . ., ck\ и такого, что II с II ф О, будет
гт
с* А (Т) с = Jc*BB*c d t — J [с*В]2 dt >>0.
оо
|
Следовательно, существует по крайней мере |
одна точка |
16 |
[0, Т], где будет с*В (t) =h 0. Это и показывает, |
что функции |
Ьъ |
. . ., Ьп линейно независимы при t ^ 0. |
|
|
Пусть в более общем случае граничные условия заданы в виде |
|
|
г |
|
|
\dG (Щ Y (®) = Н. |
(6.14) |
|
о |
|
Элементы матрицы [G = m X п) будем считать вещественными функциями ограниченной вариации. Требуется найти условия, при выполнении которых существует непрерывное решение х — = х (t) системы (6.11), удовлетворяющее (6.14), а также найти представление любого управления и = и (t), при котором такое решение существует.
68
Положим
т
А (Т) = j dG (О) У (fl). |
(6.15) |
о |
|
Здесь У (t) означает фундаментальную систему решений для урав нений (6.2).
Обозначим через Z постоянную матрицу, столбцы которой представляют собой набор всех линейно-независимых решений уравнения
|
А* (0) Z = |
0. |
(6.16) |
Если |
Теорема |
17 |
|
матрица |
|
|
|
|
г |
|
|
В = |
Z * j A (t) У -1 (0 Q (t) Q* (0 ( Y - 1 (0)* А * |
(t) dt Z (6.17) |
|
|
о |
|
|
неособая, то при любом выборе векторной функции Fh вектора Н существует решение системы (6.1), удовлетворяющее условию (6.14). При этом каждое программное управление и = и \t) пред ставляется в форме
и (0 = Q* (У -1) A*Zc + v, |
(6.18) |
где с есть вектор, определяемый единственным образом. Вектор ная функция v суммируема с квадратом и удовлетворяет условию ортогональности
т
(t) У -1 (0 Q {t) v (t) dt = 0, |
(6.19) |
о
Теорема 17*
Если матрица (6.17) неособая и матрица А (0) имеет ранг п, то каждому управлению (6.18) отвечает единственное решение х =
=х (t), удовлетворяющее условию (6.14). Если же ранг матрицы
А(0) есть к < п, то каждому управлению (6.18) отвечает семейство решений системы (6.1), удовлетворяющее условию (6.14) и завися щее от п — k произвольных постоянных.
Теорема 18
Если матрица (6.17) особая и имеет ранг /, то решение системы (6.1) удовлетворяющее (6.14), существует тогда и только тогда, когда ранг матрицы В ' также будет равен /. В этом случае каждое допустимое управление представляется в форме (6.18), однако вектор 1 зависит от п — k — / произвольных постоянных. Каждому такому управлению при k = п отвечает единственное решение х =
69