![](/user_photo/_userpic.png)
книги из ГПНТБ / Кориков А.М. Математические методы планирования эксперимента учеб. пособие
.pdf- г о -
•
ния процесса, например, существующей аппаратурой, техноло
гией, организацией. В реакторе, изготовленном из некоторого материала, температуру нельзя поднять выше температуры пла
вления этого материала.
Эксперимент обычно начинается в условиях, когда объ ект уже подвергался некоторым исследованиям. № можем ис
пользовать априорную информацию для получения представле
ния о параметре оптимизации, о факторах, |
о наилучших |
у с л о |
|
виях ведения процесса и характере |
поверхности отклика. |
||
Итак, выбор экспериментальной |
области |
факторного |
про |
странства связан с тщательным анализом априорной информа
ции. |
|
|
В ы б о р о с н о в н о г о |
у р о в н я |
* ' |
осуществляется в зависимости |
от поставленной |
задачи. |
•Если ищется оптимальный режим протекания процесса, то в ка честве основного (нулевого*) уровня принимаются наилучшие
условия, определенные из аналита априорной информации.
*Этим условиям соответствует комбинация (или несколько ком бинаций) уровней факторов. Каждая комбинация является мно гомерной точкой в факторном пространстве, Ее можно рассма
тривать как походную точку для. .построения плана эксперимен та. Построение плана эксперимента сводится к выбору зкспе- 'риме'нтальных точек, симметричных относительно нулевого уровня.
Блок-схема принятия решений при выборе основного уров ня приведена на рис.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 зависимость от числа факторов в полином неполного выс шего порядка кроме линейных членов входят выражения, харак теризующие эффекты взаимодействия. Например, при трех факто-