Планы первого порядка

Планы первого порядка [4] позволяют находить линейные уравнения регрессии (1) и нелинейные уравнения (2) с членами, учитывающими эффекты взаимодействия факторов:

ŷ = ; (1)

ŷ = = (2)

=

+ b ₁₂₃ x ₁ x ₂ x ₃ + b ₁₂₄ x ₁ x ₂ x ₄+...+ b _(k-2)(k-1)k x _k-2 x _k-1 x _k.

Для удобства программирования расчетов в состав уравнения регрессии вводят фиктивную переменную х₀ = +1 во всех опытах эксперимента:

ŷ = ;

ŷ =

Для РАМПЭ наибольшее распространение получили двухуровневые (m_j = m = 2) ортогональные D-оптимальные планы первого порядка типа
2^{( k-a )}. При таких планах все факторы в кодированном виде могут иметь только два значения (x_j = +1 и x_j = -1). Тип плана обозначает формулу для расчета числа его опытов без их повторения (N): N = m^(k-a) = 2^(k-a), где k > a, а равно 0,1, 2, 3,...

При а = 0 план типа 2^{( k-a )} является планом ПФЭ типа 2 ^k, а при а > 0 – планом ДФЭ.

Планы, отвечающие условиям ортогональности, позволяют любой коэффициент уравнения регрессии рассчитывать по одной формуле

,

где i – номер опыта в плане эксперимента; b_d – коэффициент, учитывающий эффект факторов x_d; y_i – свойства объекта, измеренные при проведении соответствующего опыта; N – число опытов в эксперименте.

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

,

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

Условием ортогональности плана эксперимента является выполнение условия

= 0 при u ¹ j и u, j равных 0, 1, 2,..., k.

Для D-оптимальных планов должны выполняться следующие условия:

при j равном 1, 2,..., k;

= N при j равном 0, 1, 2,..., k.

Выбор плана эксперимента начинается с расчета необходимого числа опытов (N_необх.) или его задания (N_зад.). При этом при расчете выбранного числа опытов в эксперименте (N) для двухуровневых планов первого порядка должны выполняться следующие соотношения:

N ³ N_необх.; (3)

N_необх. ³ k + 1; (4)

N_необх. ³ L + 1; (5)

N_зад. ³ N, (6)

где L – общее число коэффициентов в выбранном семействе полиномов (число отрезков ряда Тейлора).

При расчете N_{необх .} задаются видом полинома (типом и числом коэффициентов уравнения регрессии L), а при задании числа опытов определяют вид семейства полиномов, в котором возможно найти уравнение регрессии для данного числа опытов и числа изменяемых факторов в эксперименте:

k _max = N_зад. – 1; (7)

L _max = N_зад. – 1. (8)

Пример. Мы решили исследовать влияние на свойство y четырех факторов x_j (k = 4) и описать их зависимость уравнением регрессии в виде следующего нелинейного полинома (L = 11):

b ₁₄ x ₁ x ₄ + b ₂₃ x ₂ x ₃ + b ₂₄ x ₂ x ₄ + + b ₃₄ x ₃ x ₄.

Тогда совместное выполнение соотношений (4,5) даст:

N_необх. ³ k + 1 ³ 4+1 ³ 5;

N_необх. ³ L + 1 ³ 11+1 ³ 12;

N ≥ N_необх. ³ 12.

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

Первоначально выберем план из группы планов первого порядка типа 2^(k-а). Для четырех факторов (k = 4) возможна реализация плана полного факторного эксперимента типа 2^k (N = 2⁴ = 16) и двух планов дробного факторного эксперимента типа 2^k-a: полуреплика ПФЭ типа 2^k-1 (N = 2^4-1 = = 8) и четверть реплики ПФЭ типа 2^k-2 (N = 2^4-2 = 4).

Очевидно, что нашему условию (N ³ 12) отвечает только план первого порядка типа 2^k с N = 16. ПланПФЭ типа 2⁴ позволяет получить наиболее точные оценки коэффициентов уравнения регрессии. При ПФЭ все выборочные коэффициенты уравнения регрессии являются достаточно точными, "несмешанными" оценками соответствующих генеральных коэффициентов b_d» b_d.

