ПОИСК ГЛОБАЛЬНОГО ОПТИМУМА ДЛЯ ЗАДАЧ ОПТИМАЛЬНОГО ПРОЕКТИРОВАНИЯ СИСТЕМ ИЛИ ОПРЕДЕЛЕНИЯ ОПТИМАЛЬНЫХ ЗАКОНОВ УПРАВЛЕНИЯ.

РЕГРЕССИОННЫЙ АНАЛИЗ



При решении многих инженерных задач возникает необходимость в установлении связи между k аргументами x1 ,x2 , ..., xk и зависящей от них величиной y . Возможны следующие основные типы взаимосвязей между переменными величинами.
    1. Зависимость между неслучайными величинами. В этом случае зависимая переменная y вполне определенно задается независимыми переменными x1 ,x2 , ..., xk. Этот тип взаимосвязи представляет собой обычную функциональную зависимость между неслучайными величинами.
    2. Зависимость между случайными величинами. Между случайными величинами может существовать функциональная связь. Но между случайными величинами может существовать и стохастическая связь, которая проявляется в том, что одна из случайных величин реагирует на изменение другой изменениями своего закона распределения.
    3. Зависимость случайной величины y от неслучайных переменных x1 ,x2 , ..., xk. Природа этой взаимосвязи может быть двоякой: а) измерение зависимой переменной y связано с некоторыми ошибками измерения, а переменные x1 ,x2 , ..., xk измеряются без ошибок или эти ошибки пренебрежимо малы по сравнению с ошибкой измерения зависимой переменной; б) значения переменной y зависят не только от соответствующих значений x1 ,x2 , ..., xk но и от ряда неконтролируемых факторов, поэтому при каждом сочетании значений x1 ,x2 , ..., xk зависимая переменная y подвержена колебаниям случайного характера. Часто возникает необходимость в установлении связи между случайной величиной y и неслучайными переменными x1 ,x2 , ..., xk принимающими в каждой серии опытов определенные значения.
Величина y является случайной величиной, имеющей нормальное распределение, с переменным центром распределения М[у], изменяющимся при изменении сочетания значений факторов x1 ,x2 , ..., xk. Спучайная величина y имеет постоянную дисперсию , т. е. дисперсию, не зависящую от x1 ,x2 , ..., xk.
  Таким образом, М[у] является функцией x1 ,x2 , ..., xk, т. с на каждое изменение неслучайных величин x1 ,x2 , ..., xk случайная величина y реагирует изменением своего математического ожидания М[у]. Выражение М[у]=f(x1 ,x2 , ..., xk) называют уравнением регрессии математического ожидания случайной величины y по неслучайным величинам x1 ,x2 , ..., xk.
  Тип функции М[у]=f(x1 ,x2 , ..., xk) может быть линейным или криволинейным.
  В основе регрессионного анализа лежат следующие предположения.
    1. При каждом сочетании значений x1 ,x2 , ..., xk величина y имеет нормальное распределение.
    2. Дисперсия σ2 теоретического распределения случайной величины y постоянна.
    3. Тип функции М[у]=f(x1 ,x2 , ..., xk) известен.
    4. Независимые переменные x1 ,x2 , ..., xk измеряются с пренебрежимо малой ошибкой по сравнению с ошибкой в определении y.
    5. Переменные x1 ,x2 , ..., xk линейно независимы.
  При исследовании различных технологических процессов часто приходится иметь дело со случайной величиной y, не подчиняющейся нормальному распределению.
  В этом случае всегда можно подобрать такую функцию преобразования, которая позволит перейти от y к новой случайной величине q=f(y), распределенной приблизительно нормально. Многие асимметричные распределения удается аппроксимировать нормальным законом, перейдя от случайной величины y к случайной величине q=ln y.
  Если производить многократные повторные наблюдения над величиной y при принятом наборе значений x1 ,x2 , ..., xk, то согласно второму предположению получим дисперсию σ2 случайной величины y, которая не будет зависеть от значений x1 ,x2 , ..., xk. Практически это означает, что при многократных повторных наблюдениях над величиной у при первом, втором, третьем и т. д. наборах значений x1 ,x2 , ..., xk полученные выборочные дисперсии S12 ,S22 , ..., Sm2, являющиеся оценками дисперсии σ2 теоретического распределения, должны быть однородны. Однородность оценок дисперсии не всегда выполняется в условиях реального эксперимента. Если в этом случае удается найти функциональную зависимость , то оказывается возможным предложить такое преобразование случайной величины, которое позволит получить однородные оценки дисперсии.
  При обработке результатов многофакторного эксперимента функцию М [y] = f(xl, х2, ..., хk) обычно представляют полиномом вида

         
