Министерство образования РФ
Волгоградский государственный технический университет
Кафедра «САПР и ПК»
Курсовая работа
по моделированию
на тему:
«Математическое моделирование при активном эксперименте»
Выполнил: студент
II
-го курса
группы ИВТ-263
Б
********
Ю.В.
Проверил: Фоменков С.А.
Волгоград 2001
Основные положения теории планирования эксперимента.
Различают физическое и математическое моделирование. При физическом моделировании исследование объекта происходит при его воспроизведении в ином масштабе. Здесь возможен количественный перенос результатов эксперимента с модели на оригинал. Однако для анализа сложных объектов и процессов, каковыми являются большинство электронных схем, конструкций и технологических процессов производства радиоэлектронной техники, приборостроения, машиностроения и других промышленных отраслей, применение физического моделирования затруднительно, поскольку приходится использовать большое число критериев и ограничений, которые могут быть несовместимы, а зачастую и невыполнимы.
Математическое моделирование является методом качественного или количественного описания объектов или процессов, при этом реальный объект, процесс или явление упрощается, схематизируется и описывается определенным уравнением. В большинстве случаев математическая модель представляет собой уравнение регрессии, то есть геометрическое место точек математических ожиданий условных распределений целевой функции. Простейшим примером такой модели является уравнение парной корреляции, где на целевую функцию воздействует один фактор. На практике в реальном производстве на целевую функцию воздействуют много факторов и искомое уравнение регрессии становится многомерным.
Существует много методов отыскания уравнения регрессии, которые можно условно разделить на два класса: методы активного и методы пассивного эксперимента. Под активным экспериментом
будем понимать эксперимент, предварительный план которого составлен так, чтобы получить максимальную информацию о целевой функции при минимальной ее дисперсии и проведении минимального числа опытов (эффективный план). Такой план (например, полный факторный эксперимент) требует искусственного одновременного варьирования всеми факторами в довольно широких пределах. Методы активного эксперимента довольно хорошо разработаны в специальном разделе математической статистики, который называется "Теория планирования эксперимента".
Под математической теорией планирования эксперимента
будем понимать науку о способах составления экономных экспериментальных планов, которые позволяют извлекать наибольшее количество информации об объекте, о способах проведения эксперимента, о способах обработки экспериментальных данных, о способах использования полученных результатов для оптимизации исследуемых объектов (например, технологических процессов производства массовой продукции). Математический аппарат теории планирования эксперимента построен на сочетании методов математической статистики и методов решения экстремальных задач.
В настоящее время выделяют два основных направления теории планирования эксперимента:
планирование экстремальных экспериментов;
планирование экспериментов по выявлению механизма явлений.
В этой курсовой работе описываются в основном методы первого направления.
Любое экспериментальное исследование содержит три этапа:
этап постановки задачи;
этап планирования и проведения эксперимента;
анализ и интерпретация результатов.
Главной трудностью на этапе постановки задачи является переход с языка специальности на язык планирования эксперимента, на язык математики.
Построение математической модели технологического процесса в зависимости от поставленной задачи может преследовать следующие цели: минимизировать расход материала на единицу выпускаемой продукции при сохранении качества, произвести замену дорогостоящих материалов на более дешевые или дефицитных на распространение; сократить время обработки в целом или на отдельных операциях, перевести отдельные режимы в некритические зоны, снизить трудовые затраты на единицу продукции и т.п.; улучшить частные показатели и общее количество готовой продукции, повысить однородность продукции, улучшить показатели надежности и т.п.; увеличить надежность и быстродействие управления, увеличить эффективность контроля качества, создать условия для автоматизации процесса управления и т.п.
Прежде всего, необходимо выбрать зависимую переменную Y, которую впредь будем называть целевой функцией или параметром оптимизации, за который принимают один из показателей качества продукции либо по каждой технологической операции отдельно, либо по всему технологическому процессу сразу. Параметр оптимизации должен соответствовать следующим требованиям:
параметр должен измеряться при любом изменении (комбинации) режимов технологического процесса;
параметр должен быть статистически эффективным, то есть измеряться с наибольшей точностью;
параметр должен быть информационным, то есть всесторонне характеризовать технологический процесс (операцию);
параметр должен иметь физический смысл, то есть должна быть возможность достижения полезных результатов при соответствующих условиях процесса;
параметр должен быть однозначным, то есть должно минимизироваться или максимизироваться только одно свойство изделия.
За фактор принимают контролируемую величину объекта (изделия, процесса, операции), то есть величину, характеризующую то или иное свойство объекта или режим технологического оборудования. Эта величина, числовое значение которой измеряется в пределах (границах) изменения, должна влиять на параметр оптимизации.
При определении величин количественных оценок во внимание должны приниматься только те факторы, которые имеют четкий метрологический смысл (возможность измерения фактора с определенной точностью).
Описание исследуемого объекта нельзя получить в виде точной формулы функции, справедливой во всем диапазоне существования аргументов. Оно может быть лишь приближенным и на небольшом участке в окрестностях выбранной базовой точки. Аппроксимация искомой математической зависимости представляет собой некоторый полином - отрезок ряда Тейлора, в который разлагается неизвестная зависимость:
| (1) |
где:
В силу наличия неуправляемых и даже неконтролируемых входных переменных Xi изменение величины Y носит случайный характер, а потому уравнение (1) не дает нам точной связи между входом и выходом объекта и является лишь условным математическим ожиданием случайной величины Y, т.е. уравнением регрессии.
Чтобы отыскать коэффициенты уравнения регрессии по результатам экспериментов в N точках факторного пространства (что является типичной задачей регрессионного анализа), необходимо выполнение следующих предпосылок:
Результаты наблюдений Y 1 , Y 2 ,...,Y n выходной величины в N точках факторного пространства представляют собой независимые, нормально распределенные случайные величины.
Выборочные дисперсии опытов однородны, т.е. статистически неразличимы. Это требование означает независимость выборочной дисперсии от местоположения точки факторного пространства, в которой проводится конкретный опыт (ротатабельность).
Независимые переменные X 1 , X 2 ,...,X n измеряются с ошибкой много меньшей, чем величина возможного отклонения выходного параметра Y под влиянием неучтенных факторов.
Тогда задача отыскания коэффициентов уравнения регрессии сводится к решению системы так называемых нормальных уравнений:
| (2) |
где Y g - экспериментальные значения выходного параметра, полученные в g-й точке факторного пространства;- значение выходного параметра, найденные по уравнению регрессии в тех же точках; d - количество членов в уравнении регрессии.
Выражение (2) является основным критерием проверки правильности найденного уравнения регрессии.
Чтобы система нормальных уравнений, которая может быть представлена в виде матрицы, имела единственное решение, необходимо, чтобы матрица была невырожденной, т.е. чтобы вектор – столбцы были линейно – независимы. Чтобы величины коэффициентов уравнения регрессии не зависели от числа членов матрицы, нужно на нее наложить дополнительное условие ортогональности вектор столбцов.
1. Полный факторный эксперимент
Полным факторным экспериментом (ПФЭ) называется эксперимент, реализующий все возможные повторяющиеся комбинации уровней независимых переменных, каждая из которых принудительно варьируется на двух уровнях (Рис 1)
| Рис. 1 Схема перехода в относительные координаты |
Число этих комбинаций N=2 n определяет тип планирования.
Для гарантированного получения единственного решения системы нормальных уравнений необходимо иметь ортогональную матрицу планирования, что невозможно обеспечить в абсолютной системе единиц факторов X i , то есть тогда, когда факторы именованные (например, трудно представить 17 километров ортогональными к 12 килограммам). Поэтому необходимо провести предварительное преобразование каждого фактора - его перевод в систему относительных координат. Такое преобразование легко сделать с помощью переноса начала координат в базовую точку X * и выбора единицы отсчета D X i по каждой координате X i .
| (3) |
Это дает возможность легко построить ортогональную матрицу планирования и значительно облегчает дальнейшие расчеты, так как в этом случае верхние и нижние уровни варьирования X iв и X iн в относительных единицах будут равны соответственно x iв = +1 и x iн = -1.
Шаг варьирования по каждой переменной выбирается таким, чтобы приращение величины выходного параметра Y к базовому значению Y * при реализации шага можно было выделить на фоне "шума" при небольшом числе параллельных опытов. Если нет никаких указаний на величину шага D X i , то в первом приближении можно выбрать D X i = 0,15X * i , т.е. принять за шаг 15%-ное отклонение от базового уровня X * i . Такой шаг дает достаточную гарантию того, что фактор X i вызовет заметную реакцию Y, если связь между ними существует.
Матрица планирования должна отвечать следующим условиям:
1. Ортогональность
2. Условие нормированости
3. Симметричность относительно центра экстремума
4. Ротатабельность, т.е. координаты точек факторного пространства в матрице планирования подстраиваются так, что точность предсказания значения параметра оптимизации одинакова на равных расстояниях от центра эксперимента (базовой точки) и не зависит от направления.
Матрица планирования составляется по следующим правилам:
1. Каждая g-я строка матрицы представляет собой набор координат точки g , в которой производится эксперимент;
2. Поскольку переменные x gi принимают лишь значения +1 и -1, то все остальные переменные могут принимать те же значения, что позволяет в целях упрощения записывать в таблицу вместо +1 и -1 их знаки + и -;
3. Первая строка 1 выбирается так, чтобы управляемые переменные находились на нижнем уровне, т.е. x i1 = -1. Последующие строки при составлении матрицы планирования набираются по правилу: при построчном переборе всех вариантов частота смены знака управляемых переменных для каждой последующей переменной вдвое меньше, чем для предыдущей (см. табл. 1)
| Таблица 1 | ||||
| Матрица планирования трехфакторного эксперимента | ||||
| g | x 1 | x 2 | x 3 | |
| 1 | - | - | - | |
| 2 | + | - | - | |
| 3 | - | + | - | |
| 4 | + | + | - | |
| 5 | - | - | + | |
| 6 | + | - | + | |
| 7 | - | + | + | |
| 8 | + | + | + | |
Следует отметить, что суть матрицы не изменится, если первая строка 1 будет выбрана так, чтобы управляемые переменные находились на верхнем уровне, т.е. x i1 = +1.
Матрицы планирования любого другого типа, например, 2 4 , 2 5 и т.д. могут быть получены описаным выше способом.
Поскольку изменение выходной величины Y носит случайный характер, необходимо в каждой точке g (т.е. в точке с координатами, записаными в g-й строке) проводить m параллельных опытов и результаты наблюдений Y 1g ,Y 2g ,...,Y mg усреднять
| (4) |
Величина m может быть любой, но не меньше m=3. Тогда эксперимент делится на m серий опытов, в каждой из которых полностью реализуется матрица планирования (т.е. эксперимент проводится в N=2 n точках факторного пространства).
Одним из важнейших положений современной теории планирования эксперимента является рандомизация
. План эксперимента составляется так, чтобы рандомизировать, т.е. сделать случайными те систематически действующие факторы, которые трудно поддаются учету и контролю, для того, чтобы рассматривать их как случайные величины и учитывать статистически.
Перед реализацией плана на объекте необходимо произвести рандомизацию - с помощью таблицы равномерно распределенных случайных чисел (табл.П.6) определить последовательность реализации матрицы планирования в каждой из m серий опытов. Для этого в качестве начала выбирается любое число из табл.П.6 и записывается в столбец k 1 из табл.2 на место g=1. Остальные места этого столбца заполняют числа от 1 до N, следующие по порядку из табл.П.6 за выбранным начальным. Следует обращать внимание на то, чтобы числа в столбцах табл.2 не повторялись дважды. Пусть, например, при g=4 k 14 =8, это значит, что в первой серии испытаний точка 4 реализуется восьмой по порядку.
Аналогично рандомизируются испытания в каждой из оставшихся серий экспериментов; порядок реализации записывается в столбцах k 2 ,k 3 ,...,k m . Результаты эксперимента в каждой из серий испытаний записываются в столбцах Y 1 ,Y 2 ,...,Y m .
Проверка воспроизводимости - это проверка на выполнение второй предпосылки регрессионного анализа об однородности выборочных дисперсий S 2 g . Задача состоит в проверке гипотезы о равенстве дисперсий s 2 {Y 1 } =s 2 {Y 2 } =...s 2 {Y N } при экспериментах соответственно в точках 1 ,2 ,...,g ,...,N .
Оценки дисперсий находятся по формуле
| (5) |
Так как все дисперсии получены по выборкам одинакового объема m, то число степеней свободы для всех дисперсий одинаково и равно
| v 1 = m-1 | (6) |
Для проверки гипотезы об однородности оценок дисперсий следует пользоваться критерием Кохрена, который основан на законе распределения отношения максимальной эмперической дисперсии к сумме всех дисперсий, т.е.
| (7) |
Если вычисленное значение критерия G окажется меньше табличного значения G кр , найденного для q%-ного уровня значимости, v зн = v 2 = N - числа степеней свободы знаменателя (например для q=5%; v числ = 3 - 1 = 2; v зн =8, G кр = 0,5157, см. табл.П.5), то гипотеза об однородности дисперсий принимается. При этом всю группу дисперсий S 2 g можно считать оценкой S 2 {Y} одной и той же генеральной дисперсии воспроизводимости s 2 {Y}, откуда
| (8) |
Если проверка на воспроизводимость дала отрицательный результат, то остается признать либо невоспроизводимость эксперимента относительно управляемых переменных вследствие наличия флуктуаций неуправляемых и неконтролируемых переменных, создающих на выходе большой уровень "шума", либо наличие грубого промаха в строке, откуда взята дисперсия max{S 2 g }. В первом случае следует увеличить число параллельных опытов, во втором - найти грубый промах и заменить его на результат доброкачественного измерения при соответствующей комбинации факторов. Если это по каким-то причинам невозможно, то, чтобы не нарушать предпосылки использования критерия Кохрена, на место грубого промаха следует поместить среднюю арифметическую величину g данной строки.
Следует также отметить, что критерий Кохрена можно применять не к любой группе выборок, а только к группе выборок одинакового объема, что как раз и имеет место при полном факторном эксперименте.
Легко заметить, что исходный план (табл.1) содержит много больше строк, чем столбцев и, следовательно, из результатов эксперимента согласно условию решения нормальных уравнений (2) можно получить дополнительную информацию, т.е. расширить модель. Безусловно, это относится к средней арифметической всего эксперимента, т.е. к отклику в базовой точке b 0 , для расчета которого можно ввести фиктивную переменную x од = +1 для всех строк. Оставшиеся свободными столбцы можно использовать для нахождения оценок коэффициентов при парных взаимодействиях и т.п. При этом соответствующие величины x i x j , x i x j x l получаются простым перемножением соответствующих столбцов исходного плана.
Тогда математическая модель объекта, получающаяся в результате ПФЭ может быть представлена в виде
| Y = b 0 + b 0 x 1 + b n x n + b 12 x 1 x 2 + b (n-1) x 1 x 2 + b 123 x 1 x 2 x 3 + b 123...n x 1 x 2 x 3 x 3 | (9) |
Однако вследствие того, что из ограниченного числа опытов нельзя получить точные значения коэффициентов b i , а только их независимые оценки b i , вся математическая модель становится оценочной
| = b 0 + b 1 x 1 +...+ b n x n + b 12 x 1 x 2 + b 1...n x 1 ...x n | (10) |
Пример матрицы планирования, принцыпа ее реализации и последующей обработки экспериментальных данных приведен в табл.2 на базе трехфакторного эксперимента. В разделе "Матрица планирования эксперимента" включены не только относительные переменные x i , сочетание которых и является собственно настоящей матрицей планирования, ни и их парные и тройные взаимодействия, знание которых необходимо лишь на этапе обработки экспериментальных данных.
| Таблица 2 |
| Матрица планирования ПФЭ типа N=2 3 и обработка его результатов |
номер строки g |
Порядок реализации опытов | Матрица планирования эксперимента | Результаты эксперимента | Первичная обра ботка результатов | Проверка адекватности | |||||||||||||||||
| l | z 0
|
z 1
|
z 2
|
z 3
|
z 4
|
z 5
|
z 6
|
z 7
|
Y g1
|
... | Y gl
|
... | Y gm
|
g | S 2 g | g | (g -g ) 2 | |||||
| k 1 | ... | k l | ... | k m
|
x 0 | x 1 | x 2 | x 3 | x 1 x 2
|
x 1
x 3 |
x 2
x 3 |
x 1
x 2 x 3 |
||||||||||
| 1 | 1 | ... | 6 | ... | 8 | + | - | - | - | + | + | + | - | |||||||||
| 2 | 7 | ... | 5 | ... | 4 | + | + | - | - | - | - | + | + | |||||||||
| 3 | 3 | ... | 7 | ... | 6 | + | - | + | - | - | + | - | + | |||||||||
| 4 | 8 | ... | 2 | ... | 7 | + | + | + | - | + | - | - | - | |||||||||
| 5 | 6 | ... | 3 | ... | 2 | + | - | - | + | + | - | - | + | |||||||||
| 6 | 4 | ... | 4 | ... | 1 | + | + | - | + | - | + | - | - | |||||||||
| 7 | 2 | ... | 1 | ... | 5 | + | - | + | + | - | - | + | - | |||||||||
| 8 | 5 | ... | 8 | ... | 3 | + | + | + | + | + | + | + | + | |||||||||
Для удобства расчетов и представления формул каждый столбец может быть представлен в виде новой переменной Z ig . Тогда оценки коэффициентов уравнения регрессии легко найти по формуле
| (11) |
Легко заметить, что матрица планирования является ортогональной с линейно независимыми вектор-столбцами; отсюда следует диагональность матрицы нормальной системы уравнений, а следовательно, и взаимная независимость оценок коэффициентов уравнения регрессии.
Необходимо отметить, что получаемая модель не дает членов типа x 2 ii и, таким образом, является неполной. В большинстве случаев это не отражается на качестве модели, так как чаще всего b ii =0. Однако в случаях, когда b ii ¹ 0, модель становится неточной (неадекватной), тогда следует от ПФЭ переходить к другим принципам планирования (как правило, это случается в окрестностях частного или глобального экстремума целевой функции).
После определения оценок коэффициентов регрессии необходимо проверить гипотезу о значимости коэффициентов b i . Лучше всего это сделать в виде нуль-гипотезы, т.е. гипотезы о равенстве b i = 0. Если она подтвердилась, то коэффициент b i следует признать статистически незначимым и отбросить из искомой модели; если гипотеза не подтвердилась, то соответствующий коэффициент b i следует признать значимым и включить в модель.
Проверка гипотезы проводится с помощью t - критерия Стъюдента, который при проверка нуль-гипотезы формируется в виде
| (12) |
где S 2 {b i }- дисперсия ошибки определения коэффициента b i . При полном и дробном факторном планировании для всех i
| (13) |
Если вычисленная величина параметра t i превышает табличное значение t кр , найденное для q%-ного уровня значимости и v з =N(m-1) числа степеней свободы (например для q = 5%; v з = 16; t кр = 2,199, см.табл.П.2) то нуль-гипотеза отвергается и коэффициент считается незначимым и его следует отбросить, не включая в искомую модель.
Статистическая незначимость коэффициента b i может быть обусловлена следующими причинами:
уровень базового режима * близок к точке частного экстремума по переменной X i или по произведению переменных;
шаг варьирования D X i выбран малым;
данная переменная (или произведение переменных) не имеет функциональной связи с выходным параметром Y;
велика ошибка эксперимента вследствие наличия неуправляемых и неконтролируемых переменных.
Поскольку ортогональное планирование позволяет определять доверительные границы для каждого из коэффициентов регрессии в отдельности, то, если какой-либо из коэффициентов окажется незначимым, он может быть отброшен без пересчета всех остальных. После этого математическая модель объекта составляется в виде уравнения связи выходного параметра Y и переменных x i , включающего только значимые коэффициенты.
Чтобы проверить гипотезу об адекватности представления результатов эксперимента найденному уравнению связи (иными словами, чтобы проверить, насколько найденное уравнение соответствует экспериментальным результатам), достаточно оценить отклонение выходной величины Y g , предсказанное уравнением регрессии, от результатов экспериментов g в точках факторного пространства.
Рассеяние результатов эксперимента вблизи уравнения связи, аппроксимирующего искомую функциональную зависимость, можно охарактеризовать с помощью дисперсии неадекватности s 2 ад , оценка которой S 2 ад находится по формуле
| (14) |
с числом степеней свободы v ад = N-d, где d - число членов аппроксимирующего полинома.
Проверка адекватности состоит в выяснении соотношения между дисперсией неадекватности s 2 ад и дисперсией воспроизводимости s 2 {Y}. Если s 2 ад не превышает дисперсии опыта, то полученная математическая модель адекватно представляет результаты эксперимента, если же s 2 ад > s 2 {Y}, то описание считается неадекватным объекту.
Проверка гипотезы об адекватности проводится с использованием F-критерия Фишера.
Критерий Фишера позволяет проверить нуль-гипотезу о равенстве двух генеральных дисперсий s 2 ад и s 2 {Y}. В связи с тем, что самих генеральных дисперсий мы не знаем, F-критерий формируется как отношение
| (15) |
Если вычисленное по формуле (15) значение критерия F меньше табличного F кр , найденного для q%-ного уровня значимости, v числ = v ад = v 4 = N-d числа степеней свободы числителя и v зн = v з = N(m-1) числа степеней свободы знаменателя, то нуль-гипотеза принимается. В противном случае она отвергается и описание (модель) признается неадекватным объекту. Некоторые значения F кр (q=5%;v 4 ;v з ) приведены в табл.П.4
В ходе работы может возникнуть ситуация, когда выборочная дисперсия неадекватности S 2 ад не превосходит оценки дисперсии воспроизводимости S 2 {Y} (т.е. когда S 2 ад £ S 2 {Y}). Тогда соотношение (15) будет равно F £ 1 и неравенство F < F кр выполняется для любого числа степеней свободы v 4 и v 3 , т.е. гипотеза s 2 ад £s 2 {Y} не противоречит выборочным данным и математическая модель адекватно представляет объект.
Проверка адекватности возможна только при v ад = v 4 > 0. Число вариантов варьирования плана ПФЭ равно числу оцениваемых коэффициентов регрессии уравнения связи (N = d). Следовательно, не остается степеней свободы (v ад = 0) для проверки нуль-гипотезы об адекватности представления экспериментальных данных выбранной формой аппроксимирующего полинома. Если же некоторые коэффициенты регрессии оказались незначимыми или ими можно пренебречь в силу их малости, то число членов проверяемого уравнения в этом случае будет меньше числа вариантов варьирования (d<N), и одна или несколько степеней свободы (v ад >0) останется для проверки гипотезы адекватности.
Если гипотеза адекватности отвергается, то модель признается неадекватной экспериментальным данным. Неадекватность модели не означает ее неправильности!
Неадекватность модели может означать, что не весь перечень влияющих факторов был принят во внимание, или что необходимо перейти к более сложной форме уравнения связи, или выбрать другой шаг варьирования по одному или нескольким факторам и т.п. Однако все достижения неадекватной модели: отсев незначимых факторов, оценка дисперсии эксперимента и др. остаются в силе.
Пример 1. Методом ПФЭ найти математическую модель процесса напыления резисторов.
После консультации с экспертами и некоторых предварительных исследований было определено, что на величину сопротивления напыляемых резисторов могут оказывать влияние следующие факторы:
Состояние испарителя - "чистое", т.е. порошок для напыления сыпется в стакан испарителя впервые после промывки его сторон, или "грязное", т.е. порошок сыпется в испаритель, в котором осталось некоторое его количество от предыдущего цикла напыления; обозначим этот фактор как x 1 , причем величина x 1 = +1соответствует "чистому", а величина x 1 = -1 соответствует "грязному" состоянию испарителя;
Температура подогрева подложки x 2 , причем x 2 = +1 соответствует верхней допустимой по техпроцессу температуре, а x 2 = -1 - нижней;
Температура испарителя x 3 , причем x 3 = +1 соответствует верхней допустимой по техпроцессу температуре, а х 3 = -1
План эксперимента, его пятикратная реализация с учетом рандомизации и первичная обработка результатов представлена в таблице.
номер строки g |
Циклы | z 0 | z 1
|
z 2
|
z 3
|
z 4
|
z 5
|
z 6
|
z 7
|
Результаты, кOм | Обработка | Адекватность | ||||||||||
| g | S 2 g | g | (g -g ) 2 | |||||||||||||||||||
| k 1 | k 2 | k 3 | k 4 | k 5 | x 0 | x 1 | x 2 | x 3 | x 1 x 2 | x 1 x 3 | x 2
x 3 |
x 1
x 2 x 3 |
Y g1 | Y g2
|
Y g3
|
Y g4
|
Y g5
|
|||||
| 1 | 4 | 2 | 3 | 6 | 8 | + | - | - | - | + | + | + | - | 11,4 | 10,5 | 13,8 | 14,0 | 12,1 | 12,36 | 2,303 | 12,10 | 0,0676 |
| 2 | 3 | 3 | 6 | 2 | 5 | + | + | - | - | - | - | + | + | 18,1 | 17,4 | 15,2 | 16,8 | 19,2 | 17,34 | 2,228 | 17,08 | 0,0676 |
| 3 | 8 | 6 | 2 | 4 | 1 | + | - | + | - | - | + | - | + | 10,8 | 9,3 | 11,6 | 12,1 | 9,8 | 10,72 | 1,387 | 10,98 | 0,0676 |
| 4 | 6 | 1 | 7 | 1 | 6 | + | + | + | - | + | - | - | - | 18,8 | 29,6 | 22,0 | 22,8 | 20,7 | 21,38 | 2,752 | 21,64 | 0,0676 |
| 5 | 5 | 8 | 1 | 3 | 4 | + | - | - | + | + | - | - | + | 12,9 | 12,8 | 13,6 | 15,2 | 14,0 | 13,70 | 0,950 | 13,98 | 0,0784 |
| 6 | 2 | 5 | 5 | 7 | 2 | + | + | - | + | - | + | - | - | 12,0 | 11,6 | 14,2 | 13,4 | 12,5 | 12,74 | 1,118 | 13,00 | 0,0676 |
| 7 | 1 | 7 | 4 | 8 | 7 | + | - | + | + | - | - | + | - | 15,1 | 14,8 | 16,8 | 18,1 | 17,0 | 16,36 | 1,913 | 16,10 | 0,0676 |
| 8 | 7 | 4 | 8 | 5 | 3 | + | + | + | + | + | + | + | + | 13,5 | 11,9 | 14,3 | 17,0 | 16,2 | 14,58 | 4,227 | 14,32 | 0,0676 |
| å | 119,18 | 16,878 | - | 0,5410 | ||||||||||||||||||
При первичной обработке результатов экспериментов пользуемся формулами (4) и (5), а затем проверяем воспроизводимость опытов по (7)
Таким образом, подтверждена воспроизводимость опытов (отсутствие в данных грубых промахов), что позволяет, в свою очередь, найти среднюю дисперсию строчных выборок (дисперсию опытов) по (8)
| C | v 3 = 8·(5-1) = 32 степенями свободы |
Оценки коэффициентов уравнения регрессии ищутся по формуле (11)
и т.д. Аналогично находим b 3 = -0,55; b 12 = +0,61; b 13 = -2,30; b 23 = +0,26; b 123 = -0,81
Проверяем значимость оценок коэффициентов по критерию Стьюдента по формуле (12), предварительно найдя дисперсию оценок по формуле (13)
| ; | |||||||
| Тогда | ; | ; | |||||
| далее аналогично | t 12 = 2,602 | ; | t 13 = 9,812 | ; | t 23 = 1,109 | ; | t 123 = 3,455 |
Табличное значение критерия t i (табл.П.2) t кр (5%;v 3 =32) = 2,046, поэтому все найденные оценки коэффициентов, кроме b 23 , признаются значимыми и должны войти в модель
= 14,90 + 1,61x 1 + 0,86x 2 -0,55x 3 + 0,61x 1 x 2 -2,30x 1 x 3 - 0,81x 1 x 2 x 3
Для определения дисперсии адекватности по формуле (14) необходимо сначала найти числовые значения модели g для каждой g-ой строки матрицы планирования, а затем подсчитать сумму квадратов разностей между модельным значением и средним арифметическим g той же строки
Тогда критерий Фишера (15) дает
что доказывает адекватность найденной модели. Ее можно использовать для управления технологическим процессом испытания резисторов
2. Дробный факторный эксперимент
Полный факторный эксперимент целесообразно использовать при сравнительно небольшом числе независимых факторов (обычно не больше 5), в противном случае число вариантов варьирования N = 2 n становится непомерно большим и реализация эксперимента затрудняется. В то же время в большинстве практических задач взаимодействия внешних порядков, начиная с третьего (а то и второго), отсутствуют или пренебрежимо малы, вследствие чего излишне много степеней свободы остается на проверку гипотезы адекватности. Если заранее пренебречь взаимодействиями высших порядков, то имеется возможность получить математическую модель при меньшем числу опытов, реализовав не весь план ДФЭ, а только его часть (дробную реплику).
Эксперимент, реализующий часть (дробную реплику) полного факторного эксперимента, называется дробным факторным экспериментом
(ДФЭ). ДФЭ позволяет получить приближение искомой функциональной зависимости Y = f(X 1 ,...,X n ) в некоторой небольшой окрестности точки базового режима при минимуме опытов.
Так, для решения трехфакторной задачи можно ограничиться четырьмя вариантами (N = 4), если в планировании ПФЭ типа 2 2 произведение x 1 x 2 приравнять к третьей независимой переменной x 3 . Такое планирование, представленное матрицей табл 3, позволяет оценить свободный член b 0 и три коэффициента регрессии при линейных членах b 1 ,b 2 ,b 3 (из четырех опытов нельзя получить более четырех коэффициентов).
| Таблица 3 | ||||||||
| Полуреплика от ПФЭ типа 2 3 (планирование типа 2 3-1 ) | ||||||||
| g
|
z
0 |
z
1 |
z
2 |
z
3 |
z
4 |
z
5 |
z
6 |
z
7 |
| x
0 |
x
1 |
x
2 |
x
3 |
x
1 x 2 |
x
1 x 3 |
x
2 x 3 |
x
1 x 2 x 3 |
|
| 1
|
+ | - | - | + | + | - | - | + |
| 2
|
+ | + | - | - | - | - | + | + |
| 3
|
+ | - | + | - | - | + | - | + |
| 4
|
+ | + | + | + | + | + | + | + |
Применение ДФЭ всегда связано со смешиванием, т.е. совместной оценкой нескольких коэффициентов уравнения связи. В нашем примере, если коэффициенты регрессии b ij при парных произведениях отличны от нуля, то каждый из найденных коэффициентов будет оценкой двух теоретических коэффициентов:
b 0 ® b 0 + b 123 ; b 2 ® b 2 + b 13 ;
b 1 ® b 1 + b 23 ; b 3 ® b 3 + b 12 .
Действительно, указанные коэффициенты в таком планировании не могут быть найдены раздельно, поскольку столбцы матрицы для линейных членов и парных произведений совпадают (полностью скоррелированы). Рассмотренный план ДФЭ представляет половину плана ДФЭ типа 2 3 и называется "полурепликой" от ПФЭ типа 2 3 или планированием типа N = 2 3-1 .
При большом числе переменных можно построить дробные реплики высокой степени дробности (1/4, 1/8, 1/16 и т.д.). Дробная реплика обозначается через 2 n-p , если p переменных приравнены к соответствующим произведениям переменных.
Для правильного планирования ДФЭ необходимо использовать все полученные ранее сведения об объекте теоретического и интуитивного характера и выделить из них те переменные и произведения переменных, влияние которых на процесс минимально. При этом смешивание нужно производить так, чтобы основные оценки b 0 ,b 1 ,...,b n были смешаны с взаимодействиями, о которых заранее известно, что они не оказывают влияния на объект. Следовательно, произвольное разбиение матрицы планирования 2 3 на две части выделения полуреплики типа 2 3-1 недопустимо.
Генерирующее соотношение
служит для построения дробной реплики. Так, в рассмотренном планировании 2 3-1 мы задавали полуреплику типа 2 3 с помощью генерирующего соотношения x 3 = x 1 x 2 .
Определяющим контрастом
(ОК) называется соотношение, задающее элемент первого столбца матрицы планирования для фиктивной переменной (все они равны 1). Выражение ОК в нашем примере получается умножением левой и правой частей приведенного генерирующего соотношения на его левую часть x 3
1 = x 1 x 2 x 3 ,
так как всегда x 2 ig = 1.
Знание ОК позволяет определить всю систему совместных оценок не изучая матрицу планирования ДФЭ. Соотношения, задающие эти оценки, можно найти, последовательно перемножив независимые переменные на ОК
x 1 = x 2 x 3 ; x 2 = x 1 x 3 ; x 3 = x 1 x 2 .
Отсюда легко находим смешиваемые теоретические коэффициенты регрессии и их оценки
b 1 ® b 1 + b 23 ; b 2 ® b 2 + b 13 ; b 3 ® b 3 + b 12 .
Разрешающая способность дробных реплик определяется генерирующими соотношениями. Она тем выше, чем выше порядок взаимодействий, с которыми смешаны линейные коэффициенты, и увеличивается с ростом числа независимых переменных.
Для четверти реплики в пятифакторном планировании типа 25-2 могут быть заданы, например генерирующее соотношение
x 4 = x 1 x 2 x 3 ; x 5 = x 1 x 2
заранее полагая, что b 123 = b 12 = 0, т.е. что пара x 1 x 2 и тройка x 1 x 2 x 3 не дает значимого эффекта взаимодействия. Определяющими контрастами для этой четверть-реплики согласно вышеприведенным правилам будут соотношения
1 = x 1 x 2 x 3 x 4 ; 1 = x 1 x 2 x 5 .
Если у дробной реплики имеются два и более определяющих контраста, их необходимо перемножить между собой, используя все возможные комбинации. В случае четвертьреплики получается одна комбинация
1 = x 3 x 4 x 5
Обобщающий определяющий контраст
, построенный на основе всех полученных определяющих контрастов, полностью характеризует разрешающую способность реплик высокой степени дробности
1 = x 1 x 2 x 3 x 4 = x 1 x 2 x 5 = x 3 x 4 x 5 .
Совместные оценки здесь будут определяться соответствиями
x 0 = x 1 x 2 x 3 x 4 = x 1 x 2 x 5 = x 3 x 4 x 5 ; x 1 = x 2 x 3 x 4 = x 2 x 5 = x 1 x 3 x 4 x 5 ; x 2 = x 1 x 3 x 4 = x 1 x 5 = x 2 x 3 x 4 x 5 ; x 3 = x 1 x 2 x 4 = x 1 x 2 x 3 x 5 =x 4 x 5 ; x 4 = x 1 x 2 x 3 = x 1 x 2 x 4 x 5 =x 3 x 5 ; x 5 = x 1 x 2 x 3 x 4 x 5 = x 1 x 2 = x 3 x 4 ; x 1 x 3 = x 2 x 4 = x 2 x 3 x 5 = x 1 x 4 x 5 ; x 2 x 3 = x 1 x 4 = x 1 x 3 x 5 =x 2 x 4 x 5 ;
Эти соотношения позволяют установить, оценкой каких теоретических коэффициентов является тот или иной коэффициент регрессии, полученный при обработке результатов эксперимента
b 0 = b 0 + b 1234 + b 125 + b 345 ; b 1 = b 1 + b 234 + b 25 + b 1345 ; b 2 = b 2 + b 134 + b 15 + b 2345 ; b 3 = b 3 + b 124 + b 1235 + b 45 ; b 4 = b 4 + b 123 + b 1245 + b 35 ; b 5 = b 5 + b 12345 + b 12 + b 34 ; b 13 = b 13 + b 24 + b 235 + b 145 ; b 23 = b 23 + b 14 + b 135 + b 245 ;
Разрешающая способность этой четверти реплики невысокая, так как все линейные коэффициенты смешаны с парными взаимодействиями. Матрица планирования такой четверти реплики представлена в табл.4.
Следует иметь в виду, что ДФЭ всегда можно дополнить до ПФЭ, реализовав недостающие дробные реплики.
Вся дальнейшая работа по реализации матрицы планирования ДФЭ, проверке воспроизводимости полученных результатов, определению оценок коэффициентов регрессии и их значимости, проверке адекватности полученной математической модели не отличается от соответствующих процедур в ПФЭ.
Таблица 4
Четверть реплики от ПФЭ типа 2 5 (планирование типа 2 5-2 )
| g
|
z
0 |
z
1 |
z
2 |
z
3 |
z
4 |
z
5 |
z
6 |
z
7 |
z
8 |
z
9 |
z
10 |
z
11 |
z
12 |
z
13 |
z
14 |
z
15 |
z
16 |
z
17 |
z
18 |
z
19 |
z
20 |
z
21 |
z
22 |
z
23 |
Z
24 |
z
25 |
z
26 |
z
27 |
z
28 |
z
29 |
z
30 |
z
31 |
| x
0 |
x
1 |
x
2 |
x
3 |
x
4 |
x
5 |
x
1 x 2 |
x
1 x 3 |
x
1 x 4 |
x
1 x 5 |
x
2 x 3 |
x
2 x 4 |
x
2 x 5 |
x
3 x 4 |
x
3 x 5 |
x
4 x 5 |
x
1 x 2 x 3 |
x
1 x 2 x 4 |
x
1 x 2 x 5 |
x
1 x 3 x 4 |
x
1 x 3 x 5 |
x
1 x 4 x 5 |
x
2 x 3 x 4 |
x
2 x 3 x 5 |
x
2 x 4 x 5 |
x
3 x 4 x 5 |
x
1 x 2 x 3 x 4 |
x
1 x 2 x 3 x 5 |
x
1 x 2 x 4 x 5 |
x
1 x 3 x 4 x 5 |
x
2 x 3 x 4 x 5 |
x
1 x 2 x 3 x 4 x 5 |
|
| 1
|
+ | - | - | - | - | + | + | + | + | - | + | + | - | + | - | - | - | - | + | - | + | + | - | + | + | + | + | - | - | - | - | + |
| 2
|
+ | + | - | - | + | - | - | - | + | - | + | - | + | - | + | - | + | - | + | - | + | - | + | - | + | + | + | - | + | + | - | - |
| 3
|
+ | - | + | - | + | - | - | + | - | + | - | + | - | - | + | - | + | - | + | + | - | + | - | + | - | + | + | - | + | - | + | - |
| 4
|
+ | + | + | - | - | + | + | - | - | + | - | - | + | + | - | - | - | - | + | + | - | - | + | - | - | + | + | - | - | + | + | + |
| 5
|
+ | - | - | + | + | + | + | - | - | - | - | - | - | + | + | + | + | + | + | - | - | - | - | - | - | + | + | + | + | - | - | + |
| 6
|
+ | + | - | + | - | - | - | + | - | - | - | + | + | - | - | + | - | + | + | - | - | + | + | + | - | + | + | + | - | + | - | - |
| 7
|
+ | - | + | + | - | - | - | - | + | + | + | - | - | - | - | + | - | + | + | + | + | - | - | - | + | + | + | + | - | - | + | - |
| 8
|
+ | + | + | + | + | + | + | + | + | + | + | + | + | + | + | + | + | + | + | + | + | + | + | + | + | + | + | + | + | + | + | + |
Пример 2. Методом ДФЭ найти математическую модель процесса напыления резисторов.
Воспользуемся результатами Примера 1 и положим в качестве генерирующего соотношения равениство x 1 = x 2 x 3 (т.к. b 23 = 0). Тогда матрица планирования и результаты эксперимента (опуская промежуточные данные) будут выглядеть так
| g | x 0 | x 1 | x 2 | x 3 | S 2 g | (-) 2 | ||
| 1 | + | + | - | - | 17,34 | 2,228 | 17,41 | 0,0049 |
| 2 | + | - | + | - | 10,72 | 1,387 | 10,77 | 0,0025 |
| 3 | + | - | - | + | 13,70 | 0,950 | 13,65 | 0,0025 |
| 4 | + | + | + | + | 14,58 | 4,227 | 14,53 | 0,0025 |
Проверим воспроизводимость опытов
откуда следует, что результаты опытов получены правильно, дисперсия строчных выборок равна S 2 {y} = 8,792 / 4 = 2,198 с числом степеней свободы v 3 = 4·4 = 16.
Оценки коэффициентов уравнения регрессии
| ; |
аналогично b 2 = -1,44; b 3 = 0,05.
Проверка значимости полученных оценок начинается с определения их СКО
откуда
| ; | ; |
Табличные значения критерия t кр (5%;16) = 2,131, следовательно, модель найдена в виде
= 14,09 + 1,88x 1 - 1,44x 2 .
Проверка адекватности модели дает
| , откуда | , |
т.е. модель признается адекватной экспериментальным данным.
Сравнение моделей примера 1 и примера 2 показывает, что они имеют совершенно разный вид, а по некоторым факторам - противоположные по смыслу оценки коэффициентов. Отсюда можно сделать несколько общих выводов и рекомендаций (без подробного обоснования), пригодных для использования в рамках теории планирования экспериментов:
по одним и тем же экспериментальным данным можно построить несколько математических моделей, каждая из которых будет адекватна для своего набора оценок коэффициентов регрессии;
из всех моделей наилучшей признается та, у которой меньше членов и меньше критерий Фишера (или, если угодно, меньше дисперсия адекватности);
при большом числе факторов работу по математическому моделированию следует начинать с ДФЭ возможно большей дробности. Если модель получилась неадекватной, ее всегда можно достроить до следующей реплики вплоть до ПФЭ. Это сэкономит количество опытов, время, затраты и т.п.
Заключение .
Применение описанных выше методов математического моделирования полностью оправдало себя в условиях с небольшим числом факторов. Но при очень большом числе факторов и привлечение их к составлению математического описания исследуемого объекта методами ПФЭ или ДФЭ может потребовать увеличения объема экспериментальной работы, что редко может выполняться из-за экономических, технологических и прочих ограничений. Таким образом, возникает необходимость в предварительном отсеивании несущественных и выделении тех факторов процесса, которые оказывают наиболее заметное влияние на целевую функцию. Другим существенным затруднением для применения ПФЭ или ДФЭ в производственных условиях является метод получения оценок коэффициентов регрессии. Оценки вида (11) считаются оптимальными в смысле эффективности (минимума дисперсии), поскольку их вычисление базируется на методе наименьших квадратов, однако предварительным условием такой оптимальности являются требования независимости факторов, ортогональности и симметричности плана эксперимента, а также требование равенства дисперсий условных распределений плотности вероятности f(y/x k ). В свою очередь симметричность плана требует равного количества наблюдений, соответствующих положительным и отрицательным значениям k-го фактора.
На практике в производственных условиях требования симметричности плана и равенства дисперсий условных распределений плотности вероятности f(y/x k ) эксперимента, как правило, нарушаются, особенно в случаях, когда исследователь пытается построить модель по результатам, зафиксированными для случайной системы комбинаций производственных факторов. При этом всегда имеется выбор: либо нарушить одно из требований факторного анализа, либо потерять часть информации, пытаясь выбрать из нее только то, что согласуется с правилами ведения ПФЭ (ДФЭ).