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

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

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

- г о -

ния процесса, например, существующей аппаратурой, техноло­

гией, организацией. В реакторе, изготовленном из некоторого материала, температуру нельзя поднять выше температуры пла­

вления этого материала.

Эксперимент обычно начинается в условиях, когда объ­ ект уже подвергался некоторым исследованиям. № можем ис ­

пользовать априорную информацию для получения представле­

ния о параметре оптимизации, о факторах,

о наилучших

у с л о ­

виях ведения процесса и характере

поверхности отклика.

Итак, выбор экспериментальной

области

факторного

про­

странства связан с тщательным анализом априорной информа­

ции.

 

 

В ы б о р о с н о в н о г о

у р о в н я

* '

осуществляется в зависимости

от поставленной

задачи.

Если ищется оптимальный режим протекания процесса, то в ка­ честве основного (нулевого*) уровня принимаются наилучшие

условия, определенные из аналита априорной информации.

*Этим условиям соответствует комбинация (или несколько ком­ бинаций) уровней факторов. Каждая комбинация является мно­ гомерной точкой в факторном пространстве, Ее можно рассма­

тривать как походную точку для. .построения плана эксперимен­ та. Построение плана эксперимента сводится к выбору зкспе- 'риме'нтальных точек, симметричных относительно нулевого уровня.

Блок-схема принятия решений при выборе основного уров­ ня приведена на рис.1.1.'

 

 

 

Выбор

 

 

 

 

 

 

 

 

основного

уровня

 

 

 

 

Известная

 

Известная

 

 

Известно

Известна

подобласть,!

 

наилучшая

 

 

раилучшая

 

 

 

несколько

 

 

 

 

точка и

 

 

в которой

процесс

точка

 

 

 

[Наилучших

 

область

 

 

протекает

достаточ­

 

 

 

определения]

 

 

точек

но хорошо

 

 

Точка л е ­

Точка лежит

Имеются специаль-| Ни одной из точек

 

жит

внутри

на границе

ные

соображения

нельзя отдать

пред­

 

области

по выбору одной

почтение

 

 

области

 

 

 

 

из

точек

 

 

 

 

 

 

 

 

 

 

 

 

 

Цулевой

 

 

 

X

Выбираете^

Выбирается]

Точка

прини-

Выбирается)

Выбирается;

Ставится

уровень

(маетоя

за

наилучшая

случайная

несколько]

центр под-1

случайная

переме­

основной jypo-f

точка

точка

планов

области

точка в

щается

 

 

 

 

 

для

 

подобласти

вень

внутрь

 

 

 

 

 

 

 

1рачных

 

 

 

 

 

области

 

 

 

 

 

 

 

 

 

 

 

течек

 

 

 

 

 

 

 

 

 

 

 

 

Рис. I . I

- г г -

%

После того как нулевой уровень выбран,, переходам к с л е ­ дующему пату - выбору интервалов изменения (варьирования) факторов.

В ы б о р

и н т е р в а л о в

и з м е н е н и я

Очевидно, что для получения линейной модели нужно каж­

дый фактор варьировать на двух уровнях,

тогда задача состо ­

ит в том, чтобы для каждого фактора выбрать два уровня, на

которых он будет варьироваться в эксперименте. Интервалом

изменения факторов называется некоторое число (свое для каждого фактора), прибавление которого к основному уровню

дает верхний,

а вычитание - нижний уровни фактора. Другими

*"

• —

 

словами, интервал изменения - это

расстояние, на координат­

ной ос"» между основным и верхним

(или нижним) уровнем.

Для упрощения записи

условий эксперимента и

обработки

экспериментальных данных

масштабы по осям выбираются так,

чтобы верхний уровень соответствовал + 1 , нижний

- I , а о с ­

новной - нулю. Для факторов с непрерывной областью опреде­ ления это всегда можно сделать с помощью преобразования

где х кодированное значение фактора; натуральное значение фактора; натуральное значение основного уровня; интервал изменения;

ь- номер фактора.

Для качественных факторов, имеющих два уровня, один

- 2 3

-

уровень обозначается + 1 , а другой

- I ; порядок уровней не

имеет значения.

 

На выбор интервалов изменения накладываются естествен­

ные ограничения сверху и снизу. Интервал изменения не мо­ нет быть меньше той ошибки, с которой экспериментатор фик­ сирует уровень фактора. С другой стороны, интервал измене­ ния не может быть настолько большим, чтобы верхний или ниж­ ний уровни оказались за пределами области определения. Вну­ три этих ограничений обычно еще остается значительная не­ определенность выбора, которая устраняется с помощью инту­

итивных решений.

 

 

§ 1-2. П о л н ы й

ф а к т о р н ы й

э к с п е ­

 

р и м е н т

 

Полным факторным экспериментом называется

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

в котором реализуются все возможные сочетания уровней фак­ торов.

Если число факторов равно К , а число уровней каждого фактора равно ж-, то имеем полный факторный эксперимент типа v n , \ Рассмотрим подробнее полный факторный эксперимент типа 2 .

Нетрудно написать все сочетания уровней в эксперимен­ те с двумя факторами. Напомним, что в планировании экспериыента используются кодированные значения факторов: +1 и - I (часто для простоты записи единицы опускаются). Усло­ вия эксперимента можно записать в виде таблицы, где стро- *ки соответствуют различным опытам, а столбцы - значениям

факторов .• Такие таблицы

называются

м а т р и ц а м и

п л а н и р о в а н и я

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

В дальнейшем для удобства математических выкладок ш

в матрицу планирования

введем столбец

фиктивной

переменной

х 0 , которая принимает

во всех опытах

значение

+ 1 . Матрицу

планирования называют также матрицей независимых перемен­

ных. "Будем обозначать ее

через X .

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

эксперимента 2^ приведена в табл.

I . I .

 

 

Таблица I . I

*опыта

I

+

I

-

I

-

I

 

2

+

I

+

I

-

I

 

3

+

I

-

I

+

I

 

4

+

I

+

I

+

I

I

 

 

 

 

I

 

 

То, что записано в этой таблице, можно изобразить г е о - -метрически. Найдем в области определения фан~~ров точку, соответствующую основному уровню, и проведем через нее но­ вые оси координат, параллельные осям натуральных значений факторов. Далее выберем масштабы по новым осям так, чтобы интервал изменения для каждого фактора равнялся единице. Тогда условия проведения опытов будут соответствовать вер - пгггяям квадрата, центром которого является основной уровень, а каждая сторона параллельна одной из осей координат и рав­ на двум интервалам (см . рис . 1 . 2) .

Рис.1.2

1

Номера вершин квадрата соответствуют номерам опытов в мат­ рице планирования. Площадь, ограниченная квадратом, назы­

вается о б л а с т ь ю э к с п е р и м е н т а . Иног­

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

значений и,.

Запись матрицы планирования, особенно для многих фак­

торов, громоздка. Для ее сокращения удобно ввести услов­ ные буквенные обозначения строк. Это делается следующим образом. Порядковый номер фактора ставится в соответствие

строчной букве латинского алфавита: JCt а , Л ^ -

"6

и

т.д. Если теперь для строки матрицы планирования

выписать

латинские буквы только для факторов, находящихся на верх­ них уровнях, то условия опыта будут заданы однозначно. Опыт со всеми факторами на нижних уровнях условимся обоз-

начать ( I ) . Матрица хеширования вместе с принятыми буквен­ ными обозначениями приведена Б табл . 1 . 2 .

' $аблица 1.2

JS опыта \

 

" Г

 

 

X .