(1)

или полиномом второй степени

         
(2)

где
     , , - ... — коэффициенты уравнения регрессии.

Задачи регрессионного анализа следующие:
1) проверка гипотезы однородности дисперсий S2i
2) нахождение оценок b0 , b1, b2 .. ., bkk параметров , , ,...,.
3) определение доверительных интервалов для истинных зна чоний , , ,...,
4) проверка гипотезы адекватности.



НЕКОМПОЗИЦИОННЫЕ ПЛАНЫ И ИХ ПРИМЕНЕНИЕ


Некоторые некомпозиционные планы второго порядка

Во многих случаях целью исследования является получение математического описания изучаемого процесса. Часто из-за сложности процесса или малого объема информации неизвестную зависимость исследуемой величины &xi от k; независимых факторов представляют полиномом вида:

         
(1)

В этом случае необходимо определить коэффициенты уравнения и оценить их значимость. По результатам опытов можно определить только выборочные коэффициенты регрессии b0 , bi , bil , bii..., которые являются лишь оценками теоретических коэффициентов β0, βi, βil, βii... полинома (1). Уравнение регрессии, полученное по результатам опытов,имеет вид:

         
(2)

где y - выборочная оценка функции отклика ξ .

В ситуациях, когда априорная информация о порядке полинома отсутствует, математическую модель исследуемого процесса подбирают, начиная с простейшего линейного уравнения, последовательно увеличивая степень полинома до получения адекватной модели. Процесс получения математической модели в указанных ситуациях осуществляется следующим образом. Вначале реализуется полный факторный эксперимент 2k или эксперимент, представленный дробной репликой 2k-p , где p - число эффектов взаимодействия, замененных новыми переменными. По результатам опытов, выполненных согласно этим планам, находят коэффициенты линейного уравнения регрессии. Если это уравнение окажется неадекватным, то находят коэффициенты регрессии при эффектах взаимодействия факторов. Если уравнение регрессии с учетом взаимодействий факторов окажется также неадекватным, то выполненные ранее опыты дополняют опытами в «звездных» точках с плечом α и опытами в центре плана, число которых равно n0. Число опытов в «звездных» точках равно 2k . По результатам опытов, выполненных согласно плану 2k или 2k-p и дополнительным опытам в «звездных» точках и в центре плана, оценивают коэффициенты полинома второго порядка. Следует отметить, что исследуемый процесс часто удается описать полиномом второго порядка. В случае неадекватности полинома второго порядка переходят к планированию третьего порядка и описывают исследуемый процесс полиномом третьей степени. Если на основе априорной информации известно, что исследуемый процесс можно с достаточной точностью описать полиномом второго порядка, то для получения модели некомпозиционные планы в ряде случаев будут рациональнее центральных композиционных планов второго порядка.

При исследовании процессов и систем с двумя факторами рациональным является план типа правильного шестиугольника (рис. 1) с числом центральных точек n0>1;. На рисунке указаны номера опытов, а в скобках — уровни факторов в этих опытах. Обычно n0 принимают равным 4. В этом случае рассматриваемый план предусматривает проведение 10 опытов, из которых шесть выполняются при уровнях факторов, указанных в вершинах шестиугольника, и четыре опыта — при уровнях факторов, соответствующих центру плана. Этот план является ротатабельным и по числу опытов более экономичным, чем соответствующий ротата-бельный план второго порядка, требующий для своей реализации постановки 13 опытов: полный факторный эксперимент типа 22 содержит 4 опыта, в «звездных» точках необходимо провести 2k=2*2=4 опыта и 5 опытов в центре плана.
Другое достоинство описываемого плана состоит в том, что для фактора x1 он требует использования пяти уровней (+1; +0,5; 0; —0,5; —1), а для фактора x2 — всего трех уровней (+0,866; 0; —0,866). Центральный композиционный ротатабельный план второго порядка предусматривает использование обоих факторов на пяти уровнях (+1,414; +1; 0; —1; —1,414). На практике число уровней часто оказывается ограниченным, поэтому уменьшение числа уровней факторов представляет собой большое достоинство плана. Следует также отметить, что смена уровней факторов в процессе исследования усложняет и удорожает эксперимент.
Матрица плана, изображенного на рис. 1, представлена табл. 1. По результатам опытов, выполненных согласно этому плану, можно определить коэффициенты уравнения:

         
(3)

Таблица 1  Матрица планирования
Номер опыта x0 x1 x2 x1x2 x12 x22 y
1 +1 +1 0 0 +1 0 y1
2 +1 -1 0 0 +1 0 y2
3 +1 +0.5 +0.866 +0.433 +0.25 +0.75 y3
4 +1 +0.5 -0.866 -0.433 +0.25 +0.75 y4
5 +1 -0.5 +0.866 -0.433 +0.25 +0.75 y5
6 +1 -0.5 -0.866 +0.433 +0.25 +0.75 y6
7 +1 0 0 0 0 0 y7
8 +1 0 0 0 0 0 y8
9 +1 0 0 0 0 0 y9
10 +1 0 0 0 0 0 y10

Эти коэффициенты находят с помощью метода наименьших квадратов или по приведенным, ниже формулам:

         
;(4)

         
;(5)

         
; (6)

         
; (7)

         
; (8)

         
, (9)

где y0u — значение функции отклика в u-ом опыте в центре плана;
       x1j , x2j — кодированные значения факторов в j-ом опыте;
       yj — значения функции отклика в j-ом опыте.

Дисперсию Sy2 воспроизводимости эксперимента определяют по результатам опытов в центре плана, используя формулу:

         
, (10)

где — среднее арифметическое значение функции отклика, полученное по результатам n0 опытов в центре плана.



Дисперсии коэффициентов регрессии вычисляются по выражениям

         
; (11)


         
; (12)


         
; (13)


         
. (14)


Рассматриваемый план не является полностью ортогональным: коэффициенты b0 , b11 , b22 коррелированы. Ковариации, характеризующие статистические связи между коэффициентами b0 , b11 , b22 можно вычислить по выражениям:

         
; (15)


         
.(16)


Боксом и Бенкиным [2] для числа факторов от трех до семи рассмотрен ценный в практическом отношении класс некомпозиционных планов второго порядка. Эти планы представляют собой определенные выборки строк из полного факторного эксперимента типа 3k . В указанных планах каждая переменная варьируется всего на трех уровнях: +1, 0, —1, в то время как центральные композиционные ротатабельные планы второго порядка предусматривают использование каждого фактора на пяти уровнях. Смена уровней в процессе экспериментирования усложняет эксперимент и увеличивает его стоимость. Использование некомпозиционных планов, предусматривающих всего три уровня варьирования факторов, упрощает и удешевляет проведение эксперимента. Некомпозиционные планы характеризуются наличием в строках матрицы планирования большого числа нулей, в результате чего существенно упрощается вычисление коэффициентов модели. Кроме этого, некомпозиционные планы для 3, 4, 6 и 7 факторов требуют постановки меньшего числа опытов по сравнению с соответствующими ротатабельными центральными композиционными планами второго порядка. Числа опытов, предусмотренные некомпозиционными планами и центральными композиционными ротатабельными планами второго порядка, приведены ниже.

