Под числом степеней свободы в статистике понимают разность между числом опытов и количеством коэффициентов модели (уравнения регрессии), вычисленных по результатам опытов.
Так, например, при двух факторах число N опытов равно четырем, а для определения коэффициентов уравнения регрессии достаточно результатов трех опытов. Таким образом, в рассматриваемом случае, число степеней свободы, равное единице, может быть использовано для проверки адекватности модели. Величина и знак коэффициента регрессии указывают на вклад данного фактора в общий результат при переходе с нулевого на верхний или нижний уровень фактора. Линейным называют эффект,характеризующий линейную зависимость параметра оптимизации от соответствующих факторов.
Эффектом взаимодействияназывают эффект, характеризующий совместное влияние нескольких факторов на параметр оптимизации. Полный факторный эксперимент позволяет количественно оценить линейные эффекты и все эффекты взаимодействия. Для полного факторного эксперимента типа 22 линейное уравнение регрессии с учетом эффектов взаимодействия можно представить выражением .
При увеличении числа факторов количество возможных сочетаний уровней быстро возрастает, поэтому возникает необходимость в некоторых приемах построения матриц. Рассмотрим два наиболее простых приема.
Первый прием основан на правиле чередования знаков. В первом столбце (x1)знаки чередуются поочередно, во втором они чередуются через 2, в третьем – через 4, в четвертом – через 8, в пятом – через 16 и т. д. по степеням двойки.
Второй прием основан на последовательном достраивании матрицы. Для этого при добавлении нового фактора необходимо повторить комбинации уровней исходного плана сначала при значении нового фактора на верхнем уровне, а затем на нижнем.
Матрица планирования– это таблица, в которой представлены верхние и нижние уровни факторов в их кодовом обозначении, номера опытов, а также параметры оптимизации в кодовом обозначении (параметров оптимизации может быть один или несколько). Кодирование параметров оптимизации проще, чем кодирование факторов. Оно выполняется простой заменой натурального значения параметра на Y, например, вместо скорости подачи S вводим кодовое обозначение Y1. Каждая строка матрицы планирования представляет собой один опыт и содержит номер опыта, верхний, нижний уровни факторов и значение параметра оптимизации, полученное в результате проведения опыта. Содержание любой горизонтальной строки в матрице планирования представляет собой условия проведения одного опыта. Так, например, строка 3 в матрице, представленной в табл. (2.2), свидетельствует о том, что третий номер опыта следует проводить на нижнем уровне фактора X1, на верхнем уровне фактора Х2, а полученное значение параметра оптимизации будет представлять у3 (его находят измерением при проведении третьего опыта).Проведя все опыты, условия которых представлены в строках матрицы, получают результат многофакторного эксперимента, который в дальнейшем подлежит статистической обработке. В подавляющем большинстве технологические процессы механической и физико-технической обработки изделий машиностроения характеризуются не одним, а сразу несколькими параметрами (например, обработанная поверхность характеризуется шероховатостью, волнистостью, точностью размера, формы и т.д.).
Таблица 2.2
Матрица планирования для эксперимента N = 22 =4
№ опыта | Факторы | Параметр оптимизации | |
x1 | x2 | y | |
– | – | y1 | |
+ | – | y2 | |
– | + | y3 | |
+ | + | y4 |
Буквенное выражение содержания матрицы:
а
в
ав
Таблица 2.3
Матрица планирования для эксперимента N = 23 =8
№ опыта | Факторы | Параметр оптимизации | ||
x1 | x2 | x3 | y | |
– | – | – | y1 | |
– | – | + | y2 | |
– | + | – | y3 | |
– | + | + | y4 | |
+ | – | – | y5 | |
+ | – | + | y6 | |
+ | + | – | y7 | |
+ | + | + | y8 |
При этом, что очень важно, названные параметры могут зависеть от одних и тех же факторов.
В этих условиях при проведении опытов согласно матрице планирования целесообразно измерять не только один, но и другие интересующие экспериментатора параметры.
Такой подход позволит значительно сэкономить время на проведение экспериментов, поскольку данные опытов по другим параметрам будут использованы при обработке результатов с целью получения уравнений регрессии, связывающих независимые факторы с интересующими нас параметрами.
Описанная процедура применима и к другим матрицам планирования, содержащим 2 и более факторов (табл. 2.3 и 2.4), но эти матрицы должны обладать нижеследующими свойствами [1]:
- симметричность матрицы планирования эксперимента относительно центра эксперимента: алгебраическая сумма элементов вектора-столбца каждого фактора равна 0, т.е.
(2.8)
i-номер опыта; k-номер фактора; N – число опытов в матрице планирования.
- нормированность: сумма квадратов, элементов каждого столбца равна числу опытов
. (2.9)
- ортогональность матрицы планирования: для двух любых неравных друг другу факторов сумма скалярного произведения двух вектор-столбцов равна 0, то есть
. (2.10)
- рототабельность: точка в матрице планирования подбирается так, что точность предсказаний параметра оптимизации одинакова на равных расстояниях и не зависит от направления.
При увеличении числа факторов от 2 до к матрицу планирования N=22 =4 можно последовательно достраивать в соответствии с числом факторов, которое выбрал экспериментатор при этом он получит ту матрицу, которая требуется.
Такая процедура достраивания представлена в табл. 2.4.
В табл. 2.2 и 2.4 в крайнем правом столбце указаны строчные буквы латинского алфавита, совокупность которых представляет собой содержание матрицы планирования в более компактном буквенном выражении.
Подготовку, проведение многофакторных экспериментов можно выполнять в последовательности: обоснование эксперимента; выбор параметров оптимизации; выбор факторов, расчет уровней, их кодирование; выбор плана, матрицы планирования; проведение эксперимента; обработка результатов эксперимента; проверка уравнения регрессии на адекватность; анализ результатов; формулировка выводов и научно-обоснованных рекомендаций.
Таблица 2.4
Матрица планирования при увеличении числа факторов
№ | x0 | x1 | x2 | x3 | x4 | x5 | Буквенное выражение содержания матрицы |
+ | + | + | + | + | + | abcde | |
+ | – | + | + | + | + | bcde | |
+ | + | – | + | + | + | acde | |
+ | – | – | + | + | + | cde | |
+ | + | + | – | + | + | abde | |
+ | – | + | – | + | + | bde | |
+ | + | – | – | + | + | ade | |
+ | – | – | – | + | + | de | |
+ | + | + | + | – | + | abce | |
+ | – | + | + | – | + | bce | |
+ | + | – | + | – | + | ace | |
+ | – | – | + | – | + | ce | |
+ | + | + | – | – | + | abe | |
+ | – | + | – | – | + | be | |
+ | + | – | – | – | + | ae | |
+ | – | – | – | – | + | e | |
+ | + | + | + | + | – | abcd | |
+ | – | + | + | + | – | bcd | |
+ | + | – | + | + | – | acd | |
+ | – | – | + | + | – | cd | |
+ | + | + | – | + | – | abd | |
+ | – | + | – | + | – | bd | |
+ | + | – | – | + | – | ad | |
+ | – | – | – | + | – | d | |
+ | + | + | + | – | – | abc | |
+ | – | + | + | – | – | bc | |
+ | + | – | + | – | – | ac | |
+ | – | – | + | – | – | c | |
+ | + | + | – | – | – | ab | |
+ | – | + | – | – | – | b | |
+ | + | – | – | – | – | a | |
+ | – | – | – | – | – | (1) |
Проведение опытов
После выбора плана эксперимента, параметров оптимизации, определения уровней и интервалов варьирования факторов переходят к эксперименту.
Каждая строка матрицы – это условие опыта. Для исключения систематических ошибок рекомендуется опыты, предусмотренные матрицей, проводить в случайной последовательности, а не в той последовательности, которая приведена в колонке «номер опыта» матрицы планирования. Эта процедура носит название рандомизации эксперимента во времени.
Случайную последовательность опытов можно определить, воспользовавшись таблицей случайных чисел (табл. 2.5) или выполнив компьютерную рандомизацию. Если требуется, например, провести в случайной последовательности восемь опытов матрицы, то из случайного места таблицы (2.5) последовательно выписываем числа, лежащие в интервале от 1 до 8, при этом отбрасываем уже выписанные числа и числа, больше восьми.
Так, например, начиная с числа 87 (1-я строка табл. 2.5) [3], получаем следующую последовательность реализации опытов: 72831456 вместо последовательности: 345678.
Таблица 2.5
Фрагмент таблицы случайных чисел
Для компенсации влияния случайных погрешностей каждый опыт рекомендуется повторить п раз.
Дата добавления: 2019-10-16; просмотров: 692;