Файл: Иоффе, А. Д. Теория экстремальных задач [учеб. пособие].pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 15.10.2024

Просмотров: 132

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

436 ГЛ. 10. ПРИЛОЖ ЕНИЕ ТЕОРИИ К ЗАДАЧАМ

того как возбуждение будет снято) его энергия достигла бы заданной величины за минимальное время. Сейчас мы построим полный синтез задачи при любой величине О < е < 1 п затем проведем сравнение точного решения

п приближенного решения

(см. Евтушенко [1]).

Положив

х =

у, приведем задачу к стандартному

виду задачи оптимального управления:

 

 

ty -> inf;

х =

у,

у = — (1— еи)х (0 < в < 1),

х (0) =

*0,

х (0) =

у0, х2(7,) + х2(7,) = 1,

 

 

^

1 .

Применим принцип максимума Понтрягнна в гамильто­ новой форме (теорема 1 § 2.4). Функция Н имеет вид

Н — ру — <7(1— е«) х

и, следовательно, сопряженная система такова:

Р =

(1 — е«) <7,

<7 —

Р-

Для р и <7 выполнены условия

трансверсальности

р М= —и* (О,

qiti) =

— v-y(U),

пли

q + (1 — ей) <7 =

0,

|

 

<7(/|) =

~ р х (7,),

q (7,) =

рх (7,). J

Видно, что функции x(t) и <7(7) удовлетворяют одному и тому же уравнению.

Принцип максимума дает такое выражение для оп­

тимального управления:

 

 

 

 

 

|

1,

если

qx >

0,

(2)

\

0,

если

qx <

0,

 

или, короче, u(t) = Q(qx),

где

0(Я) — функция Хеви­

сайда.

 

 

 

 

 

Найдем единственное решение написанных соотно­

шений.

Ясно,

что р Ф 0, ибо

иначе

из

(1) мы получили

бы, что

q =

0, а значит, и

р = q ^

0, чего не может

оыть в силу принципа максимума.

Значит, можно счи­

тать, что р =

1.

 

 

 


 

§ 10.3.

ВОЗБУЖДЕНИЕ

ОСЦИЛЛЯТОРА

437

Далее,

легко

понять,

что

u(ti — 0 ) = 1 ,

ибо иначе

энергия Е достигалась

бы за

время меньшее, чем t{

(ибо, если и = 0,

то энергия

не меняется). Следова­

тельно, из

(1) и (2)

получаем,

что

 

0 < q {tt) х (*,) = — л: (t,) х (/,).

Таким образом, точка (x(/i), x (t,) ) лежит во второй или четвертой четвертях фазовой плоскости. Для опре­ деленности допустим, что вектор (at(/i) , х(^])) лежит в четвертой четверти:

х((,) = cos а,

х(tj) = sin а,

— л/2 < а < 0.

Тогда, в

силу (1)

(напомним,

что ц = 1 )

 

 

Я (ti) = —х (/,) =

cos (а + л/2),

 

 

q (/,) = х (tx) = sin (а +

я/2).

 

Получилось,

что вектор

(q(ti), q{t\)) повернут

относи­

тельно вектора

(x(i{), x{t\))

на угол л/2 против часовой

стрелки.

теперь,

что

в

невозбужденном состоянии

Отметим

(когда и =

0)

точки

(x{t),x(t))

и (q(t),q(t))

совер­

шают движение по кругам в фазовой плоскости, а в возбужденном состоянии — по эллипсам х2-{-х2/( 1— е) = = а, в обоих случаях движение совершается по часо­ вой стрелке.

Будем решать нашу задачу «с конца», более того,

оставим нашу систему в возбужденном состоянии еще на

некоторое время (после момента /i) до тех пор, пока х

не станет равным нулю. Обозначим этот

момент через

Т > t\. Точка (х(Т), х{Т)) расположена

на вертикаль­

ной осп. Допустим, что в момент времени Т вектор

iq(T), q(T)) занял

такое положение:

 

q (Т) =

pcosco,

 

<НТ) = psinw,

 

со =

со

(а).

(Впоследствии мы

вычислим

со (а ), и окажется, что

о) (а) > 0 .) Теперь будем двигаться «в обратную сторону».


43 8

ГЛ. 10. ПРИЛОЖ ЕНИЕ ТЕОРИИ К ЗАДАЧАМ

Сначала

(при Т t <

р) наш

осциллятор возбуж­

ден, а

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

и точка

(х,х), и точка (q,q)

совершают движение по эллипсам. Первым вертикаль­

ной оси достигнет вектор (q,q),

а вектор

(х,х)

к этому

моменту будет

образовывать с

осью х

угол,

равный

(— 1)<в(а). В

этот момент

произошло последнее

переключение, до которого жесткость равнялась единице и, следовательно, оба вектора совершали равномерное круговое движение. При этом вращении (но при движе­ нии в обратном направлении!) первым достигнет верти­ кальной оси вектор (х,х), и это будет момент предпо­ следнего переключения. В этот момент вектор (q, q) бу­ дет образовывать с горизонтальной осью угол ш(а) + я и т. д.

Мы приходим к тому, что при своем движении жест­ кость переключается тогда и только тогда, когда фазо­

вая траектория пересекает либо вертикальную ось

тогда жесткость

меняется с

1 — е на

1),

либо прямую

х = (— tgco(a))A:

 

тогда

жесткость

меняется

с

1

на

1 — е).

 

вычислить

кривую

 

ш(а).

Вблизи

точки

 

Остается

 

 

t — ti функция x(t) имеет вид

 

 

 

 

 

 

 

 

 

 

 