{Буквенные

 

 

 

 

 

 

 

 

 

 

 

j обозначения

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j строк

 

+ I

. - 1

 

- 1

1

( D '

 

+

I

+

i

 

- 1

 

а

 

 

 

1>

 

+

I

-

1

 

+ 1

 

 

 

 

 

 

+

I

+ 1

I

+ 1

1

а/о

 

 

 

 

 

 

 

 

 

Вели для двух факторов все возможные комбинации уровней

легко найти прямым перебором

(или

просто

запомнить), то

:,

ростом числа факторов возникает необходимость в некотором

приеме построения матриц. Рассмотрим три приема, основанные на переходе от матриц меньшей размерности к матрицам большей размерности.

Л р в н й п р и е - м : записать исходный план для

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

внями фактора.

 

 

 

 

Пусть

исходный план записан в табл . 1 . 2 . Проверьте, что

ма­

трица планирования для трехфакторной задачи,

построенная

с

помощью первого

приема, будет иметь следующий

вид:

 

( I ) ;

a ; i ;

аЛ; С ;

а с - , 1>с; aic.

 

 

Рассмотрим

в т о р о й

п р и е м . Для

этого введем

 

- 27 -

правило перемножения столбцов матрицы. При построчном пе -

ремножекин двух столбцов матрицы произведение единиц с од­

ноименными знаками дает + 1 , а с разноименными - I . Восполь­

зовавшись этим правилом, получим для случая, который мы

рассматриваем, вектор-столбец ацэс^. Далее повторим еще раз

исходный план, а у столбца произведений знаки поменяем на

обратные. Вы можете убедиться, что таким образом построе­

на матрица в таб.1.3.

 

 

 

 

 

 

 

 

 

 

 

 

 

Таблица 1.3

» опыта;

 

 

.,

J _ . .

 

!

 

 

— !

 

 

i

¥

 

 

t

 

 

x 1

 

 

 

 

 

 

 

 

 

 

J

 

 

i

 

 

1

!

+ I

1

-

1

I

-

I

 

+

I

!

*

i

 

 

j

2

i

+ I

 

- 1

 

+ I

- I

 

 

\

Hi

3

i

+ I

i

 

+ 1

 

- I

 

- I

 

 

i

I

Ъ

4-

j

+ I

i

 

+ 1

 

+ I

i

+ I

i

 

5

j

+

I

 

-

1

 

-

I

-

I

 

j

 

 

i

 

 

6

j

+

I

 

-

1

 

+

I

+

J.

 

 

 

 

 

t

%

i

 

 

 

7

'•

+ I

 

+ 1

 

- I

i

+ I

 

 

i

 

i

 

 

 

 

I

+

I

 

+

I

 

+

I

i

-

I

 

 

 

 

;

4'

8

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T p

e

T

И Й

n

 

p и e м

основан

на

iгравиле

чередова-

кия знаков. В первом столбце знаки меняются поочередно, во

втором столбце они чередуются через два, в третьем - через

4, в четвертом - через 8 и т.д. по степеням двойки.

§ 1-3. С в о й с т в а

п о л н о г о

ф а к т о р ­

н о г о

э к с п е р и м е н т а - т и п а 2 ^

Здесь мы выясним, какими общими свойствами обладают ма­

трицы планирования полного факторного эксперимента типа 2\

-г в -

независдао от числа факторов. Говоря о свойствах матриц,

имеем в виду те из них, которые определяют качество мо­

дели. Ведь эксперимент и планируется для того, чтобы полу­

чить модель, обладающую некоторыми оптимальными свойстзаш. Два свойства следуют непосредственно из построения мат­

рицы» Первое из них - с и м м е т р и ч н о с т ь относи­

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

зом: алгебраическая

с у ю »

элементов

гэктор-столбпа

г'ддого

сактота

равна тлю,

или

 

 

 

 

 

Л

 

 

 

 

 

 

 

И * . . « О ,

 

 

( 1 Л )

 

j -

i

*

 

 

 

 

где I -

номер фактора,

1 =

1 , 2 , . . . ,

К. , X

- число

опытов.

Зторсе свойство

-

тал

называемое

у с л о в и е

к ор-

м z р о в к и: сумма квадратов элетлехдов

какого

стслб -

ца равна числу опытов,

или

 

 

 

Это следствие

того, что значения факторов в матрвде зада­

ются +1 и -

I .

Мы рассмотрели свойства отдельных столбцов матриц-

планирования. Теперь остановился на свойстве совокупное^;:

столбцов..

Сумма почленных произведений любых двух

вектор- .

столбцов

матрицы равна нулю,

ш

 

 

Е. a c - k j о , с/u,

t,u, = i , 2

k.

( i . 3 )

Это важное ово2":ство называется

о р т о г о н а л ь н о -

 

 

 

- 2 9

-

 

 

 

с т

ь ю

матрицы планирования. Отсюда планы полного

фактор­

ного

эксперимента называют

о р т

о - т о н а л ь н ы м и .

 

 

Последнее, четвертое свойство

называется р о т

а т

а

-

б е л ь

л о с т ь ю , т . е .

точки в

матрице планирования

под-

бкраются так, что точность предсказания значений параметра оптимизации одинакоьа на равных расстояниях от центра экс­

перимента л не зависит от направления. Доказательство этого

утверждения приведем нике в § 3 - 2 . Планы, для которых имеет

место свойство ротатабельности, называют

р о т а т а б е -

л ь

н ы м и .

 

 

 

В дальнейшем т убедимся,

что свойство

ортогональности

дает

возможность значительно

снизить вычислительные т р у с о ­

сти, Еозшжаюпяе при расчете коэффициентов регрессии, г на­ личие ротатгбельноста влечет равномерность распределения дисперсии предсказанных значений параметра оптимизации в изу­ чавши области йакторкого пространства. В этом смысле полный'

факторный эксперимент типа

2 является

оптимальный при по­

строении линейной модели.

 

 

 

§" 1-4. Э ф ф е к т

ф а к т о р а

и

э ф ф е к т ы

 

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

При варьировании факторами на двух уровнях можно полу­

чить математическую

модель

объекта в виде полинома первого

и неполного высшего

порядка.

 

 

3 зависимость от числа факторов в полином неполного выс­ шего порядка кроме линейных членов входят выражения, харак­ теризующие эффекты взаимодействия. Например, при трех факто-

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