
- •ВВЕДЕНИЕ
- •1. ЭКСПЕРИМЕНТ КАК ПРЕДМЕТ ИССЛЕДОВАНИЯ
- •1.1. Понятие эксперимента
- •1.2. Виды экспериментальных исследований
- •1.3. Математическая модель объекта исследования
- •1.4. Факторы
- •1.5. Параметр оптимизации
- •Вопросы и задания для самоподготовки
- •2. ПРОСТЫЕ СРАВНИВАЮЩИЕ ЭКСПЕРИМЕНТЫ
- •2.1. Предварительная обработка экспериментальных данных
- •2.2. Статистические гипотезы
- •2.5. Сравнение двух рядов наблюдений
- •2.6. Сравнение двух дисперсий
- •2.7. Проверка гипотезы о законе распределения
- •2.7.1. Общие сведения
- •2.7.2. Пример проверки гипотезы о нормальном законе распределения экспериментальных данных
- •Вопросы и задания для самоподготовки
- •3. АНАЛИЗ РЕЗУЛЬТАТОВ ПАССИВНОГО ЭКСПЕРИМЕНТА
- •3.1. Элементы дисперсионного анализа
- •3.1.1. Общие сведения
- •3.1.2. Пример применения однофакторного дисперсионного анализа
- •3.2. Характеристика видов связей между рядами наблюдений
- •3.4. Определение тесноты связи между случайными величинами
- •3.5. Парная линейная корреляция
- •3.6. Статистическое изучение корреляционной связи
- •3.6.2. Исключение из массива первичной информации промахов
- •3.6.4. Измерение степени тесноты связи, оценка ее существенности
- •3.8. Линейная множественная регрессия
- •Вопросы и задания для самоподготовки
- •4. МНОГОФАКТОРНЫЕ ЭКСПЕРИМЕНТЫ
- •4.1. Общие сведения
- •4.2. Полный факторный эксперимент
- •4.2.1. Кодирование факторов
- •4.2.2. Матрицы планирования эксперимента
- •4.2.3. Рандомизация опытов
- •4.2.4. Проведение эксперимента
- •4.2.5. Проверка однородности дисперсии параллельных опытов, воспроизводимости эксперимента
- •4.2.6. Расчет коэффициентов регрессии, проверка их значимости
- •4.2.7. Проверка адекватности модели
- •4.2.8. Пример применения планов первого порядка полного факторного эксперимента
- •4.3. Дробный факторный эксперимент
- •4.4. Применение плана первого порядка дробного факторного эксперимента для исследования технических устройств
- •4.4.1. Теоретическая модель усилителя
- •4.5. Планы второго порядка
- •4.5.1. Ортогональные планы
- •4.5.2. Ротатабельные планы второго порядка
- •5.1. Метод покоординатной оптимизации
- •5.3. Симплексный метод планирования
- •ПРИЛОЖЕНИЯ

где j – номер столбца в матрице планирования; i – номер строки. Дисперсии коэффициентов уравнения регрессии следующие:
|
Sbj2 Sвосп2 |
n |
Sbjj/ 2 |
Sвосп2 |
n |
|
/ji2 , |
|
|||||
|
x2ji , |
x |
|
||||||||||
|
|
|
|
i 1 |
|
|
|
|
i 1 |
|
|
|
|
2 |
2 |
n |
|
2 |
2 |
|
2 |
|
k |
2 |
2 |
(4.67) |
|
|
|
|
|
|
|
||||||||
Sbju |
Sвосп |
xjixui , |
Sb0 |
Sb/ 0 |
xjjSb/ jj . |
||||||||
|
|
i 1 |
|
|
|
|
|
|
j 1 |
|
|
|
Коэффициенты уравнения регрессии, получаемые с помощью ортогональных планов второго порядка, определяются с разной точностью, ортогональные планы первого порядка обеспечивают одинаковую точность коэффициентов. Проверяют значимость коэффициентов по t –критерию и адекватность уравнения по критерию Фишера.
сии, постоянна для всех точек, находящихсяИна равном расстоянии от центра эксперимента. Экспериментатору заранее неизвестно, где на-
4.5.2. Ротатабельные планы второго порядка
Ротатабельным называют планирование, для которого диспер-
ходится та часть поверхности откликаД, которая представляет для него особый интерес, поэтому следует стремиться к тому, чтобы количест-
сия параметра оптимизации ~y , предсказанного уравнением регрес-
во для всех равноотстоящихботАцентра эксперимента точек. Действительно, удаление от центра точек 5-8 в 2 раз меньше, чем удаление точек 1-4 и, следовательно, коэффициенты уравнения регрессии оп-
во информации, содержащееся в уравнении регрессии, было одинако-
ределяются с разл чнойид сперсией. Бокс и Хантер предложили ротатабельные планы 2-го порядка. Для того чтобы композиционный план
был ротатабельным, величину звёздного плеча выбирают из условия |
|||||||
С |
|
k 1 |
|
|
|||
|
k |
|
|
|
|
||
2 |
4 |
при k 5 |
и 2 |
4 |
|
при k 5. |
(4.68) |
|
|
|
Значения звёздных плеч и числа точек в центре ротатабельных планов приведены в табл. 4.22.
Матрица планирования для ротатабельного плана представлена табл. 4.23.
88

