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

Метод указ по МО

.pdf
Скачиваний:
7
Добавлен:
08.06.2015
Размер:
575.04 Кб
Скачать

-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