Другие методы построения моделей 2-го порядка

№6
Дисциплина

плана, впервые разработал Г. Данциг в 1947 г.
Этот метод

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

больше, чем размерность факторного пространства:

выбираются интервалы варьирования по каждому из значимых факторов, необходимые для

определения ребра симплекса. Область проведения поиска должна включать в себя не менее 8 симплексов. Например, имеем два параметра: температура в интервале от 100 до 200 градусов и давление от 1 атмосферы до 5 атмосфер. Для данного факторного пространства одной кодированной единице по температуре будет соответствовать около 10 градусов, а по давлению – около 0,4 атмосферы;
2. Задается значение ребра симплекса в единицах варьирования соответствующих переменных (проводится их кодирование). Ребро симплекса обычно принимается равным единице. Факторное пространство изменяется в кодированных значениях от 0 до минимум 8 единиц. По каждому параметру кодированной единице пространства может соответствовать свой размер;

из осей факторного пространства (а, б), или в промежутке между

ними (с). Два первых варианта упрощают расчеты координат точек, но ведут поиск вдоль одного из параметров;
4. Рассчитываются координаты точек симплекса сначала в кодированных значениях, потом переводятся в натуральные и реализуется исследования на объекте. Например, для указанной ранее температуры имеем ось в кодированных значениях от 0 до 10, что соответствует натуральным значениям от 100 до 200 и через пропорцию находим температуру в натуральных единицах;
5. Симплекс перемещается по факторному пространству за счет замены вершины симплекса с наихудшем значением по сравнению с остальными вершинами на вершину, являющуюся зеркальным отображением отброшенной вершины.
6. Расчет завершается после того, как симплекс сделает полный оборот вокруг наилучшей вершины.

а

б

с

для построения равносторонней фигуры, например, для 2-х параметрической модели строим

равносторонний треугольник с углами, равными 60 градусов. Первая точка задается произвольно, вторая точка лежит вдоль одной из координат факторного пространства на расстоянии 1 и третья точка вычисляется через тригонометрические функции Sin(60) и Cos(60), что дает (0,866; 0,5);
Нахождение зеркальной точки выполняется тоже геометрически через определение центра тяжести оставшейся фигурки (для n=2 – это отрезок линии), который вычисляется как среднее значение координат её точек. Для построения новой точки к координатам центра тяжести добавляются расстояния от него до худшей точки с учетом знаков (направлений движения симплекса).

расчету коэффициентов уравнения регрессии в виде полинома второго порядка со

взаимосвязями, мы имеем семь экспериментальных точек и уравнение с шестью неизвестными. Используя метод наименьших квадратов подбираем наилучшие коэффициенты для описания полученных данных.

Данные подход кажется более простым, но имеет основной недостаток – трудно выбрать правильный размер симплекса, небольшой размер приведет к медленному поиску области экстремума и малых размеров области для хорошего описания экстремума из-за точности используемой методики Большой размер симплекса может сразу же увести поиск за пределы факторного пространства.

смесевых систем не отличается от задачи факторного эксперимента второго порядка

(ЦКП). Однако здесь надо учитывать условие:



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

является правильным n-1 симплексом. Тогда факторное пространство можно представить в

n-1 мерном пространстве в виде, например треугольника для n=3.



Планирование на симплексах осуществляется равномерным разбросом экспериментальных точек. Получаются {n, m} -решетки, где n — число компонентов смеси; m — порядок полинома.

0
Трехпараметрические модели

+ β3X3
1
2
3
(1)

= β1X1 + β2X2 + β2X3+ β123Х1Х2Х3
1
2
3
4
(2)

решетки
Y = β1X1 + β2X2 + β3X3 + β1123Х12Х2Х3

+ β1223Х1Х22Х3 + β1233Х1Х2Х32

1

2

3

4

5

6

(3)

+ β3X3+ β12Х1Х2 + β13Х1Х3 + β23Х2Х3
1
2
3
4
5
6
(4)

β1X1 + β2X2 + β3X3 + β12Х1Х2 + β13Х1Х3 +

β23Х2Х3 + β123Х1Х2Х3

1

2

3

4

5

6

7

(5)

Y = β1X1 + β2X2 + β3X3 + β12Х1Х2 +

β13Х1Х3 + β23Х2Х3 + β1123Х12Х2Х3 + β1223Х1Х22Х3 + β1233Х1Х2Х32

1

2

3

4

5

6

7

8

9

(6)

β3X3 + β12Х1Х2 + β13Х1Х3 + β23Х2Х3 + β123Х1Х2Х3 +

γ12Х1Х2(Х1 - Х2) + γ13Х1Х3(Х1 – Х3) + γ23Х2Х3(Х2 – Х3)

1

2

3

4

5

6

7

8

9

10

(7)

β3X3 + β12Х1Х2 + β13Х1Х3 + β23Х2Х3 + + γ12Х1Х2(Х1

- Х2) + γ13Х1Х3(Х1 – Х3) + γ23Х2Х3(Х2 – Х3) + δ12X1X2(X1-X2)2 + δ13X1X3(X1-X3)2 + δ23X2X3(X2-X3)2 + β1123Х12Х2Х3 + β1223Х1Х22Х3 + β1233Х1Х2Х32

1

2

3

4

11

12

7

5

9

15

8

10

6

14

13

(8)

по разным путям от простых к сложным.
При отсутствии экстремумов

функции в факторном пространстве усложнение моделей идет по цепочке:

1

2

3

4

5

6

6

5

4

4

Модель 1
Точка для проверка
Модель 2
Точки для проверки
Модель 3
Точки для проверки
Модель 4

подбор модели можно проводить по следующей цепочке:
1
2
3
5
6
7
8
9
10
4
Модель 2
Точки для

проверка
Модель 4
Точки для проверки
Модель 5
Точки для проверки
Модель 6
Точки для проверки
Модель 7
Точки для проверки
Модель 8

11

12

13

14

16

15

17

18

19

20

21

22

опытов и неизвестных коэффициентов в уравнении, поэтому их расчеты доведены

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

а, например для второй модели остальные коэффициенты



По такой же схеме идут расчеты коэффициентов и в других моделях.

(метод наименьших квадратов и «Поиск решения»), которые были рассмотрены ранее.

В этом случае надо заботиться о количестве проведенных опытов и числе неизвестных коэффициентов (N > L) и выбирать уравнение, которое обеспечит свое существование только внутри факторного пространства. Последнее достигается введением дополнительного условия X1 + X2 + X3 =1 или произвести замену последнего фактора через предыдущие X3 = 1 – (X1 + X2)
Для решения данных задач и построения графиков имеется программа в среде MS Excel, которая позволяет выбрать модель, опередить таблицу экспериментальных данных, провести эксперимент и выполнить расчет с построением графика.
Порядок выполнения расчета показан далее

результаты
Получаем ответ
Настраиваем график