
Метод указ по МО
.pdf
-52 -
x(k) x(k 1) L 1 2(x(k 1) )13 b(b(x(k 1) )1 c(x(k 1) )2 ) p / 2 ,
2a(x |
(k 1) )3 |
c(b(x |
(k 1) ) c(x(k 1) ) |
|
|
, |
k 1, 2, , |
2 |
) q / 2 |
||||||
|
2 |
|
1 |
|
|
|
где | x(0) | r и
L192 max(1, a2 )r2 24 max(| b |4,| c |4 ) .
Аформула метода Ньютона принимает вид
x(k) x(k 1) 18a(x(k 1) ) |
2 (x(k 1) )2 3c2 (x(k 1) )2 |
3ab2 (x(k 1) )2 |
1 |
|||||||||
|
|
|
1 |
2 |
|
|
|
1 |
|
2 |
|
|
|
|
|
6a(x(k 1) ) |
2 c2 |
|
bc |
|
|
|
|||
|
|
|
bc |
2 |
|
|
|
|
|
|
|
|
|
|
|
6(x(k 1) )2 b2 |
|
|
|||||||
|
|
|
|
|
|
|
|
|
1 |
|
|
|
2(x(k 1) )3 |
b(b(x(k 1) ) c(x(k 1) ) |
2 |
) p / 2 |
|
|
|||||||
|
1 |
|
|
|
1 |
|
|
|
, |
k 1, 2, . |
|
|
2a(x(k 1) )3 c(b(x(k 1) ) c(x(k 1) ) |
|
) q / 2 |
|
|||||||||
|
2 |
|
|
|
||||||||
|
|
2 |
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
Задание № 7 |
|
|
|
|
|
|
|
Для заданной функции |
f , |
зависящей от двух переменных x1, x2 , |
найти |
|||||||||
точку минимума с точностью |
0,05, применяя |
|
|
|
а) метод спуска;
б) метод Ньютона.
Варианты
1)f (x1, x2 ) x14 2x24 (3x1 2x2 )2 2x1 4x2
2)f (x1, x2 ) x14 x24 (2x1 3x2 )2 4x1 2x2
3)f (x1, x2 ) x14 3x24 (x1 2x2 )2 2x1 2x2
4)f (x1, x2 ) x14 4x24 (3x1 2x2 )2 6x1 2x2
5)f (x1, x2 ) x14 5x24 (x1 x2 )2 4x1 2x2
6)f (x1, x2 ) x14 x24 (4x1 3x2 )2 2x1 6x2
- 53 -
7)f (x1, x2 ) x14 6x24 (5x1 x2 )2 4x1 2x2
8)f (x1, x2 ) x14 5x24 (4x1 3x2 )2 8x1 4x2
9)f (x1, x2 ) x14 4x24 (3x1 2x2 )2 6x1 2x2
10)f (x1, x2 ) x14 3x24 (2x1 3x2 )2 4x1 4x2
11)f (x1, x2 ) x14 2x24 (3x1 4x2 )2 6x1 2x2
12)f (x1, x2 ) x14 x24 (3x1 x2 )2 8x1 2x2
13)f (x1, x2 ) x14 2x24 (3x1 2x2 )2 2x1 4x2
14)f (x1, x2 ) x14 x24 (2x1 7x2 )2 4x1 2x2
15)f (x1, x2 ) x14 3x24 (x1 2x2 )2 2x1 2x2
16)f (x1, x2 ) x14 4x24 (3x1 8x2 )2 6x1 2x2
17)f (x1, x2 ) x14 5x24 (2x1 9x2 )2 4x1 2x2
18)f (x1, x2 ) x14 x24 (4x1 3x2 )2 2x1 6x2
19)f (x1, x2 ) x14 6x24 (5x1 x2 )2 4x1 2x2
20)f (x1, x2 ) x14 5x24 (4x1 3x2 )2 8x1 4x2
21)f (x1, x2 ) x14 4x24 (7x1 3x2 )2 6x1 2x2
22)f (x1, x2 ) x14 3x24 (2x1 3x2 )2 4x1 4x2
23)f (x1, x2 ) x14 2x24 (3x1 5x2 )2 6x1 2x2
24)f (x1, x2 ) x14 x24 (3x1 7x2 )2 8x1 2x2
25)f (x1, x2 ) x14 3x24 (5x1 2x2 )2 2x1 2x2
- 54 -
§8. Простейшая задача вариационного исчисления
Простейшая задача вариационного исчисления состоит в нахождении минимума или максимума функционала
f (x) F (t, x(t), x (t)) dt
на множестве гладких функций x(t), удовлетворяющих условиям
x( ) A, x( ) B .
Решение задачи, если оно существует, будет решением уравнения Эйлера
d F |
|
F |
|
||
|
|
|
|
||
dt x |
(t, x, x ) |
x |
(t, x, x ) 0. |
||
|
|
Поэтому уравнение Эйлера называют необходимым условием экстремума.
Всякое решение x(t) уравнения Эйлера, удовлетворяющее краевым условиям
x( ) A, x( ) B , называют экстремалью простейшей вариационной
задачи. Таким образом, решение простейшей вариационной задачи следует искать среди экстремалей.
Экстремаль x0 (t) будет решением простейшей вариационной задачи,
если удовлетворяет двум условиям – усиленному условию Лежандра и условию Якоби; их называют достаточными условиями экстремума.
1. Усиленное условие Лежандра: для задачи на минимум должно
выполняться неравенство
|
2F |
(t, x0 |
(t), x0 |
(t)) 0, |
t [ , ], |
|
x 2 |
||||||
|
|
|
|
а для задачи на максимум должно выполняться неравенство
|
2F |
(t, x0 |
(t), x0 |
(t)) 0, |
t [ , ], |
|
x 2 |
||||||
|
|
|
|
2. Условие Якоби: существует решение u(t) уравнения Якоби
- 55 -
|
2F |
d 2F |
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
2 (t, x0 (t), x0 (t)) |
|
|
|
|
(t, x0 (t), x0 |
(t)) |
|||
|
x |
dt x x |
u(t) |
|||||||
|
|
|
|
|
|
|
d |
|
2F |
|
|
|
|
|
|
|
|
|
|
(t, x0 (t), x0 |
(t))u (t) |
|
0 , |
dt |
|
x |
2 |
|
||||
|
|
|
|
|
|
|
такое, что u( ) 0 и u(t) 0 при t .
Пример 8.1. Найдем решение задачи
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
10 |
|
|
|
|
|
||||
|
|
|
|
f (x) |
|
( sin t |
|
|
|
|
|
|
dt min , |
||||||||||||||
|
|
|
|
|
2x(t) 4x (t) ) |
|
|
||||||||||||||||||||
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x(0) 2 , x( ) 4. |
|
|
|
|
|
|
|
||||||||||
Составим уравнение Эйлера и находим экстремаль: |
|
|
|
|
|
||||||||||||||||||||||
|
|
|
|
|
F |
|
|
|
|
|
|
|
|
|
|
|
|
|
10 |
|
|
||||||
|
|
|
|
|
(t, x, x ) ( sin t 2x 4x ) |
, |
|
|
|||||||||||||||||||
|
|
|
|
|
F |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
9 |
, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
x |
(t, x, x ) 40( sin t 2x |
4x ) |
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
F |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
9 |
, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
x |
|
(t, x, x ) 20( sin t 2x |
4x ) |
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 ( sin t |
|
|
|
|
|
9 |
|
|
|
|
|
|
9 |
0 - уравнение Эйлера, |
|||||||||||||
2x 4x ) |
( sin t 2x 4x ) |
|
|||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
9 |
|
|
|
1 |
|
|
|
|
||
|
|
|
|
|
|
|
|
|
( sin t |
2x 4x ) |
|
|
, |
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
( sin t |
|
|
|
|
9 |
|
|
2 |
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
2x 4x ) |
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
sin t 2x 4x C1e t /18 , |
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
x |
1 |
x |
1 |
sin t C1e t /18 , |
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
2 |
|
|
4 |
|
|
|
|
|
|
|
|
|
|
|
|
|||
x(t) |
9 |
C e |
t /18 C |
2 |
et / 2 |
|
1 |
(sin t 2cost) - общее решение уравнения |
|||||||||||||||||||
|
|
|
|||||||||||||||||||||||||
20 |
1 |
|
|
|
|
|
20 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Эйлера,
- 56 -
x(0) 2 , |
|
9 |
|
C C |
2 |
|
|
1 |
2, |
|
|
|
|
||||||||
20 |
|
10 |
|
|
|
|
|||||||||||||||
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|||||||
x( ) 4, |
|
9 |
|
C e /18 |
|
C |
e / 2 |
1 |
4, |
||||||||||||
20 |
|
|
|||||||||||||||||||
|
|
|
|
1 |
|
|
|
|
|
2 |
|
|
10 |
|
|||||||
|
|
|
|
|
|
|
|
C10 2,82 , |
|
C20 0,63, |
|||||||||||
x |
|
(t) |
9 |
C e t /18 C |
|
|
|
et / 2 |
|
1 |
(sin t 2cost) - экстремаль. |
||||||||||
|
|
|
|
|
20 |
||||||||||||||||
|
0 |
20 |
10 |
|
|
|
|
20 |
|
|
|
|
|
|
|
Проверим, выполняются ли усиленное условие Лежандра и условие Якоби:
2F |
(t, x0 (t), x0 (t)) 1440( sin t 2x 4x0 (t))8 |
1440C108 e 4t / 9 |
0 , |
|
x 2 |
||||
|
|
|
усиленное условие Лежандра на минимум выполняется, |
|
|
|||||||||||||||||
|
|
2F |
(t, x0 (t), x0 |
(t)) 360( sin t 2x 4x0 (t))8 |
360C108 e 4t / 9 , |
||||||||||||||
|
|
x2 |
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
2F |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
720C108 e 4t / 9 , |
|||
|
|
(t, x0 (t), x0 (t)) 720( sin t 2x 4x0 (t))8 |
|||||||||||||||||
|
x x |
||||||||||||||||||
|
|
|
|
|
4 |
|
|
1 |
u(t) 0 - уравнение Якоби, |
||||||||||
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
u (t) |
9 |
u (t) |
36 |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
u(t) D e t / 2 D et /18 - общее решение уравнения Якоби, |
|||||||||||||||||
|
|
|
|
1 |
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
u(0) 0, |
|
D1 D2 |
0, |
D2 D1, |
|||||||||||
|
|
|
|
|
|
|
u(t) D (e t / 2 |
et /18). |
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
Положим D1 1. Тогда |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
u(t) et /18 e t / 2 0 |
при t (0, ]. |
||||||||||||||
Следовательно, экстремаль x0 (t) является точкой минимума: |
|||||||||||||||||||
|
|
|
|
|
|
10 |
5t / 9 |
|
10 9 |
|
5 / 9 |
|
|||||||
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
fmin |
f (x0 ) |
|
|
C10 e |
|
|
|
|
dt C10 |
|
(1 e |
|
|
) 47137,95. |
|||
|
|
|
|
|
|
|
5 |
|
|
||||||||||
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|

