bbdesign

Бокс-Бенкен- проект

Синтаксис

dBB = bbdesign(n)
[dBB,blocks] = bbdesign(n)
[...] = bbdesign(n,param,val)

Описание

dBB = bbdesign(n) генерирует проект Box-Bennken для n факторы. n должно быть целым числом 3 или больший. Выходная матрица dBB m -by- n, где m количество запусков в проекте. Каждая строка представляет один запуск с настройками для всех факторов, представленных в столбцах. Значения факторов нормированы так, что точки куба берут значения между -1 и 1.

[dBB,blocks] = bbdesign(n) запрашивает заблокированный проект. Область выхода blocks является вектором m -by-1 чисел блоков для каждого запуска. Блоки указывают запуски, которые должны быть измерены при сходных условиях, чтобы минимизировать эффект межблочных различий на оценки параметров.

[...] = bbdesign(n,param,val) задает одну или несколько необязательных пар параметр/значение для проекта. В следующей таблице перечислены допустимые пары параметр/значение.

ПараметрОписаниеЗначения
'center'

Количество центральных точек.

Целое число. Значение по умолчанию зависит от n.

'blocksize'

Максимальное число точек на блок.

Целое число. Значение по умолчанию является Inf.

Примеры

Следующее создает проект Box-Bennken с 3 факторами:

dBB = bbdesign(3)
dBB =
    -1    -1     0
    -1     1     0
     1    -1     0
     1     1     0
    -1     0    -1
    -1     0     1
     1     0    -1
     1     0     1
     0    -1    -1
     0    -1     1
     0     1    -1
     0     1     1
     0     0     0
     0     0     0
     0     0     0

Центральная точка запускается 3 раза, чтобы позволить более равномерную оценку отклонения предсказания во всем проектном пространстве.

Визуализируйте проект следующим образом:

plot3(dBB(:,1),dBB(:,2),dBB(:,3),'ro',...
         'MarkerFaceColor','b')
X = [1 -1 -1 -1 1 -1 -1 -1 1 1 -1 -1; ...
     1 1 1 -1 1 1 1 -1 1 1 -1 -1];
Y = [-1 -1 1 -1 -1 -1 1 -1 1 -1 1 -1; ...
     1 -1 1 1 1 -1 1 1 1 -1 1 -1];
Z = [1 1 1 1 -1 -1 -1 -1 -1 -1 -1 -1; ...
     1 1 1 1 -1 -1 -1 -1 1 1 1 1];
line(X,Y,Z,'Color','b')
axis square equal

См. также

Представлено до R2006a
Для просмотра документации необходимо авторизоваться на сайте