Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Автоматизация научных исследований..pdf
Скачиваний:
17
Добавлен:
15.11.2022
Размер:
6.27 Mб
Скачать

Дисперсия o?(jr) в этом случае зависит не только от расстояния

между точкой х и центром плана; следовательно, план X не является ротатабельным.

План X называется ненасыщенным, если N > к +1, и насыщен­ ным, если N = к + ]. Здесь N - число экспериментов в плане; +1) - число оцениваемых коэффициентов.

6.3. Общие требования к плану эксперимента. О критериях планирования эксперимента

В главе 5 мы при определении оценок параметров модели иссле­ дуемого объекта исходим из того, что результаты наблюдений полу­ чены в некотором заданном множестве точек. На практике, однако, часто оказывается возможным свободно выбирать условия проведе­ ния опытов в пределах некоторых границ. Выбор числа и условий проведения экспериментов, обеспечивающих получение наилучшего в определенном смысле результата исследования, составляет цель планирования эксперимента. Разработан ряд критериев оптимально­ сти планов эксперимента, важнейшие из которых будут рассмотрены ниже. Следует отметить, что оптимальный выбор плана эксперимента существенным образом зависит от конкретных особенностей иссле­ дуемого объекта, таких как вид его модели, стоимость отдельных опытов, области варьирования независимых переменных и т.п.

Одной из важнейших характеристик плана, влияющей, с одной стороны, на стоимость и длительность исследования, а с другой - на точность результатов, является число экспериментов N

Заметим, что план с минимально возможным числом эксперимен­ тов N = к + \ (насыщенный план) не позволяет проверить адекват­ ность модели. Поэтому обычно выбирают N > к +1, где к +1 - число оцениваемых параметров модели исследуемого объекта.

Важное значение для оценки качества плана эксперимента имеет вид информационной матрицы плана

М = FTF.

Матрица М должна быть невырожденной, г.е. |м| Ф0, где |м| - оп­ ределитель матрицы М . Только в этом случае формула (6.2) имеет

единственное решение вида

 

a = (FTF y 'F TY

(6.28)

Перечислим теперь некоторые критерии планирования экспери­ мента, используемые в практических исследованиях.

1. Критерий ортогональности плана. Этот критерий требует та­ кого выбора плана X для оценки коэффициентов модели заданного вида, при котором информационная матрица плана диагональная.

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

(FTF)~' является диагональной, и, следовательно, ковариации оценок коэффициентов равны нулю. Это значит, например, что замена нулем любого коэффициента в уравнении модели не изменит значений оце­ нок остальных коэффициентов. Такое свойство ортогональных планов оказывается очень полезным, когда точный вид модели исследуемого объекта неизвестен и исследователь использует экспериментальные данные для отбора переменных, существенно влияющих на выходную величину.

2. Критерий ротатабельности. Этот критерий требует такого расположения экспериментальных точек в области планирования Qx,

при котором дисперсия <т?(х) оценки значений зависимой переменной

в точке х зависит только от расстояния от этой точки до центра пла­ на. Такой критерий хорошо согласуется с требованием равнозначно­ сти (с точки зрения точности оценки зависимой переменной) всех на­ правлений от центра плана.

Названные критерии обеспечивают некоторые полезные и удобные свойства оценок коэффициентов, однако они никак не связаны с требо­ ванием максимальной точности построения модели. Критерии, приво­ димые ниже, обеспечивают оптимальность планов эксперимента с точки зрения точности оценки параметров модели или зависимой переменной.

3. Критерий А-оптималыюсти. Этот критерий требует такого выбора плана, при котором матрица С = (FTF y ] имеет минимальный след (т.е. сумма диагональных элементов матрицы С минимальна). Поскольку в соответствии с (5.40) диагональный элемент сп матрицы

С пропорционален дисперсии оценки z-го коэффициента, то критерий А-оптимальности, по существу, требует минимизации средней дис­ персии оценок коэффициентов модели исследуемого объекта.

4. Критерий Д-оптимальности. Этот критерий требует такого

расположения точек в области Qt , при котором определитель матри­

цы С = (F'F)~' минимален (или определитель матрицы М = F1F

максимален). £>-оптимальный план минимизирует объем эллипсоида рассеяния оценок коэффициентов модели исследуемого объекта.

5. Критерий G-оптимальности. Критерий G-оптимальности тре­ бует такого расположения точек в области Qt , при котором достигает­ ся наименьшая величина максимальной дисперсии оценки зависимой переменной в области Qv В отличие от критериев А- и D-опти- мальности, связанных с точностью нахождения коэффициентов моде­ ли, критерий G-оптимальности требует максимальной точности оценки зависимой переменной.

Все перечисленные критерии связаны с предположением, что вид модели исследуемого объекта известен. Однако на практике часто возникает такая ситуация, когда исследователь не знает истинного вида модели. В таком случае эксперимент обычно сначала планирует­ ся, исходя из простейшего предположения о линейности модели от­ носительно варьируемых переменных. После проведения опытов про­ веряется адекватность линейной модели. Если модель адекватна, экс­ перимент заканчивают. В противном случае переходят к построению модели более высокого порядка (например, к построению квадратич­ ной модели) и так далее до получения адекватной модели.

На практике часто полезно стремиться к тому, чтобы один и тот же план удовлетворял одновременно ряду критериев, например, был бы D-оптимальным и ортогональным и т.п. В общем случае такого сочетания свойств не наблюдается, однако иногда в одном плане со­ четается ряд полезных свойств. Например, легко убедиться в справед­ ливости следующего утверждения: ортогональный план для линейной модели (5.3) одновременно является и ротатабельным, если

(FrF)~' = С = g l, g - константа, I - единичная матрица.

В самом деле, согласно (6.20) дисперсия оценки зависимой пере­ менной в точке плана х для таких планов определяется как

4 .0 = / Т|W fJ f y■'fix )* 2 = go2(i + i x?)