Конспект (1062950), страница 7
Текст из файла (страница 7)
(55)
К оэффициенты bj определяются в соответствии с выражением:
(56)
А дисперсии коэффициентов равны
(57)
Для перехода к виду (51) необходимо пересчитать коэффициент b0:
(58)
При применении ЦКОП получение идентичной информации (т.е. с одинаковой точностью предсказания выходной величины) во всех направлениях исследуемого пространства невозможно, т.к. информационные поверхности не являются сферами.
Центральный композиционный рототабельный план (ЦКРП).
Более удачным является рототабельное планирование, при котором информационная поверхность приближается к сферической. Это достигается тем, что выбирая удаленные от центра плана «звездные» точки на осях координат, они дополняются информацией из центра плана, равноточной во всех направлениях. Удельный вес этой информации в общем объеме увеличивается, что достигается увеличением числа опытов в центре плана. Таким образом, в ЦКРП число опытов в центре плана зависит от числа учитываемых в эксперименте факторов. Это увеличивает количество опытов, но дает возможность получать равноточную информацию. Для сокращения количества опытов можно отказаться от постановки параллельных опытов для оценки воспроизводимости, которая в этом случае может быть оценена по экспериментам в центре плана.
Чтобы композиционный план был рототабельным, величина «звездного» плеча выбирается из условий:
при k<5
при k>=5. (59)
Подсчитанные значения звездного плеча и число точек в центре плана приведены в табл. 9.
Таблица 9
k | 2 | 3 | 4 | 5 | 6 | 7 |
| 1,414 | 1,682 | 2,00 | 2,00 | 2.38 | 2,83 |
mo | 5 | 6 | 7 | 8 | 9 | 14 |
Матрица планирования для ЦКРП приведена в таблице 10. Столбцы, соответствующие взаимодействию факторов, в ней отсутствуют. Из сравнения матрицы ЦКРП с матрицей для ПФЭ типа 23 видно, что значения, приведенные в этих столбцах до восьмого опыта, были бы идентичны. Начиная с девятого опыта значения, соответствующие взаимодействию факторов, равны нулю и не влияют на оценку значимости соответствующего коэффициента в полиноме вида (34). Поэтому приводить эти столбцы в матрице ЦКРП необязательно.
Таблица 10
Матрица центрального композиционного рототабельного плана
Номер опыта | Xo | X1 | X2 | X3 | X12 | X22 | X32 | Yj |
1 | +1 | —1 | —1 | —1 | +1 | +1 | +1 | Y1 |
2 | +1 | +1 | —1 | —1 | +1 | +1 | +1 | ... |
3 | +1 | —1 | +1 | —1 | +1 | +1 | +1 | ... |
4 | +1 | +1 | +1 | —1 | +1 | +1 | +1 | ... |
5 | +1 | —1 | —1 | +1 | +1 | +1 | +1 | ... |
.6 | +1 | +1 | —1 | +1 | +1 | +1 | +1 | ... |
7 | +1 | —1 | +1 | +1 | +1 | +1 | +1 | ... |
8 | +1 | +1 | +1 | +1 | +1 | +1 | +1 | ... |
9 | +1 | — | 0 | 0 | | 0 | 0 | ... |
10 | +1 | + | 0 | 0 | | 0 | 0 | ... |
11 | +1 | 0 | — | 0 | 0 | | 0 | ... |
12 | +1 | 0 | + | 0 | 0 | | 0 | ... |
13 | +1 | 0 | 0 | — | 0 | 0 | | ... |
14 | +1 | 0 | 0 | + | 0 | 0 | | ... |
15 | +1 | 0 | 0 | 0 | 0 | 0 | 0 | Y15 |
16 | +1 | 0 | 0 | 0 | 0 | 0 | 0 | ... |
17 | +1 | 0 | 0 | 0 | 0 | 0 | 0 | ... |
18 | +1 | 0 | 0 | 0 | 0 | 0 | 0 | ... |
19 | +1 | 0 | 0 | 0 | 0 | 0 | 0 | ... |
20 | +1 | 0 | 0 | 0 | 0 | 0 | 0 | Y20 |
Можно заметить, что матрица ЦКРП не соответствует условиям ортогональности для столбцов с квадратичными членами полинома. Поэтому оценка коэффициентов полинома 2-го порядка не будет являться независимой, т.е. коэффициенты bjj коррелированы между собой и со свободным членом b0. Но этот недостаток ЦКРП компенсируется более высокой точностью определения выходного параметра во всех направлениях на одинаковом расстоянии от центра.
Коэффициенты и их дисперсии определяются из условий:
(60)
(61)
(62)
(63)
(64)
(65)
(66)
(67)
Значения констант, входящих в приведенные выше выражения, приведены в таблице 11.
Таблица 11
Вычисление коэффициентов регрессии при ротатабельном
планировании для k<=7
Число независимых переменных k | | a1 | a2 | a3 | a4 | a5 | a6 | a7 |
2 | 1,412 | 0,2 | 0,1 | 0,125 | 0,25 | 0,1251 | 0,0187 | 0,1 |
3 | 1,682 | 0,1663 | 0,0568 | 0,0732 | 0 125 | 0,0625 | 0,0069 | 0,0568 |
4 | 2,00 | 0,1428 | 0,0357 | 0,0417 | 0,0625 | 0,0312 | 0,0037 | 0,0357 |
5 | 2,378 | 0,0988 | 0,0191 | 0,0231 | 0,0312 | 0,0156 | 0,0015 | 0,0191 |
6 | 2,828 | 0,0725 | 0,0098 | 0,0125 | 0,0156 | 0,0078 | 0,0005 | 0,0098 |
7 | 3,333 | 0., 0398 | 0,0052 | 0,0066 | 0,0078 | 0,0039 | 0,0002 | 0,0052 |
Дисперсию воспроизводимости, как было сказано выше, определяют по опытам в центре плана. Дисперсия адекватности определяется по выражению, включающему в себя величины т.н. остаточной дисперсии и дисперсии воспроизводимости и их степени свободы: