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

Материалы по дисциплине «Методы оптимизации» (90

..pdf
Скачиваний:
14
Добавлен:
15.11.2022
Размер:
505.45 Кб
Скачать

Так как кривая должна проходить через начало координат, то нужно положить C2 = 0 (при u = 0 ). Видно, что искомая кривая

– брахистохрона есть циклоида. Постоянная C1 должна быть най-

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

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

Постановка задачи. Пусть в прозрачной среде с переменной оптической плотностью даны две точки A и B . Требуется определить траекторию луча света, идущего от A к B .

Решение задачи основано на так называемом принципе Фер-

ма: из всех кривых, соединяющих A и B ,

траектория луча света

есть линия, по которой свет проходит из A в B в кратчайшее вре-

мя.

 

 

 

 

 

Ограничиваясь плоским случаем, примем за плоскость рас-

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

света

плоскость

XOY

и

пусть

A = (x0 , y0 ), B

= (x1 , y1 ), а

x0 x x1

некоторая кривая, со-

единяющая эти точки.

Обозначим через v (x , y) скорость света в точке (x , y) . Пусть

y

=

y (x) есть уравнение кривой, по которой движется точка из A

в

 

 

 

 

 

 

 

 

dS

 

1 +

y2

B . Скорость движения точки равна

v

=

 

=

 

 

dx , где

dt

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

– элемент времени. Отсюда dt =

 

1

+ y2

dx . Интегрируя, по-

 

v (x , y)

лучим время T , потребное для покрытия луча света от точки A до

точки B по кривой y = y (x):

 

 

 

 

 

 

 

 

 

 

dy 2

 

 

 

 

 

 

 

 

 

 

 

x1

1 +

 

 

 

 

 

 

 

 

 

 

 

T

=

dx

dx .

 

 

 

 

 

 

 

 

v (x ,

y)

 

 

 

 

 

 

 

 

 

x0

 

 

 

 

 

 

 

 

 

Согласно принципу Ферма задача определения траектории луча света сводится к определению линии, для которой интеграл T (функционал T ) принимает наименьшее значение: T min .

Рассмотрим частный случай,

когда v (x , y) = y . Тогда инте-

грант функционала зависит лишь от

y и

y, поэтому уравнение

Эйлера есть уравнение энергии L

yLy

= C .

Имеем

 

 

 

 

 

 

 

y

y

 

1 + y2

=

C .

 

 

y 1 +

y2

y

 

 

 

 

 

 

 

 

41

Экстремалями функционала являются окружности в полу-

 

 

плоскости y

>

0 с центрами на оси OX .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В самом деле

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

2

(1

+

y

2

)

=

C ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

1

+

y2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C y 1 +

y2

 

= − 1 ,

 

C 2 y2 (1 + y2 ) = 1 ,

 

y2

 

=

1 C 2 y2

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C 2 y2

 

 

 

 

 

 

 

 

dy

=

1

C 2 y2

 

 

y dy

 

 

 

 

 

=

dx

 

1

 

 

d (1 C 2 y2 )

=

dx ,

 

 

 

 

 

 

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

dx

 

C y

 

 

 

 

1

C

2

y

2

 

C

2C

1 C

2

y

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

1

C 2 y2

 

= x +

 

C1 ,

 

 

1

 

 

(1

C 2 y2 )

= (x + C1 )2 .

 

 

 

 

 

 

 

 

C

 

 

 

C 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Окончательно имеем следующий вид семейства экстремалей:

(x + C1 )2

+ y2 =

1

.

 

Задача 7.

 

C 2

Известно, что существует дважды непрерывно-

дифференцируемая вектор-функция (x1 (t) , x2 (t)), которая удовле-

творяет

 

 

 

 

начальным

условиям

x1 (1) = 1 , x1 (2) =

2 , x2 (1) = 0 , x2 (2)

= 1 и доставляет минимум

функционалу

 

 

 

 

 

2

2

+

2

)dt .

 

