|
Раздел "Математика\Statistics Toolbox"
Список функций Statistics Toolbox
\ \ Функции планирования эксперимента
|
|
Генерация матрицы центрального композиционного плана |
Синтаксис
D = ccdesign(nfactors)
D = ccdesign(nfactors,'pname1',pvalue1,'pname2',pvalue2,...)
[D,blk] = ccdesign(...)
Описание
D = ccdesign(nfactors) функция предназначена для генерации матрицы D центрального композиционного плана для nfactors факторов. Размерность матрицы плана эксперимента D равна nfactors x nfactors. Количество опытов в матрице равно числу факторов. Значения факторов нормализованы таким образом, чтобы вершины куба плана были равны -1 и 1.
[D,blk] = ccdesign(nfactors) функция возвращает вектор номеров блоков blk опытов в плане эксперимента. Блоки позволяют сгруппировать опыты таким образом, чтобы минимизировать влияние изменений внешних условий на значения оцениваемых параметров.
[...] = ccdesign(nfactors,'pname1',pvalue1,'pname2',pvalue2,...) такой вариант синтаксиса позволяет задать следующие дополнительные параметры 'pname' и их значения pvalue:
| 'center' |
Число центральных точек в плане: |
| Целое число |
Определяет число центральных точек в плане |
| 'uniform' |
Число центральных точек выбирается для обеспечения равной точности в области проведения эксперимента |
| 'orthogonal' |
Число центральных точек выбирается для обеспечения ортогональности плана. Принимается как значение по умолчанию. |
| 'fraction' |
Вид дробной реплики от полного факторного плана. Задается в виде степени 1/2. Например: |
| 0 |
Полный план |
| 1 |
1/2 реплика |
| 2 |
1/4 реплика |
| 'type' |
Вид центрального композиционного плана. Возможные значения 'inscribed', 'circumscribed', 'faced'. |
| 'blocksize' |
Максимальное количество точек в блоке плана |
Дополнительные аргументы задаются в виде пары: 'название параметра', значение параметра.
Примеры использования функции генерации матрицы центрального композиционного плана
Генерация матрицы центрального композиционного плана для 3 факторов
>> D = ccdesign(3)
D =
-1.0000 -1.0000 -1.0000
-1.0000 -1.0000 1.0000
-1.0000 1.0000 -1.0000
-1.0000 1.0000 1.0000
1.0000 -1.0000 -1.0000
1.0000 -1.0000 1.0000
1.0000 1.0000 -1.0000
1.0000 1.0000 1.0000
-1.6818 0 0
1.6818 0 0
0 -1.6818 0
0 1.6818 0
0 0 -1.6818
0 0 1.6818
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
Генерация матрицы центрального композиционного плана для 5 факторов разделенного на блоки
>> [D,blk] = ccdesign(5)
| D= |
blk= |
-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
-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
1 1 1 -1 -1
1 1 1 1 1
-2 0 0 0 0
2 0 0 0 0
0 -2 0 0 0
0 2 0 0 0
0 0 -2 0 0
0 0 2 0 0
0 0 0 -2 0
0 0 0 2 0
0 0 0 0 -2
0 0 0 0 2
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 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
2
2
2
2
2
2
2
2
2
2
1
1
1
1
1
2
2
2
2
2
|
Генерация матрицы 1/2 дробной реплики центрального композиционного плана для 5 факторов разделенного на блоки с числом центральных точек в плане равном 7. Максимальное количество точек в блоке плана - 20.
>> [D,blk] = ccdesign(5, 'center', 7, 'fraction', 1, 'blocksize', 20)
D =
-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
-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
1 1 1 -1 -1
1 1 1 1 1
-2 0 0 0 0
2 0 0 0 0
0 -2 0 0 0
0 2 0 0 0
0 0 -2 0 0
0 0 2 0 0
0 0 0 -2 0
0 0 0 2 0
0 0 0 0 -2
0 0 0 0 2
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 0 0 0 0
blk =
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
2
2
2
2
2
2
2
2
2
2
1
1
1
2
2
2
2
Генерация матрицы 1/2 дробной реплики центрального композиционного плана для 5 факторов разделенного на блоки с числом центральных точек расположенных с целью обеспечения равной точности в области проведения эксперимента. Вид центрального композиционного плана - 'inscribed'.
>> [D,blk] = ccdesign(5, 'center', 'uniform', 'fraction', 1, 'type', 'inscribed')
D =
-0.5000 -0.5000 -0.5000 -0.5000 0.5000
-0.5000 -0.5000 -0.5000 0.5000 -0.5000
-0.5000 -0.5000 0.5000 -0.5000 -0.5000
-0.5000 -0.5000 0.5000 0.5000 0.5000
-0.5000 0.5000 -0.5000 -0.5000 -0.5000
-0.5000 0.5000 -0.5000 0.5000 0.5000
-0.5000 0.5000 0.5000 -0.5000 0.5000
-0.5000 0.5000 0.5000 0.5000 -0.5000
0.5000 -0.5000 -0.5000 -0.5000 -0.5000
0.5000 -0.5000 -0.5000 0.5000 0.5000
0.5000 -0.5000 0.5000 -0.5000 0.5000
0.5000 -0.5000 0.5000 0.5000 -0.5000
0.5000 0.5000 -0.5000 -0.5000 0.5000
0.5000 0.5000 -0.5000 0.5000 -0.5000
0.5000 0.5000 0.5000 -0.5000 -0.5000
0.5000 0.5000 0.5000 0.5000 0.5000
-1.0000 0 0 0 0
1.0000 0 0 0 0
0 -1.0000 0 0 0
0 1.0000 0 0 0
0 0 -1.0000 0 0
0 0 1.0000 0 0
0 0 0 -1.0000 0
0 0 0 1.0000 0
0 0 0 0 -1.0000
0 0 0 0 1.0000
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
blk =
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
2
2
2
2
2
2
2
2
2
2
1
1
1
2
2
2
Графическое представление 3-х видов центрального композиционного плана для 3 факторов
>> nfactors=3;
>> D = ccdesign(nfactors, 'type', 'inscribed');
>> x= D(:,1);
>> y= D (:,2);
>> z= D (:,3);
>> subplot(1,3,1)
>> plot3(x,y,z,'o')
>> grid on
>> D1 = ccdesign(nfactors, 'type', 'circumscribed');
>> x1= D1(:,1);
>> y1= D1 (:,2);
>> z1= D1 (:,3);
>> subplot(1,3,2)
>> plot3(x1,y1,z1,'dr')
>> grid on
>> D2 = ccdesign(nfactors, 'type', 'faced');
>> x2= D2(:,1);
>> y2= D2 (:,2);
>> z2= D2 (:,3);
>> subplot(1,3,3)
>> plot3(x2,y2,z2,'+r')
>> grid on
\ \
|
