следовательно,
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). Следовательно,
любая интегральная кривая, начинающаяся в точке (хг, . . ., хп), удовлетворяющей (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 |
возрастала, за исключением того случая, когда (xlt . . ., х„) является оптималь ной точкой функции f 0.
Итак, процесс интегрирования системы (32.42) дает возможность минимизи ровать функцию f о (х) независимо от наличия у нее стационарных точек.
Аналогичный процесс минимизации можно указать для отыскания оптималь ного значения функции /„ ПРИ наличии ограничений вида (32.32). Для этого рассмотрим систему уравнений:
I
4 r = - gvadlf о
Система (32.43) должна иметь функции /у своими интегралами. Это дает воз можность определить величины р ь . . ., pv единственным образом из уравнений
S dXs д// _ о (/= 1, .... О-
с—1 dt dxs
Из этой системы находим
где М = (рь . . ., щ).
Покажем, что система (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 |
|
Обозначим через Z постоянную матрицу, столбцы которой пред ставляют собой набор всех линейно-независимых решений урав нения
Теорема 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 . |
Приложение 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.
Теорема 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 заданы при t £ [0, |
Т], вещественны |
и непрерывны. |
Порядок системы (1) п. |
функцию |
|
Требуется найти векторную |
|
х = х ( / ) , |
(2) |
непрерывно дифференцируемую, заданную при /£ [О, Г], удовле творяющую системе (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
найдем |
о |
|
|
А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)
Теорема 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) |