Полный факторный эксперимент (ПФЭ)


ПФЭ – называется эксперимент, реализующий все возможные повторяющиеся комбинации уровней независимых переменных, каждая из которых принудительно варьируется на этих двух уровнях. Число этих комбинаций clip_image002 определяет тип планирования.

g

clip_image004

clip_image006

clip_image008

1

-1

-1

-1

2

+1

-1

-1

3

-1

+1

-1

4

+1

+1

-1

5

-1

-1

+1

6

+1

-1

+1

7

-1

+1

+1

8

+1

+1

+1

граф пфэclip_image012

Требования:

1. Симметричность относительно центра эксперимента — clip_image014

2. Условие нормировки clip_image016

3. Ортогональность clip_image018

4. Ротатабельность, т.е. координаты точек в матрице планирования подбираются таким образом, что точность предсказания значения параметра стабилизации одинакова на равных расстояниях от центра эксперимента и не зависит от направления.

g

Порядок проведения эксперимента.

Планирование

Результаты

Обр. рез-тов

clip_image020

clip_image022

clip_image024

clip_image026

clip_image028

clip_image030

clip_image032

clip_image034

1…

L

…m

clip_image036

clip_image038

1…

l

…m

clip_image040

clip_image042

clip_image044

clip_image046

clip_image042[1]clip_image044[1]

clip_image042[2]clip_image046[1]

clip_image044[2]clip_image046[2]

clip_image042[3]clip_image044[3]clip_image046[3]

1

7

6

3

+

+

+

+

Y11

Y1l

Y1m

clip_image055

clip_image057

2

1

3

1

+

+

+

+

Y21

Y2l

Y2m

clip_image059

clip_image061

3

3

7

6

+

+

+

+

Y31

Y3l

Y3m

clip_image063

clip_image065

4

6

1

2

+

+

+

+

Y41

Y4l

Y4m

clip_image067

clip_image069

5

2

2

4

+

+

+

+

Y51

Y5l

Y5m

clip_image071

clip_image073

6

5

8

5

+

+

+

+

Y61

Y6l

Y6m

clip_image075

clip_image077

7

4

5

8

+

+

+

+

Y71

Y7l

Y7m

clip_image079

clip_image081

8

8

4

7

+

+

+

+

+

+

+

+

Y81

Y8l

Y8m

clip_image083

clip_image085

ПФЭПолучается неполную неточную модель, в которой нет ни одного слагаемого при clip_image089, но это допустимо т.к. квадратичное взаимодействие будет в точке экстремума.

При реализации плана необходимо:

1. В каждой точке этого плана необходимо получить выборку.

2. Надо воспользоваться рандомизацией, т.к. существует влияние неучтённых факторов.

Для каждой точки: clip_image091; clip_image093. Для сравнения дисперсий можно воспользоваться критерием Кохрена или Бартлетта. Дисперсию эксперимента нужно считать по формуле: а) clip_image095 — для Кохрена; б) Средневзвешенная дисперсия по Бартлетту.

Если дисперсии неразличимы, то, стало быть, есть грубый промах: 1. Выбросить грубый промах; 2. На его место становится среднее арифметическое.

Теперь находим коэффициенты модели, проверяем их на значимость, проверяем адекватность модели.

clip_image097;clip_image099;clip_image101 т.к. clip_image103

введём фиктивную переменную clip_image105, тогда clip_image107

Получаем исходные коэффициенты модели, т.е. нашли не истинные коэффициенты, а их оценки. Сравним их с НУЛЁМ – проверим на значимость.

clip_image109.

Если неравенство справедливо, то данный коэффициент признаём незначимым, а если нет, то оставляем его в нашей модели. Количество значимых коэффициентов модели обозначим за d.

clip_image111 — Это наша модель. clip_image113 — Предсказанное значение выходной величины. Мерой отклонения экспериментальных данных от предсказанных служит дисперсия адекватности.

clip_image115, где N>d. Если N = d, то мы не можем оценить адекватность нашей модели.

Для оценки адекватности модели пользуются критерием Фишера: clip_image117. Если неравенство выполняется, то модель адекватна, если нет, то неадекватна. Неадекватность модели ещё не говорит о том, что модель неверна. Возможные причины неадекватности модели:

1. Постулируемый вид модели не подходит для описания данного объекта.

2. Существует фактор, оказывающий сильное влияние на выходную величину, которым экспериментатор не управлял.

3. Мы находимся на вершине параболы, в этом случае большое влияние оказывают квадраты факторов, которые в нашей модели отсутствуют.

4. Шаг варьирования факторов слишком мал (слишком велик). Рекомендуется брать 10-15% от базовой точки.

5. Слишком большой шум эксперимента.

Для обработки результатов измерений, которые характеризуются статистически различными дисперсиями строчных выборок, может быть применён взвешенный метод наименьших квадратов: clip_image119, Тi имеет одинаковые дисперсии во всех точках факторного эксперимента.

Пример ПФЭ

Пусть влияют факторы:

X1 – Чистый (+1) или Грязный (-1) стаканчик.

X2 – Температура подогрева подложки 3600(-1) и 3800(+1)

X3 – Температура испарителя

Установка по напылению.

Полный факторный эксперимент (ПФЭ)

Составляем план ПФЭ:

g

X0

X1

X2

X3

X1X2…

clip_image036[1]

clip_image038[1]

clip_image125

1

+

…………

12.6

1.97

0.25

2

+

+

…………

17.3

2.33

0.16

3

+

+

…………

10.8

1.65

4

+

+

+

…………

20.0

2.86

5

+

+

…………

13.8

1.35

6

+

+

+

…………

12.4

1.44

7

+

+

+

…………

16.3

3.25

8

+

+

+

+

…………

14.4

1.71

clip_image127 < clip_image129

Если нер-во выполняется, то самая большая дисперсия не отличается от других, следовательно, можем определить дисперсию опыта. Дисперсия один из самых грубых параметров статистики по способу их нахождения.

Ряд clip_image038[2] тоже ряд С.В., но числа в квадрате и он подчиняется распределению c2 Пирсона.

clip_image131; clip_image133

Определим коэффициенты модели: clip_image135

clip_image137; clip_image139; clip_image141;clip_image143clip_image145; clip_image147; clip_image149; clip_image151; clip_image153.

Ни один коэфф. Кроме b0 не может быть больше других в 10 раз. Если не один коэфф. Нельзя отбросить то мы потеряли один (несколько) значимых факторов. У нас не осталось степеней свободы. Значимость коэффициентов определяется критерием Стьюдента.

clip_image155 Т.к. clip_image157. Упростим и получим clip_image159, а clip_image161. Все числа, которые больше табличного значения должны войти в модель, а остальные отбрасываем.

Получили следующую модель: clip_image163. Число членов уравнения –d.Находим дисперсию адекватности:

clip_image165(m элементов содержит выборка для каждого clip_image167). Когда 1-ца в числителе – жесткая оценка, Когда m – мягкая оценка- верхняя и нижняя границы существования модели. Для оценки адекватности модели пользуются критерием Фишера: clip_image169. Т. е. Условие не выполняется, и модель неадекватно описывает объект.

Неадекватность модели не означает её неправильность, всё сделано правильно, но что-то недоделано.

1. Неполный список факторов.

2. Слишком много выброшенных коэффициентов.

3. Противоречие в самой модели.

Если внутри модели есть противоречия, то они могут возникнуть от слишком точных измерений. Необходимо либо переделать эксперимент, либо обработать его более грубым методом. Таким методом является дробный факторный эксперимент.

Загрузка...