Математическое моделирование при активном эксперименте (стр. 3 из 6)

Таблица 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 1x 3	x 2x 3	x 1x 2x 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 может быть обусловлена следующими причинами:

1. уровень базового режима
   * близок к точке частного экстремума по переменной X i или по произведению переменных;
2. шаг варьирования D X i выбран малым;
3. данная переменная (или произведение переменных) не имеет функциональной связи с выходным параметром Y;
4. велика ошибка эксперимента вследствие наличия неуправляемых и неконтролируемых переменных.

Поскольку ортогональное планирование позволяет определять доверительные границы для каждого из коэффициентов регрессии в отдельности, то, если какой-либо из коэффициентов окажется незначимым, он может быть отброшен без пересчета всех остальных. После этого математическая модель объекта составляется в виде уравнения связи выходного параметра 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. Методом ПФЭ найти математическую модель процесса напыления резисторов.

После консультации с экспертами и некоторых предварительных исследований было определено, что на величину сопротивления напыляемых резисторов могут оказывать влияние следующие факторы: