Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Зубов В.И. Математические методы исследования систем автоматического регулирования

.pdf
Скачиваний:
32
Добавлен:
25.10.2023
Размер:
12.11 Mб
Скачать

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

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

Соседние файлы в папке книги из ГПНТБ