Таблица 4.22
Значения звёздных плеч в ротатабельных планах второго порядка
Параметр |
|
|
Значения параметров при числе независимых факторов |
|||||||||||||||||||
плана |
|
|
2 |
|
|
3 |
|
4 |
|
5 |
5 |
|
6 |
|
|
6 |
|
7 |
7 |
|||
Ядро плана |
22 |
|
|
23 |
|
24 |
|
25 |
25-1 |
|
26 |
|
|
26-1 |
|
27 |
27-1 |
|||||
Звёздное |
1,414 |
|
1,682 |
2,000 |
|
2,378 |
2,000 |
2,828 |
|
2,378 |
|
3,333 |
2,828 |
|||||||||
плечо |
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Число |
точек |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
в |
центре |
5 |
|
|
6 |
|
7 |
|
10 |
6 |
|
15 |
|
|
9 |
|
21 |
14 |
||||
плана n0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
Ротатабельный план второго порядка |
|
|
|
|
Таблица 4.23 |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Номер опыта |
|
|
|
|
|
|
|
|
Факторы |
|
|
|
|
|
|
|
|
Результат |
||||
|
|
|
|
|
|
X0 |
|
X1 |
|
X2 |
X1X2 |
|
X12 |
|
|
X22 |
|
|
yi |
|||
Ядро плана |
|
1 |
|
|
+1 |
|
-1 |
-1 |
+1 |
|
+1 |
|
|
|
+1 |
|
|
y1 |
||||
|
|
|
2 |
|
|
+1 |
|
+1 |
-1 |
-1 |
|
+1 |
|
|
|
+1 |
|
|
y2 |
|||
|
|
|
3 |
|
|
+1 |
|
-1 |
+1 |
-1 |
|
+1 |
|
|
|
+1 |
|
|
y3 |
|||
|
|
|
|
|
|
|
|
|
|
И |
|
|
|
|
|
|||||||
|
|
|
4 |
|
|
+1 |
|
+1 |
+1 |
+1 +1 |
|
+1 |
|
|
y4 |
|||||||
Звёздные |
|
5 |
|
|
+1 |
+1,414 |
0 |
0 |
|
+2 |
|
|
|
0 |
|
|
|
y5 |
||||
точки |
|
|
6 |
|
|
+1 |
-1,414 |
0 |
0 |
|
+2 |
|
|
|
0 |
|
|
|
y6 |
|||
|
|
|
|
|
|
|
|
|
|
|
|
Д |
|
|
|
|
|
|
|
|
|
|
|
|
|
7 |
|
|
+1 |
|
0 |
+1,414 |
0 |
|
0 |
|
|
|
+2 |
|
|
y7 |
|||
|
|
|
8 |
|
|
+1 |
|
0 |
-1,414 |
0 |
|
0 |
|
|
|
+2 |
|
|
y8 |
|||
|
|
|
|
|
|
|
|
|
А |
|
|
|
|
|
|
|
|
|
|
|
||
Центр |
|
|
9 |
|
|
+1 |
|
0 |
0 |
0 |
|
0 |
|
|
|
0 |
|
|
|
y9 |
||
плана |
|
|
10 |
|
|
+1 |
|
0 |
0 |
0 |
|
0 |
|
|
|
0 |
|
|
|
y10 |
||
|
|
|
|
|
|
|
|
б |
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
11 |
|
|
+1 |
|
0 0 |
0 |
|
0 |
|
|
|
0 |
|
|
|
y11 |
|||
|
|
|
12 |
|
|
+1 |
|
0 |
0 |
0 |
|
0 |
|
|
|
0 |
|
|
|
y12 |
||
|
|
|
|
|
|
и |
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
13 |
|
|
+1 0 |
0 |
0 |
|
0 |
|
|
|
0 |
|
|
|
y13 |
УчитываяСспецифический характер ротатабельного плана в общем виде, можно получить формулы для расчёта коэффициентов уравнения регрессии и их дисперсий:
|
|
|
|
|
|
|
|
A |
|
|
k |
|
|
|
|
|
|
|
|
b |
|
|
|
2 2 k |
2 0y 2 c jjy ; |
||
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
0 |
|
n |
|
|
j 1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
bj c/n jy ; |
|
||||
|
|
|
A |
|
k 2 k jjy c |
k |
|
||||||
b |
jj |
|
|
|
c2 |
2 1 jjy 2 c 0y ; |
|||||||
|
|||||||||||||
|
|
|
n |
|
|
|
|
|
|
|
j 1 |
|
|
|
|
|
|
|
|
|
|
bju |
c2 |
juy ; |
|
||
|
|
|
|
|
|
|
|
n |
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
(4.69)
(4.70)
(4.71)
(4.72)
89

