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

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

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

следовательно,

i

lj gradfj.

 

grad,/о = grad/0 — £

(32.37)

/=!

 

Числа Ku . . %i определяются

из условия

ортогональности

вектора (32.37) к векторам grad/, . Отсюда получаем систему I алгебраических уравнений для определения Ях, . . ., А,,:

(grad fj, grad/0) = £ ^-(gradf/f grad/,)

(/ =

1, . . 0 -

(32.38)

1=1

 

 

 

Без ограничения общности можно считать, что векторьучгаб /у,

(/ = 1, . . ., I) линейно-независимы в

точке

(хх, . . .,

хп), ибо

в противном случае из них можно было бы выбрать независимую систему. Матрицу системы (32.38) обозначим через С; ее опреде­

литель не обращается в

нуль. Пусть

R — матрица, столбцами

которой

являются

векторы

grad f x, . . .,

grad ft. Тогда вектор

А = (A.j,

. . . . Xf)

можно

 

представить

в

форме

 

 

А

=

C-'R* grad /о.

(32.39)

Знак * означает транспонирование матрицы. Из (32.39) вытекает,

что

величины

А^,

 

К;

являются

функциями

переменных

Х х, .

.

х п.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Теорема 77

 

 

 

 

 

 

 

 

Если выполняются условия:

 

 

 

 

 

 

 

 

 

1)

функция

f 0

при

условии (32.32) имеет оптимальное значе­

ние в некоторой

точке

(х(0), . . .,

х ^у,

 

 

 

 

 

 

 

2)

вектор grad, f 0

обращается

в

нуль

при

условии

(32.32)

 

 

 

 

 

 

. . ., лф°->); |grad, f0\ > а

 

 

 

П

х\~у

лишь в одной точке (х<°>,

> 0 п р и

^

+ оо и точка (х{0),

. .

 

х п(0>) удовлетворяет

(32.32),

то любая

интегральная

кривая

xs =

xs (t)

(s

 

.,

п)

системы

 

 

 

dxs

___ d/„

+

__ 2

 

 

. . . , п ) ,

 

 

(32.40)

 

 

dt

dxs

t <

s=

 

 

 

 

 

 

 

 

 

/= 1

 

 

 

 

 

 

 

 

 

такая, что xs =

xs (s =

1,

. . ., n) при t = 0

и (xx,

. . .,

xn) удовле­

творяет

(32.32),

обладает свойством

xs (t) -> xf'1 (s

=

1,

. .

n)

при

t

+oo.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Вообще говоря, интегральная кривая может достигать опти­ мальной точки в конечное время. Здесь предполагается также, что правые части системы (32.40) удовлетворяют условию Лип­ шица.

Доказательство. Каждая из функций fj, определяемая фор­ мулой (32.32), является интегралом системы (32.40). Следовательно,

3 2 0

любая интегральная кривая, начинающаяся в точке (хг, . . ., хп), удовлетворяющей (32.32), во все время движения будет удовле­ творять (32.32). Вычислим полную производную от / 0 в силу системы (32.40). Тогда получим

— ferad /о. grad, /о) = — (grad /0 — A, grad, f0) =

— (grad< /о)2-

Предположим, что любое движение системы (32.40) продолжимо при t £ [0, оо). Возьмем е > 0 и покажем, что существует такое Т (е), что при t > Г (е) будет

S = = 1

Доказательство этого утверждения осуществляется точно так же, как и в теореме 75, с той лишь разницей, что числа К и б опреде­ ляются с учетом (32.32). Следовательно, теорему 77 можно счи­ тать доказанной.

З а м е ч а н и е . Если

ограничения имеют

вид неравенств

/;

0 (/ = 1...........I),

(32.41)

то, повышая размерность исходного пространства, ограничения вида (32.41) можно свести к ограничениям вида (32.32), полагая

 

?/ = // + 4 +;-(/ =

1. •••• о-

 

 

 

Переменные хп+/ (/ = 1, . . ., /) выбираются так,

чтобы было // = 0.

 

 

П р и м е р .

Пусть функция /„ (xj...........хп) принимает оптимальное значе­

ние в точке ( х |° \

. . ., х*,0))