^(0 =

^ c o s (V l

— е (t — t^ +

y),

 

 

(3)

откуда,

используя

равенства

x(t^) =

cos a,

я (г^) =

sin а,

сразу получается,

что

 

 

 

 

 

 

 

 

 

 

 

 

/ Г

1

-

е

 

sin Y —

s i n

 

a

COS у — c o s

а

(4)

 

/

 

 

 

 

A /

l

 

-

е ’

~Л~

 

Вблизи

точки t =

tl

функция

q{t)

имеет

вид

q(t) =

=

£ c o s (]/l

— е (t— ^) +

б). Учитывая

равенства q(t,) =

=

— sin a, q (t^ =

cos a,

будем

иметь:

 

 

 

 

 

В

/ 1

е

s i n 2 a

 

 

c o s

a

, COS 6

 

s i n

a

(5)

/

1

-

,

sin 6 =

 

 

 

 

 

В

1•

 

s

 

 

 

В / l —е

 

 

 

 

Момент времени t2, когда происходит последнее пере­ ключение, определяется из условия q (/2) = 0, q (t2) > 0, т. е.

/ 1 — e(f2 — /,) + б = — я/2.


 

§ 10.3. ВОЗБУЖ ДЕНИЕ

ОСЦИЛЛЯТОРА

439

Отсюда

и из (3), (4), (5) находим, что

 

х (t2) = A sir; (Y — б );

 

 

 

В У 1

 

 

X (t2) — А У

1— 8 COS (y — 6) =

 

6

 

8 sin a cos а

е sm а cos а

( )

 

 

 

ВУ\ — ъ

)А18 s i n 2 а

j

Формулы (6) определяют точку последнего переклю­ чения в зависимости от а, т. е. от положения фазовой точки в момент t\. При изменении а от 0 до —п/2 точка {x(t2), x(t2) ) описывает на фазовой плоскости кривую /, концы которой лежат на оси абсцисс. Из формул (6) находим, что

t g ы ( а ) =

х (t2)

е sin а cos а.

x(t2)

 

Проведем сравнение точного и приближенного реше­ ний. Обозначим xjx через tgy. Тогда оптимальное управление задается так:

и = 9 (tg со (а) tg у) = 0 (е sin а cos а tg у).

В то же время, как это следует из работы Евтушенко [1], приближенная формула для оптимального управле­ ния такова: й = 0 (—tg у).

Комментарий к гл. 10. Вариационный принцип Ферма был вы­ двинут в 1669 г., принцип Гюйгенса — в 1690 г. (Гюйгенс [1]). Оп­ тико-механическая аналогия привела И. Бернулли к решению задачи о брахистохроне. Серия задач, включающая задачи 1—3 § 10.1, подробно исследована в курсе анализа Гурса.

Задачи на быстродействие для объектов, описываемых диффе­ ренциальными включениями, изучались Болтянским [3] и Благо­

датских [1].

О теореме Хелли см. Данцер, Грюибаум, Кли [1]. Доказатель­ ства теоремы, опирающиеся на теорему об очистке, содержатся в ра­ ботах Карлина [1] и Пшеничного [4]. Решение задачи Юнга, сходное с нашим, содержится в работе Джона [1].

Задача из п. 10.3.1 рассматривалась Евтушенко [1]. Решение, изложенное у нас, получено Барыкиным.


ЗАДАЧИ

Раздел открывается циклом старинных задач. Далее следуют

разнообразные задачи на экстремум из

алгебры

и геометрии.

Часто их решают без применения

теории.

Однако применение тео­

рии почти во всех случаях дает

решение

не более

сложное, чем

то, которое получается с помощью искусственных приемов. В за­ ключительной части раздела приведены стандартные задачи из ва­ риационного исчисления и оптимального управления. Стандартные задачи конечномерного анализа и линейного программирования мы не приводим. Их можно почерпнуть из многих задачников, см., на­ пример, Демидович [I] и Заславский [1]. Упомянем еще два сбор­ ника экстремальных задач из геометрии: Зетель [1] и Шклярский, Ченцов, Яглом [1].

Задачи группируются по определенным темам. В каждой теме есть задачи, снабженные решением; в конце раздела есть краткие указания к решениям некоторых задач. Решение большинства за­ дач этого раздела получается применением принципа Лагранжа.

1. К л а с с и ч е с к а я и з о п е р и м е т р и ч е с к а я з а д а ч а . Среди кривых на плоскости, имеющих заданную длину, найти ту, которая ограничивает наибольшую площадь.

2. На стороне ВС треугольника АВС найти точку Е так, чтобы параллелограмм ADEF, у которого точки D и F лежат соответ­ ственно на сторонах АВ и АС, имел наибольшую площадь (Евклид).

3. Среди шаровых сегментов, имеющих заданную площадь по­ верхности, найти сегмент наибольшего объема (Архимед).

4.Найти кратчайшее расстояние от заданной точки на плоско­ сти до заданного конического сечения (Аполлоний).

5.Среди всех «-угольников, имеющих заданный периметр, найти «-угольник наибольшей площади (Зенодор).

6.На данной прямой найти точку С так, чтобы сумма расстоя­

ний от С до точек А и В была минимальной (Герои).

7.Разделить число 8 на две части так, чтобы произведение их произведения на разность было минимальным (Тарталья).

8.Найти наклонную прямую, двигаясь по которой без трения под воздействием силы тяжести, материальная точка достигнет за­ данной линии за кратчайшее время (Галилей).

9.Среди цилиндров, вписанных в шар, найти цилиндр с макси­ мальным объемом (Кеплер).

10. Найти прямоугольный

треугольник наибольшей площади,

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

равна заданному числу (Ферма).

11. Среди конусов, вписанных в шар, найти конус с максималь­ ной боковой поверхностью (Лопиталь).