|
|
|
|
|
|
|
|
2 |
|
|
2A 2 |
k 2 2 |
|
|
|
|||||||||
|
|
|
|
|
|
|
|
Sb0 |
|
|
|
|
|
|
|
|
|
Sвосп |
; |
|
(4.73) |
|||
|
|
|
|
|
|
|
|
|
n |
|
|
|
|
|||||||||||
|
|
|
|
|
|
2 |
|
|
A k 1 k 1 c2 |
2 |
|
|
||||||||||||
|
|
|
|
|
|
Sbjj |
|
|
|
|
|
|
|
|
|
|
|
|
|
Sвосп ; |
(4.74) |
|||
|
|
|
|
|
|
|
|
|
n |
|
|
c2 |
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
2 |
|
|
2 |
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
Sbju |
|
|
|
|
Sвосп , |
|
|
(4.75) |
|||||
|
|
|
|
|
|
|
|
|
|
|
n |
|
|
|||||||||||
|
n |
|
|
|
|
|
|
|
n |
|
|
|
|
|
|
n |
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
jy |
|
||||||||
где 0y X0i yi ; |
|
juy X ji Xui yi ; |
|
|
X ji yi |
; |
||||||||||||||||||
|
i 1 |
|
|
|
|
|
|
i 1 |
|
|
|
|
|
|
|
|
|
i 1 |
|
|||||
|
n |
|
|
|
|
|
|
n |
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
jjy X 2 y |
|
; |
|
c |
|
|
; |
|
|
A |
|
|
; |
|||||||||||
i |
|
n |
|
|
|
2 k 2 k |
||||||||||||||||||
|
ji |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
i 1 |
|
|
|
|
|
X 2ji |
|
|
|
|
|
|
|
||||||||||
|
nk |
|
|
|
|
|
|
i 1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
; |
|
|
|
n |
n n . |
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
k 2 n1 |
|
|
|
1 |
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
Матрица ротатабельного планирования оказывается неортого- |
|||||||||||||||||||||||
нальной, так как |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
n |
|
X 2 |
0; |
|
|
n |
|
X |
|
|
0; |
j u. |
(4.76) |
|||||||
|
|
|
|
X |
0i |
|
|
X 2 |
|
2 |
|
|||||||||||||
|
|
|
|
i 1 |
ui |
|
|
|
|
|
ji |
|
|
ui |
|
|
И |
|
||||||
|
|
|
|
|
|
|
|
|
|
i 1 |
|
|
|
|
|
|
|
|
|
|
|
|
||
|
Следовательно, если какой-либо из квадратичных эффектов ока- |
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
Д |
|
|
|
||||||||
зался незначимым, то после его исключения коэффициенты уравне- |
||||||||||||||||||||||||
ния регрессии необходимо пересчитать заново. |
|
|
|
|||||||||||||||||||||
|
При использован |
|
ротатаАельных планов второго порядка дис- |
персию воспроизводимостибможно определить по опытам в центре плана. В связи с эт м при проверке адекватности уравнения регрессии, полученного по ротатабельному плану второго порядка, поступают следующим образом:
находят остаточную сумму квадратов
2 |
|
|
n |
~ |
2 |
|
|
(4.77) |
|
С S1 |
|
yi yi |
|
|
|
||||
|
|
i 1 |
|
k 2 k 1 |
|
|
|||
с числом степеней свободы d.f |
1 |
n |
; |
|
|||||
|
|
||||||||
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
по опытам в центре плана определяют дисперсию воспроиз- |
|||||||||
водимости |
|
n0 |
|
|
|
|
|
|
|
2 |
|
~ |
2 |
|
(4.78) |
||||
S2 |
|
y0i y0i |
|
i1
счислом степеней свободы d.f2 n0 1;
находят сумму квадратов, характеризующих неадекватность
90

