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

книги из ГПНТБ / Иоффе, А. Д. Теория экстремальных задач [учеб. пособие]

.pdf
Скачиваний:
61
Добавлен:
22.10.2023
Размер:
15.51 Mб
Скачать

Г л а в а 10

ПРИЛОЖЕНИЕ ТЕОРИИ К РЕШЕНИЮ ЗАДАЧ

§ 10.1. Задачи геометрической оптики

На примере задач геометрической оптики легко про­ следить те геометрические и физические идеи, которые лежат в основании формализма Гамильтона — Якоби — Веллмана в вариационном исчислении и теории опти­ мального управления, а также метода возмущений в теории экстремальных задач. Это связано прежде всего с тем, что математическая формулировка задач геомет­ рической оптики необычайно проста и не заслоняет их физической природы. Поэтому каждый результат допу­ скает физическое истолкование и, наоборот, из есте­ ственных физических соображений легко угадывается форма необходимых математических теорем. С другой стороны, многие интересные задачи можно интерпрети­ ровать как задачи геометрической оптики, — таковы, на­ пример, задачи оптимального быстродействия для авто­ номных объектов.

Большая часть рассуждений в начале параграфа не

претендует

на строгость. Их цель — показать, как, от­

талкиваясь от простых и интуитивно ясных соображе­

ний, можно прийти к содержательным математическим

результатам.

Ферма. Пусть имеются прозрачная

10.1.1.

Принцип

неоднородная и, для простоты, изотропная среда и

источник света, расположенный в точке

x0 = (x,J, х

xj).

Обозначим

через v(x)

скорость распространения

света

в точке х =

(х1, х2, х3).

Спрашивается,

по какому

пути

будет идти световой луч из точки х0 до точки х,? Отве­ чая на этот вопрос, Ферма выдвинул следующий вариа­ ционный принцип, который был первым принципом та­ кого рода для физической проблемы. Согласно прин­ ципу Ферма, свет избирает такую траекторию, чтобы

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

время его движения по ней

*i

 

= .f v f e -

(1)

Хо

 

было не больше, чем время движения по любой кривой,

соединяющей точки х0 и xi. В дальнейшем этот принцип получил уточнение: траектория распространения света не обязательно доставляет минимум функционалу Т, но она обязана быть стационарной точкой этого функцио­ нала. В такой форме принцип оказывается верным и для неизотропных сред, где абсолютная величина скорости зависит от ее направления.

В общем случае неизотропной среды ее оптические свойства в точке х описываются множеством V(x), ко­ торое представляет собой совокупность всех возможных векторов скоростей света в точке х. Само множество ^(л:) при этом меняется от точки к точке (неоднород­ ность среды). Это означает, иначе говоря, что задано

многозначное отображение х —*V(x). При этом траек­ тория движения светового луча от точки х0 до точки хЛ в соответствии с принципом Ферма может рассматри­ ваться как решение следующей задачи на оптимальное быстродействие:

Г —> inf;

xe=V(x),

х (t0) =

х0, х (^) = *i.

(2)

Нижнюю грань

в задаче

(2)

обозначим Т(хи х0).

 

Приведем несколько примеров. Задача о брахистохроне, сфор­

мулированная во

введении, сводится к минимизации функционала

 

(*i. Vi)

У 1+'■J—

 

 

 

 

 

dx.

У> 0.

 

 

(*о. №>)

Vv

 

 

 

 

 

 

 

 

 

Мы видим, что это задача плоской геометрической оптики, где ско­

рость света задается равенством v (х, у) = Y У ■

Известно, что плоскость Лобачевского допускает реализацию в

полуплоскости

у > 0,

если элемент длины задать формулой

ds2 = y~2(dx2+

dy2).

Полученное так риманово многообразие на­

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

<*|. !Л>

I

У

 

V i + t f 2 dx -> inf.

(*о, Уо)

§ 10.1. ЗАДАЧИ ГЕОМЕТРИЧЕСКОЙ ОПТИКИ

42 3

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

оптики с

v (х, У) = У-

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

V(X, у) = у~1-

Рассмотрим задачу оптимального управления автономным объ­

ектом:

 

 

 

Т -> inf;

 

 

 

 

х =

(фх, и), и е U,

 

 

 

х ( 0 )

=х0, х ( Т ) = Х \ .

 

 

Ее

можно редуцировать

к

задаче

(2), если

положить Р(л:) =

= ср(*.U ) .

 

 

 

 

 

 

В задачах геометрической

оптики

удобно

считать,

что все множества V(x)

выпуклы и содержат нуль. Пер­

вое

предположение,

хотя и кажется поначалу

слишком

