На втором шаге выяснилось, что максимальное значение целевой функции больше значения на первом шаге (-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. Получилось, что наиболее точным является метод Гаусса-Зейделя в силу его высокой помехозащищённости в смысле выбора рабочей точки. Вторым по точности является метод Бокса-Уилсона, потому что у него тоже высокая помехозащищенность (составляющая градиента оценивается по четырём точкам), но точность нахождения экстремума сильно зависит от интервала варьирования и рабочего шага. Третий по точности – метод случайного поиска, четвёртый – метод дифференцирования. Следует так же отметить, что наиболее точным получился метод Гаусса-Зейделя случайно, из-за вида целевой функции. Так как эффективность каждого из рассмотренных методов зависит от конкретных условий, в частности от формы поверхности отклика: плавности (узкие гребни, овраги и т.д.). Точность нахождения экстремума зависит от правильного подбора метода поиска, наиболее подходящего для конкретных условий (форма поверхности отклика, уровень шумов и др.).
Уважаемый посетитель!
Чтобы распечатать файл, скачайте его (в формате Word).
Ссылка на скачивание - внизу страницы.