Определение оптимальных параметров и экстремумов целевой функции методами Бокса-Уилсона и случайного поиска, страница 3

На втором шаге выяснилось, что максимальное значение целевой функции больше значения на первом шаге (-2576,85 и -2577,42  соответственно), а величины коэффициентов регрессии приблизились к нулю. Поэтому необходимо прекратить дальнейшее движение к экстремуму. И так, глобальный экстремум (максимум) для данной целевой функции равен y=-2576,85, при координатах x11= -606,759, x21= -498,5.

Метод случайного поиска.

Исходные данные:

1. начальная точка  

2. рабочий шаг .

Таблица случайных чисел имеет вид:

Таблица 3.

Таблица случайных чисел

9

23

29

20

19

6

20

9

23

26

5

23

25

9

10

4

12

8

4

22

24

1

21

8

2

26

1

18

1

12

Порядок выполнения работы.

Выбирается начальная точка у0.

Проводиться опыт в пробной точке, отстоящей от нулевой на расстояние r. На окружности точка, в которой проводится пробный опыт, определяется случайным образом.

Пусть r=10.

Берется таблица случайных чисел, выбираются точки со значением, меньшим 10: 28; 17; 8; 22; 6.

Число 8 принимаем в качестве размера вектора x, его направление определяем в зависимости от числа, которое стоит перед числом 8. если предыдущее число четное, то направление вектора положительное, если нечетное – отрицательное.

; .

Знак вектора x2 определяется от последующего числа: если последующее число четное, то знак x2 принимается отрицательным, если нечетное – положительным.

Затем определяется значение функции в точке у1 и сравнивается со значением у0. если у0>у1, то движение будет осуществляться в направлении, противоположном у1 на рабочий шаг r, который больше, чем r (r>r).

В у2 снова находится значение случайной точки у3, r остается тем же.

Значения целевой функции у от значений векторов  и  представлены в таблице 4.

Таблице 4.

Значения целевой функции при различных величинах векторов

Величины векторов

Координаты точки

Значение функции

x10

x20

у

-6

6

-494

-490,305

-4741,76

9

-9

-503

-479,864

-3929,31

5

-5

-508

-469,617

-3413,6

-9

9

-499

-460,078

-3341,58

4

4

-503

-450,28

-2996,12

8

8

-511

-440,688

-2713,95

4

4

-515

-430,89

-2599,24

1

-1

-516

-420,841

-2580,73

-8

8

-508

-411,249

-2581,27

2

2

-510

-401,349

-2660,7

1

1

-511

-391,4

-2806,3

1

1

-512

-381,45

-3021,64

По таблице определили максимум целевой функции у=-2580,73 при значениях координат x1=-516, x2=-420,841.

Вывод:   В ходе выполнения работы были найдены экстремумы целевой функции (максимумы)  методом Бокса–Уилсона и методом случайного поиска. Для первого метода y=-2576,85, при координатах x11=-606,759, x21=-498,5; а для второго          у=-2580,73 при значениях координат x1=-516, x2=-420,841. В шестой работе были найдены экстремумы этой же функции методом дифференцирования y=-2595, при координатах x11=-560, x21=-460 и методом Гаусса-Зейделя    y=-2575, при координатах x11=-610, x21=-500. Получилось, что наиболее точным является метод Гаусса-Зейделя в силу его высокой помехозащищённости в смысле выбора рабочей точки. Вторым по точности является метод Бокса-Уилсона, потому что у него тоже высокая помехозащищенность (составляющая градиента оценивается по четырём точкам), но точность нахождения экстремума сильно зависит от интервала варьирования и рабочего шага. Третий по точности – метод случайного поиска, четвёртый – метод дифференцирования. Следует так же отметить, что наиболее точным получился метод Гаусса-Зейделя случайно, из-за вида целевой функции. Так как эффективность каждого из рассмотренных методов зависит от конкретных условий, в частности от формы поверхности отклика: плавности (узкие гребни, овраги и т.д.). Точность нахождения экстремума зависит от правильного подбора метода поиска, наиболее подходящего для конкретных условий (форма поверхности отклика, уровень шумов и др.).