Для построения ортогонального и D-оптимального плана ПФЭ типа 2⁴ воспользуемся одним из распространенных приемов, заключающемся в следующем.

Делается заготовка плана в виде таблицы (план-матрица эксперимента), в которой предусматривается не менее N строк и (L_max +1) столбцов для х_d.

В первый столбец таблицы заносят номера строк, соответствующие номерам опытов. Во второй столбец – кодированные значения фиктивного фактора х₀ (во всех строках плана х₀ = + 1). В третий столбец – кодированные значения первого фактора (х₁) в виде последовательного чередования значений (+1) и (–1). В следующем, четвертом столбце, для х₂ выбранная комбинация чередований в предыдущем столбце знаков (+1) и (–1) удваивается, например: после двух знаков (+1) следуют два знака (–1). По аналогичному принципу удвоения комбинации чередования знаков предыдущего столбца заполняются и последующие столбцы для всех оставшихся факторов x_d.

Столбцы для оценки эффектов взаимодействия факторов (х₁₂, х₁₃, х₂₃ и др.) заполняются числами, полученными путем перемножения знаков для соответствующих факторов в соответствующих строках таблицы (х₁₂ = х₁х₂, х₁₃ = х₁х₃ и т.д.).

Правильность составления плана проверяется по выполнению условия его D-оптимальности .

Построенный по этому приему план приведен в табл. 14.

Нетрудно проверить, что данный план является ортогональным и D-оптимальным.

План с натуральными значениями факторов X_j строится исходя из плана с кодированными значениями путем замены знаков (+1) и (–1) на соответствующие им натуральные значения для данного фактора.

Довольно часто на практике приходится задаваться не видом полинома, а числом опытов из-за дефицита ресурсов для проведения эксперимента (времени, средств и др.). В этом случае выбор плана эксперимента начинают с расчета параметров полинома, которые возможно определить при N_зад..

Допустим, что N_зад. = 10, k = 4, и уравнение регрессии необходимо получить в следующем виде (L = 5):

ŷ = b₀ +b₁x₁ +b₂x₂ +b₃x₃ +b₄x₄.

Таблица 14

План ПФЭ типа 2⁴

Номер опыта, i	Кодированные значения факторов xd	y
х0	х1	х2	х3	х4	х12	х13	х14	х23	х24	х34
1	+1	+1	+1	+1	+1	+1	+1	+1	+1	+1	+1
2	+1	-1	+1	+1	+1	-1	-1	-1	+1	+1	+1
3	+1	+1	-1	+1	+1	-1	+1	+1	-1	-1	+1
4	+1	-1	-1	+1	+1	+1	-1	-1	-1	-1	+1
5	+1	+1	+1	-1	+1	+1	-1	+1	-1	+1	-1
6	+1	-1	+1	-1	+1	-1	+1	-1	-1	+1	-1
7	+1	+1	-1	-1	+1	-1	-1	+1	+1	-1	-1
8	+1	-1	-1	-1	+1	+1	+1	-1	+1	-1	-1
9	+1	+1	+1	+1	-1	+1	+1	-1	+1	-1	-1
10	+1	-1	+1	+1	-1	-1	-1	+1	+1	-1	-1
11	+1	+1	-1	+1	-1	-1	+1	-1	-1	+1	-1
12	+1	-1	-1	+1	-1	+1	-1	+1	-1	+1	-1
13	+1	+1	+1	-1	-1	+1	-1	-1	-1	-1	+1
14	+1	-1	+1	-1	-1	-1	+1	+1	-1	-1	+1
15	+1	+1	-1	-1	-1	-1	-1	-1	+1	+1	+1
16	+1	-1	-1	-1	-1	+1	+1	+1	+1	+1	+1

Оценим возможности двухуровнего плана первого порядка с N = 10:

k _max = N_зад. – 1= 10–1= 9;

L _max = N_зад. – 1= 10–1= 9.