Число факторов 3 4 5 6 7
Некомпозиционный план (выборка из плана типа 3k) 15 27 46 54 62
Центральный композиционный ротатабельный план второго порядка 20 31 32 58 92

Схема некомпозиционного плана второго порядка для трех факторов изображена на рис.2. Этот план предусматривает проведение 15 опытов.
Номера опытов указаны на схеме. В центре плана, т. е. при нахождении всех факторов на нулевых уровнях (x1= x2= x3= 0 ), предусмотрено проведение трех опытов (рис. 2, опыты 5; 10; 15). Матрица плана, изображенного на рис. 2, представлена в табл. 2.

Таблица 2 Матрица некомпозиционного плана второго порядка для трех факторов (выборка из плана эксперемента 33)
Номер опыта x0 x1 x2 x3 x1x2 x1x3 x2x3 x12 x22 x32 y
1 +1 +1 +1 0 +1 0 0 +1 +1 0 y1
2 +1 +1 -1 0 -1 0 0 +1 +1 0 y2
3 +1 -1 +1 0 -1 0 0 +1 +1 0 y3
4 +1 -1 -1 0 +1 0 0 +1 +1 0 y4
5 +1 0 0 0 0 0 0 0 0 0 y5
6 +1 +1 0 +1 0 +1 0 +1 0 +1 y6
7 +1 +1 0 -1 0 -1 0 +1 0 +1 y7
8 +1 -1 0 +1 0 -1 0 +1 0 +1 y8
9 +1 -1 0 -1 0 +1 0 +1 0 +1 y9
10 +1 0 0 0 0 0 0 0 0 0 y10
11 +1 0 +1 +1 0 0 +1 0 +1 +1 y11
12 +1 0 +1 -1 0 0 -1 0 +1 +1 y12
13 +1 0 -1 +1 0 0 -1 0 +1 +1 y13
14 +1 0 -1 -1 0 0 +1 0 +1 +1 y14
15 +1 0 0 0 0 0 0 0 0 0 y15

Для четырех факторов матрица планирования приведена в табл. 3, для пяти факторов — в табл. 4, для шести — в табл. 6, для семи — в табл. 7.

Таблица 3 Матрица некомпозиционног плана второго порядка для четырех факторов
(выборка из плана эксперимента 34)
Номер опыта x0 x1 x2 x3 x4
1
2
3
4
5
6
7
8
9
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
-1
-1
0
0
0
0
0
+1
-1
+1
-1
0
0
0
0
0
0
0
0
0
+1
+1
-1
-1
0
0
0
0
0
+1
-1
+1
-1
0
10
11
12
13
14
15
16
17
18
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
-1
-1
0
0
0
0
0
0
0
0
0
+1
+1
-1
-1
0
0
0
0
0
+1
-1
+1
-1
0
+1
-1
+1
-1
0
0
0
0
0
19
20
21
22
23
24
25
26
27
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
-1
-1
0
0
0
0
0
0
0
0
0
+1
+1
-1
-1
0
+1
-1
+1
-1
0
0
0
0
0
0
0
0
0
+1
-1
+1
-1
0

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











         
(17)

Коэффициенты уравнения (17) могут быть определены с помощью метода наименьших квадратов. Боксом и Бенкиным [2] для определения коэффициентов этого уравнения получены формулы, которые для трех, четырех, пяти и семи факторов имеют вид:

         
; (18)


         
; (19)


         
; (20)


         
, (21)


где n0 — число опытов в центре плана;
       u — номер параллельного опыта в центре плана;
       y0u — значение функции отклика в u-ом опыте;
       N — число опытов в матрице планирования; j — номер опыта в матрице планирования;
       i , l — номера факторов;
       xijxlj — кодированные значения i-го и l-го факторов в j-м опыте;
       yj — значение функции отклика в j-м опыте;
       k — число факторов;
       A , B , C , D , p — константы, зависящие от числа факторов.

Таблица 4 Матрица некомпозиционногo плана второго порядка для пяти факторов
(выборка из плана эксперимента 35)
Номер опыта x0 x1 x2 x3 x4 x5 Номер опыта x0 x1 x2 x3 x4 x5
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
-1
-1
0
0
0
0
0
0
0
0
+1
+1
-1
-1
0
0
0
0
0
0
0
+1
-1
+1
-1
0
0
0
0
+1
+1
-1
-1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
+1
+1
-1
-1
0
0
0
0
+1
-1
+1
-1
0
0
0
0
0
0
0
0
0
0
0
+1
-1
+1
-1
0
0
0
0
0
0
0
0
+1
+1
-1
-1
0
0
0
0
0
0
0
0
0
0
0
+1
-1
+1
-1
0
0
0
0
+1
-1
+1
-1
0
0
0
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
0
0
0
0
+1
+1
-1
-1
0
0
0
0
+1
+1
-1
-1
0
0
0
0
0
0
0
+1
+1
-1
-1
0
0
0
0
0
0
0
0
0
0
0
0
+1
+1
-1
-1
0
0
0
+1
-1
+1
-1
0
0
0
0
+1
+1
-1
-1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
+1
-1
+1
-1
0
0
0
0
0
0
0
0
+1
-1
+1
-1
0
0
0
0
0
0
0
0
0
0
0
+1
-1
+1
-1
+1
-1
+1
-1
0
0
0
0
0
0
0

Для трех, четырех, пяти и семи факторов дисперсии Sb02 , Sbi2 , Sbil2 , Sbii2 коэффициентов регрессии определяют по формулам:

         
; (22)



         
; (23)



         
; (24)



         
, (25)



где Sy2 —дисперсия воспроизводимости эксперимента;
       B1— константа, зависящая от числа факторов.

Значения констант A , B , C , D , p и n0 для 3, 4, 5 и 7 факторов приведены в табл. 5.

Таблица 5  Значения констант в формулах для вычисления коэффициентов регрессий и дисперсий
Число факторов A B B1 C D p n0
3 1/8 1/4 13/48 -1/16 1/4 2 3
4 1/12 1/8 3/16 -1/48 1/4 2 3
5 1/16 1/12 11/96 -1/96 1/4 2 6
7 1/24 1/16 4/54 -1/144 1/8 3 6



Таблица 6 Матрица некомпозиционногo плана второго порядка для шести факторов
(выборка из плана эксперимента 36)
Номер опыта x0 x1 x2 x3 x4 x5 x6 Номер опыта x0 x1 x2 x3 x4 x5 x6
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
-1
-1
-1
-1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
+1
+1
-1
-1
+1
+1
-1
-1
+1
+1
+1
+1
-1
-1
-1
-1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
+1
+1
-1
-1
+1
+1
-1
-1
+1
+1
+1
+1
-1
-1
-1
-1
0
0
0
+1
-1
+1
-1
+1
-1
+1
-1
0
0
0
0
0
0
0
0
+1
+1
-1
-1
+1
+1
-1
-1
0
0
0
0
0
0
0
0
0
0
0
+1
-1
+1
-1
+1
-1
+1
-1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
+1
-1
+1
-1
+1
-1
+1
-1
0
0
0
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
-1
-1
-1
-1
0
0
0
0
0
0
0
0
+1
+1
+1
+1
-1
-1
-1
-1
0
0
0
0
0
0
0
0
0
0
0
+1
+1
+1
+1
-1
-1
-1
-1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
+1
+1
-1
-1
+1
+1
-1
-1
0
0
0
+1
+1
-1
-1
+1
+1
-1
-1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
+1
-1
+1
-1
+1
-1
+1
-1
+1
+1
-1
-1
+1
+1
-1
-1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
+1
-1
+1
-1
+1
-1
+1
-1
+1
-1
+1
-1
+1
-1
+1
-1
0
0
0