сильным, на самом деле не является ограничительным. Подобно тому, как в теореме Боголюбова первоначаль­ ный интегрант можно заменять выпуклым по производ­ ной, так и в задачах на быстродействие можно, не изменяя значения задачи, «овыпуклять» множество воз­ можных скоростей, заменяя его выпуклой оболочкой.

10.1.2. Принцип Гюйгенса и уравнение Гамильтона — Якоби. Для задач геометрической оптики этот принцип можно записать в очень естественной и красивой форме.

Пусть А(х)

и В( х) — многозначные отображения из

Rn

в R". Через А ° В обозначим, как и в § 8.1, суперпози­

цию отображений Л и В:

 

 

 

 

 

 

 

 

( А о В ) ( х ) =

U

А(1).

 

 

 

 

 

 

 

1 ^ В

(х)

 

 

 

 

Рассмотрим

однопараметрическое

семейство множеств

Тогда

 

£Гt (*0) =

es R" |Т (х\ *0) < /} .

 

 

t + s (*о) = ( Г

 

W

= (

T s ° Г

 

(х0).

(3)

f

t ° 3 ~ s )

t)

Соотношение

(3) и есть принцип Гюйгенса.

 

 

Действительно, вложение

{ Т s ° ГГt)

(х0) с= &~t+s{x0)

означает,

что, если

точки

г свет

может

достичь

из

точки х за время s, а точки х из точки х0— за время t, то точка z достижима из точки х0 по крайней мере за

время

t +

s. С

другой

стороны,

вложение

5r'i+s(x0)c г

сд ( 3 ~ s ° & ~ t ) { x q )

означает, что,

если

точки z

можно до­

стичь

из

х0 за

время

/_+ s,

то

на

траектории можно

42 4

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

отметить точку х, в которую свет придет за время /. Но тогда из х в z можно дойти за время, не большее s. Физическое содержание принципа Гюйгенса состоит в том, что он позволяет описывать процесс распростране­ ния света при помощи волнового фронта — границы множества х(х0) .

Пусть теперь

Т(х; х0) = inf [t > 0 <= Т ,(х0)}.

Функция Т(х,х0), называемая обычно оптической дли­ ной пути от точки Xq до х, является, очевидно, S-функ­ цией задачи (2), если в качестве возмущения рассмат­ ривать точки, в которые свет должен попадать в конце траектории. Покажем сейчас, как с помощью принципа Гюйгенса и простых (хотя и не абсолютно строгих) рассуждений можно получить уравнение Гамильтона — Якоби для задачи (2).

Допустим, что оптическая длина Т(х; х0) является дифференцируемой функцией переменного х. (Разу­ меется, так хорошо дело обстоит далеко не всегда. На­ пример, в задаче о наименьшей поверхности вращения функция Т не является всюду дифференцируемой.) Предположим далее, что любая точка х некоторой об­

ласти

U cz R?!

соединяется с точкой х0

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

траекторией светового луча (исходящего

из х0), кото­

рую обозначим х (- ,х 0).

Рассмотрим

точку х* е U.

Ве­

дущая

в нее

траектория

светового

луча

есть х * (-,х 0),

а оптическая

длина пути

из х0 в х*

есть

7'(х*;х0) =

Г*.

За последующее время At свет, двигаясь по траектории х*(-,Хо), пройдет путь х*(Г*; х0) At = у* At (если пре­ небречь малыми более высокого порядка). При этом получается равенство

Т (х, + vt At; ха) = Т (х„; х,) + At = T t + At.

Но если бы в точке х* свет избрал другую скорость V, то по принципу Гюйгенса и по определению оптической длины получилось бы, что

Т (х, +

v At; х0) Т (х„; х0) +

At.

В итоге получаем уравнение

 

 

At = sup (х, -f

At; х0) — Т (х,;

х0)) =

 

PEEV (ДЩ

=

sup

(Г (х,; х0) I у)) м

 

 

 