и не имеет других максимумов или минимумов, однако

имеет стационарные точки,

где grad / 0

обращается в нуль. Для отыскания опти­

мальной точки рассмотрим систему

 

 

 

 

 

 

 

 

 

 

ж

т

' 8 -

1 .........................» > ■

( 3

2 4 2 >

Эта система

имеет интеграл

 

 

 

 

 

 

 

 

 

п

 

 

 

 

 

 

 

 

 

s=i its 4"

{ х 1'

•> х п )‘

 

 

 

Все ее точки покоя располагаются на гиперплоскости

 

 

 

 

Ул. = У% = • • • = Уп = 0.

 

 

 

Возьмем некоторую точку (х1; . . ., х„,

у х...........уп),

причем (уи .

. .,

уп)

выберем так, чтобы вдоль

интегральной кривой

(xs (t),

ys (t)) (s = 1,

. . . , n),

проходящей через эту точку при / =

0, функции возрастали. Так как на любой

интегральной кривой интеграл сохраняет постоянное значение, то вдоль указан­ ной кривой функция f0 будет монотонно убывать. Следует иметь в виду, что числа

(г/j, . . ., рп)

не всегда можно выбрать так, чтобы функции y2s (t) возрастали при

 

П

возрастании

Однако всегда их можно выбрать так, чтобы функция ^ i/s W

 

s=l

2 1 В . И . З у б о

321

возрастала, за исключением того случая, когда (xlt . . ., х„) является оптималь­ ной точкой функции f 0.

Итак, процесс интегрирования системы (32.42) дает возможность минимизи­ ровать функцию f о (х) независимо от наличия у нее стационарных точек.

Аналогичный процесс минимизации можно указать для отыскания оптималь­ ного значения функции /„ ПРИ наличии ограничений вида (32.32). Для этого рассмотрим систему уравнений:

I

/ = 1

(32.43)

4 r = - gvadlf о

Система (32.43) должна иметь функции /у своими интегралами. Это дает воз­ можность определить величины р ь . . ., pv единственным образом из уравнений

S dXs д// _ о (/= 1, .... О-

с—1 dt dxs

Из этой системы находим

М = — C -xR*y,

(32.44)

где М = (рь . . ., щ).

Покажем, что система (32.43) имеет также интеграл

S r f + 2/о.

S— 1

Действительно, полная производная этого выражения в силу системы (32.43) имеет вид

~

Е

у\ + 2?оj =

— 2gradJ /0у + 2 grad* / 0 (у +

ЯМ) =

= 2Rh.y +

2 grad* f0RM =

2 (RC^R* grad f0y — grad* f0RC~1R*y) = 0,

так как матрица С симметрична.

 

 

 

 

 

Возьмем точку (х,, . . ., хп;

y lt . .

., г/п) так,

чтобы были

удовлетворены ус­

ловия (32.32) и чтобы вдоль интегральной кривой

 

 

 

 

(Xs = xs (0;

Ув =

Уз (0) (S =

1 ...........«),

 

проходящей через эту точку при t

 

0, функция

п

y 2s (t) возрастала при воз-

=

^

растании t. Тогда

условия (32.32)

 

 

S = 1

во все

время движения,

будут выполнены

а функция /о будет убывать вдоль указанной интегральной кривой. Этот процесс вычислений дает возможность построить сколь угодно точно оптимальную точку функции f0.

Приложение 1

ПОСТРОЕНИЕ ПРОГРАММНЫХ ДВИЖЕНИЙ В ЛИНЕЙНЫХ УПРАВЛЯЕМЫХ СИСТЕМАХ

Настоящее приложение содержит условие существования про­ граммного движения в линейной управляемой системе, подчи­ няющегося линейным граничным условиям общего вида. Наряду с этим дано полное описание множества допустимых управлений.

Рассмотрим систему дифференциальных уравнений

 

 

i =

/>(/) х +

Q (0 u + F (О-

Требуется найти условия, при выполнении которых существует

непрерывное решение

 

 

 

 

 

 

х =

х (0

(2)

системы

(1),

удовлетворяющее

граничному условию

 

 

 

т

 

 

 

 

 

J[dG(0)]x(0) = H,

(3)

 

 

 

о

 

 

а также требуется найти представление любого

уравнения и =

= и (/),

при

котором

такое решение системы

(1) существует.

Далее все такие управления будем называть допустимыми, если

к тому же будет u (t) £

Ь 2 [0, 74; будем предполагать, что

эле­

менты матриц Р и Q, а также компоненты вектора F являются

вещественными, непрерывными функциями, заданными при

t 6

G [О, Г]. Решение (2) будем разыскивать на том же промежутке.

Элементы матрицы G будем считать вещественными функциями

ограниченной вариации.

Вектор х имеет размерность п, вектор

и — размерность г, а матрица G имеет m строк.

 

Обозначим через Y (t) фундаментальную систему решений для

линейных однородных

уравнений

 

x = P { t ) x .

Положим

т

 

Л ( 0 - f[dG(0)]K(0).

(4)

i

 

21*

323

Обозначим через Z постоянную матрицу, столбцы которой пред­ ставляют собой набор всех линейно-независимых решений урав­ нения

А* (0) Z = 0.

(5)

Теорема 1

Если матрица

т

В + Z* J А (0 У-1 (/) Q (t) Q* (t) [Y-1 (t ) f A* (t) dtZ (6)

о

неособая, то при любом выборе векторной функции F и вектора Н существует решение (2) системы (1), удовлетворяющее условию (3). При этом каждое допустимое управление представимо в форме

 

и (/) == Q* (У-1)*А*С А- V,

(7)

где

С — вектор, определяемый единственным

образом.

Векторная функция V суммируема с квадратом и удовлетво­

ряет

условию ортогональности

(8)

 

т

 

J A (t) У-1 (t) Q{t)V(t) dt = 0.

о

Теорема 2

Если матрица (6) неособая и матрица А (0) имеет ранг п, то каждому управлению (7) отвечает единственное решение (2), удовлетворяющее условию (3). Если же ранг матрицы А (0) есть k < п, то каждому управлению (7) отвечает семейство реше­ ний (2), удовлетворяющее условию (3), зависящее от п k произ­ вольных постоянных.

Теорема 3

Если матрица (6) особая и имеет ранг /, то решение (2), удовле­ творяющее условию (3), существует тогда и только тогда, когда ранг матрицы В' также будет равен I. В этом случае каждое допустимое управление представимо в форме (7), однако вектор С зависит от п k I произвольных постоянных. Каждому та­ кому управлению при k = п отвечает единственное решение (2), при k < п — семейство решений, зависящее от п k произволь­ ных постоянных.

Здесь через В' обозначена матрица, получаемая из В приписы­

ванием дополнительного столбца

 

/ г

\

Bm. k+1- Z* [ \ A (0

(0 F (0 dt - H j .

32 4

Приложение 2

О РАЗРЕШИМОСТИ ЗАДАЧИ ПЕРЕХОДА С ПОМОЩЬЮ ИМПУЛЬСНЫХ УПРАВЛЕНИЙ

Пусть задана система

 

 

 

 

 

 

 

 

 

 

 

 

 

x = P x -l-Q u ,

 

 

 

(1)

где

х — n-мерный

вектор

фазового состояния

системы;

и — г-

мерный вектор

управления;

Р

и

Q — вещественные, непрерыв­

ные

и ограниченные

матрицы,

элементы

которых заданы

при

t

0.

 

 

 

 

 

 

 

 

 

 

 

 

Будем говорить, что управление и в системе (1) импульсное,

если промежуток его задания

[0,

Т ] можно разбить точками 0 =

= tx <

t2 < . . . <

4

= Т

так,

что в каждом

из промежутков

[tj,

tj+1)

(/ =

1, . .

 

k

1)

u

сохраняет

постоянное значение.

При

этом и =

0 при

четном

у

и

и = и [у] при

нечетном

у,

где

и [у] — постоянный вектор пространства Ег. Вектор и [у] условно

можно называть величиной

импульса, а промежуток

[tj, tj+1)

длительностью

импульса.

релейно-импульсных управлений.

Рассмотрим

еще случай

Не ограничивая общности,

можно управление и

определить

как релейно-импульсное, если это управление является импульс­ ным и если компоненты векторов и [у] могут принимать лишь

два значения:

-±_ 1. Обозначим через

У фундаментальную систему

решений

для

линейной

системы

 

 

 

 

 

 

 

 

 

х =

Рх

 

 

 

(2)

и положим В =

Y 1Q.

 

 

 

 

 

 

 

 

 

 

 

Теорема 1

 

 

 

Чтобы существовало импульсное управление, переводящее

систему

(1)

из любого

начального

состояния х = х0

при t

= 0

в конечное состояние х — 0 при t =

Т, необходимо и достаточно,

чтобы на

промежутке

[0, 74

существовали

такие

точки

т1(

т2, . . . .

tj,

что

среди

векторов, образующих

столбцы матриц

В (ts) (s

=

1, . .

., /), имеется

п линейно-независимых.

 

З а м е ч а н и е .

Можно показать,

что точки Tt, . . .,

тj существуют тогда

и только тогда, когда строки матрицы линейно-независимы в промежутке (0, Т),

что, в свою очередь, эквивалентно неособенное™ матрицы Jт ВВ* dt, где В*

0

транспонированная по отношению к В матрица. Это показывает, что условие теоремы 1 может быть выражено аналитически через матрицы Р и Q.

325

Теорема 2

Чтобы существовало релейно-импульсное управление, пере­ водящее систему (1) из любого начального состояния х = х 0 при t — 0, xjj < h2 в состояние х = 0 при t = Т, необходимо и до­

статочно, чтобы было выполнено условие теоремы 1. Здесь h — достаточно малая положительная величина.

З а м е ч а н и е . Если положение равновесия системы (2) асимптотически устойчиво по Ляпунову, то при выполнении условий теоремы (1) при помощи релейно-импульсного управления эту систему за конечное время можно перевести

в состояние х = 0

из любого начального состояния х = х0 при t =

0.

Рассмотрим

теперь

квазилинейную систему

 

 

х =

Рх-(- Pu + pF (/, х, и),

(3)

где р, — малый параметр; F — вещественная непрерывная, п- мерная вектор-функция, непрерывно дифференцируемая по отно­ шению к компонентам вектора х.

Теорема 3

Если выполнены условия теоремы 1, то можно указать такие два числа / г > 0 и р 0 > 0 , что система (3) может быть переведена из любого начального состояния х = х0 при t = 0, х^ < h2 в со­

стояние х = 0 при t = Т для любого выбранного |р | < р 0.

П р и л о ж е н и е 3

К ТЕОРИИ СУЩЕСТВОВАНИЯ РЕШЕНИЙ КРАЕВЫХ ЗАДАЧ ДЛЯ СИСТЕМ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ

Рассмотрим систему линейных дифференциальных уравнений

х = .Рх + Q.

(1)

Будем считать, что элементы квадратной матрицы Р и компоненты

вектора Q заданы при t £ [0,

Т], вещественны

и непрерывны.

Порядок системы (1) п.

функцию

 

Требуется найти векторную

 

х = х ( / ) ,

(2)

326

непрерывно дифференцируемую, заданную при /£ [О, Г], удовле­ творяющую системе (1) и совокупности краевых условий

 

 

т

 

 

 

 

 

 

 

j [dG(т)] х (т) =

R .

 

 

(3)

 

 

о

 

 

 

 

 

Матрица G имеет m строк и п столбцов, вещественна,

задана при

t

[О, Т ]

и элементы ее

являются

функциями

ограниченной

вариации на этом сегменте;

R является /и-мерным вещественным

вектором.

 

 

 

 

 

 

 

Наряду с системой (1) рассмотрим систему однородных урав­

нений

 

х = Рх

 

 

 

(4)

 

 

 

 

 

 

и

фундаментальную систему решений

Y =

Y (t),

т

Y = Е при

t

= 0, где

Е — единичная

матрица.

Если

через

обозначить

текущий момент времени и через

t — некоторый фиксированный

момент, то по формуле Коши найдем

 

 

X

 

х (т) = Y (т) Y-1(0 х (0 +

} К.(т) Y-1(6) Q (0) dQ.

(5)

 

t

 

Предположим, что решение (2) системы (1), удовлетворяющее условию (3), существует. Тогда это решение представимо в форме

(5). Подставляя (5) в (3), находим

т

\{dG(x)]Y(x) dxY-1 (t)x(t) +

о

Тт

+ \ [dG (т)] { Y (т) Y ' 1(0) Q(0)d0 = R.

(6)

о

t

 

Сделаем в левой части (6) следующие преобразования. Разобьем интеграл от t до т на два:

и затем в интеграле Jт dG (т) JX сделаем замену порядка интегри-

оо

рования. Далее, полагая

г

Л(т)=|[Л?(е)]К(0),

X

327

найдем

о

 

 

А0

 

 

У 1 (х) Q (т) dx +

 

 

 

(т) У"1(т) Q (г) dr — R.

(7)

Здесь

через z обозначен вектор У-1 (t) х (t), а через

А 0— ма­

трица

А (0).

 

 

Обозначим через Л матрицу, столбцы которой являются

линейно-независимыми

решениями системы уравнений

 

 

 

л ; л = о .

(8)

Система (7) может рассматриваться как система линейных алгебраических уравнений, служащая для определения вектора г. Известно, что эта система имеет решение тогда и только тогда, когда будет выполнена система равенств

т

 

{ А*А (т) У-1 (т) Q (т) dx = A*R.

(9)

о

 

Теорема 1

Для того чтобы существовало решение (2) системы (1), удовле­ творяющее условиям (3), необходимо и достаточно, чтобы было выполнено условие (9).

Рассмотрим систему нелинейных дифференциальных уравне­

ний

 

 

 

 

 

 

х — Рх -(- F -f pG (t, х, р).

(10)

Здесь

векторная

функция

G предполагается заданной

при

t 6 [0,

Т\, х 6 Еп,

| (х | < р 1(

вещественной непрерывной

и не­

прерывно дифференцируемой относительно компонент вектора х и параметра р.

Предположим теперь,

что m = п и что матрица А 0 неособая.

Заменив в (7) Q векторной функцией F + pG, получим

т

1(т) [F (т) -f pG (т, х, р)] dx +

х = У (t) j У

о

 

I

+ у (О^о"1|R — | Л (т) Y~X(т) [F (т) + pG (т, х, р)] dx \ . (11)

328

Теорема 2

Если квадратная матрица Л 0 порядка п является неособой, то существует такое р 0, что при любом | р | < р 0 система (10) имеет единственное непрерывное решение (2), удовлетворяющее усло­ вию (3).

Нелинейный оператор, стоящий в правой части (11) в про­

странстве непрерывных векторных

функций х (/), заданных

при

t 6 (0,

Т\,

удовлетворяет условию сжатия при

любом достаточно

малом

| р |

< РоСледовательно,

существует

единственное

не­

прерывное решение х = х (t) системы (11); дифференцированием этой системы можно показать, что это решение непременно удовле­ творяет также дифференциальным уравнениям (10). Подставляя (11) в (3), найдем, что это решение будет также удовлетворять усло­ вию (3).

Пусть по-прежнему m = п и матрица Л 0 имеет ранг (п — k),

где k > 0. Тогда, подставляя в (7) и (9) вместо Q вектор

F + pG,

находим

 

 

 

х ( 0

т

 

 

= П О YJ1(т) [F (т) + pG уг, х, р)] dx +

 

 

о

 

 

+ У (9Л 1с + у (9Л2 R

А (т) У-4 (т) (F (т) +

 

 

+ pG (Т, х, р)) dx ,

( 12)

где А ! и А 2 — некоторые постоянные матрицы, определяемые единственным образом на матрице А 0; С — ^-мерный вектор с про­ извольными вещественными компонентами си . . ., ck.

Из (9) получим

т

 

} А*Л (т) у -1 (т) [F (т) + pG (т, х, Р)1 dx == A*R.

(13)

о

 

Далее будем исходить из того, что для системы (10) при р = 0 существует решение (2), удовлетворяющее условию (3). Тогда условие (13) примет вид

т

 

j А*Л (т) У '1 (т) G (т, х, р) dx = 0.

(14)

329

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