Можно сформулировать три причины, по которым используются расщепленные участки:
–непрактичность или невозможность дифференциации (выделения) основных экспериментальных единиц для некоторых факторов;
–необходимость включения дополнительных факторов в уже проводящийся эксперимент;
–желание более точно оценить эффекты факторов частей объектов и их взаимодействий с факторами основных объектов при некоторой потере точности на факторах основных объектов.
При этом неправильным будет использование рандомизированной блоковой схемы для одного из нескольких факторов с последующим расщеплением каждого объекта для уровней второго фактора, дальнейшим расщеплением этих полученных подобъектов для уровней третьего фактора и т. д. Правильным будет использование одной и той же экспериментальной единицы, но при этом рассмотрение подобъектов в качестве объектов, а основных объектов – в качестве блоков, и смешивание полученных структур АВС.
При однофакторном эксперименте искомая математическая модель функции y = f(x) может быть найдена лишь в результате совместной обработки всех полученных значений x и y. Экспериментальные точки занимают некоторую полосу вокруг искомой кривой y = f (x), ограниченную штриховы-
y
y = f(x)
x 
Рис. 3.1. Графическое изображение однофакторной зависимости (сплошная линия) и ее полосы неопределенности (штриховые линии)
ми линиями (рис. 3.1).
Отклонения экспериментальных точек от номинальной кривой могут быть вызваны разными причинами: например, неполнотой модели и учитываемых факторов, случайным характером самих исследуемых процессов и т. д. При этом невозможно разделить отклонения, вызванные погрешностями измерения x, от погрешностей, вызванных неточностью измерения y. Поэтому описанием по-
16
грешности исходных данных может быть лишь указание ширины полосы их разброса вокруг найденной кривой y = f (x). При этом полоса разброса не обязательно будет иметь постоянную ширину по всей своей длине и одинаковую конфигурацию (штриховые линии на рис. 3.1). Вопрос о форме полосы погрешностей должен анализироваться в каждом конкретном случае, так как одна и та же кривая с различной степенью приближения может быть аппроксимирована различными функциями (полиномом степени n, параболой, гиперболой, эллипсом и т. д.).
Основной помехой при установлении вида исследуемой зависимости является случайный разброс экспериментальных данных. Если случайный разброс исходных данных практически отсутствует (диффузность экспериментальных данных мала), то кривая просто проводится через эти точки
(рис. 3.2).
y |
y |
y |
|
y = f(x) |
|
|
|
|
|
|
|
x |
|
|
|
x |
|
|
|
x |
|
Рис. 3.2. Аппроксимирующая кривая |
|
Рис. 3.3. Экспериментальные данные |
|||||
|
при малой диффузности |
||||||
|
|
с большой диффузностью |
|||||
|
экспериментальных данных |
|
|||||
|
|
|
|
||||
При этом может оказаться, что некоторые точки все-таки не лежат на этой кривой, поэтому их следует рассматривать как промахи.
Если диффузность экспериментальных данных значительна (рис. 3.3), то для обработки таких данных привлекаются различные статистические методы установления вида однофакторных зависимостей.
Плавной линией обводится контур (рис. 3.4), включающий все экспериментальные точки (штриховые линии), а искомая кривая будет осевой линией контура (сплошная кривая). Если при этом для сохранения плавности границ контура какие-то точки необходимо оставить за пределами границ, то их
17
y |
|
y |
|
y = f(x) |
следует рассматривать как промахи или |
||
|
|
||||||
|
|
|
аномально большую случайную по- |
||||
|
|
|
|
||||
|
|
|
|||||
|
|
|
|
|
|
|
грешность. |
|
|
|
|
|
|
|
Однако при большом рассеянии |
|
|
|
|
|
|
|
результатов эксперимента форма кон- |
|
|
|
|
|
|
|
тура может иметь случайные, бессмыс- |
|
|
|
|
|
|
|
|
|
|
|
|
x |
|
ленные очертания. В этом случае огра- |
|
|
|
|
|
|
|
|
ничиваются лишь установлением уров- |
|
|
|
|
x |
|||
|
Рис. 3.4. Метод обведения контура |
ня и наклона искомой зависимости, по- |
|||||
|
|
|
|
|
|
|
|
лагая ее прямой линией, проходящей по центру обведенной контуром полосы точек.
Данный метод применяется при очень большой диффузности экспериментальных точек (рис. 3.5).
y |
|
|
|
|
|
Обведенное контуром поле точек |
|
y |
|
|
|||
|
|
|
||||
|
|
y = f(x) |
делят на несколько частей и в каждой |
|||
|
|
|
|
|
||
|
|
|
|
|
||
|
|
|
|
|
|
из них находят медианный центр (пе- |
|
|
|
|
|
|
|
|
|
|
|
|
|
ресечение горизонтали и вертикали, |
|
|
|
|
|
|
слева и справа и ниже и выше которых |
|
|
|
|
|
|
оказывается равное количество точек), |
|
|
|
|
|
x |
а затем через эти медианные центры |
|
|
|
|
|
проводят плавную кривую. Если пред- |
|
|
|
|
|
|
x |
|
|
|
Рис. 3.5. Метод медианных центров |
полагается, что зависимость может |
|||
быть прямой линией, то достаточно разделить поле на две области; если не прямой – от трех до пяти областей. Достоинством данного метода является его нечувствительность к промахам.
Рассматриваемую функцию y = f (x) полагают состоящей из двух слага-
емых: y = f0 (x)+ f1(x). |
|
|
|
|
|
|
|
||||||
|
|
На рис. 3.6 участок I представляется в виде y = y0 + f1(x). Затем выде- |
|||||||||||
ляется функция |
f1(x) как |
f1(x) = y − y0 , и этот остаток полагается равным |
|||||||||||
f |
1 |
(x) = f |
2 |
(x)+ f |
3 |
(x). Участок II |
полагают равным f |
2 |
(x) = a x2 . Теперь |
||||
|
|
|
|
|
|
|
|
|
1 |
||||
можно построить остаток f |
3 |
(x) = y − y |
0 |
−a x2 и найти его форму, и т. д. |
|||||||||
|
|
|
|
|
|
|
|
|
1 |
|
|
||
18
Полученный многочлен и будет ап- |
y |
|
y |
|
|
||||
проксимирующей функцией. При этом |
|
|
||
|
|
II |
||
всегда аппроксимирующая функция будет |
|
|
||
определенным приближением реальной |
|
|
|
|
функции, следовательно, будет присут- |
|
|
I |
|
ствовать погрешность, называемая по- |
y0 |
|
|
|
грешностью адекватности. |
|
|
|
|
|
|
x |
||
Погрешность адекватности моде- |
|
|
|
|
|
|
|
x |
|
|
|
|
||
ли – погрешность в описании данного явления, возникающая из-за недостаточного соответствия аппроксимирующей
функции всем особенностям формы экспериментальной кривой.
Для достижения наилучшего соответствия модели описываемому явлению обычно приходится усложнять модель, что приводит к противоречию между компактностью модели и точностью описания экспериментальных данных. Рациональным будет прекращение усложнения модели, когда она еще относительно проста и при этом имеет удовлетворительную погрешность адекватности. Ориентиром здесь может служить примерное равенство погрешности адекватности принятой модели и ширины полосы ее неопределенности вследствие случайного разброса исходных данных, т. е. принимаемая модель не противоречит данному полю экспериментальных данных.
После того как вид графической зависимости установлен, встает задача выбора ее математического описания. Обычно используют наиболее компактные выражения, представляющие собой элементарные функции (степенную, показательную, дробно-рациональную).
Чтобы проверить, является ли данная функция степенной вида y = axn , необходимо прологарифмировать правую и левую части и найти величины lg x и lg y . Это будет прямая в координатах lg y −lg x (рис. 3.7).
Если поле экспериментальных точек группируется относительно прямой линии, то можно говорить о том, что модель в виде степенной функции не противоречит данному полю точек и выбор данной функции может быть принят окончательно. Тогда точка пересечения прямой с осью lg y будет
определять величину a, а наклон прямой к оси lg x – величину n.
19
lg y |
|
lg y |
|
|
|
Если же экспериментальные точки не |
||||
|
|
|
лежат вдоль прямой, то модель в виде сте- |
|||||||
|
|
|
|
|
||||||
|
|
|
|
|
пенной функции не является адекватной, и |
|||||
|
|
|
|
|
необходимо перейти к проверке, является |
|||||
|
|
|
|
|
ли |
данная |
функция |
показательной |
||
|
|
|
|
lg x |
(y = aebx ). Для проверки соответствия по- |
|||||
|
|
|
|
казательной функции экспериментальным |
||||||
|
|
|
|
lg x |
данным используют полулогарифмический |
|||||
|
Рис. 3.7. Степенная функция |
масштаб, т. е. строят график в координатах |
||||||||
в логарифмических координатах |
||||||||||
|
|
(lg y = lg a +n lg x) |
|
ln y − x : |
ln y =ln a +bx . Если эксперимен- |
|||||
тальные точки укладываются на прямую линию, то можно говорить о соот- |
||||||||||
ветствии искомой модели показательной функции, и точка пересечения пря- |
||||||||||
мой с осью ln y будет определять параметр a, |
а наклон прямой к оси х – па- |
|||||||||
раметр b. |
|
|
|
|
|
|
|
|
||
|
|
Чтобы проверить, является ли искомая функция дробно-рациональной |
||||||||
|
= |
axm |
|
|
|
|
1 |
− x (рис. 3.8). |
|
|
y |
|
, строят график в координатах |
|
|
||||||
|
|
b + cx |
n |
|
|
|
y |
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
3 |
1/y |
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
c/a |
|
|
|
|
|
|
|
|
|
– с |
0 |
|
|
b |
|
x |
|
|
|
|
Рис. 3.8. Дробно-рациональная функция в координатах 1/y – x: |
|||||||
|
|
|
1 – для функции y = a/x; 2 – y = a/(c + x); 3 – y = a/(b – x) |
|
||||||
Замена координат x на 1/х и y на 1/y допустима лишь в том случае, если сдвиг по этой координате отсутствует. Если имеется сдвиг, то уравнение, например, гиперболы y = c + a/x в координатах 1/y – 1/х не будет прямой. В этом случае используют метод последовательных приближений.
Например, предполагают, что y = b ax− cx + d ; задают ряд возможных значений b; вычисляют значения 1/(x – b) и останавливаются на том значе-
20