v(BV(*„)

 

§ IO.f. ЗАДАЧИ ГЕОМЕТРИЧЕСКОЙ ОПТИКИ

425

или, если провести это рассуждение для произвольной точки х е U,

sup (Т'(х-,х0) |и) = 1.

Если вспомнить определение опорной функции множе­ ства А:

s (л И ) = sup (ц 11),

то станет ясно, что получилось следующее уравнение в частных производных для оптической длины:

s(T' {х\ х0) |V ( x ) ) = 1.

(4)

Это и есть уравнение

Гамильтона Якоби

для задач

геометрической оптики.

 

 

 

 

Приведем пример. Если среда изотропна, т. е. если

 

V (х) — S (0, v (*)),

где

5 (|. г) = {х 11 л: — 11 =

г},

ТО

 

 

г I ч I

 

S (n I s (0, г )) =

 

и уравнение Гамильтона — Якоби приобретает вид

 

у ^т ГдТ/

((х\; халг)0 )\2_

 

(5)

'•?*\

дх1

}

V

I

 

 

 

 

10.1.3. Принцип максимума. Рассмотрим задачу опти­ мального быстродействия автономным объектом:

7’ —>inf; х = ср(х, и), u<=U, х(0) = х0, х{Т) = хи

и обозначим, как и раньше, V(х) = ср(х, U). Пусть (**(•), «*(• ))— оптимальный управляемый процесс в этой задаче. Тогда в силу принципа максимума най­ дется такое нетривиальное решение р(() сопряженного уравнения

р =

— Ф* (*. (0,

«. 00) Р,

(6)

что

и* (0)) = max (/) |<р {X, (/), и)).

(60

(р (t) (X. (0,

Положим

u^U

 

 

 

 

о (р, х) =

S I V (х)) =

sup |<р (х, и)).

 

 

 

u^U

 

42 6

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

Если

множество U компактно, то в силу теоремы 3 из

§ 4.4

функция х —>о(р, х)

регулярно локально выпукла

и соотношения (6) и (6')

можно переписать в виде

р <= дха (р (t), х, (0),

хе дра (р (t), х, (0).

Это естественная форма принципа максимума для задач геометрической оптики, записанных в форме (2).

10.1.4. Задачи.

З а д а ч а 1. Г е о д е з и ч е с к и е на п о л у п л о с к о ­

ст и Пу а н к а р е .

Всоответствии со сказанным ранее, нужно решить

такую задачу:

(*,. г/.)

,---------

Г

V l + y '2- dx-*mi\ у > 0.

J

У

(*о. Уо)

 

Интегрант здесь не зависит от х, и следовательно, урав­ нение Эйлера допускает интеграл энергии (см. § 2.2)

у V I у' = С,

т. е. приводится к такому дифференциальному уравне­ нию первого порядка:

у dy

dx.

У С2 - у2

 

Его интегрирование дает следующее выражение для экстремалей:

р2 + (х + С,)2 = С2.

Таким образом, экстремалями являются полуокружно­ сти с центром, лежащим на оси у — 0. Через любую пару точек (х0, у0) и {хх, у х), х0 фх\, проходит един­ ственная экстремаль, которую многими разными спосо­ бами можно включить в поле экстремалей, покрываю­ щее всю полуплоскость у > 0. Например, можно взять пучок полуокружностей, имеющих общую точку на оси у = 0. Интегрант в этой задаче, как и в любой изотроп­ ной задаче геометрической оптики, квазирегулярен. По­ этому из теоремы 3' § 7.4 следует, что найденная экст­ ремаль является решением поставленной задачи.

§ 10.1. ЗАДАЧИ ГЕОМЕТРИЧЕСКОЙ ОПТИКИ

427

З а д а ч а 2. Б р а х и с т о х р о н а .

Не раз говорилось уже, что задача о брахистохроне допускает такую формализацию:

(*.. гл) ,--------- -

\

(*о. Уо)

Снова можно написать интеграл энергии

y j l ! L = d X. V а2 —у

Это уравнение легче всего проинтегрировать с помощью подстановки

у = a2sin2 (t/2).

При этом решение находится в параметрической форме

* =

-тг it — sin 0 + с,

У =

- у

(1 — cost).

Это — циклоиды.

Как

и в первой задаче, можно

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

сечения

с осью у — 0.) Так что и в этой задаче в силу

теоремы

3' § 7.4

экстремали являются решениями.

З а д а ч а 3.

М и н и м а л ь н а я п о в е р х н о с т ь

в р а щ е н и я .

 

Эта задача записывается так:

 

(хи г/,)

________

 

J

у V\ + у'2 dx-> in f; у ^ О .

 

(х„. у„)

 

Мы ее интерпретировали выше как задачу плоской гео­ метрической оптики со скоростью, обратно пропорцио­ нальной у. Если в предыдущих двух задачах скорость

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

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

y = c V l + y ' 2

имеет решение

У — с ch ( -7 + d)

Продолжим исследование в симметричном случае, ког­

да х0 =

хи у0 = у\. Тогда решения уравнения Эйлера

приобретают вид

y = c c b j .

Таким образом, экстремалями являются цепные линии. Легко понять, что их совокупность двукратно покрывает внутренность угла, обра­ зованного прямыми, про­ ходящими через начало координат и касатель­ ными к кривой у — due (рис. 13). Вне этого угла провести допустимую цеп­ ную линию не удается.

Если

точки (*о, Уо) и

(хи У\)

лежат на сторо­

нах угла, то через них проходит одна цепная ли­ ния, если внутри — то две, одна нз которых ка­ сается сторон угла, а дру­

гая — нет. Вторую нз этих экстремалей легко окружить полем, и следовательно, по теореме 3' § 7.4 она дает локальный минимум. На первой же экстремали не вы­ полнено условие Якоби, ибо близкие к ней экстремали пересекаются с ней в точках, близких к тем, в которых она касается сторон угла.

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

§ 10.1. ЗАДАЧИ ГЕОМЕТРИЧЕСКОЙ ОПТИКИ

429

шем задачу 3 как вариационную задачу в параметриче

ской

форме

 

 

 

 

 

1

 

 

 

 

 

[

У V х 2+

у2 dx

ini;

 

 

о

 

 

 

 

 

х (0) = х0, у( 0) = Уо,

Х(1) =

х 1,

у{\) = у и

где

т — некоторый

параметр, а х и

у — производные

по т. Система уравнений Эйлера в этой задаче выгля­ дит так:

 

d

УУ

 

Y x 2 +

ft

=

о,

 

dx

 

 

V * 2 +

У2

 

 

 

 

 

 

d I

ух

= 0.

 

 

 

 

dx V V х 2 + у2

 

 

 

 

 

 

 

 

Рассмотрим

интеграл

второго из

этих

уравнений

 

 

 

1/Х

С.

 

 

 

 

 

 

=

 

 

 

 

 

\ х 2+ у2

 

 

 

 

Если с ф 0,

то

у ф 0

и,

значит,

у >

0,

х ф 0 и, следо­

вательно, х имеет постоянный знак. Поэтому все сво­ дится к рассмотренному уже случаю. Предположим те­ перь, что с — 0. Тогда на всякой экстремали в каждой точке либо у = 0, либо х — 0. Легко видеть, что такая экстремаль состоит из трех участков, на первом из ко­

торых

х =

0, а у

меняется

от у0 до нуля, на

втором

у = 0 ,

а х

меняется от х0 до x h и на третьем х — 0 , а у

меняется от нуля до у0.

экстремали двух

н

Таким

образом,

имеются

видов:

цепные линии и ломаные экстремали, состоящие из двух вертикальных кусков и участка х0 sg х ^ х и у = 0. «Оптическая длина пути» по экстремалям первого вида, соединяющим точки (—х, у) и (х, у), равна, как не­

трудно подсчитать, у2Ф(х/у) =

с2(р(х/с), где

I

 

<p(g)==2jch 2xdx,

y = c c h ~ ,

о

 

а по экстремалям второго вида у2. Эти функции при­ нимают одинаковые значения на сторонах угла y ^ k 0\x\, целиком лежащего внутри угла, образован­ ного касательными к кривой у = chx, проходящими че­ рез начало координат.

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

Функция

S (х, у) = min (у2Ф (х/у), у2)

оказывается 5-функцией нашей задачи. Мы видим, что она является ломаной, т. е. не непрерывно дифферен­ цируемой.

З а д а ч а 4.

О б о б щ е н н а я

и з о п е р и м е т р и-

ч е с к а я з а да ча .

 

 

Рассмотрим

такую задачу:

 

 

V2

[ {xv yu)dt~* sup;

х =

и, ij = v,

x{Q) =

6

у(0) = у{Т),

 

(9)

x(T),

Т — фиксировано,

 

 

{и, t i ) e 4 c

R2.

 

Она несколько отличается от задач геометрической оптики, ибо в (9) — другой функционал, но вместе с тем имеет с ними много общего. Функционалом в (9) яв­ ляется площадь. Множество А можно интерпретировать как совокупность локальных скоростей материальной точки. Таким образом задача (9) состоит в отыскании такой замкнутой кривой, ограничивающей наибольшую

площадь, которую

материальная точка

может обойти

за заданное время.

Если А — единичный

круг, то полу­

чаем классическую

изопериметрическую

задачу. Если

А — сдвиг единичного крута на постоянную величину, то получается задача Чаплыгина о самолете, который дол­ жен за заданное время облететь наибольшую площадь при наличии ветра и т. д.

Задача (9) есть задача оптимального управления и ее можно решать с помощью принципа максимума в га­ мильтоновой форме (теорема 1 § 2.4). Составим функ­ цию Понтрягина

Н = ри + qv — (Я0/2) (xv уи)

 

и сопряженную систему

 

р = X0v/2 = Х0у/2; — q = Я0ц/2 = Хах/2.

(10)

Если Хо = 0, то р и q постоянны и в силу принципа максимума скорости и и v были бы тоже постоянны, но

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