bbdesign

Проект поля-Behnken

Синтаксис

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

Описание

dBB = bbdesign(n) генерирует проект Поля-Behnken для 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.

Примеры

Следующее создает проект Поля-Behnken с 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
Для просмотра документации необходимо авторизоваться на сайте