QalOGUGtk0
.pdf3r3 cos3 32 r cos 32 r3 cos3 32 r cos 92 r3 cos3 0 – верно. Ответ: u r, 32 r cos 12 r3 cos3 .
Литература
1.Тихонов А.Н., Самарский А.А. Уравнения математической физики. М., МГУ – Наука, 2004.
2.Араманович И.Г., Левин В.И. Уравнения математической физики. М.,
Наука, 1964.
3.Вуколов Э. А., Ефимов А.В., Земсков В.Н., Каракулин А.Ф., Лесин В.В., Поспелов А.С., Терещенко А.М. Сборник задач по математике для ВТУЗОВ. Методы оптимизации. Уравнения в частных производных. Интегральные уравнения. М., Наука, 1990.
4.Кузнецов Л.А. Сборник заданий по высшей математике. Типовые расчеты. Санкт-Петербург – Москва – Краснодар, Лань, 2005
191
Мартынов О.М.
Задачи для подготовки к междисциплинарному экзамену по дисциплине «Методы оптимизации» (специальность «Прикладная математика
и информатика»)
Вариационными задачами называются задачи о поиске экстремума функционалов.
Определение 1. Переменная величина v называется функционалом, зависящим от функции y x , если каждой функции y x из некоторого класса
функций M соответствует значение v , то есть имеет место соответствие: |
|||
функции y x соответствует число |
v . |
|
|
Обозначение: v v y x . |
|
|
|
|
|
|
|
Класс M функций (кривых), на которых определен функционал, называет- |
|||
ся его областью определения. |
|
|
|
Будем полагать, что функционал |
v y x , определен на элементах y x |
||
|
|
|
|
линейного нормированного пространства функций, в котором каждому элементу y x поставлено в соответствие действительное число y , называемое нормой элемента, при этом выполняются следующие условия:
1) y 0 и y =0 тогда и только тогда, когда y 0 (0 – нулевой элемент);
2) y y ;
3) y z y z
для любых элементов y, z , принадлежащих пространству, и любого действительного числа .
Далее, в основном, будем рассматривать пространства C0 и C1 . Пространство C 0 x0 , x1 состоит из непрерывных функций (кривых) y x , опре-
деленных на отрезке x0 , |
x1 . В пространстве C 0 x0 , x1 норма вводится |
|||||||||||||||||
следующим образом: |
|
y |
|
0 |
max |
|
|
y x |
|
. |
|
|
|
|||||
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
x x0 ;x1 |
|
|
|
|
|
|
|
|
|
|
|
Пусть y x C0 x0 , |
x1 |
|
|
|
|
|
|
|
|
|
|
|||||||
и 0 |
|
|
– произвольное число. |
|
||||||||||||||
Определение 2. - окрестностью нулевого порядка кривой |
y x называ- |
|||||||||||||||||
ется совокупность кривых y x C0 x0 , x1 , такая, что |
|
|||||||||||||||||
|
|
y y |
|
0 max |
|
|
y x y x |
|
. |
|
||||||||
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
x x0 ;x1 |
|
|
|
|
|
|
|
192
Пространство C1 x0 , x1 состоит |
из непрерывных функций (кривых) |
y x , определенных на отрезке x0 , |
x1 и имеющих на этом отрезке не- |
прерывную производную. В пространстве C1 x0 , x1 норма вводится следующим образом:
|
|
|
|
|
|
|
y |
|
1 |
max |
|
y x |
|
|
max |
|
y ' x |
|
. |
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
x x0 ;x1 |
|
|
|
|
|
|
x x0 ;x1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
Пусть y x C1 x0 , x1 |
и 0 |
|
– произвольное число. |
|
|
|
|||||||||||||||||||
Определение 3. |
- окрестностью первого порядка кривой |
y x называ- |
|||||||||||||||||||||||
ется совокупность кривых y x C1 x0 , |
x1 , такая, что |
|
|
|
|||||||||||||||||||||
|
|
y y |
|
1 |
max |
|
y x y x |
|
|
max |
|
y ' x y ' |
x |
|
. |
||||||||||
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
x x0 ;x1 |
|
|
|
|
|
|
|
x x0 ;x1 |
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Аналогично вводится норма в пространстве C m x0 , x1 функций, имеющих непрерывные производные до порядка m включительно, то есть
|
|
|
|
m |
|
|
y p x |
|
y |
|
|
|
m max |
|
|
. |
|
|
|
|
||||||
|
|
|
|
p 0 x x0 ;x1 |
|
|
|
|
|
|
|
|
|
|
|
Определение 4. Кривые y x , на которых сравниваются значения функционала, называются допустимыми кривыми или кривыми сравнения. Определение 5. Обозначим через y x допустимую кривую, на которой
функционал достигает экстремума, а через |
y x произвольную допусти- |
мую кривую. Разность y x y x y x |
называется вариацией кривой |
y x . |
|
Вариация y x есть функция x и принадлежит тому же функциональному пространству, что и функция y x . Используя вариацию y x , можно
представить любую допустимую кривую y x в виде
y x y x y x .
Используется так же и другая запись
y x y x y x ,
где y x – фиксированная функция, а – числовой параметр. Очевидно,
что при 0 |
справедливо равенство y x y x . |
|
|
|
|
|
|
|
||||||||
Определение |
6. |
Если приращение |
|
функционала |
v v y x y |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
v y x |
|
можно |
|
представить |
|
в |
|
|
|
|
виде |
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
v L y x , |
y |
y x , |
y max |
|
y |
|
, где L y x , |
y |
– линейный |
|||||||
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
по отношению к y функционал и |
y x , y 0 |
при max |
|
y |
|
0 , то |
||||||||||
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
193 |
|
|
|
|
|
|
|
|
|
|
|
Теорема |
1 (Необходимое условие экстремума). Если функционал |
|
v y x , |
имеющий вариацию, достигает максимума или минимума на |
|
|
|
|
кривой y y x , где y x |
– внутренняя точка области определения функ- |
||||||||||||||||
ционала, то при y x y x |
вариация функционала равна нулю: v 0 . |
||||||||||||||||
Исследуем на экстремум функционал |
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
x1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
v y |
|
x |
|
|
F |
|
x, y |
|
x |
|
, y ' |
|
x |
|
dx , |
(1) |
|
|
|
|
|
|
x0 |
|
|
|
|
|
|
|
|
|
|
|
|
причем граничные точки допустимых кривых закреплены: |
y x0 y0 и |
||||||||||||||||
y x1 y1 . Будем считать, |
что функция |
F x, |
y, y ' имеет непрерывные |
частные производные до второго порядка включительно. Эта задача называется простейшей задачей вариационного исчисления.
Теорема 2. Для того, чтобы функционал (1) достигал на функции y x (из
некоторого класса M функций) слабого экстремума, необходимо, чтобы эта функция удовлетворяла уравнению Эйлера (впервые им было опубликовано в 1744 году, Леонард Эйлер (1707 – 1783) – великий математик, большую часть своей жизни провел в России, по происхождению швейцарец):
Fy dxd Fy ' 0 .
Замечание 2. Уравнение Fy dxd Fy ' 0 в развернутом виде записывается
так: Fy Fxy ' Fyy ' y ' Fy ' y ' y '' 0 . Интегральные кривые уравнения Эйлера y y x, C1, C2 называются экстремалями. Только на экстремалях
может достигаться экстремум функционала (1). Для нахождения кривой, реализующей экстремум функционала (1), надо проинтегрировать уравнение Эйлера и определить обе произвольные постоянные, входящие в общее решение этого уравнения, из условий на границе y x0 y0 , y x1 y1 .
Только на удовлетворяющих этим условиям экстремалях может реализовываться экстремум функционала.
Замечание 3. Так как всякий сильный экстремум функционала является и слабым, то теорема 2 дает необходимое условие и сильного экстремума функционала (1).
Далее рассмотрим два обобщения простейшей задачи вариационного ис-
числения. Первым |
из |
них |
является задача на экстремум функционала |
|||||||
v y x , зависящего от производных высших порядков функции y x : |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x1 |
|
|
x, y, y ', ... , y n |
|
|
|
|
|
|
F |
dx , |
(2) |
|||||
|
v y |
x |
|
|
|
|||||
|
|
|
|
x0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
195 |
|
|
|
y |
|
C1x C2 C1 exdx C3x C4 |
u C1x C2 |
C1, |
|
du C1dx |
|
||
|
x |
|
|
x |
|
||||
|
|
|
dx, |
v e |
|
|
|||
|
|
|
dv e |
|
|
|
|
C1x C2 C1 ex C1exdx C3x C4 C1x C2 2C1 ex C3x C4 .
Для определения констант C1, C2, C3, C4 воспользуемся граничными условиями:
C2 2C1 C4 0 |
|
|
C1 1 |
|
||||
|
C1 C3 1 |
|
|
|
2 |
|
||
C2 |
|
|
C2 |
. |
||||
|
|
C |
e C C |
|
e |
|
0 |
|
C |
2 |
4 |
C |
|
||||
|
1 |
3 |
|
3 |
|
|
||
|
|
|
2e |
|
|
|
0 |
|
C2e C3 |
|
|
C4 |
|
Следовательно, искомая экстремаль имеет вид y xex .
Ответ: y xex .
Другим обобщением простейшей задачи вариационного исчисления является задача об экстремуме функционала, зависящего от нескольких функций:
|
x , y |
|
x , ... , y |
|
x |
|
x1 |
|
x, y |
x , ... , y |
|
x , y ' |
x , ... , y ' |
x dx , |
|||
v y |
2 |
n |
|
|
F |
n |
|||||||||||
1 |
|
|
|
|
|
|
|
1 |
|
|
1 |
n |
|
||||
|
|
|
|
|
|
|
|
x0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(5) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
где функция |
|
|
|
yn , |
' |
, |
|
' |
имеет непрерывные частные про- |
||||||||
F x, y1, ... , |
y1 |
... , yn |
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
изводные до второго порядка включительно по всем своим аргументам и yk x C1 x0 , x1 , k 1, ... , n.
Граничные условия в этой задаче имеют вид
yk x0 yk0 , yk x1 yk1 , k 1, ... , n . |
(6) |
Теорема 2 для данного случая обобщается следующим образом. |
|
Теорема 4. Для того, чтобы набор функций y1 x , ... , yn x C1 x0 , x1
доставлял слабый экстремум функционалу (5), необходимо, чтобы эти функции удовлетворяли системе дифференциальных уравнений Эйлера
Fyk dxd Fyk' 0, k 1, ... , n .
Пример 2. Найти функции y1 x и y2 x C1 x0 , x1 , на которых может достигаться экстремум функционала v y1, y2 при указанных граничных
197
|
|
|
|
1 |
y1y2 dx, |
|
|
|
|
|
|
условиях: |
v y1, y2 |
y1' y2' |
y1 0 y2 0 1, |
y1 |
1 e, |
y2 1 |
|||||
|
|
|
|
0 |
|
|
|
|
|
|
|
|
1 |
. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
e |
|
|
|
|
|
|
|
|
|
|
Решение. |
, y1' , y2' y1' |
|
|
|
|
|
|
|
Имеем |
||
F x, y1, y2 |
y2' y1 y2 |
Fy1 y2 |
, Fy2 |
y1, Fy ' |
y2' |
, |
Fy ' y1' |
. Сле- |
|||
|
|
|
|
|
|
|
1 |
|
|
2 |
|
довательно, система уравнений Эйлера имеет вид:
Fy |
|
d |
|
F |
|
|
0 |
|
|
|
|
|
|
|
|
|
|||
dx |
|
' |
|
|
|
'' |
0 |
|
|
|
|||||||||
|
1 |
|
|
y1 |
|
|
y2 |
|
y2 |
y1 |
C1e x C2ex , y2 |
C3e x C4ex . |
|||||||
|
|
|
|
d |
|
|
|
|
|
|
|
y '' |
0 |
||||||
|
|
|
|
|
|
|
|
|
|
y |
|
|
|
||||||
Fy2 |
|
|
|
Fy2' |
0 |
|
1 |
|
|
1 |
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
dx |
|
|
|
|
|
i 1,2,3,4 используем граничные условия: |
|||||||||
Для определения Ci |
|||||||||||||||||||
C1 C2 1 |
|
|
|
C 0 |
|
|
|
|
|||||||||||
C3 C4 1 |
|
|
|
|
1 |
|
|
|
|
|
|||||||||
|
|
|
C |
2 |
1 |
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
C e 1 |
|
C |
e e |
|
C 1 . |
|
|
|
|||||||||||
|
1 |
|
|
|
2 |
|
|
|
|
|
|
3 |
|
|
|
|
|
||
|
C e |
1 |
|
C |
|
e e |
1 |
|
C |
4 |
0 |
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
3 |
|
|
|
4 |
|
|
|
|
|
|
|
|
|
|
|
|
|||
Следовательно, искомые функции имеют вид: y ex , |
y e x . |
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
2 |
Ответ: y ex , y |
e x . |
|
|
|
|
|
|||||||||||||
|
|
|
|
|
1 |
|
|
|
2 |
|
|
|
|
|
|
|
|
|
Рассмотрим теперь вариационные задачи на условный экстремум. Задачи вариационного исчисления, в которых на искомые функции накладываются, кроме граничных условий, дополнительные ограничения, называются задачами на условный экстремум.
Рассмотрим следующую задачу об экстремуме функционала (5), зависящего от нескольких функций, с граничными условиями (6), при дополнительных ограничениях, заданных уравнениями связи
x, y , ... , y |
n |
, y ' |
, ... , y ' |
0, |
i 1, ... , m, |
m n . |
(7) |
|
i |
1 |
1 |
n |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Эта задача вариационного исчисления называется задачей Лагранжа (Жо- зеф-Луи Лагранж (1736 – 1813) – французский математик и механик). Введем функцию Лагранжа рассматриваемой задачи
L |
x, y , ... , y , y', |
... , y' |
, |
, ... , |
F |
x, y , ... |
, |
y |
n |
, y ' |
, ... , y ' |
|
||||
|
|
1 |
n |
1 |
n 1 |
m |
|
|
1 |
|
|
1 |
n |
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
m |
|
|
|
|
y1', ... , yn' |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
yn , |
, |
|
|
(8) |
|
||||
|
|
|
|
i x i x, y1, ... , |
|
|
|
|
||||||||
|
|
|
|
i 1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
где i x C1 x0 , x1 |
– произвольные функции (множители Лагранжа). |
|||||||||||||||
|
|
|
|
|
|
|
198 |
|
|
|
|
|
|
|
|
|
При решении задачи Лагранжа используется следующее необходимое условие экстремума функционала (5).
Теорема 5. Если функции y1 x , ... , yn x доставляют слабый экстремум функционалу (5) при условиях (6), (7), то существуют множители Лагран-
жа i x , i 1, ... , m, |
при которых эти функции удовлетворяют системе |
|||||||
уравнений Эйлера |
|
|
|
|
|
|
|
|
Lyk |
|
d |
L |
' 0, |
k 1, ... , n, |
(9) |
||
dx |
||||||||
|
|
|
yk |
|
|
|
||
|
|
|
|
|
x |
|
|
|
записанной для функционала u y1, |
... , yn 1 |
L dx . |
|
|||||
|
|
|
|
|
x0 |
|
|
С помощью теоремы 5 решение задачи об условном экстремуме функцио- |
|
нала v y1, y2 |
, ... , yn сводится к исследованию экстремума функционала |
u y1, ... , yn |
без дополнительных условий (7). |
При использовании теоремы 5 для решения задачи Лагранжа искомые
функции yk x , k 1, ... , n , |
и множители Лагранжа |
i x , i 1, ... , m, |
|||||||||
определяются из системы n m уравнений (9) и (7). |
|
|
|
||||||||
Пример 3. Найти функции |
y1 x и |
y2 x , на которых может достигаться |
|||||||||
экстремум функционала v y1, |
y2 в следующей задаче Лагранжа: |
||||||||||
1 |
|
|
|
|
|
|
|
|
|
|
|
v y1, y2 y1'2 y2' 2 dx, |
y1 |
0 2, |
|
y2 0 0, |
y1 1 2ch1, |
y2 1 2sh1, |
|||||
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
y ' y |
2 |
0 . |
|
|
|
|
|
||
|
|
1 |
|
|
|
|
|
|
|
|
|
Решение. Составим функцию Лагранжа |
|
|
|
|
|
||||||
L y1'2 y2' 2 y1' |
y2 , x . |
|
|
|
|||||||
Найдем ее производные: L |
0, L |
y ' |
2 y ' , |
L |
, |
L |
y ' |
2 y ' . |
|||
|
y1 |
|
|
|
1 |
|
y2 |
|
2 |
||
|
|
|
1 |
|
|
|
|
|
|
2 |
|
Составим систему из уравнений Эйлера и уравнения связи для определе- |
||||||||||
ния функций y1 |
x , |
y2 x и x : |
|
|
||||||
' 2 y '' |
0 |
2 y '' ' 0 |
2 y '' ' 0 |
|
||||||
|
|
1 |
|
|
1 |
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
0 y2''' y2' 0 |
2 y2'' |
0 2 y2'' 0 |
2 y2''' ' 0 2 y2''' 2 y2' |
||||||||
|
|
|
|
|
|
|
|
|
|
|
y ' y |
2 |
|
0 |
y ' |
y |
2 |
0 |
y '' y ' |
|
|
1 |
|
|
1 |
|
|
1 |
2 |
|
||
Пусть |
y ' z z '' z 0 . Характеристическое уравнение для этого обык- |
|||||||||
|
|
2 |
|
|
|
|
|
|
|
|
новенного |
дифференциального |
уравнения имеет вид r 2 1 0 |
||||||||
|
|
|
|
|
|
|
|
|
199 |
|