- •Планирование и проведение эксперимента при производстве электронно-вычислительных средств
- •Введение
- •Модуль 1 «методология математического моделирования» практическое занятие № 1
- •1.1 Экспериментальный анализ случайной величины
- •4 Определение оценок математического ожидания , дисперсии и среднего квадратического отклонения
- •1.2 Проверка статистических гипотез
- •1.3 Решение типовых примеров
- •1.4 Задачи для решения
- •1.5 Контрольные вопросы
- •Практическое занятие № 2
- •2.1 Метод ранговой корреляции
- •2.2 Решение типового примера
- •2.3 Задачи для решения
- •2.4 Контрольные вопросы
- •Практическое занятие № 3
- •3.1 Однофакторный дисперсионный анализ
- •3.2 Решение типового примера
- •3.3 Задачи для решения
- •3.4 Контрольные вопросы
- •Практическое занятие № 4
- •4.1 Двухфакторный дисперсионный анализ
- •4.2 Трехфакторный дисперсионный анализ
- •4.3 Решение типового примера
- •4.4 Задачи для решения
- •4.5 Контрольные вопросы
- •Практическое занятие № 5
- •5. Методы насыщенных и сверхнасыщенных планов
- •5.1 Метод насыщенных планов
- •5.2 Метод сверхнасыщенных планов (метод случайного баланса)
- •5.3 Решение типового примера
- •5.4 Задачи для решения
- •5.5 Контрольные вопросы
- •Модуль 2 «активный эксперимент» практическое занятие № 6
- •6. Полный факторный эксперимент
- •6.1 Планирование эксперимента
- •6.2 Проведение эксперимента
- •6.3 Обработка и анализ результатов эксперимента
- •6.4 Решение типового примера
- •6.5 Задачи для решения
- •6.6 Контрольные вопросы
- •Практическое занятие № 7
- •7.1 Дробный факторный эксперимент
- •7.2 Решение типового примера
- •7.3 Задачи для решения
- •7.4 Контрольные вопросы
- •Практическое занятие № 8
- •8 Центральные композиционные планы
- •8.1 Центральный композиционный ортогональный план (цкоп)
- •8.2 Центральный композиционный рототабельный план (цкрп)
- •8.3 Решение типового примера
- •Задачи для решения
- •8.5 Контрольные вопросы
- •Модуль 3 «пассивный эксперимент» практическое занятие №9
- •9.1 Метод регрессионного анализа
- •9.2 Решение типового примера
- •9.3 Задачи для решения
- •9.4 Контрольные вопросы
- •Модуль 4 «методы оптимизации» практическое занятие № 10
- •10.1 Методы оптимизации
- •10.2 Метод Гаусса-Зайделя
- •10.3 Метод случайного поиска
- •10.4 Градиентные методы
- •10.5 Метод крутого восхождения (метод Бокса-Уилсона)
- •10.6 Симплексный метод
- •10.7 Решение типового примера
- •10.8 Задачи для решения
- •10.9 Контрольные вопросы
- •Литература
- •Приложение а
- •Статистические таблицы
4.5 Контрольные вопросы
1. Какого типа практические задачи обычно решают методом дисперсионного анализа?
2. Как математически формулируется задача однофакторного дисперсионного анализа?
3. В чем заключается основная идея метода дисперсионного анализа?
4. Каким образом производится оценивание существенности влияния фактора в однофакторном дисперсионном анализе?
5. Как производится оценивание влияния двух факторов и их взаимодействий в двухфакторном дисперсионном анализе?
Практическое занятие № 5
5. Методы насыщенных и сверхнасыщенных планов
При исследовании сложных процессов исследователю приходится иметь дело с большим количеством факторов, которые способны оказать влияние на функцию отклика исследуемого процесса. Для первоначального построение «грубой модели» процесса желательно оставить только те факторы, которые оказывают сравнительно существенное влияние на функцию отклика, отбросив на первом этапе факторы, оказывающие незначительное влияние. Это можно сделать с помощью насыщенных и сверхнасыщенных планов.
5.1 Метод насыщенных планов
Насыщенные планы – планы, для которых число степеней свободы равно
N–k=1, (5.1)
то есть число вариантов условий проведения эксперимента (число номеров опытов) должно быть на единицу больше число рассматриваемых факторов.
Необходимым условием применения насыщенных планов является отсутствие влияния эффекта взаимодействия факторов на функцию отклика исследуемого процесса. Соблюдение этого условия основано на предпосылке, что на выходной параметр исследуемого процесса оказывают влияние лишь линейные эффекты и не влияют взаимодействия факторов.
При этом используют дробные реплики ПФЭ, стремясь к тому, чтобы все экспериментальные данные, полученные при N условий проведения эксперимента, были бы использованы для оценки коэффициентов при соответствующих переменных.
Если предполагается, что на функцию отклика исследуемого процесса способны оказывать влияние 15 факторов, то для отсеивания несущественных или оказывающих незначительное влияние факторов может быть использован ДФЭ типа 215-11 с числом различных условий эксперимента (минимальным числом опытов) N=16. Условие (5.1) в этом случае выполняется, так как N–k=16-15=1.
Число опытов N=16 предусматривает применение ПФЭ типа 24. Полином первого порядка в этом случае имеет следующий вид:
Y=b0+b1X1+b2X2+b3X3+b4X4+b12X1X2+b13X1X3+
+b14X1X4+b23X2X3+b24X2X4+b34X3X4+b123X1X2X3+ (5.2)
+b124X1X2X4+b134X1X3X4+b234X2X3X4+b1234X1X2X3X4.
Из приведенного полинома 1-го порядка (5.2) видно, что имеется 15 коэффициентов (без учета коэффициента b0). Поэтому, заменяя все члены полинома (5.2), учитывающие эффект влияния взаимодействия ранее выбранных четырех из пятнадцати рассматриваемых факторов, на одиннадцать оставшихся, получаем полином 1-го порядка с пятнадцатью факторами:
Y=b0+b1X1+b2X2+b3X3+b4X4+b5X5+b6X6+b7X7+b8X8+
+b9X9+b10X10+b11X11+b12X12+b13X13+b14X14+b15X15. (5.3)
В (5.3) имеем дело не с ПФЭ типа 24, а с ДФЭ типа 215-11, на основании которого можно оценить все пятнадцать коэффициентов b1, b2, b3,…, b15.
X5=X1X2X3X4; X10=X1X2;
X6=X1X2X3; X11=X1X3;
X7=X1X3X4; X12=X1X4;
X8=X1X2X4; X13=X2X3;
X9=X2X3X4; X14=X2X4;
X15=X3X4.
Проведя соответствующую замену в матрице ПФЭ типа 24 при использовании значений рассматриваемых в эксперименте 15-ти факторов, получим матрицу ДФЭ типа 215-11 (таблица 5.1). После проведения экспериментов производится вычисление коэффициентов по формуле (5.10).
Таблица 5.1 – Матрица насыщенного планирования
Номер опыта |
X0б |
X1б |
X2б |
X3б |
X4б |
X5б |
X6б |
X7б |
X8б |
X9б |
X10б |
X11б |
X12б |
X13б |
X14б |
X15б |
|
1 |
+ |
– |
– |
– |
– |
+ |
– |
– |
– |
– |
+ |
+ |
+ |
+ |
+ |
+ |
|
2 |
+ |
+ |
– |
– |
– |
– |
+ |
+ |
+ |
– |
– |
– |
– |
+ |
+ |
+ |
|
3 |
+ |
– |
+ |
– |
– |
– |
+ |
+ |
+ |
+ |
– |
+ |
+ |
– |
– |
+ |
|
4 |
+ |
+ |
+ |
– |
– |
+ |
– |
– |
+ |
+ |
+ |
– |
– |
– |
– |
+ |
|
5 |
+ |
– |
– |
+ |
– |
– |
+ |
– |
+ |
+ |
+ |
– |
+ |
– |
+ |
– |
|
6 |
+ |
+ |
– |
+ |
– |
+ |
– |
+ |
– |
+ |
– |
+ |
– |
– |
+ |
– |
|
7 |
+ |
– |
+ |
+ |
– |
+ |
– |
+ |
+ |
– |
– |
– |
+ |
+ |
– |
– |
|
8 |
+ |
+ |
+ |
+ |
– |
– |
+ |
– |
– |
– |
+ |
+ |
– |
+ |
– |
– |
|
9 |
+ |
– |
– |
– |
+ |
– |
– |
+ |
+ |
+ |
+ |
+ |
– |
+ |
– |
– |
|
10 |
+ |
+ |
– |
– |
+ |
+ |
+ |
– |
– |
+ |
– |
– |
+ |
+ |
– |
– |
|
11 |
+ |
– |
+ |
– |
+ |
+ |
+ |
– |
+ |
– |
– |
+ |
– |
– |
+ |
– |
|
12 |
+ |
+ |
+ |
– |
+ |
– |
– |
+ |
– |
– |
+ |
– |
+ |
– |
+ |
– |
|
13 |
+ |
– |
– |
+ |
+ |
+ |
+ |
+ |
– |
– |
+ |
– |
– |
– |
– |
+ |
|
14 |
+ |
+ |
– |
+ |
+ |
– |
– |
– |
+ |
– |
– |
+ |
+ |
– |
– |
+ |
|
15 |
+ |
– |
+ |
+ |
+ |
– |
– |
– |
– |
+ |
– |
– |
– |
+ |
+ |
+ |
|
16 |
+ |
+ |
+ |
+ |
+ |
+ |
+ |
+ |
+ |
+ |
+ |
+ |
+ |
+ |
+ |
+ |
|
Факторы, при которых коэффициенты в результате проведенной оценки по критерию Стьюдента оказались незначимыми, отбрасываются. На первых этапах исследования, когда создается «грубая» модель исследуемого процесса, допускается отсеивание несущественных факторов, исходя из значений полученных коэффициентов.
Если рассматривать процесс с числом факторов k=17, то число опытов ПФЭ типа 24 будет недостаточным. Ближайшее минимальное число опытов можно получить с помощью ПФЭ типа 25, которое составляет N=32. Число опытов в данном случае значительно превышает число учитываемых в эксперименте факторов, но облегчается замена эффектов взаимодействия на линейные эффекты. Все линейные эффекты могут быть введены в план вместо эффектов взаимодействия более высокого порядка, чем парные (по сравнению с k=15), а следовательно, менее значимыми с точки зрения их влияния на функцию отклика. Действительно,
X6=X1X2X3X4X5; X7=X1X2X3X4; X8=X1X3X4X5;
X9=X1X2X4X5; X10=X2X3X4X5; X11=X1X2X3;
X12=X1X3X4; X13=X1X4X5; X14=X1X3X5;
X15=X2X3X4; X16=X2X4X5; X17=X2X3X5.
Однако объем экспериментальной работы в данном случае увеличится не пропорционально увеличению числа рассматриваемых факторов, в отличие от предыдущего случая.
При k=9; 17; 33 и т.д. использование дробных реплик от ПФЭ ведет к значительному увеличению числа опытов соответственно N=16; 32; 64 и т.д. Для того, чтобы увеличить насыщенность планов, разработаны ортогональные планы с N=12; 20; 24; 36 и т.д. Однако применение метода насыщенных планов для исследования сложных процессов ограничено, так как эффект влияния взаимодействия факторов на выходной параметр может быть значительным.
