3.7. Метод симплекс-планирования

Симплекс - это выпуклая геометрическая фигура с числом вершин на единицу больше, чем число факторов. На плоскости симплекс имеет вид треугольника, в трёхмерном пространстве - это тетраэдр и в k-мерном пространстве фигура, имеющая k+1 вершину (фиг.3.5).

Стратегия симплекс-планирования определяется правилами, знание которых позволяет правильно записать исходный симплекс, найти координаты следующих точек при движении и перейти от симплекс-плана к планированию эксперимента, выполняемого для описания почти стационарной области.

Центр симплексов помещают в начало координат, а расстояние между его вершинами (сторону симплекса) принимают равным единице. Ориентируют симплекс в факторном пространстве таким образом, чтобы вершина лежала на оси , а остальные вершины располагались симметрично относительно координатных осей. В этом случае координаты вершин симплексов можно определить, зная радиусы вписанной ( ) и описанной ( ) гипосфер.

При симплекс-планировании движение в область оптимума в факторном пространстве осуществляется перемещением симплекса путём последовательного отбрасывания его вершин, имеющих минимальное значение параметра оптимизации, и поворота симплекса вокруг грани, противоположной отброшенной вершине. Например, если в точках 1,2,3 поставлены эксперименты (фиг.3.6) и наихудшее значение параметра оптимизации оказалось в точке 1, то треугольник поворачивают вокруг стороны 2-3, в результате чего получают точку 4, в которой ставят очередной эксперимент. Из точек 2,3 и 4 выбирают наихудшую и снова делают поворот вокруг стороны, противолежащей наихудшей точке. Процедуру повторяют до тех пор, пока симплекс не начнёт вращаться вокруг точки оптимума.

Если поверхность отклика близка к линейной, то движение из центра каждого симплекса через грань, противоположную вершине с минимальным значением параметра оптимизации, совпадает по направлению с движением по градиенту, рассчитанному по результатам опытов, выполненных в вершинах симплексов.

Первая задача симплекс-планирования состоит в построении матрицы исходного симплекса, в которую записывают координаты экспериментальных точек. Матрицу исходного симплекса строят путём последовательного перехода от одного симплекса к другому, увеличивая размерность симплексов на единицу. Координаты вершин симплексов можно определить из следующей матрицы:

.

Величины радиусов для вписанной ( ) и описанной ( ) сфер - мерного симплекса определяются формулами:

.

Результаты определения координат вершин симплексов для k=7 приведены в табл.3.14.

Матрицей пользоваться удобнее, если её привести к такому виду, чтобы она соответствовала симплексу с радиусом сферы, равным единице. Для этого все элементы числовой матрицы в табл.3.14 необходимо разделить на величину максимального элемента для симплекса определённой размерности.

Для симплекса, характеризующего два фактора , все элементы матрицы следует разделить на 0,578; при k=3 на 0,612, при k=4 на 0,632 и т.д. Деление на величину максимального элемента приводит к тому, что суммы квадратов элементов для всех столбцов становятся равными. Это условие, а также ортогональность матрицы и симметричное расположение независимых переменных относительно центра эксперимента обеспечивают оптимальность симплекс-планирования. Преобразованные оптимальные матрицы для двух и четырёх факторов приведены в табл.3.15 и 3.16.

Таблица 3.14

Численные значения координат вершин симплекса для k=7

№ опыта								отклик
1	-0,500	-0,289	-0,204	-0,158	-0,129	-0,109	-0,0945
2	0,500	-0,289	-0,204	-0,158	-0,129	-0,109	-0,0945
3	0	0,578	-0,204	-0,158	-0,129	-0,109	-0,0945
4	0	0	0,612	-0,158	-0,129	-0,109	-0,0945
5	0	0	0	0,632	-0,129	-0,109	-0,0945
6	0	0	0	0	0,645	-0,109	-0,0945
7	0	0	0	0	0	0,654	-0,0945
8	0	0	0	0	0	0	0,662

Таблица 3.15

Матрица симплекс-планирования для двух факторов

№ опытов в исходном симплексе	Факторы		Отклик
1	-0,865	-0,500
2	+0,865	-0,500
3	0	1

Таблица 3.16

Матрица симплекс-планирования для четырёх факторов

№ опытов в исходном симплексе					отклик
1	-0,790	-0,457	-0,323	-0,250
2	0,790	-0,457	-0,323	-0,250
3	0	0,914	-0,323	-0,250
4	0	0	0,969	-0,250
5	0	0	0	1,0

После построения матрицы исходного симплекса решается вторая задача симплекс-планирования, а именно, движение в область оптимума.

Перемещение симплекса в факторном пространстве осуществляется путём его поворота вокруг грани, противоположной вершине с минимальным значением параметра оптимизации. Координаты новой вершины следующие:

,

где - координаты новой точки; - среднее из координат всех точек симплекса, кроме "плохой"; - координата "плохой" точки. После построения нового симплекса и проведения опытов в найденной точке вновь решается вопрос о том, какая из точек может быть исключена. Последовательное отображение различных величин симплекса приводит к тому, что центр симплекса перемещается по некоторой ломанной линии. Траектория его движения зависит от порядка, в котором отображаются вершины симплекса. Центр симплекса движется с максимальной скоростью, если порядок отображения вершин соответствует периодически повторяющемуся ряду всех вершин (фиг.3.7). Незакономерная последовательность отображения вершин (некоторые вершины отображаются чаще остальных) замедляет движение центра симплекса.

Фиг.3.7. Движение по поверхности отклика при симплекс-планировании:

a) без ошибок; b) с колебаниями; с) при зацикливании.

При перемещении симплекса в факторном пространстве может возникнуть ситуация, при которой в новой точке величина параметра оптимизации будет минимальной (как в предшествующей точке). При отбрасывании этой точки происходит возврат к исходному симплексу, а далее симплексы начинают колебаться относительно одной грани. такое возможно, если вершины симплекса располагаются на гребне поверхности отклика. В этом случае рекомендуется уменьшить длину грани симплекса и эксперименты повторить.

При движении симплекса в факторном пространстве может наступить его вращение вокруг одной из вершин - зацикливание симплекса (фиг.3.7). Это свидетельствует о том, что в результате движения достигнута область оптимума. Следует иметь в виду, что причиной зацикливания симплекса может быть и ошибка эксперимента.

Симплекс-планирование является одним из эффективных методов поиска оптимальной области. Метод обладает многими достоинствами: прост для вычислений; не предъявляется особых требований к точности эксперимента (наихудшую ситуацию можно определить приближённо); обеспечивается высокая эффективность при минимальном числе экспериментов. Симплекс-планирование относится к методам, которые можно успешно использовать в промышленных условиях.