В каждом из N опытов (номер опыта и = 1+N) оценка дисперсии единичного результата равна
где ти – число повторностей и-го опыта.
Средневзвешенная оценка дисперсии единичного результата рассчитывается по всем оценкам дисперсии единичного результата опытов:
а) при различных ти
где - число степеней свободы средневзвешенной оценки дисперсии; ти– 1 = fu– «вес» соответствующей и-ой оценки дисперсии, равный числу степеней свободы fu;
б) прити = т = const
где N(m-1)=f– число степеней свободы средневзвешенной оценки дисперсии.
Прежде чем пользоваться соотношениями (28) и (29) для расчета средневзвешенных уточненных оценок дисперсии (чем больше число степеней свободы, тем более точной будет оценка дисперсии), надо доказать однородность исходных оценок дисперсии.
Определение «однородные» в статистике означает «являющиеся оценкой одного и того же параметра» (в данном случае – дисперсии σ 2).
Если измеряемая случайная величина уикраспределена по нормальному закону во всем исследуемом диапазоне, то независимо от значений
идисперсия σ не будет изменять своей величины и оценки этой дисперсии должны быть однородными. Однородность этих оценок проявляется в том, что они могут отличаться друг от друга лишь незначительно, в пределах, зависящих от принятой вероятности и объема экспериментальных данных.Если ти = т и f = const, то однородность оценок дисперсий можно проанализировать при помощи критерия Кохрена Gkp. Вычисляют отношение максимальной дисперсии S2(yuk)maxк сумме всех дисперсий
и сравнивают это отношение с величиной критерия Кохрена Gkp (P; f; N). Если G < Gkp, то оценки однородны.
Таблица значений критерия Кохрена в зависимости от числа степеней свободы числителя fu, числа сравниваемых дисперсий N и принятого уровня значимости q= 1 – Р дана в приложении.
Если число повторностей в опытах различно (flt ≠ const), однородность оценок дисперсии можно проанализировать с помощью критерия Фишера FТ.Для этого из N оценок дисперсии выбирают 2: максимальную S2(yuk)max и минимальную S2(yuk)min. Если вычисленное значение Fих отношения меньше Ft,
то все N оценок дисперсии будут однородны.
Значения критерия Фишера FTданы в приложении в зависимости от принятого уровня значимости qи числа степеней свободы f1иf2 оценок S2(yuk)max и S2(yuk)min соответственно.
Если оценки дисперсии непосредственно измеряемого параметра у оказались неоднородными, т.е. оценками различных дисперсий, то средневзвешенная оценка не может быть рассчитана. И кроме того, величины укуже нельзя считать подчиняющимися нормальному закону, при котором дисперсия может быть лишь одной и неизменной при любом у.
Причиной нарушения нормального закона распределения может быть наличие оставшихся грубых ошибок (анализ на грубые ошибки либо не проводился, либо проведен недостаточно тщательно).
Другой причиной может быть наличие активного фактора, ошибочно отнесенного исследователем к неактивным и не снабженного системой стабилизации. Поскольку условия изменились, этот фактор стал значимо влиять на процесс.
9. Планирование и обработка результатов однофакторных экспериментов
9.1 Формализация экспериментальных данных методом наименьших квадратов
Влияние какого-либо фактора на выход процесса может быть выражено зависимостью у = f(C). Если конкретному значению Сисоответствует единственное значение уи, то такая зависимость называется функциональной. Эту зависимость получают путем строгих логических доказательств, не нуждающихся в опытной проверке. Например, площадь квадрата ω может быть представлена функциональной зависимостью от размера стороны квадрата а: ω = а2.
Если уиостается неизменным в то время как Сиизменяется, то у не зависит от С. Например, угол при вершине квадрата равный π/2, не зависит от размера стороны аи.
Если для оценки величин уии Сииспользуются данные наблюдений, величины случайные, то функциональная зависимость между ними существовать не может.
Измерив отдельно сторону а и площадь ω квадрата, можно убедиться, что полученные результаты не могут быть представлены с абсолютной точностью зависимостью ω = а 2.
К формализации экспериментальных данных, т.е. построению по ним описывающей процесс зависимости, исследователь прибегает, когда не может составить эвристическую (детерминированную) математическую модель из-за недостаточного понимания механизма процесса или его чрезмерной сложности.
Полученная в результате формализации экспериментальных данных эмпирическая математическая модель имеет меньшую ценность, чем отражающая механизм процесса эвристическая математическая модель, которая может предсказать поведение объекта за пределами изученного диапазона изменения переменных.
Приступая к эксперименту с целью получения эмпирической математической модели, исследователь должен определить необходимый объем опытных данных с учетом количества принятых к исследованию факторов, воспроизводимости процесса, предполагаемой структуры модели и обеспечения возможности проверки адекватности уравнения.
Если по результатам эксперимента, состоящего из двух опытов, получено линейное однофакторное уравнение у = b0+ b1С, то построенная по этому уравнению прямая обязательно пройдет через эти экспериментальные точки. Следовательно, для того чтобы проверить, насколько хорошо эта зависимость описывает данный процесс, надо поставить опыт хотя бы еще в одной точке. Этот дополнительный опыт дает возможность осуществить корректную процедуру проверки пригодности уравнения. Однако проверку обычно проводят не по одной дополнительной точке, которая не участвовала в определении коэффициентов уравнения, а по всем экспериментальным точкам, число которых (N) должно превышать число коэффициентов уравнения (N')
Так как N > N', решение такой системы требует специального подхода.
9.2 Симметричный и равномерный план однофакторного эксперимента
Задача в значительной степени упростится, если при планировании эксперимента, можно будет обеспечить условие:
ΣCu=0 (1)
При натуральной размерности факторов выполнить условие ΣCu=0 невозможно, т. к. в этом случае величина фактора должна иметь как положительные значения, так и отрицательные.
Если же точку отсчета величины фактора перенести в середину диапазона изменения фактора (центр эксперимента)
то появляется возможность удовлетворить условию
в виде , где С'u=Сu – С0.Для равномерного плана Сu – С(u-1) = λ = const,
где λ – интервал варьирования фактора.
Условие
может быть выполнено, если для обозначения величины фактора использовать безразмерные выражения:xu =
,отсюда легко увидеть, что условие
эквивалентно условию и такие планы называют симметричными.При составлении плана диапазон фактора ориентировочно ограничивают величинами Сmin и Сmax, назначенными после изучения литературы по теме исследования. От опыта к опыту предусматривают такое изменение величины фактора, которое позволило бы достоверно уловить имеющимися в распоряжении исследователя приборами изменение выхода процесса
.С учетом величины λ и диапазона (Сmax– Cmin) определяют число опытов, округляя его до нечетного N:
Затем определяют величины факторов в каждом из N опытов и уточняют исследуемый диапазон фактора СN – С1:
= ,где хu – безразмерное выражение фактора, аналогичное полученному по соотношению
Для расчета коэффициентов уравнения используем формулу:
,множители аjuи знаменатель lj берем из приложения.
Число опытов эксперимента может быть четным или нечетным, и, как правило, должно быть больше числа коэффициентов N' уравнения.
Чем больше разность (N – N'), тем с большей точностью можно получить оценки коэффициентов данного уравнения и тем в большей степени эти оценки будут освобождены от влияния случайных неуточненных факторов.