При шести факторах формулы для вычисления коэффициентов регрессии и их дисперсий имеют вид:

         
; (26)



         
; (27)



         
      для b14 ; b25 ; b36;(28)


         
       для остальных bil;(29)


         
; (30)



         
; (31)



         
; (32)



         
; (33)



         
; (34)



         
; (35)



         
; (36)



         
; (37)



         
       для S2b14 ; S2b25 ; S2b36;(38)


         
       для остальных Sbil2;(39)


         
.(40)



Таблица 7 Матрица некомпозиционного плана второго порядка для семи факторов
(выборка из плана эксперимента 37)
Номер опыта x0 x1 x2 x3 x4 x5 x6 x7 Номер опыта x0 x1 x2 x3 x4 x5 x6 x7
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
0
0
0
0
0
0
0
0
+1
+1
+1
+1
-1
-1
-1
-1
0
0
0
0
0
0
0
0
+1
+1
+1
+1
-1
-1
-1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
+1
+1
+1
+1
-1
-1
-1
-1
+1
+1
-1
-1
+1
+1
-1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
+1
+1
+1
+1
-1
-1
-1
-1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
+1
-1
+1
-1
+1
-1
+1
+1
+1
-1
-1
+1
+1
-1
-1
0
0
0
0
0
0
0
0
+1
+1
-1
-1
+1
+1
-1
-1
0
0
0
0
0
0
0
+1
-1
+1
-1
+1
-1
+1
-1
+1
+1
-1
-1
+1
+1
-1
-1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
+1
-1
+1
-1
+1
-1
+1
-1
+1
-1
+1
-1
+1
-1
+1
-1
0
0
0
0
0
0
0
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
-1
0
0
0
0
0
0
0
0
+1
+1
+1
+1
-1
-1
-1
-1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
-1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
+1
+1
+1
+1
-1
-1
-1
-1
0
0
0
0
0
0
0
+1
+1
+1
+1
-1
-1
-1
-1
+1
+1
-1
-1
+1
+1
-1
-1
+1
+1
-1
-1
+1
+1
-1
-1
0
0
0
0
0
0
-1
+1
+1
-1
-1
+1
+1
-1
-1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
+1
-1
+1
-1
+1
-1
+1
-1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
+1
-1
+1
-1
+1
-1
+1
-1
0
0
0
0
0
0
0
+1
-1
+1
-1
+1
-1
+1
-1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0

Дисперсию Sy2 воспроизводимости эксперимента определяют по результатам n0 опытов в центре плана:

         
,(41)

где             ,        .

Адекватность полученной модели проверяют с помощью F - критерия Фишера табл. 8

         
.(42)

Дисперсию Sад2 адекватности вычисляют по формуле

         
,(43)

где SR —сумма квадратов отклонений эмпирических значений yj функции отклика от ее значений , вычисленных по модели, во всех точках      плана;
       f — число степеней свободы,  f = N - k' - (n0 - 1) ;
       k' — число коэффициентов аппроксимирующего полинома.

Если найденное значение критерия Fp меньше табличного при принятом уровне значимости и соответствующих числах степеней свободы, то гипотеза адекватности полученной модели принимается.
Рассмотренные Боксом и Бенкиным некомпозиционные планы для числа факторов от трех до семи имеют высокую степень ортогональности, а именно: только свободный член b0 и коэффициенты bii при квадратичных членах коррелированы друг с другом.
При четырех и семи факторах указанные планы являются ротатабельными, а при другом числе факторов эти планы являются почти ротатабельными.