- 57 -
Задание № 8
Найти решение простейшей вариационной задачи
|
|
m |
|
|
f (x) |
dt min ( max ) , |
|||
M ( a(t) bx(t) cx (t) ) |
|
|||
|
|
|
|
x( ) A, x( ) B ,
где числа M , m, , , b, c , A, B и функция a(t) задаются, m - целое и четное, m 2, , | c | 1.
|
|
|
|
Варианты |
1) |
M 1, m 2, 0, 1, b 1, c 2 , |
|||
|
A 0, B 1, |
a(t) 2t 1 |
||
2) |
M 1, m 4, 1, 2 , b 2, c 3, |
|||
|
A 1, B 0, |
a(t) 3t2 2 |
||
3) |
M 1, m 6, 0, , b 3, c 2 , |
|||
|
A 0, B 1, |
a(t) sin t |
||
4) |
M 1, m 4, 0, / 2, b 3, c 3, |
|||
|
A 1, B 0, |
a(t) cos t |
||
5) |
M 2 , m 2, 0, 3, b 1, c 2 , |
|||
|
A 0, B 1, |
a(t) t2 1 |
||
6) |
M 3, m 6, 1, 2 , b 1, c 2, |
|||
|
A 1, B 1, |
a(t) t 2 |
||
7) |
M 4, m 2, 2, 0, b 1, c 2, |
|||
|
A 0, B 1, |
a(t) |
|
|
|
t |
|||
8) |
M 6, m 6, 1, 3, b 1, c 3, |
|||
|
A 1, B 1, a(t) 4t2 1 |