Из данных равенств следует, что план с 10 опытами позволяет решить нашу задачу (k = 4 и L = 5), так как при десяти опытах можно оценить в эксперименте одновременное влияние до девяти факторов (k ≤ k _max ) по уравнению регрессии в виде полинома с девятью коэффициентами (L ≤ L _max ).

При N_зад. = 10 из планов первого порядка возможна реализация только плана ДФЭ, так как план ПФЭ для четырех факторов насчитывает 16 опытов.

Воспользуемся методом дробных реплик для построения планов ДФЭ типа 2^(k-a). Наиболее близким по числу опытов к N_зад. = 10 является полуреплика (1/2 часть) плана ПФЭ, т.е. план ДФЭ типа 2^(4-1) с числом опытов N = 8. Проверка показывает, что план типа 2^(4-1) пригоден для решения поставленной задачи, так как выполняются следующие соотношения:

N ³ k + 1 ³ 4+1 ³ 5 (8 > 5);

N ³ L + 1 ³ 5+1 ³ 6 (8 > 6).

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

Построение планов ДФЭ начинают по тому же приему, что и при построении планов ПФЭ для числа факторов, равных разности (k–а). В нашем случае k–a = 4–1 = 3.

Поэтому построим первоначально заготовку плана ДФЭ типа 2^(4-1) в виде плана ПФЭ типа 2³, предусмотрев в нем не менее 6 колонок для x_d (табл. 15).

Таблица 15

Заготовка плана ДФЭ типа 2^(4-1)

Но- мер опыта i	Кодированные значения факторов xd	y
х0	х1	х2	х3	х4	х12	...
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

При заполнении столбца для фактора х₄ принцип удвоения чередований уровней, применяемый для построения планов ПФЭ, не подходит, так как его использование в данном столбце даст только знаки (+1) и такой план не будет являться ортогональным. Если же для заполнения столбца х₄ воспользоваться произведением двух и более других факторов в одной строке плана (так называемым генерирующим соотношением), то тогда план будет и ортогональным и D-оптимальным.

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

Составим и проанализируем все возможные варианты генерирующего соотношения для х₄:

х₄ = х₁х₂ (I); x₄ = x₁x₃ (II); x₄ = x₂x₃ (III); x₄ = x₁x₂x₃ (IV).

Для данных генерирующих соотношений рассчитаем определяющие контрасты умножением левой и правой частей соответствующего генерирующего соотношения на х₄:

; ; ; .

Так как х₄ = ± 1, то и определяющие контрасты можно выразить равенствами

1 = х₁х₂х₄; 1 = х₁х₃х₄; 1 = x₂×x₃×x₄; 1 = x₁×x₂×x₃×x₄.

Перемножив левые и правые части определяющих контрастов на каждый фактор, можно определить в плане столбцы с одинаковым порядком чередования знаков (+1) и (-1), например для фактора х₁:

; ;

; .

Эти равенства показывают, что при генерирующих соотношениях I-IV выборочный коэффициент b₁ будет служить оценкой влияния на y не только фактора x₁, но и других:

b₁» b₁ + b₂₄; b₁» b₁ + b₃₄; b₁» b₁ + b₁₂₃₄; b₁» b₁ + b₂₃₄.

Эффекты взаимодействия трех и более факторов обычно близки к нулю и ими можно пренебрегать [4]:

b₁» b₁ + b₂₄; b₁» b₁ + b₃₄; b₁» b₁; b₁» b₁.

Поэтому можно считать, что несмешанные оценки эффекта влияния фактора x₁ на свойство y могут быть получены при реализации плана ДФЭ с генерирующими соотношениями III и IV для фактора х₄. Результаты проверки на "смешиваемость" остальных эффектов приведены в табл. 16.

Данные табл. 16 показывают, что при любом генерирующем соотношении точными (несмешанными) будут пять коэффициентов. Для генерирующего соотношения I точными будут все коэффициенты, оценивающие эффект фактора х₃, при II – эффект фактора х₂, при III – эффект фактора х₁, а при IV – линейные эффекты всех факторов.

Таблица16

Параметры проверки разрешающей силы дробной реплики типа 2^(4-1)

Параметр	Выражение для определения параметра
Генерирующее соотношение	x4 = x1x2	x4 = x1x3	x4 = x2x3	x4 = x1x2x3
Определяющий контраст	1 = x1x2x4	1 = x1x3x4	1 = x2x3x4	1 = x1x2x3x4
Оценки точности («смешанности»*) коэффициентов уравнения регрессии	b0» b0+b124 b1» b1 + b24 b2» b2 + b14 b3» b3 + b1234 b4» b4 + b12 b12» b12+b4 b13» b13+b234 b14» b14+b2 b23» b23+b134 b24» b24+b1 b34» b34+b123	b0» b0+b134 b1» b1 + b34 b2» b2 + b1234 b3» b3 + b14 b4» b4 + b13 b12» b12+b234 b13» b13+b4 b14» b14+b3 b23» b23+b124 b24» b24+b123 b34» b34+b1	b0» b0+b234 b1» b1 + b1234 b2» b2 + b34 b3» b3 + b24 b4» b4 + b23 b12» b12+b134 b13» b13+b124 b14» b14+b123 b23» b23+b4 b24» b24+b3 b34» b34+b2	b0» b0+b1234 b1» b1 + b234 b2» b2 + b134 b3» b3 + b124 b4» b4 + b123 b12» b12+b34 b13» b13+b24 b14» b14+b23 b23» b23+b14 b24» b24+b13 b34» b34+b12
* - жирным шрифтом выделены несмешанные коэффициенты

Так как по заданию нам необходимо получить уравнение регрессии с линейными эффектами всех факторов, и эти эффекты должны быть наиболее точными, то выбираем генерирующее соотношение IV и в соответствии с ним заполняем колонку плана для х₄ (табл. 17).

Таблица 17

План ДФЭ типа 2^(4-1) с генерирующим соотношением x₄ = x₁x₂x₃

Но- мер опыта i	Кодированные значения факторов	y
х0	х1	х2	х3	х4	х12 и х34	х13 и х24	х14 и х23
1	+1	+1	+1	+1	+1	+1	+1	+1
2	+1	-1	+1	+1	-1	-1	-1	+1
3	+1	+1	-1	+1	-1	-1	+1	-1
4	+1	-1	-1	+1	+1	+1	-1	-1
5	+1	+1	+1	-1	-1	+1	-1	-1
6	+1	-1	+1	-1	+1	-1	+1	-1
7	+1	+1	-1	-1	+1	-1	-1	+1
8	+1	-1	-1	-1	-1	+1	+1	+1
9	+1	0	0	0	0	0	0	0
10	+1	0	0	0	0	0	0	0

Этот план с N = 8 является ортогональным и D-оптимальным. По данному плану есть возможность оценить еще 6 эффектов парного влияния факторов, однако, как показывают данные табл. 16, расчеты приведут к получению смешанных коэффициентов уравнения регрессии (т.е. неточно отражающих парное влияние соответствующих факторов), так как комбинации знаков (изменения значений факторов в опытах) совпадают у х₁₄ и х₂₃, х₁₃ и х₂₄ и др.

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

В соответствии с общепринятыми рекомендациями запланируем опыты для определения при нулевых кодированных значениях всех исследуемых факторов, т.е. в центре области изменения факторов (табл. 17).

В качестве планов первого порядка для проведения РАМПЭ можно использовать не только дробные реплики ПФЭ, но и некоторые другие планы ДФЭ, например планы Плакетта-Бермана [4].

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

Алгоритмы расчетов при РАМПЭ по планам первого порядка зависят от наличия повторений опытов. Познакомьтесь с ними самостоятельно [4].

После реализации плана эксперимента первого порядка довольно часто найденное уравнение регрессии оказывается неадекватным. В этом случае обычно переходят к выполнению РАМПЭ для поиска уравнения регрессии в семействе полиномов второго порядка по результатам специально спланированных экспериментов.

Дата добавления: 2014-01-07