Планирование эксперимента ортогонально к дискретному дрейфу


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

Разбиение проводится с учётом требования ортогональности вектор-столбцов планирования как между собой, так и к вектору дрейфа, в пределах каждого блока. В частности, при исследованиях с помощью планов типа 2n для получения полиномиальной модели независимо от дрейфа обычную матрицу ПФЭ разбивают на p блоков

N = n1 + n2 + …+ np ,

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

clip_image002.

Такому произведению можно удовлетворить, используя в качестве блоков дробные реплики от ПФЭ так, чтобы блоковый эффект был смешан только с эффектами высших порядков.

Допустим, что изучается процесс с тремя неизвестными переменными и исследователя интересуют не только линейные эффекты, но и парные взаимодействия, при этом имеет место дрейф, обусловленный неоднородностью сырья. Предполагается провести ПФЭ типа 23 с реализацией первых четырёх опытов на одной партии сырья, а оставшиеся четыре – на другой.

Процедура разбиения матрицы планирования 23 на блоки состоит в следующем:

1. Первый блок строится как матрица планирования ДФЭ типа 23-1 с генерирующим соотношением, например, x3= x1x2;

2. Второй блок строится как матрица планирования ДФЭ типа 23-1 с генерирующим соотношением, например, x3= —x1x2;

3. Вводится четвёртая независимая переменная, характеризующая межблоковый дрейф. Эта переменная приравнивается тройному взаимодействию, которым почти всегда можно пренебречь Pб=x1x2x3.

Таблица 5.5

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

g

z0

z1

z2

z3

z4

z5

z6

z7

Y’g

x0

x1

x2

x3

x1x2

x1x3

x2x3

x1x2x3

1

+

+

+

+

Y1

I

блок

2

+

+

+

+

Y2

3

+

+

+

+

Y3

4

+

+

+

+

+

+

+

+

Y4

5

+

+

+

+

Y5+DYt

II

блок

6

+

+

+

+

Y6+DYt

7

+

+

+

+

Y7+DYt

8

+

+

+

+

Y8+DYt

Очевидно, что по самому принципу построения вектор-столбцы каждого блока ортогональны; две дробные реплики образуют матрицу ПФЭ типа 23. Ортогональность столбцов варьируемых факторов и дрейфа позволяет получить раздельные оценки коэффициентов регрессии (матрица планирования приведена в табл.5.5). Эксперимент в точках факторного пространства, принадлежащих первому блоку, даёт возможность получить значения выходной величины Y1, Y2, Y3, Y4, реализация второго блока имеет информацию о величинах Y5 + DYt; Y6 + DYt; Y7 + DYt; Y8 + DYt, где DYt – изменение выхода, обусловленное сменой сырья (т.е. дискретным дрейфом).

Легко убедится, что при подсчёте коэффициентов уравнения составляющие дрейфа DYt взаимно умножаются и коэффициенты определяются в «чистом виде», как если бы дрейфа не было. Например:

clip_image004

и так далее, кроме b0.

Аналогично проводится разбиение на два блока матрицы 24, при этом строятся две полуреплики типа 24-1 с генерирующим соотношением x4 = x1x2x3 и x4 = —x1x2x3, а блоковая составляющая равна Pб = x1x2x3x4.

Полный факторный эксперимент типа 24 может быть разбит и на 4 блока, каждый из которых строится как четверть реплика типа 24-2 с генерирующими соотношениями:

1) x3 = x1x2; x4 = x1x2x3;

2) x3 = —x1x2; x4 = x1x2x3;

3) x3 = x1x2; x4 = —x1x2x3;

4) x3 = —x1x2; x4= —x1x2x3.

Реализация матрицы планирования осуществляется по блокам. Так, например, в ПФЭ типа 23 сначала реализуются четыре точки первого блока на первой партии сырья, а затем – четыре точки второго блока на второй партии сырья. Планирование в условиях дискретного дрейфа позволяет проводить дублирующие серии экспериментов. Рандомизация порядка проведения опытов в каждой из серий проделывается поблочно, т.е. для ПФЭ типа 23 по таблице случайных чисел находят сначала числа от 1 до 4 для реализации первого блока, а затем числа от 5 до 8 для реализации второго.

Проверка воспроизводимости производится обычным способом по критерию Кохрена.

Коэффициенты уравнения регрессии находятся по обычным формулам (5.11).

Коэффициент межблокового дрейфа (равный для ПФЭ типа 23 clip_image006) опре-

деляется по формуле

clip_image008. (5.17)

Проверка значимости всех коэффициентов проводится с помощью tкритерия Стьюдента.

Математическое описание объекта будет иметь вид:

clip_image010. (5.18)

В том случае, если количество значимых коэффициентов аппроксимирующего уравнения будут меньше, чем число экспериментов N=2n, проводится проверка адекватности модели с использованием F-критерия Фишера.

Загрузка...