S2 |
S2 |
S2 |
(4.79) |
3 |
1 |
2 |
|
счислом степеней свободы d.f3 d.f1 d.f2 ;
проверяют адекватность по F-критерию:
S2 |
/d.f |
3 |
|
|
|
F |
3 |
|
, |
(4.80) |
|
|
/d.f2 |
||||
S22 |
|
|
уравнение адекватно, если F FТ .
Если модель неадекватна, следует повторить эксперименты на меньшем интервале варьирования факторов или перенести центр плана в другую точку факторного пространства. В тех случаях, когда адекватность модели по-прежнему не достигается, рекомендуется перейти к планам третьего порядка.
1.Как зависит число опытов отДвида Ипринимаемой математической модели?
2.Чем можно объяснитьАширокое распространение полиномиальных моделей?
3.Дайте определениебполного факторного эксперимента.
4.Что характеризуют -коэффициенты?
5.Перечислитеиэтапы планирования и реализации полного факторного эксперимента.
6.Что называютСкод рованием факторов? Зачем его проводят?
7.Геометрическое представление планов типа 2k . Постройте
9.Перечислите свойства матрицы планирования полного факторного эксперимента.
10.Что называют рандомизацией опытов? Зачем ее проводят?
11.Какие опыты называют параллельными?
12.Как и для чего проводится проверка однородности дисперсии параллельных опытов?
13.Что означает понятие воспроизводимости эксперимента?
14.Как оценить ошибку эксперимента?
15.Какой метод применяется при расчете коэффициентов уравнения регрессии? Запишите формулу расчета b-коэффициентов.
91

16. Что называют взаимодействием факторов и как оно учитывается при планировании полного факторного эксперимента?
17. Что называют взаимодействием первого, второго, третьего и т.д. порядка? Как определяется число возможных взаимодействий
факторов? |
|
|
|
18. |
Какиесуществуютспособыпроверкизначимостиb-коэффициентов? |
||
19. |
Чем может быть обусловлена незначимость коэффициентов |
||
уравнения регрессии? |
|
|
|
20. |
Как и для чего проводится проверка адекватности уравнения |
||
регрессии? |
|
|
|
21. |
Что называют дробным факторным экспериментом? |
||
22. |
Дайте определение дробной реплики полного факторного |
||
эксперимента. |
|
|
|
23. |
Порядок планирования дробного факторного эксперимента. |
||
24. |
Какие планы называют насыщенными? |
||
25. |
Явление смешивания оценок -коэффициентов в дробном |
||
26. |
Что называют генерирующим соотношением и определяю- |
||
щим контрастом? |
|
И |
|
|
|
||
27. |
Ортогональные и ротатабельные планы второго порядка. |
||
28. |
|
|
Д |
Определение звёздных плеч и количества опытов в центре |
|||
планов второго порядка. |
|
|
|
29. |
Определен е коэфф циентовАуравнения регрессии и провер- |
||
ка их значимости в ортогональных и ротатабельных планах. |
|||
30. |
|
б |
|
Оценка адекватности модели, построенной с помощью пла- |
|||
нов второго порядка. |
|
|
|
|
и |
|
|
5. ПЛАНИРОВАНИЕ ЭКСПЕРИМЕНТОВ ПРИ ПОИСКЕ |
|||
|
ОПТИМАЛЬНЫХ УСЛОВИЙ |
||
|
С |
|
|
Во многих случаях инженерной практики перед исследователем возникает задача не только выявления характера связи между двумя или несколькими рядами наблюдений, но и нахождения таких численных значений факторов, при которых отклик достигает своего экстремального значения (максимума или минимума). Эксперимент, решающий эту задачу, называется экстремальным. В этом случае задача сводится к оптимизационной и формулируется следующим образом: требуется определить такие координаты экстремальной точки
92