- 58 -
9)M 5, m 4, 0, 1, b 1, c 3,
A 0, B 2, a(t) 3t2 1
10)M 3, m 8, 2, 2 , b 1, c 3,
A 0, B 1, |
a(t) et 1 |
11)M 2 , m 10, 0, 4 , b 1, c 2,
A 2, B 1, a(t) 2t4
12)M 4 , m 6, 0, 2 , b 1, c 2 ,
A 0, B 3, a(t) 1 sin 2t
13)M 5, m 2, 0, , b 1, c 2,
A 2, B 0, a(t) 4 3cos t
14)M 1, m 2, 0, 1, b 1, c 2 ,
|
A 0, B 1, |
a(t) 3t2 2 |
|
|
15) |
M 1, m 4, 1, 2 , |
b 2, c 3, |
||
|
A 1, B 0, |
a(t) 2t 1 |
|
|
16) |
M 2 , |
m 6, |
0, , |
b 3, c 2 , |
|
A 0, |
B 1, |
a(t) 1 cos t |
|
17)M 1, m 2, 0, / 2, b 3, c 3,
A 1, B 0, a(t) 1 2sin t
18)M 2, m 4, 0, 3, b 1, c 2 ,
|
A 0, B 1, |
a(t) t2 1 |
||
19) |
M 5, m 8, |
1, 2 , b 1, c 2, |
||
|
A 1, B 1, |
a(t) t 2 |
||
20) |
M 4 , m 4, 2, 0, b 1, c 2, |
|||
|
A 0, B 1, |
a(t) 2 |
|
|
|
t |
- 59 -
21) M 6, m 2, 1, 3, b 1, c 3,
A 1, B 1, a(t) 4t2 t
22)M 5, m 2, 0, 1, b 1, c 3,
A 0, B 2, a(t) 3t2 1
23)M 3, m 10, 2, 2 , b 1, c 3,
A 0, B 1, |
a(t) e t 1 |
24)M 2 , m 6, 0, 4 , b 1, c 2,
A 2, B 1, a(t) 5t2 3
25)M 4, m 8, 0, 2 , b 1, c 2 ,
A 0, B 3, a(t) 2 cos 3t
- 60 -
Библиографический список
1. Корнеенко, В. П. Методы оптимизации: учебник / В. П.
Корнеенко. – М.: Высшая школа, 2007. – 664 с.
2. Пантелеев, А.В. Методы оптимизации в примерах и задачах:
учебное пособие / А. В. Пантелеев, Т. А. Летова. – М.: Высшая школа,
2002. – 554 с.
3.Пантелеев, А. В. Вариационное исчисление в примерах и задачах. / А. В. Пантелеев. – М.: МАИ, 2000. – 228 с.
4.Васильев, Ф. П. Методы оптимизации / Ф. П. Васильев. – М.:
Наука, 2002. – 824 с.
- 61 -
Содержание
|
Стр. |
§1. Методы решения задач линейного программирования |
4 |
§2. Матричные игры |
14 |
§3. Транспортная задача |
19 |
§4. Метод выпуклого программирования |
29 |
§5. Метод динамического программирования |
33 |
§6. Максимальный поток в транспортной сети |
41 |
§7. Численные методы оптимизации |
47 |
§8. Простейшая задача вариационного исчисления |
53 |
Библиографический список |
59 |