Рандомизация


Перед реализацией плана на объекте необходимо произвести рандомизацию – с помощью таблицы равномерно распределённых случайных чисел (табл. П.6) определить последовательность реализации матрицы планирования в каждой из m серий опытов. Для этого в качестве начала выбирается любое число из табл. П.6 и записывается в столбец k1 из табл. 5.2 на место g=1. Остальные места этого столбца заполняют числа от 1 до N, следующие по порядку из табл.П.6 за выбранным начальным. Следует обращать внимание на то, чтобы числа в столбцах табл. 5.2 не повторялись дважды. Пусть, например, при g=4 k4=8, это значит, что в первой серии испытаний точка clip_image002реализуется восьмой по порядку.

Аналогично рандомизируются испытания в каждой из оставшихся серий экспериментов; порядок реализации записывается в столбцах k2, k3, …, km. Результаты эксперимента в каждой из серий испытаний записываются в столбцах Y1, Y2, …, Ym.

Проверка воспроизводимости – это проверка выполнения второй предпосылки регрессионного анализа об однородности выборочных дисперсий clip_image004. Задача состоит в проверке гипотезы о равенстве дисперсий s2{Y1}=s2{Y2}=…=s2{YN} при экспериментах соответственно в точках clip_image006. Оценки дисперсий находятся по формуле

clip_image008. (5.5)

Так как все дисперсии получены по выборкам одинакового объёма m, то число степеней свободы для всех дисперсий одинаково и равно

clip_image010. (5.6)

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

clip_image012 (5.7)

Если вычисленное значение критерия G окажется меньше табличного значения Gкр, найденного для q%-го уровня значимости, nзн=n2=N – числа степеней свободы знаменателя (например, для q=5%; nчисл=3-1=2; nзн=8, Gкр=0,5157, см. табл. П.5), то гипотеза об однородности дисперсий принимается. При этом всю группу дисперсий clip_image004[1] можно считать оценкой S2{Y} одной и той же генеральной дисперсии воспроизводимости s2{Y}, откуда

clip_image015. (5.8)

Если проверка на воспроизводимость дала отрицательный результат, то остаётся признать либо невоспроизводимость эксперимента относительно управляемых переменных вследствие наличия флуктуаций неуправляемых и не контролируемых переменных, создающих на выходе большой уровень «шума», либо наличие грубого промаха в строке, откуда взята дисперсия clip_image017. В первом случае следует увеличить число параллельных опытов, во втором – найти грубый промах и заменить его на результат доброкачественного измерения при соответствующей комбинации факторов. Если это по каким-то причинам невозможно, то, чтобы не нарушать предпосылки использования критерия Кохрена, на место грубого промаха следует поместить среднюю арифметическую величину clip_image019данной строки.

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

Легко заметить, что исходный план (табл. 5.1) содержит много больше строк, чем столбцов и, следовательно, из результатов эксперимента согласно условию решения нормальных уравнений (5.2) можно получить дополнительную информацию, то есть расширить модель. Безусловно, это относится к средней арифметической всего эксперимента, т.е. к отклику в базовой точке b0, для расчёта которого можно ввести фиктивную переменную xg=+1 для всех строк. Оставшиеся свободными столбцы можно использовать для нахождения оценок коэффициентов при парных взаимодействиях и т.п. При этом соответствующие величины xixj, xixjx1 получаются простым перемножением соответствующих столбцов исходного плана. Тогда математическая модель объекта, получающаяся в результате ПФЭ может быть представлена в виде

Y=b0+b0x1+…+bnxn+b12x1x2+…+b(n-1)xn-ixn+b123x1x2x3+…+b123…nx1x2x3xn . (5.9)

Однако вследствие того, что из ограниченного числа опытов нельзя получить точные значения коэффициентов bi, а только их независимые оценки bi, вся математическая модель становится оценочной

clip_image021. (5.10)

Таблица 5.2

Матрица планирования ПФЭ типа N=23 и обработка его результатов

Номер строки

g

Порядок реализации опытов

l

Матрица планирования

эксперимента

Результаты

эксперимента

Первичная обработка результатов

Проверка адекватности

z0

z1

z2

z3

z4

z5

z6

z7

Ygi

Ygl

Ygm

clip_image023

clip_image025

clip_image027

clip_image029

k1

kl

km

x0

x1

x2

x3

x1x2

x1x3

x2x3

x1x2x3

1

1

6

8

+

+

+

+

2

7

5

4

+

+

+

+

3

3

7

6

+

+

+

+

4

8

2

7

+

+

+

+

5

6

3

2

+

+

+

+

6

4

4

1

+

+

+

+

7

2

1

5

+

+

+

+

8

3

8

3

+

+

+

+

+

+

+

+

Загрузка...