книги из ГПНТБ / Кориков А.М. Математические методы планирования эксперимента учеб. пособие
.pdf•пах может быть три двойных и один |
тройной эшфект Ь „ х , х , |
||
Поясним понятая "эффект фактора" и "эффект эзаикодей- |
|||
с т в ж " . Вклад фактора в величину |
параметра оптимизации при |
||
переходе от нижнего уровня к верхнему уровню |
называется |
||
э ф ф е к т о м |
ф а к т о р а |
(иногда его |
называют л и |
нейным,основным или главным эффектом). Он численно равен
удвоенное коэффициенту полинома. Для качественных факто
ров, изменяемых на двух уровнях, основной |
уровень |
не име |
||||
ет физического сшсла . Поэтому понятие "эффект фактора" |
||||||
является здесь естественным. Если еффект одного фактора |
||||||
зависит от уровня, |
на котором находится другой фактор, тс |
|||||
говорят, |
что имеет |
место |
э ф ф е к т |
в з а и м о д е й |
||
с т в и я |
двух факторов |
(иногда |
говорят |
двойной |
или пар |
|
ный эффект взаимодействия |
(эффект |
взаимодействия |
первого |
порядка). Также выясняется смысл тройного эффекта взаимо действия (эффект взаимодействия второго порядка) и т . д . Численно эффект взаимодействия равен удвоенному соответст вующему коэффициенту полинома.
§ 1-5» Д р о б н ы е р е п л и к и
Число опытов в полном факторном эксперименте превыша ет число коэффициентов линейной модели, причем тем больше, чем больше факторов. Разность между числом опытов и числом коэффициентов во многих случаях оказывается очень велика, и возникает естественное желание сократить число необходи мых опытов. Заманчиво сократить кх число за счет той ян-
- 3 1 -
гТорькцш, которая несущественна при построении линейных моделей.
Итак, нам заранее известно, что объект описывается ли нейным уравнением (предполагается, что эффекты, взаимодей ствия отсутствуют) или нас интеррсуют только линейные чле ны. Возникает вопрос: как построить ортогональный план, по зволяющим определить коэффициенты линейного уравнения, ко торый содержит меньшее число опытов, чем полный факторный эксперимент?
Попробуем построить такой план для случая h, = 4 . План полного факторного эксперимента состоит из 16 точек. Лдя получения линейной зависимости требуется, как кишмум, 5 точек, а чтобы выполнялось свойстзо ортогональности, таких точек делгло быть, как минимум, шесть.. Посмотрим, существу ет ли ортогональный план, для которого число точек
бт. < 16 .
Такой план существует - это план для случая "к = 3, представ ленный в табл.1.4.
т - |
i |
|
Г" |
"T |
опыта; |
|
xz |
i |
j |
I |
! |
+ |
+ |
! |
+ |
j |
2 |
! |
+ |
_ |
i |
+ |
i |
з |
! |
- |
+ |
1 |
+ |
|
4 |
! |
- |
|
|
+ |
i |
5 |
i |
+ |
+ |
i |
- |
i |
6 |
i |
+ |
- |
! |
- |
j |
7 |
j |
- |
+ |
• i |
|
j |
8 |
i |
- |
i |
|
||
|
i |
|
"Г |
|
т |
|
f |
х 1 х |
з ! |
+ |
1 |
+ |
! |
- |
i |
.f |
j |
- |
j |
- |
! |
+ |
i |
- |
| |
I |
- |
i |
|
- |
i |
- |
i |
|
|
||
|
ij |
+ |
j |
+ |
+ |
i |
T a r t i r w ry Т , 4
"T X 1X2X3
+ |
t |
+ |
|
- |
I |
- |
|
+• |
i |
- |
|
|
|
||
- |
i |
+ |
|
- |
|||
- |
i |
||
+ |
i |
+ |
|
|
+ |
||
+ |
i |
||
|
|||
i |
|
||
|
|
ъг
3 связи с тегл, что эЖекты взапшде.2ств:хя прпяшазтся равными нулю, МОЙИО воспользоваться лобым из столбцов, ха
рактеризующих эффекты взаимодействия для четвертое перемеГ'
ной, например |
столбцом |
1 , ^ |
^ (или - |
Х , х г х ^ ) . Приняв |
|
для |
фактора х ^ столбец |
х , x t x 3 |
, получш |
план, призеденны:'! |
|
в табл . 1 . 5 . |
|
|
|
|
|
|
|
|
|
|
Таблица Х.Ъ |
й |
опыта |
|
|
|
|
|
I |
+ |
|
+ |
|
|
2 |
+ |
|
+ |
|
|
3 |
|
|
+ |
|
|
4 |
|
|
+ |
|
5+
6+
8 |
|
I |
I |
|
|
||
Этот план содержит половину опытов полного факторного |
|||
эксперимента и носит |
название |
п о л у р е п л и к и . |
|
Используют также £ р |
е п л и к |
и ' . |
J р е п л z к и и |
т . д . |
|
|
|
Символически дробные реплики записывают следующие об
разом: |
|
где К. - общее число факторов; (k-t) |
- число факторов в |
плане полного факторного эксперимента, к которому прирав нивается дробная реплика; Ь - число линейных эффектов,при равненных к эффектам взаимодействия.
- 3 5 |
- |
Б табл.I.6 приведены условные обозначения дробных реп
лик и,количество опытов. Из таблицы видно, что целесообраз
ность применения дробных реплик возрастает с ростом коли
чества факторов.
т
Количест?о!
«акторов j'. Дробная реплика
3 |
I |
т |
|
|
Q |
|
||
• |
j |
|
реплика |
от |
2 |
|
||
|
|
|
||||||
4 |
• |
- |
|
реплика |
от |
2 4 |
|
|
|
1 2 |
|
|
|
|
5 |
|
|
5 |
I |
Т |
|
|
от |
2 |
|
|
~ реплика |
|
|
||||||
|
i |
|
|
|
|
|
|
|
6 |
• |
I |
|
реплика |
от |
2 6 |
|
|
|
I |
° |
|
|
|
7 |
|
|
7 |
; |
т |
|
|
|
|
||
I |
|
|
реплика от |
2 |
|
|||
|
|
|
|
|||||
5 |
:1 т~ реплика |
от |
25 |
|
||||
6 |
I |
т |
|
реплика |
от |
2 |
6 |
|
- |
|
|
|
|||||
|
• |
4 |
|
|
|
7 |
|
|
7 |
I |
т |
|
от |
|
|
||
! |
± реплика |
2 |
|
|
||||
8 |
|
т |
реплика |
от |
2 |
8 |
|
|
|
•~ |
|
|
|
||||
9 |
1 т |
|
реплика |
от |
2 |
9 |
|
|
|
— |
|
|
|
||||
10 |
j |
- |
i реплика |
от |
2 1 |
0 |
||
|
I |
64 |
|
|
|
|
|
|
I I |
|
|
|
реплика |
от |
2 1 |
1 |
|
12 |
I |
т |
|
|
|
|
ТО |
|
|
|
|
реплика |
от |
2 Х |
" |
||
|
|
|
|
|||||
13 |
!256 |
|
|
|
j o |
|||
j 512 |
реплика от |
2 |
|
|||||
ТА |
|
|
|
|
|
|||
1102A Р е г о г а к |
|
|
|
|
||||
|
а 0 |
1 |
|
|
||||
15 |
|
|
|
реплика от 2 1 Э |
|
|
« |
|
Таблица 1.6 |
|
Услов |
Кол-во |
;Кол-во |
ное |
опытов |
jопытов для |
обоз |
|ддя дро- |
{ПОЛНОГО |
начение •бной'ре |
{факторного |
|
|
плики |
(эксперим. |
2 3 - 1
0 5 - 2
,6-3
,7-4
, 5 - 1
,6 - 2
,7-3
0 8 - 4
„9-5
,10-6
,11-7
2 13-9
,14-10 Л5-11
4 8
8 |
. 16 |
8 32
8 64
8 128
16 32
16 64
16 128
16 |
256 - |
16 512
16 1024
16 2048
-16 4096
16 8192
16 16384
16 32768
- 3 4 -
Перед постановкой эксперимента го дробным репликам не обходимо решить - каким эффектом взаимодействия можно пре небречь и к какому это приведет риску.
Поставив х ^ в рассмотренном примере на место тройного эффекта взаимодействия, получили
|
Х 4 " Х 1 |
Л 2 Х 5 |
* |
|
|
Это соотношение называют |
г е н е р и р у ю щ и м |
|
|||
с о о т н о ш е н и е м , |
так |
как оно генерирует, или |
с о |
||
здает, дробную реплику. Умножив обе части генерирующего |
с о |
||||
отношения на I |
, |
|
|
|
|
|
асц |
- х , х Л х 3 х ^ |
|
||
получим в левой |
части |
единичный столбец, который обозначил |
|||
через |
|
|
|
|
|
|
I - |
^ , Х 2 Х 4 Х ^ . |
( 1 # 4 ) |
||
Это произведение |
называют |
о п р е д е л я ю щ и м |
|
||
к о н т р а с т о м . |
Итак, определяющим контрастом назы |
||||
вается символическое обозначение произзедения столбцов, |
|
||||
элементы которого равны |
+1 или - I . |
|
Определяющий контраст'позволяет установить разрешаю щую способность дробной реплики. Разрешающую способность
дробных реплак. мы обсудим в следующем параграфе, |
а |
сейчас |
||
отметим, |
что рассмотренные |
здесь дробные реплики, |
содерка- |
|
щие 2 |
"точек, называют |
р е г у л я р н ы м и |
.'Сущест |
вуют также нерегулярные дробные реплики, содержащие количество точек, не равное 2 'brt , например jя (2 к,) .
* При разработке технологических гроцессов дробные реп-.
- о б
ЛИКИ очень широко используют на стадии крутого восхождения (см.§ 4 - 3 ) , а также при математическом описании локальной области факторного пространства с узким интервалом измене ния переменных. Здесь же заметим, что при иллюстрации идеи оптимального использования факторного пространства в табл. В.2 (см.§ В-3) наш выписана дробная реплика типа 2 3 " * 1 . Применяя дробную реплику для планирования эксперимента при взвешивании трех объектов, можно не беспокоиться,о том, что найденные величины являются совместными оценками для
линейных эффектов и эффектов взаимодействия. Здесь из физи ческих соображений ясно, что эффекты взаимодействия сущест вовать не могут.
§ 1-6. Р а е р е т а ю щ а я |
с п о с о б н о с т ь |
д р о б н ы х |
р е п л и к |
Разрешающая способность дробной реплики задается сис
темой смешивания данной реплики. Она будет максимальной, если линейные эффекты смешаны с эффектами взаимодействия наибольшего возможного порядка.
Определяющий контраст позволяет определить систему
смешивания дробной реплики. Для того чтобы определить,
какой эффект смешан с данным, нужно помногить обе части определяющего контраста на столбец, соответствующий дан
ному эффекту." |
|
|
|
|
3—т |
|
|
|
|
При построении полуреплики 2 |
существует |
всего |
две |
|
возможности: приравнять х , к + £ х |
или |
- ^ |
X j . |
Поэто- |
i
му, если I = г , 1 г х } , то, учитывая, чтох.^ = 1, для имеем
|
|
|
|
|
|
- |
об |
- |
|
|
|
|
|
|
|
2 |
|
|
|
|
|
i |
|
|
|
|
|
Х ^ - х , х г |
x s - х.г осъ |
, |
|
' |
(Т . 5) |
||
для |
х , |
находим |
|
|
|
|
|
|
|
||
|
|
|
|
X , - « , . « - * - Х . Х . , |
|
|
|
„ . „ |
|||
для |
х 3 |
имеем |
|
|
|
|
|
|
|
|
|
|
|
|
|
х » — л , л а |
а: * - ж, |
|
|
|
(1 . 7) |
||
|
Соотношения (1.5). - ( 1 . 7 ) указывают |
на равенство |
стол |
||||||||
бцов в матрице планирования, например,.столбцы для х ( и |
|||||||||||
з^*зодинаковы. Поэтому коэффициент |
$ ( |
будет |
оценивать |
||||||||
сумму |
р , |
р г 5 |
. Это |
записывается следующим |
образом: |
|
|||||
|
* i |
— |
Pi+ |
Ргг. |
^ |
Рг+ |
Р » . |
- |
, V |
Pit • |
|
Полуреплики, в которых основные эффекты смешаны с двой ными эффектами взаимодействия, носят название штанов с раз решающей способностью Ш (по наибольшему числу факторов в определяющем контрасте); Такие планы принято обозначать:
|
При выборе полуреплики |
2^~^ возможны восемь |
решений: |
||||
|
*Е? х^ = х , хг |
|
5. |
= х , х 5 ,. |
|
||
|
2 ; |
= - Х 1 |
Х г ^ |
6 . j t j , = - X , х 5 |
, |
||
|
3 . х н = х г |
х 3 , |
7 . х ч |
= х , х ь " х & , |
|||
• |
4 . Х ц = - Х г Х ь , |
8 . х „ = |
- l ^ t j . |
||||
|
Разрешающая способность |
этих.полуреплик |
различна, Так, |
реплики 1-6 имеют по три фактора в определяющем контрасте,
а-.7,8*по четыре. Умножая определяющий контраст I |
= х , х г х 4 х ^ |
последовательно на х , , х г , х 3 , х 4 , определим |
систему |
смешивания данной реплики: |
|
- ы -
Реплики, в которых нет ни одного линейного эффекта, смешан ного с другим линейным эффектом или парным взаимодействием,
а все парные ззаилодействия смешаны друг с другом, носят
название планов с разрешающей способностью* 17 (по наиболь шему числу факторов в определяющем контрасте). Они имеют
обозначение |
2 j y * . Полуреплика, заданная |
определяющим кон- |
- частом I = |
+ x , x l x i x 4 , имеет только |
четные комбинации |
букв в каждой строке. Ее южно записать следующим образок,
считай отроку |
( I ) |
четной: |
|
|
|
|
|
|
|
|||
( I ) , |
аЛ , |
bd |
, |
a d |
, |
ас, |
ccL, be, |
abed. |
|
|||
А полуреплика, |
заданная |
I |
= - |
x , x , I |
, t ^ , имеет ™олько |
|||||||
нечетные |
комбинации |
|
|
|
|
|
|
|
|
|
||
а , о, |
с , |
c t , |
aid |
, |
dud |
, а б с , |
bed. |
|
||||
Такие полурешшки |
называют |
г л а в н ы м и , |
так как они |
|||||||||
обладают |
наибольшей разрешающей |
способностью. |
|
|
||||||||
|
|
|
|
|
|
|
|
С |
-г |
|
|
|
Рассмотрим |
теперь |
полуреплику.2 |
|
. При выборе полу- |
||||||||
реплики 2^~* в распоряжении экспериментатора |
тлеется |
22 |
||||||||||
возможных |
варианта. Так, x s можно приравнять |
к одному |
из |
6 парных взаимодействий. В этом случае получим полурепли ку с разрешающей способностью Ш. Очевидно, это будеп ' ;:э' лучший выбор полуреплики. Далее, х 5 можно приравнять" к одному из четырех тройных взаимодействий. Тогда получил план с разрешающей способностью 1У, и все линейные эффе- -
)
- 3 6 -
та будут и.лешаны с тройными взаимодействиями. И, наконец,
гюлурешшг-'. может быть задана генерирующими соотношениями
х5 = х 1 х 2 а с й х 4 или X s = - x 1 x a x J x J l . Определяющие кон
трастами в |
этом |
случае будут |
I = x t х г |
J C 3 |
x 4 xs |
ш Z = |
|
= - |
Х , х г |
х 5 х 4 |
х 5 . Такие полуреплики |
нося1 ! |
название пла- , |
||
нов |
с разрешающей способностью У и обозначаются |
2 у - х . |
|||||
. |
Мы не |
станем рассматривать |
полуреплики 2 ^ ; |
2 7 - 1 и |
т. д . Такими полурепликами редко пользуются на практике.
ст
Зедь пол"решшка 2 |
требует 32 опыта, а для |
эксперимен- |
татора выгодны планы 2 |
или 2 , требующие |
соответствен |
но 16 или 8 опытов. Поэтому с ростом числа факторов возрас-
т?°.т дробность |
применяемых |
реплик. |
|
|
|
|
||
?азберемся |
в выботзе |
Т |
реплпк на примере реплики 2 |
S-? |
||||
- |
. |
|||||||
|
|
|
4 |
• |
|
|
|
|
Здесь возможны |
12 р |
зекиъ, |
ec.^i х ^ приравнять |
парному |
вза- |
|||
№*оде;?ствию, а Х у - |
тройному. Допустим, |
выбран |
следующий |
|||||
вариант: х.^= У^хь.. |
х 5 |
= |
х,лхгхъ. |
Тогда |
определяющими |
|||
контрастами являются: I |
= х 1 х д х А |
и 1-= |
х 1 Х-г^ъ x s |
• |
Если перемножить эти определяющие контрасты, то получится
трсл>е |
соотношение I = х г |
х ^ х 5 . Чтобы полностью охарак |
||
теризовать разрешающую способность |
реплики, необходимо за |
|||
писать |
о б о б щ а ю щ и й |
о п р е д е л я ю щ и й |
||
к о н т р |
а с т ; |
|
|
|
|
1 |
— X^ Х^ X ^ *** X ^ X ^ X j * * |
X f - X j ^ j ^ ^ . |
Система смешивания определяется умножением обобщающе го определяющего контраста последовательно на х 1 , х г , 1 3 и т . д .
X , — o c s x 4 mt x , x t x A x s •» x a x , x t , |
|||
X^ |
X ^ X ^ X ^ X j i * * X ^ X g |
X^ Хд X y , |
|
X & — X , X ^ — X b X s X^ X j «• K , t t I r , |
|||
|
X , |
X J X J •«• X , X j |
X} X k Xg, |
X g " X t X 4 X ^ X ( — X t X 4 — X , X t X j , |
|||
X ^ j |
- |
X j X , X f c — X t X v X f — X , Xg, |
|
X , X S |
X j X 4 I j — X , X t X 4 » X t X j . |
||
Получается довольно сложная система омешиванвн лшнея- |
ных в&фектоъ с эффектами взахмодейоиии первого, второго, третьего и четвертого порядков. Вели, например, коэффици
енты in~~ |
p u + p w |
+ |
f u |
И |
|
. |
P« + |
f W |
P»» |
оадпвтеея о» нуля, |
|||
то возникает сомнения, можно ли пренебрегать другое |
вер |
|||||
ными взаимодействиями, |
с которыми |
ffuafm |
линейные эффек |
|||
ты. Тогда следует поставить вторую серию опытов, хворав |
||||||
нужным образом другую \ реплику. |
|
|
|
|
||
|
4 |
|
|
|
|
|
При этом можно воспользоваться |
м в * |
о д*о |
• • п е |
|||
р е в а л а " . |
Смысл этого метода заключается в |
той, |
что |
вторая четверть-реплика получается ив первой путем вше^- ненкя. всех знаков матрицы на обратные. Тогда в ооЪбяавцш определяющем контрасте тройные произведения имеют знак, противоположный их знаку в первой четверть-реплике. Трой ные произведения определяют парные взаимодействия в сов местных оценках для линейных эффектов. Усредняя результа ты обеих четверть-реплик, можно получить линейные аффек ты, не смешанные с парными въааюдеЗотвняия.