Методика регрессионного анализа основана на предположении, что учтены все или, по крайней мере, все существенные факторы, иначе полученная математическая модель окажется неадекватной в изучаемом диапазоне изменения переменных. Привлечение всего множества переменных к составлению математического описания может потребовать непомерного объема экспериментальной и вычислительной работы, что зачастую невыполнимо в силу технологических, экономических и прочих ограничений. Возникает необходимость предварительного отсеивания несущественных переменных и выделения тех входных воздействий
Если число всех возможных факторов, влияющих на объект, не превышает 6 – 7, то для предварительного изучения объекта можно применить методы дробного или полного факторного эксперимента. Однако при большом числе рассматриваемых факторов методы ПФЭ и даже ДФЭ, предназначенные для тщательного изучения поверхности отклика, оказываются слишком громоздкими и трудоемкими для постановки отсеивающих опытов. В случае изучения более 8 – 10 факторов, если эксперименты недороги и если заведомо известно, что лишь немногие переменные являются существенными, следует применять метод случайного баланса (МСБ).
Важнейшей теоретической предпосылкой МСБ является априорное знание того, что из всей совокупности рассматриваемых переменных
Под "шумовым полем" обычно понимают случайные помехи
Основная идея метода заключается в том, что вместо дробных реплик, которые представляют собой систематические ортогональные выборки из ПФЭ, берутся случайные выборки. Тогда вектор-столбцы матрицы планирования можно считать не коррелированными (не связанными) или слабо коррелированными друг с другом. Совместные оценки оказываются смешанными случайным образом. Появляется возможность с высокой надежностью выделить и независимо оценить все доминирующие переменные.
Из сказанного следует, что МСБ обладает меньшей чувствительностью, чем ПФЭ и ДФЭ. Под чувствительностью метода обычно понимается способность выделять коэффициенты
2.1 Составление матрицы планирования
Построение матрицы планирования для проведения отсеивающих опытов выполняют на основе предпосылки, что исследуемые факторы должны быть смешаны случайным образом. Все линейные факторы zi (i = l, 2, ..., n) разбивают на группы, при этом стремятся заведомо взаимодействующие факторы включить в одну группу. Если же нет априорных сведений о физике процесса, то разбивку факторов по группам можно производить формально, с использованием таблицы (или программы) случайных чисел. Затем для каждой группы составляют МП на основе ПФЭ или ДФЭ.
В задании требуется исследовать 8 факторов z1, z2, z3, z4, z5, z6, z7, z8 и их взаимодействия и с помощью МСБ выделить самые существенные. Разобьем факторы на две группы: 1 – z1, z2, z3, z4, 2 - z5, z6, z7, z8.
Общая МП для МСБ составляется в виде таблицы 2.3 путем построчной стыковки таблиц 2.1 и 2.2 после рандомизации их строк с помощью таблицы случайных чисел или компьютерной программы для реализации процесса рандомизации. Общая МП в данном примере может иметь вид таблицы 2.3.
Таблица 2.1
g | z1 | z2 | z3 | z4 |
1 | -1 | -1 | -1 | -1 |
2 | 1 | -1 | -1 | -1 |
3 | -1 | 1 | -1 | -1 |
4 | 1 | 1 | -1 | -1 |
5 | -1 | -1 | 1 | -1 |
6 | 1 | -1 | 1 | -1 |
7 | -1 | 1 | 1 | -1 |
8 | 1 | 1 | 1 | -1 |
9 | -1 | -1 | -1 | 1 |
10 | 1 | -1 | -1 | 1 |
11 | -1 | 1 | -1 | 1 |
12 | 1 | 1 | -1 | 1 |
13 | -1 | -1 | 1 | 1 |
14 | 1 | -1 | 1 | 1 |
15 | -1 | 1 | 1 | 1 |
16 | 1 | 1 | 1 | 1 |
Таблица 2.2
g | z5 | z6 | z7 | z8 |
1 | -1 | -1 | -1 | -1 |
2 | 1 | -1 | -1 | -1 |
3 | -1 | 1 | -1 | -1 |
4 | 1 | 1 | -1 | -1 |
5 | -1 | -1 | 1 | -1 |
6 | 1 | -1 | 1 | -1 |
7 | -1 | 1 | 1 | -1 |
8 | 1 | 1 | 1 | -1 |
9 | -1 | -1 | -1 | 1 |
10 | 1 | -1 | -1 | 1 |
11 | -1 | 1 | -1 | 1 |
12 | 1 | 1 | -1 | 1 |
13 | -1 | -1 | 1 | 1 |
14 | 1 | -1 | 1 | 1 |
15 | -1 | 1 | 1 | 1 |
16 | 1 | 1 | 1 | 1 |
Таблица 2.3
g | g' | z1 | z2 | z3 | z4 | g'' | z5 | z6 | z7 | z8 |
1 | 4 | 1 | 1 | -1 | -1 | 2 | 1 | -1 | -1 | -1 |
2 | 3 | -1 | 1 | -1 | -1 | 16 | 1 | 1 | 1 | 1 |
3 | 2 | 1 | -1 | -1 | -1 | 5 | -1 | -1 | 1 | -1 |
4 | 5 | -1 | -1 | 1 | -1 | 4 | 1 | 1 | -1 | -1 |
5 | 9 | -1 | -1 | -1 | 1 | 11 | -1 | 1 | -1 | 1 |
6 | 11 | -1 | 1 | -1 | 1 | 3 | -1 | 1 | -1 | -1 |
7 | 15 | -1 | 1 | 1 | 1 | 10 | 1 | -1 | -1 | 1 |
8 | 7 | -1 | 1 | 1 | -1 | 6 | 1 | -1 | 1 | -1 |
9 | 16 | 1 | 1 | 1 | 1 | 9 | -1 | -1 | -1 | 1 |
10 | 8 | 1 | 1 | 1 | -1 | 1 | -1 | -1 | -1 | -1 |
11 | 1 | -1 | -1 | -1 | -1 | 14 | 1 | -1 | 1 | 1 |
12 | 13 | -1 | -1 | 1 | 1 | 7 | -1 | 1 | 1 | -1 |
13 | 6 | 1 | -1 | 1 | -1 | 13 | -1 | -1 | 1 | 1 |
14 | 10 | 1 | -1 | -1 | 1 | 12 | 1 | 1 | -1 | 1 |
15 | 14 | 1 | -1 | 1 | 1 | 8 | 1 | 1 | 1 | -1 |
16 | 12 | 1 | 1 | -1 | 1 | 15 | -1 | 1 | 1 | 1 |
Результаты МСБ анализируются с помощью диаграмм рассеяния либо с помощью выборочных ортогональных МП.
2.2 Проведение эксперимента на объекте исследования
Проведем, используя лабораторный стенд, эксперимент в районе базовой точки: х1=30, х2=40, х3=30, х4=40, х5=20, х6=40, х7=40, х8=20. Примем за интервал варьирования Δxi=10 (I=1,…,8). Число параллельных опытов m=2. Целевую функцию запишем как сумму: y=y1+y2.
Таблица 2.4
N | t,c | X1 | X2 | X3 | X4 | Y1 | N | t,c | X5 | X6 | X7 | X8 | Y2 | Y | Yср | S2 |
1 | 26,1 | 40 | 60 | 5 | 15 | 111,3 | 1 | 10,5 | 30 | 65 | 50 | 20 | 102,3 | 213,6 | ||
2 | 26,8 | 40 | 60 | 5 | 15 | 101,8 | 2 | 11,1 | 30 | 65 | 50 | 20 | 103,7 | 205,5 | 209,55 | 32,805 |
3 | 19,4 | 20 | 60 | 5 | 15 | 84,9 | 3 | 81,1 | 30 | 65 | 70 | 40 | 205,8 | 290,7 | ||
4 | 20,1 | 20 | 60 | 5 | 15 | 92,2 | 4 | 81,8 | 30 | 65 | 70 | 40 | 202,1 | 294,3 | 292,5 | 6,48 |
5 | 12,9 | 40 | 40 | 5 | 15 | 77,5 | 5 | 23,2 | 10 | 65 | 70 | 20 | 106,7 | 184,2 | ||
6 | 13,9 | 40 | 40 | 5 | 15 | 75,4 | 6 | 24,1 | 10 | 65 | 70 | 20 | 109,3 | 184,7 | 184,45 | 0,125 |
7 | 31 | 20 | 40 | 25 | 15 | 107,1 | 7 | 17,9 | 30 | 85 | 50 | 40 | 137,1 | 244,2 | ||
8 | 32 | 20 | 40 | 25 | 15 | 113 | 8 | 18,3 | 30 | 85 | 50 | 20 | 130,7 | 243,7 | 243,95 | 0,125 |
9 | 54,1 | 20 | 40 | 5 | 35 | 62,9 | 9 | 52,4 | 10 | 85 | 50 | 40 | 88,5 | 151,4 | ||
10 | 55,3 | 20 | 40 | 5 | 35 | 64,8 | 10 | 52,9 | 10 | 85 | 50 | 40 | 100,5 | 165,3 | 158,35 | 96,605 |
11 | 70,6 | 20 | 60 | 5 | 35 | 88,4 | 11 | 14,5 | 10 | 85 | 50 | 20 | 106 | 194,4 | ||
12 | 71,5 | 20 | 60 | 5 | 35 | 78,8 | 12 | 14,9 | 10 | 85 | 50 | 20 | 105,5 | 184,3 | 189,35 | 51,005 |
13 | 88,1 | 20 | 60 | 25 | 35 | 141,4 | 13 | 47,5 | 30 | 65 | 50 | 40 | 96,1 | 237,5 | ||
14 | 88,9 | 20 | 60 | 25 | 35 | 155,3 | 14 | 47,9 | 30 | 65 | 50 | 40 | 101,2 | 256,5 | 247 | 180,5 |
15 | 39,8 | 20 | 60 | 25 | 15 | 144,9 | 15 | 26,6 | 30 | 65 | 70 | 20 | 153,7 | 298,6 | ||
16 | 40,8 | 20 | 60 | 25 | 15 | 141,2 | 16 | 27,6 | 30 | 65 | 70 | 20 | 163,4 | 304,6 | 301,6 | 18 |
17 | 92,7 | 40 | 60 | 25 | 35 | 190,4 | 17 | 41,1 | 10 | 65 | 50 | 40 | 79,1 | 269,5 | ||
18 | 93,7 | 40 | 60 | 25 | 35 | 186,9 | 18 | 44,8 | 10 | 65 | 50 | 40 | 74,2 | 261,1 | 265,3 | 35,28 |
19 | 43,9 | 40 | 60 | 25 | 15 | 183,3 | 19 | 5,5 | 10 | 65 | 50 | 20 | 71 | 254,3 | ||
20 | 45,1 | 40 | 60 | 25 | 15 | 190,3 | 20 | 8,2 | 10 | 65 | 50 | 20 | 66,7 | 257 | 255,65 | 3,645 |
21 | 7 | 20 | 40 | 5 | 15 | 63,4 | 21 | 76,1 | 30 | 65 | 70 | 40 | 162,9 | 226,3 | ||
22 | 8,4 | 20 | 40 | 5 | 15 | 56,5 | 22 | 76,5 | 30 | 65 | 70 | 40 | 161 | 217,5 | 221,9 | 38,72 |
23 | 79,7 | 20 | 40 | 25 | 35 | 111,7 | 23 | 31,4 | 10 | 85 | 70 | 20 | 149,6 | 261,3 | ||
24 | 80,1 | 20 | 40 | 25 | 35 | 110 | 24 | 32,3 | 10 | 85 | 70 | 20 | 154,1 | 264,1 | 262,7 | 3,92 |
25 | 35,7 | 40 | 40 | 25 | 15 | 147 | 25 | 72,8 | 10 | 65 | 70 | 40 | 113 | 260 | ||
26 | 36,5 | 40 | 40 | 25 | 15 | 158,2 | 26 | 73,2 | 10 | 65 | 70 | 40 | 106 | 264,2 | 262,1 | 8,82 |
27 | 60,6 | 40 | 40 | 5 | 35 | 85,7 | 27 | 56,8 | 30 | 85 | 50 | 40 | 130,2 | 215,9 | ||
28 | 61,5 | 40 | 40 | 5 | 35 | 77,5 | 28 | 57,3 | 30 | 85 | 50 | 40 | 128,4 | 205,9 | 210,9 | 50 |
29 | 83,4 | 40 | 40 | 25 | 35 | 148,6 | 29 | 36,2 | 30 | 85 | 70 | 20 | 203,3 | 351,9 | ||
30 | 84,4 | 40 | 40 | 25 | 35 | 149 | 30 | 37,4 | 30 | 85 | 70 | 20 | 205,1 | 354,1 | 353 | 2,42 |
31 | 75,7 | 40 | 60 | 5 | 35 | 109,7 | 31 | 78,7 | 10 | 85 | 70 | 40 | 135,7 | 245,4 | ||
32 | 76,7 | 40 | 60 | 5 | 35 | 112,4 | 32 | 79,1 | 10 | 85 | 70 | 40 | 152 | 264,4 | 254,9 | 180,5 |
Таблица 2.5