Таблица 8 Значения F - критерия Фишера при 5%-ном уровне значимости
Число степеней свободы для меньшей дисперсии fE Число степеней свободы для большей дисперсии fi
1 2 3 4 5 6 12 24 Безкон
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
22
24
26
28
30
40
60
120
Безкон
164.4
18.5
10.1
7.7
6.6
6.0
5.5
5.3
5.1
5.0
4.8
4.8
4.7
4.6
4.5
4.5
4.5
4.4
4.4
4.4
4.3
4.3
4.2
4.2
4.2
4.1
4.0
3.9
3.8
199.5
19.2
9.6
6.9
5.8
5.1
4.7
4.5
4.3
4.1
4.0
3.9
3.8
3.7
3.7
3.6
3.6
3.6
3.5
3.5
3.4
3.4
3.4
3.3
3.3
3.2
3.2
3.1
3.0
215.7
19.2
9.3
6.6
5.4
4.8
4.4
4.1
3.9
3.7
3.6
3.5
3.4
3.3
3.3
3.2
3.2
3.2
3.1
3.1
3.1
3.0
3.0
3.0
2.9
2.9
2.8
2.7
2.6
224.6
19.3
9.1
6.4
5.2
4.5
4.1
3.8
3.6
3.5
3.4
3.3
3.2
3.1
3.1
3.0
3.0
2.9
2.9
2.9
2.8
2.8
2.7
2.7
2.7
2.6
2.5
2.5
2.4
230.2
19.3
9.0
6.3
5.1
4.4
4.0
3.7
3.5
3.3
3.2
3.1
3.0
3.0
2.9
2.9
2.8
2.8
2.7
2.7
2.7
2.6
2.6
2.6
2.5
2.5
2.4
2.3
2.2
234.0
19.3
8.9
6.2
5.0
4.3
3.9
3.6
3.4
3.2
3.1
3.0
2.9
2.9
2.8
2.7
2.7
2.7
2.6
2.6
2.6
2.5
2.5
2.4
2.4
2.3
2.3
2.2
2.1
244.9
19.4
8.7
5.9
4.7
4.0
3.6
3.3
3.1
2.9
2.8
2.7
2.6
2.5
2.5
2.4
2.4
2.3
2.3
2.3
2.2
2.2
2.2
2.1
2.1
2.0
1.9
1.8
1.8
249.0
19.4
8.6
5.8
4.5
3.8
3.4
3.1
2.9
2.7
2.6
2.5
2.4
2.3
2.3
2.2
2.2
2.1
2.1
2.1
2.0
2.0
2.0
1.9
1.9
1.8
1.7
1.6
1.5
254.3
19.5
8.5
5.6
4.4
3.7
3.2
2.9
2.7
2.5
2.4
2.3
2.2
2.1
2.1
2.0
2.0
1.9
1.9
1.8
1.8
1.7
1.7
1.7
1.6
1.5
1.4
1.3
1.0

Доверительный интервал Δb0 коэффициента b0 находится из выражения:

         
,(44)

где   t — табличное значение критерия Стьюдента табл. 9;
     .

Аналогично определяем доверительные интервалы коэффициентов bi , bil , bii:

         
;(45)


         
;(46)


         
.(47)


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

Таблица 9 Значение t - критерия Стьюдента при 5%-ном уровне значимости
Число степеней свободы Значения t-критерия Число степеней свободы Значения t-критерия Число степеней свободы Значения t-критерия
1
2
3
4
5
6
7
8
9
10
12.71
4.303
3.182
2.776
2.571
2.447
2.365
2.306
2.262
2.228
11
12
13
14
15
16
17
18
19
20
2.201
2.179
2.160
2.145
2.131
2.120
2.110
2.101
2.093
2.086
21
22
23
24
25
26
27
28
29
30
2.080
2.074
2.069
2.064
2.060
2.056
2.052
2.048
2.045
2.042