Пособие Моделирование1234 (538757), страница 6
Текст из файла (страница 6)
который содержит коэффициентов.
В рамках выбранной модели планирования в виде алгебраических полиномов строится план эксперимента путем варьирования каждого из факторов на нескольких уровнях q относительно исходной точки
, представляющей центр эксперимента.
Виды планов экспериментов. Эксперимент, в котором реализуются все возможные сочетания уровней факторов, называется полным факторным экспериментом (ПФЭ). Если выбранная модель планирования включает в себя только линейные члены полинома и их произведения, то для оценки коэффициентов модели используется план эксперимента с варьированием всех k факторов на двух уровнях, т. е. д=2. Такие планы называются планами типа 2 , где N=2 — число всех возможных испытаний.
Начальный этап планирования эксперимента для получения коэффициентов линейной модели основан на варьировании факторов на двух уровнях: нижнем и верхнем
— симметрично расположенных относительно основного уровня
,. Геометрическая интерпретация показана на рисунке а.
Так как каждый фактор принимает лишь два значения то для стандартизации и упрощения записи условий каждого испытания и обработки выборочных данных эксперимента масштабы по осям факторов выбираются так, чтобы нижний уровень соответствовал -1, верхний +1, а основной — нулю. Это легко достигается с помощью преобразования вида
где — кодированное значение i-гo фактора;
— натуральное значение фактора;
— нулевой уровень;
— интервал варьирования фактора.
Расположение точек для ПФЭ типа показано на рисунке б. Выписывая комбинации уровней факторов для каждой экспериментальной точки квадрата, получим план D полного факторного эксперимента типа
.
При этом планы можно записывать сокращенно с помощью условных буквенных обозначений строк. Для этого порядковый номер фактора ставится в соответствие строчной букве латинского алфавита: и т. д.
Таблица 3.1
Затем для каждой строки плана выписываются латинские буквы только для факторов, находящихся на верхних уровнях; испытание со всеми факторами на нижних уровнях обозначается как (1). Запись плана в буквенных обозначениях показана в последней строчке.
Геометрическая интерпретация ПФЭ приведена на рисунке, а его план ниже:
Таблица 3.2
Номер испытания | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | |
-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) | c | b | bc | a | ac | ab | abc |
Полный факторный эксперимент дает возможность определить не только коэффициенты регрессии, соответствующие линейным эффектам, но и коэффициенты регрессии, соответствующие всем эффектам взаимодействия. Эффект взаимодействия двух (или более) факторов появляется при одновременном варьировании этих факторов, когда действие каждого из них на выход зависит от уровня, на которых находятся другие факторы.
Для оценки свободного члена и определения эффектов взаимодействия
план эксперимента D расширяют до матрицы планирования X путем добавления соответствующей «фиктивной переменной»: единичного столбца
и столбцов произведений, как показано, например, для ПФЭ типа
в таблице 3.2
Таблица 3.3
Номер испытания | План ПФЭ | Реакция У | |||||||
1 | +1 | -1 | -1 | -1 | +1 | +1 | +1 | -1 | |
2 | +1 | -1 | -1 | +1 | +1 | -1 | -1 | +1 | |
3 | +1 | -1 | +1 | -1 | -1 | +1 | -1 | +1 | |
4 | +1 | -1 | +1 | +1 | -1 | -1 | +1 | -1 | |
5 | +1 | +1 | -1 | -1 | -1 | -1 | +1 | +1 | |
6 | +1 | +1 | -1 | +1 | -1 | +1 | -1 | -1 | |
7 | +1 | +1 | +1 | -1 | +1 | -1 | -1 | -1 | |
8 | +1 | +1 | +1 | +1 | +1 | +1 | +1 | +1 |
Как видно из рассмотренных планов экспериментов типов и
, количество испытаний в ПФЭ значительно превосходит число определяемых коэффициентов линейной модели плана эксперимента, т. е. ПФЭ обладает большой избыточностью и поэтому возникает проблема сокращения их количества.
Рассмотрим построение планов так называемого дробного факторного эксперимента. Пусть имеется простейший полный факторньй эксперимент типа . Используя матрицу планирования, приведенную в табл. 3.4, можно вычислить коэффициенты и представить результаты в виде уравнения
Таблица 3.4
Если в выбранных интервалах варьирования уровня процесс можно описать линейной моделью, то достаточно определить три коэффициента: и
. Таким образом, остается одна степень свободы, которую можно использовать для минимизации числа испытаний. При линейном приближении
и вектор-столбец
(табл. 3.4) можно использовать для нового фактора
. Поставим в табл. 3.4 этот фактор в скобках над взаимодействием
. В этом случае раздельных оценок, которые имели место в ПФЭ типа 2 , уже не будет и оценки смещаются следующим образом: при постулировании линейной модели все парные взаимодействия не учитывают. Таким образом, вместо восьми испытаний в полном факторном эксперименте типа
необходимо провести только четыре. Правило проведения дробного факторного эксперимента формулируется так: для сокращения числа испытаний новому фактору присваивается значение вектор-столбца матрицы, принадлежащего взаимодействию, которым можно пренебречь.
При проведении эксперимента из четырех испытаний для оценки влияния трех факторов пользуются половиной ПФЭ типа , так называемой «полурепликой». Если приравнять
и
, то можно получить вторую «полуреплику». Для обозначения дробных реплик, в которых d линейных эффектов приравнены к эффектам взаимодействия, пользуются условным обозначением
. Например, «полуреплика» от
записывается в виде
, а «четвертьреплика» —
.
Когда модель планирования анализируется методами дисперсионного анализа, применяют планы дисперсионного анализа. Если при постановке эксперимента реализуются все возможные совокупности условий, то говорят о полных классификациях дисперсионного анализа. Если проводится сокращение перебора вариантов — это неполная классификация дисперсионного анализа. Сокращение перебора может проводиться случайным образом (без ограничения на рандомизацию) или в соответствии с некоторыми правилами (с ограничениями на рандомизацию). Чаще всего в качестве таких планов используют блочные планы и планы типа латинского квадрата.