Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Кориков А.М. Математические методы планирования эксперимента учеб. пособие

.pdf
Скачиваний:
21
Добавлен:
23.10.2023
Размер:
7.9 Mб
Скачать

•пах может быть три двойных и один

тройной эшфект Ь „ х , х ,

Поясним понятая "эффект фактора" и "эффект эзаикодей-

с т в ж " . Вклад фактора в величину

параметра оптимизации при

переходе от нижнего уровня к верхнему уровню

называется

э ф ф е к т о м

ф а к т о р а

(иногда его

называют л и ­

нейным,основным или главным эффектом). Он численно равен

удвоенное коэффициенту полинома. Для качественных факто­

ров, изменяемых на двух уровнях, основной

уровень

не име­

ет физического сшсла . Поэтому понятие "эффект фактора"

является здесь естественным. Если еффект одного фактора

зависит от уровня,

на котором находится другой фактор, тс

говорят,

что имеет

место

э ф ф е к т

в з а и м о д е й ­

с т в и я

двух факторов

(иногда

говорят

двойной

или пар­

ный эффект взаимодействия

(эффект

взаимодействия

первого

порядка). Также выясняется смысл тройного эффекта взаимо­ действия (эффект взаимодействия второго порядка) и т . д . Численно эффект взаимодействия равен удвоенному соответст­ вующему коэффициенту полинома.

§ 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

И

 

.

+

f W

P»»

оадпвтеея о» нуля,

то возникает сомнения, можно ли пренебрегать другое

вер­

ными взаимодействиями,

с которыми

ffuafm

линейные эффек­

ты. Тогда следует поставить вторую серию опытов, хворав

нужным образом другую \ реплику.

 

 

 

 

 

4

 

 

 

 

 

При этом можно воспользоваться

м в *

о д*о

• • п е ­

р е в а л а " .

Смысл этого метода заключается в

той,

что

вторая четверть-реплика получается ив первой путем вше^- ненкя. всех знаков матрицы на обратные. Тогда в ооЪбяавцш определяющем контрасте тройные произведения имеют знак, противоположный их знаку в первой четверть-реплике. Трой­ ные произведения определяют парные взаимодействия в сов­ местных оценках для линейных эффектов. Усредняя результа­ ты обеих четверть-реплик, можно получить линейные аффек­ ты, не смешанные с парными въааюдеЗотвняия.

Соседние файлы в папке книги из ГПНТБ