6
исследуемая величина (величины) «y» зависит от многих факторов x1, x2,..., то есть задача отыскания математической модели
y f x1 , x2 ,...,xn .
Традиционный подход в определении вида модели заключался в проведении однофакторных экспериментов по определению сначала частных зависимостей y1 f1 x1 , y2 f 2 x2 , ... yn f n xn , когда меняется поочередно один фактор, а остальные фиксируются, а затем уже происходит обобщение математической модели.
Однако с ростом числа факторов и их уровней варьирования катастрофически возрастает число необходимых опытов. Это возрастание получило образное название «проклятия размерностей» [2]. Полное количество опытов определяется по формуле
n qk ,
где q – число уровней варьирования каждого фактора; k – количество факторов.
Тогда при |
q=5 и k=2 |
n=52 =25, |
а при k=6 (любая электротехническая задача) |
|||||
n=56=15625! |
|
|
|
|
|
|||
|
Поэтому раньше экспериментальные исследования занимали много времени. |
|||||||
Например, |
полное |
сопротивление |
в электротехнике определяется формулой |
|||||
|
|
|
|
|
|
, |
|
|
z |
r 2 x 2 |
|
, где |
r |
x L . |
На отыскание данной зависимости ушло |
||
|
|
|||||||
|
|
|
|
|
S |
|
|
|
несколько десятилетий. Другой пример – на установление закона равноускоренного движения S=at2/2. Галилей затратил более 25 лет своей жизни.
Проход через «проклятие размерностей» указал Рональд Фишер, английский статистик в 1925-1929 гг. Он предложил одновременно варьировать все факторы.
Однако этот подход не был задействован до начала 60-х годов прошлого века. Лишь в это время возникла настойчивая необходимость уменьшения затрат времени на проведение исследований. С этого момента теория планирования эксперимента (ТПЭ) начала и продолжает бурно развиваться до настоящего времени.
Анализ априорной информации
7
Перед тем, как проводить экспериментальные исследования, необходимо выявить всю имеющуюся и накопленную человечеством информацию по существу изучаемого вопроса.
Частично эта задача формализуется, а именно, большую часть информации можно выявить из информационных литературных и патентных исследований.
Однако большая часть работы по выявлению априорной (имеющейся до опыта) информации является «скрытой» от исследователя, и здесь необходимо «догадаться» использовать те или иные сведения, подчас из других областей знаний.
В ТПЭ все принимаемые решения, основанные на догадках, аналогиях, ассоциациях, различных косвенных соображениях, наводящих на мысль,
называются эвристическими [3].
Если для принятия решения нет никаких данных, и его приходится принимать на собственный страх и риск, то его называют волевым. Эти понятия являются в значительной степени субъективными, зависящими от таланта исследователя.
Важнейшими задачами при анализе априорной информации являются: а) выбор влияющих факторов и их интервалов варьирования; б) выбор вида математической модели.
При выборе влияющих факторов часто используется метод экспертных оценок (ранжирования), заключающийся в том, что различные специалисты после опроса оценивают по шкале значимости, по их мнению, существенные факторы. После анализа такой информации выбираются наиболее значимые факторы.
В качестве математической модели чаще всего используют аддитивную
(линейную) и мультипликативную модели:
yˆ b0 b1 x1 b2 x2 ... bk xk |
(1) |
или
yˆ b0 x1b1 x2b2 ...xkbk ,
где хi – факторы регрессионной модели, bi – коэффициенты регрессии, y – значение функции (отклик). Последняя модель после логарифмирования легко преобразуется в аддитивную модель.
Следует отметить, что большая часть физических явлений описывается мультипликативными моделями.
Таким образом, на этапе анализа априорной информации проводится выбор вида математической модели, выбор влияющих факторов и их интервалов
8
варьирования. Эти процедуры являются неформализуемыми, субъективными и сильно зависящими от таланта исследователя этапами работы.
Планирование эксперимента, его проведение и статистическая обработка
Вотличие от предыдущего следующие этапы хорошо разработаны для моделей 1-го и 2-го порядков. Этот этап является хорошо формализуемым, то есть для его проведения существуют методики, правила поведения и рекомендации для экспериментатора в тех или иных случаях. Для этого этапа имеются алгоритмы и пакеты прикладных программ. Построение планов эксперимента первого и второго порядка рассматриваются далее в этом разделе.
Вслучае использования аддитивной модели последующими этапами являются
[4]:
1) Составление матрицы планирования.
2) Проведение эксперимента согласно матрице плана эксперимента.
3) Статистическая обработка результатов эксперимента:
a)Поиск коэффициентов регрессии с использованием метода наименьших квадратов (МНК), суть которого заключается в нахождении коэффициентов регрессии, обеспечивающие минимум функции:
N |
y |
n |
y |
2 |
min , |
|
|
ˆ n |
|
|
n 1
где yn , yˆ n – соответственно экспериментальные значения отклика и рассчитанные по регрессионной модели.
Для линейчатой модели
|
k |
l |
|
y b0 |
bi xi |
bij xi x j |
(2) |
|
i 1 |
i, j 1;i j |
|
коэффициенты регрессии определяются по формулам:
|
N |
|
N |
|
N |
|
|
|
yn |
|
yn xin |
|
yn xin x jn |
|
|
b |
n 1 |
; b |
n 1 |
; b |
n 1 |
, |
(3) |
0 |
N n |
i |
N |
ij |
N |
|
|
|
|
|
|
|
где N – количество проведенных опытов;
n – порядковый номер опыта.
|
|
9 |
|
|
|
|
b) Проверка воспроизводимости |
|
опытов |
с помощью дисперсии |
|||
воспроизводимости: |
|
|
|
|
|
|
|
N |
Q |
yn 2 |
|
||
|
ynq |
|
||||
S 2 |
n 1q 1 |
|
|
, |
(4) |
|
|
|
|
||||
b |
N Q 1 |
|
|
|
||
|
|
|
|
|
||
где Q — число повторных опытов (одинаковое для всех опытов матрицы), q — порядковый номер повторных опытов;
Q
ynq
y |
|
|
q 1 |
; |
|
n |
Q |
||||
|
|
|
|||
|
|
|
|
c)Определение дисперсии адекватности модели:
N
yn yˆn 2
S 2 |
|
n 1 |
|
, |
|
|
|||
ад |
|
|
f |
|
|
|
|
||
где f N k 1 – число степеней свободы.
d)Определение адекватности модели с помощью критерия Фишера:
(5)
(6)
F |
Sад2 |
. |
(7) |
|
|||
|
Sb2 |
|
|
Модель адекватна (гипотеза об |
адекватности не |
отвергается), если |
|
FFт f ; N Q 1 ; , где α – выбранный уровень значимости принятия решения.
e)Проверка значимости коэффициентов регрессии с помощью t-критерия Стьюдента при выбранном уровне значимости α, если значение коэффициента
регрессии bi больше доверительного интервала |
bi, то можно считать, что гипотеза |
о значимости не отвергается, то есть при bi bi |
считают коэффициенты регрессии |
значимыми, где bi tinSb 
N . Незначимость некоторых коэффициентов регрессии может быть обусловлена близостью основного уровня для данного фактора к точке локального экстремума или малостью интервала варьирования;
f)Анализ дальнейших действий (либо переходят к другим интервалам, либо
кдругим планам эксперимента, либо останавливаются, либо обобщают локальные модели в общую по методу Брандона).
10
Рекомендации по выбору средств измерений
Повышение точности эксперимента и степени адекватности его модели может быть достигнуто не только использованием методики ТПЭ и увеличения числа повторных опытов, но и за счет использования более точной аппаратуры СИ. Однако здесь надо исходить из понятия эффективности эксперимента как отношения достигнутого эффекта к затратам при проведении эксперимента [3].
Разброс исходных данных ИД всегда складывается из 3-х составляющих: 0
– диффузности объекта измерения (невоспроизводимости или среднеквадратического отклонения воспроизводимости), М – погрешности адекватности модели и С.И . – погрешности средств измерений.
Если эти погрешности некоррелированы (статистически независимы), то
ИД 
20 2М 2С.И . .
Задача состоит в рациональном выборе средств измерений с точки зрения их точности. Для простоты рассуждений примем, что М С.И. и М 0 . Тогда возможны три случая:
a) Очень точная аппаратура С.И . 0 . При этом ИД 
20 2С.И . 0
и разброс полностью определяется диффузностью (нестабильностью параметнов) объекта. Это плохо, т.к. для усреднения результата надо провести большое число опытов, но излишне точная аппаратура требует и большего времени измерений.
|
|
|
|
|
2 |
20 |
|
|
|
|
||
Если |
С.И. 0 / 3, то ИД |
0 |
|
|
1,05 0 , а приборы можно использовать |
|||||||
9 |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
||
более |
грубые. |
Значит |
при |
|
С.И. 0 / 3 можно |
повысить эффективность |
||||||
эксперимента. |
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
||||||
|
b) При С.И . 0 |
ИД |
|
20 |
2С.И . 1,14 0 , т.е. разброс исходных данных |
|||||||
возрастает всего |
на 14% |
в сравнении со случаем |
С.И . 0 . В этом случае |
|||||||||
эффективность статистической обработки сильно возрастает с увеличением 
n . Однако надо помнить, что систематические погрешности ( М и часть С.И . ) не уменьшаются.
c) При С.И . 0 ИД 
2С.И . 20 С.И . . Если порядок ИД устраивает экспериментатора, то нет необходимости в организации многократных измерений и их статистическом усреднении.