(x1

+ x2

 

x2

 

Найти эту функцию.

Решение. Известно, что если вектор-функция (x1 (t) , x2 (t)) ми-

нимизирует функционал, то она является решением следующей системы:

L

 

 

 

d

 

 

 

L

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

0

x

 

 

 

x

 

 

 

dt

 

 

,

1

 

 

 

 

 

 

1

 

 

 

L

 

 

 

d

 

 

 

L

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

0

x

 

 

 

 

 

 

x

 

2

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

где L – интегрант функционала (подынтегральная функция). В данном случае система имеет вид:

d

 

 

 

 

 

 

.

dt (2x1 )

=

0

2x2

d

(2x2 )

= 0

dt

 

 

 

 

42

Ее решение имеет вид

′′

= 0 , x2

 

′′

= 0},

 

{x1

=

 

C1 t + C2 , x2

= C3 e

t

+ C4 e

t

}.

{x1

 

x2

 

 

 

 

Из

краевых

 

 

 

условий

 

 

 

находим

 

константы

C1 = 1 ,

C2

=

0 ,

 

C3

=

 

 

1

 

,

C4

=

 

 

e2

. Тогда решение задачи

 

 

e2

1

1

e2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

et

 

 

 

 

 

 

e2

et

 

 

 

 

 

 

 

 

 

 

имеет вид

t ,

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e

1

 

 

 

e

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Задача 8. Найти экстремали функционала

F (y), зависящего

от старших производных. Предполагается, что искомые функции являются дважды непрерывно-дифференцируемыми на отрезке [x0 , x1 ] и удовлетворяют некоторым граничным условиям. Имеем

вариационную задачу:

 

F (y)

 

 

=

 

1 (360 x2

y

 

 

 

 

y′′2 )dx

 

min ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y(0) = 0 , y(1) = 0 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y (0) = 1 , y

(1) =

 

 

 

2 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение. Согласно теореме о необходимых условиях мини-

 

 

мума данной вариационной задачи, решения задачи должны

 

 

удовлетворять

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

уравнению

 

 

 

 

Эйлера – Пуассона

 

 

L

 

 

d

 

 

L

 

 

 

 

 

 

d

2

L

 

 

 

 

Здесь интегрант

 

 

 

 

L (x , y ,

y′′).

 

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

= 0 .

L

=

 

 

 

y

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

dx2 y′′

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Уравнение Эйлера – Пуассона имеет вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

360 x

2

+

 

d 2

 

 

 

 

 

′′

)

 

 

=

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2 y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

или y(4)

=

 

180 x2 .

 

 

Последовательно интегрируя левую и пра-

 

 

вую

 

 

 

 

 

 

 

 

 

 

 

 

 

части

 

 

 

 

 

 

 

 

 

уравнения,

 

 

 

 

 

 

 

получим

 

 

y(3)

= 60 x3

+ C

 

,

y(2)

 

 

=

 

15 x4

 

+ C x + C

2

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y′ = 3x5

 

 

+

 

C

 

x2

 

 

+ C

 

x + C

 

,

y =

x6

 

+

 

C1 x3

 

+

C

2

x2

 

+ C

 

x +

C

 

 

 

 

2

 

 

 

 

6

 

 

6

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

2

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

3

 

 

4 .

Константы C1 , C2 , C3 , C4

 

 

находим из начальных условий. Оконча-

 

 

тельно экстремаль функционала F (y) имеет вид

 

 

 

 

 

 

 

 

 

 

 

y =

 

 

 

x6

 

+

 

3x3

 

 

3 x

2

+ x .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

43

13. Изопериметрическая задача – задача Дидоны

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

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

 

x

 

 

F0 (y)

= 1

L0

(x , y , y) dx min ,

 

x0

 

 

 

x

 

 

Fi (y)

= 1

Li

(x , y , y) dx = Ci (i = 1 m) ,

 

x0

 

 

y (x0 ) = y0 ,

y (x1 ) = y1 .

Для изопериметрической задачи естественным образом определяются понятия решения и локального решения. Считается также, что функции Li (x , y , y), i = 1 m и их частные производные

по y и yнепрерывны по совокупности переменных. Для реше-

ния изопериметрической задачи используется правило множителей Лагранжа, согласно которому решение исходной задачи должно удовлетворять уравнению Эйлера для интегранта

 

 

 

m

Li ,

соответствующего

функционалу

L =

L0

+

λi

 

 

 

i = 1

 

 

 

F (y)

 

 

(y) +

m

 

 

=

F0

λi Fi (y), где λi , i = 1 m – константы, не все од-

i = 1

новременно равные нулю, называемые множителями Лагранжа. Начиная решение задачи Дидоны, отметим, прежде всего, что искомая кривая должна быть выпуклой. В противном случае существовала бы прямая, имеющая с рассматриваемой областью S внутри этой кривой по крайней мере два не имеющих общих точек замкнутых отрезка AB и DE . Отразив зеркально BCD в BCD ,

44

получим кривую той же длины l , но ограничиваемая площадь будет больше (рис. 1).

Рис. 1 Рис. 2

Примем за ось OX прямую AB , Тогда площадь, ограниченная кривой y = y (x), расположенной над осью OX , выразится инте-

гралом S (y)

b

 

 

= y dx , где a

и b

– абсциссы точек A и B (рис. 2).

 

a

 

 

Поставленная изопериметрическая задача будет иметь вид:

F0 (y)

b

 

= −

y dx min ,

 

a

 

F1 (y)

b

 

=

1 + y2 dx = l ,

 

a

 

y (a)

= 0 , y (b) = 0 .

Следуя правилу множителей Лагранжа, составляем функцио-

нал

F (y) = b

[y + λ 1 + y2 ]dx .

a

 

Интегрант этого функционала не зависит от x , а лишь от y и y, поэтому имеет место частный случай уравнения Эйлера – ин-

теграл энергии yLy

L =

C :

λ y

 

y

 

+ y

λ

1 + y2 = C .

 

 

 

 

+

 

y2

1

 

 

 

 

Приведя в левой части уравнения к общему знаменателю, получим:

45

 

 

 

λ

y

2

 

+ y 1

+

y

2

 

λ (1 +

y

2

)

= C

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 + y2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

или

 

 

 

 

 

 

λ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

=

 

C +

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

+

y2

 

 

 

 

 

 

 

 

 

 

 

Положим

dϕ

 

 

y

=

 

tg ϕ .

Тогда

y = C

+

λ cos ϕ

;

y

=

 

λ sin ϕ

=

tg ϕ ,

 

 

 

 

 

 

 

 

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

=

λ cos ϕ ,

dx

=

λ cos ϕ dϕ , x =

λ sin ϕ

+

C1 . Следова-

 

 

 

dϕ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

x

=

 

λ sin ϕ

+

C1 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

=

 

λ cos ϕ

+

C .

 

 

 

 

 

 

 

 

 

 

 

 

 

В явной форме семейство экстремалей исследуемой задачи

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(x C1 )2

+ (y C)2

= λ2 .

в функционал F0 (y), получим ис-

 

 

Подставляя значения y и dx

комое значение минимума.

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

46

Оглавление

1.

Введение.........................................................................................

3

2.

Цели и задачи дисциплины........................................................

8

3.

Требования к уровню освоения содержания курса ...............

9

4.Сведения из Государственного образовательного стандарта

ирабочей программы дисциплины «Методы

 

 

оптимизации»........................................................................

9

5.

Содержание разделов дисциплины.........................................

12

6.

Список рекомендуемой литературы.......................................

15

7.

Вопросы к экзамену по курсу «Методы оптимизации» .....

16

8.

Контрольная работа № 1 по теме «Оптимизация в

 

 

 

конечномерных пространствах» .....................................

18

9.

Контрольная работа № 2 по теме «Вариационное

 

 

 

исчисление» .........................................................................

19

10.

Дисциплина «Методы оптимизации» в итоговой

 

 

 

государственной аттестации............................................

20

11.

Варианты расчетно-графической работы...........................

23

12.

Разбор решений некоторых вариационных задач .............

35

13.

Изопериметрическая задача – задача Дидоны ..................

44

47

Учебное издание

Майорова Наталия Львовна

Материалы по дисциплине «Методы оптимизации»

Методические указания

Редактор, корректор И.В. Бунакова Компьютерная верстка Е.Л. Шелеховой

Подписано в печать

2009 г. Формат 60х84/16.

Бумага тип. Усл. печ. л. . Уч.-изд. л. 1,8.

Тираж 100 экз. Заказ

.

Оригинал-макет подготовлен в редакционно-издательском отделе ЯрГУ.

Отпечатано на ризографе.

Ярославский государственный университет. 150000 Ярославль, ул. Советская, 14.

48

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]